Hello -
I am new to fipy but quite enjoying it thus far. I am running into
one problem with visualization. I would like to display the results
from solving a diffusion eqn on a 2D surface embedded in 3D, e.g.
a sphere. I was trying to get the mesh to display with MayaviViewer
but immediately run into the error below.
I assume numerix.sum used to have an index argument but no longer?
Beyond that, it looks like _getStructure is going to check the mesh
dimension and won't create the appropriate polygons for my surface mesh.
Instead it will see 3 dimensions and start looking for tetrahedra?
Any suggestions? Am I looking at the wrong viewer module all together?
regards,
charless
----------------------------
In [43]: mesh
Out[43]: GmshImporter2D()
In [44]: mesh.getDim()
Out[44]: 3
In [45]: from fipy.viewers.mayaviViewer.mayaviViewer import MayaviViewer
In [46]: viewer = MayaviViewer(phi)
---------------------------------------------------------------------------
exceptions.TypeError Traceback (most
recent call last)
/home/eecs/fowlkes/bioregistration/adjoint/<ipython console>
/usr/lib/python2.4/site-packages/fipy/viewers/mayaviViewer/mayaviViewer.py
in __init__(self, vars, limits, title)
99
100 for var in self.vars:
--> 101 self.structures.append(self._getStructure(var.getMesh()))
102
103 def _getStructure(self, mesh):
/usr/lib/python2.4/site-packages/fipy/viewers/mayaviViewer/mayaviViewer.py
in _getStructure(self, mesh)
106
107 from fipy.tools import numerix
--> 108 lengths = len(cellVertexIDs[0]) -
numerix.sum(numerix.MA.getmaskarray(cellVertexIDs), index = 1)
109
110 cellDict = {2 : [], 4: [], 6: [], 8: [], 'polygon' : []}
TypeError: sum() got an unexpected keyword argument 'index'
In [47]: python: Python/ceval.c:2531: PyEval_EvalCodeEx: Assertion
`tstate != ((void *)0)' failed.
Aborted
