Hi everyone,
thanks to you and your valuable help my code is now working nicely in 2D. I
thought the fipy structure would allow smooth transition to 3D case, as all
instructions are very general and mesh independant.
Unfortunately, I have kind of a black box error when evaluating phi (which
is the DistanceVariable of my level set problem)
to your knowledge, is there an obvious thing to modify when we go from a 2D
mesh to a 3D mesh ?
Julien
ERROR: An unexpected error occurred while tokenizing input
The following traceback may be corrupted or invalid
The error message is: ('EOF in multi-line statement', (16, 0))
---------------------------------------------------------------------------
IndexError Traceback (most recent call last)
/home/julien/Dropbox/3Dtest/fipy3Dlevelset.py in <module>()
269
--> 270 phi.calcDistanceFunction()
271
272
/usr/local/lib/python2.6/dist-packages/FiPy-2.2_dev-py2.6.egg/fipy/models/levelSet/distanceFunction/distanceVariable.pyc
in calcDistanceFunction(self, narrowBandWidth, deleteIslands)
242
243 """
--> 244 self._calcDistanceFunction(narrowBandWidth =
narrowBandWidth, deleteIslands = deleteIslands)
245 self._markFresh()
246
/usr/local/lib/python2.6/dist-packages/FiPy-2.2_dev-py2.6.egg/fipy/models/levelSet/distanceFunction/distanceVariable.pyc
in _calcDistanceFunction(self, extensionVariable, narrowBandWidth,
deleteIslands)
261 self._value = MA.where(tmp, -1, self._value)
262
--> 263 adjVals = numerix.take(self._value, cellToCellIDs)
264 adjInterfaceValues = MA.masked_array(adjVals, mask =
(adjVals * self._value) > 0)
265 dAP = self.mesh._cellToCellDistances
/usr/local/lib/python2.6/dist-packages/FiPy-2.2_dev-py2.6.egg/fipy/tools/numerix.pyc
in take(a, indices, axis, fill_value)
1020 ## `indices` is a masked array.
1021 ##
-> 1022 taken = MA.take(a, MA.filled(indices, 0), axis=axis)
1023
1024 mask = MA.getmask(indices)
/usr/lib/python2.6/dist-packages/numpy/ma/core.pyc in __call__(self, a,
*args, **params)
4844 method = getattr(narray(a, copy=False), self.__name__)
4845 try:
-> 4846 return method(*args, **params)
4847 except SystemError:
4848 return getattr(np,self.__name__).__call__(a, *args,
**params)
IndexError: index out of range for array
WARNING: Failure executing file: <fipy3Dlevelset.py>
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy]
