Hi, thanks for replies. On Mon, Dec 1, 2008 at 4:17 PM, Daniel Wheeler <[EMAIL PROTECTED]> wrote: > > Okay, I hadn't noticed that you are running "CURRENT" r2472. I went > back to "CURRENT" and it runs through at least some time steps of > "examples/cahnHilliard/inputTanh1D.py" without issue. What version of > numpy are you using?
ok, here are my versions: [EMAIL PROTECTED]:~/pywork/CURRENT> svn info Path: . URL: http://matforge.org/svn/fipy/tags/CURRENT Repository Root: http://matforge.org/svn/fipy Repository UUID: d80e17d7-ff13-0410-a124-85740d801063 Revision: 2865 Node Kind: directory Schedule: normal Last Changed Author: guyer Last Changed Rev: 2472 Last Changed Date: 2008-02-08 19:26:15 +0000 (Fri, 08 Feb 2008) ======= [EMAIL PROTECTED]:~/pywork/CURRENT> python -c "import matplotlib; print matplotlib.__version__" 0.98.3 ======= In [1]: import numpy In [2]: numpy.__version__ Out[2]: '1.2.1' In [3]: from numpy.core import ma as MA --------------------------------------------------------------------------- ImportError Traceback (most recent call last) /home/ionut/pywork/CURRENT/<ipython console> in <module>() ImportError: cannot import name ma In [4]: from numpy import ma as MA In [5]: ====== For the example which you mentioned "examples/cahnHilliard/inputTanh1D.py", it is started, I can see for two seconds the plot, and after that I receive: [EMAIL PROTECTED]:~/pywork/CURRENT> python examples/cahnHilliard/inputTanh1D.py /usr/local/lib64/python2.5/site-packages/fipy/viewers/matplotlibViewer/__init__.py:35: UserWarning: Matplotlib1DViewer efficiency is improved by setting the 'datamax' and 'datamin' keys return Matplotlib1DViewer(vars = vars, title = title, limits = limits) Traceback (most recent call last): File "examples/cahnHilliard/inputTanh1D.py", line 233, in <module> exec(fipy.tests.doctestPlus._getScript()) File "<string>", line 166, in <module> File "/usr/local/lib64/python2.5/site-packages/fipy/terms/term.py", line 121, in solve matrix, RHSvector = self.__buildMatrix(var, boundaryConditions, dt) File "/usr/local/lib64/python2.5/site-packages/fipy/terms/term.py", line 88, in __buildMatrix return self._buildMatrix(var, boundaryConditions, dt) File "/usr/local/lib64/python2.5/site-packages/fipy/terms/binaryTerm.py", line 66, in _buildMatrix termMatrix, termRHSvector = self.term2._buildMatrix(var, boundaryConditions, dt = dt) File "/usr/local/lib64/python2.5/site-packages/fipy/terms/binaryTerm.py", line 64, in _buildMatrix matrix, RHSvector = self.term1._buildMatrix(var, boundaryConditions, dt = dt) File "/usr/local/lib64/python2.5/site-packages/fipy/terms/diffusionTerm.py", line 272, in _buildMatrix self._getCoefficientMatrix(mesh, self.coeffDict['cell 1 diag']), numerix.zeros(N,'d')) File "/usr/local/lib64/python2.5/site-packages/fipy/terms/diffusionTerm.py", line 167, in _getCoefficientMatrix interiorCoeff = numerix.take(interiorCoeff, mesh._getCellFaceIDs()) File "/usr/local/lib64/python2.5/site-packages/fipy/tools/numerix.py", line 1231, in take taken = taken.filled() AttributeError: 'numpy.ndarray' object has no attribute 'filled' So, in my case, some examples are still not working, as I said. Ok, now let's add two more questions: 1. There is possible to use/specify periodic boundary conditions for simulations? 2. There is possible to use for Source Term also an integral function? Cheers, Ioan
