Dear fipy users,
I hope you can help me to resolve what goes wrong in my script/ fipy
installation
or what I do wrong. I just want to calculate the partial derivatives (vs. x,y,
and z)
of my "phi". Thank you very much for your help!
Best wishes,
Ferenc
2D case:
>>> from fipy import *
>>> grid2D = Grid2D(nx=10, ny=10, dx=1.0, dy=1.0)
>>> phi = CellVariable(name="phi", mesh=grid2D,value=0.0)
>>> phi.setValue(GaussianNoiseVariable(mesh=grid2D,mean=0.5,variance=0.00001))
>>> source=phi.getGrad()[1]
No message in case of 2D grid!
3D case:
>>> from fipy import *
>>> grid2D = Grid3D(nx=10, ny=10, nz=10, dx=1.0, dy=1.0, dz=1.0)
>>> phi = CellVariable(name="phi", mesh=grid2D,value=0.0)
>>> phi.setValue(GaussianNoiseVariable(mesh=grid2D,mean=0.5,variance=0.00001))
>>> source=phi.getGrad()[1]
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File
"/home/x_ferta/PYTHON_MODULES/lib/python2.7/site-packages/FiPy-2.1.3-py2.7.egg/fipy/variables/variable.py",
line 1335, in __getitem__
unit=self.getUnit(),
File
"/home/x_ferta/PYTHON_MODULES/lib/python2.7/site-packages/FiPy-2.1.3-py2.7.egg/fipy/variables/variable.py",
line 244, in getUnit
return self._extractUnit(self.getValue())
File
"/home/x_ferta/PYTHON_MODULES/lib/python2.7/site-packages/FiPy-2.1.3-py2.7.egg/fipy/variables/variable.py",
line 471, in getValue
value = self._calcValue()
File
"/home/x_ferta/PYTHON_MODULES/lib/python2.7/site-packages/FiPy-2.1.3-py2.7.egg/fipy/variables/gaussCellGradVariable.py",
line 90, in _calcValue
return inline._optionalInline(self._calcValueIn, self._calcValuePy, N, M,
ids, orientations, volumes)
File
"/home/x_ferta/PYTHON_MODULES/lib/python2.7/site-packages/FiPy-2.1.3-py2.7.egg/fipy/tools/inline.py",
line 18, in _optionalInline
return pythonFn(*args)
File
"/home/x_ferta/PYTHON_MODULES/lib/python2.7/site-packages/FiPy-2.1.3-py2.7.egg/fipy/variables/gaussCellGradVariable.py",
line 75, in _calcValuePy
contributions = numerix.take(self.faceGradientContributions, ids, axis=1)
File
"/home/x_ferta/PYTHON_MODULES/lib/python2.7/site-packages/FiPy-2.1.3-py2.7.egg/fipy/tools/numerix.py",
line 969, in take
taken = a.take(indices, axis=axis)
File
"/home/x_ferta/PYTHON_MODULES/lib/python2.7/site-packages/FiPy-2.1.3-py2.7.egg/fipy/variables/variable.py",
line 1339, in take
return numerix.take(self.getValue(), ids, axis)
File
"/home/x_ferta/PYTHON_MODULES/lib/python2.7/site-packages/FiPy-2.1.3-py2.7.egg/fipy/variables/variable.py",
line 471, in getValue
value = self._calcValue()
File
"/home/x_ferta/PYTHON_MODULES/lib/python2.7/site-packages/FiPy-2.1.3-py2.7.egg/fipy/variables/faceGradContributionsVariable.py",
line 45, in _calcValue
return self.mesh._getAreaProjections() *
self.var.getArithmeticFaceValue().getNumericValue()
File
"/home/x_ferta/PYTHON_MODULES/lib/python2.7/site-packages/FiPy-2.1.3-py2.7.egg/fipy/variables/variable.py",
line 608, in getNumericValue
value = self.getValue()
File
"/home/x_ferta/PYTHON_MODULES/lib/python2.7/site-packages/FiPy-2.1.3-py2.7.egg/fipy/variables/variable.py",
line 471, in getValue
value = self._calcValue()
File
"/home/x_ferta/PYTHON_MODULES/lib/python2.7/site-packages/FiPy-2.1.3-py2.7.egg/fipy/variables/cellToFaceVariable.py",
line 47, in _calcValue
alpha = self.mesh._getFaceToCellDistanceRatio()
File
"/home/x_ferta/PYTHON_MODULES/lib/python2.7/site-packages/FiPy-2.1.3-py2.7.egg/fipy/meshes/numMesh/uniformGrid3D.py",
line 367, in _getFaceToCellDistanceRatio
numerix.ravel(YZdis.swapaxes(0,2))), axis=1)
IndexError: axis 1 out of bounds [0, 1)
>>>
--
Ferenc Tasnadi <[email protected]>
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
