gmane.comp.python.fipy

Re: Solver stagnated error

franck kalala — 2008-08-28T18:43:02

I thinking about that. I was using one document in which I picked up this example. As you can in the attached file I sent to you there is non boundary conditions for phi and I think may be the problem should be there. I have check every where in the document I have, they have not mentioned initial conditions, sure they forgoted to mentioned them. Daniel according to plot in the attached file, I want to try or I guess these boundary conditions: phi=0.9 at x=0 phi=0.3 at x=1 and other values of phi (the approximations) are in between. we can try these boundaries. if this not work, I will join by the week end the author to know if there are other boundary conditions or if he just forgot it. Regards 2008/8/28 Daniel Wheeler public.gmane.org>

Re: Solver stagnated error

Daniel Wheeler — 2008-08-28T16:47:40

Sorry. What are the boundary conditions for phi? On Thu, Aug 28, 2008 at 12:35 PM, Daniel Wheeler public.gmane.org> wrote:

Re: Solver stagnated error

Daniel Wheeler — 2008-08-28T16:47:29

What should are the boundary conditions for phi? On Tue, Aug 26, 2008 at 1:58 PM, franck kalala public.gmane.org> wrote:

Re: example.diffusion.cricle

Jonathan Guyer — 2008-08-28T11:55:55

On Aug 27, 2008, at 4:53 PM, franck kalala wrote: How strange. Please apply the patch below and report back the error messages (if you're using svn, you can just "revert" later; otherwise, work on a copy). Index: examples/diffusion/circle.py =================================================================== --- examples/diffusion/circle.py(revision 2734) +++ examples/diffusion/circle.py(working copy) < at >< at > -141,14 +141,11 < at >< at > >>> viewer = None >>> if __name__ == '__main__': - ... try: - ... from fipy.viewers.gistViewer.gist2DViewer import Gist2DViewer - ... viewer = Gist2DViewer(vars=phi, - ... limits={'datamin': -1, 'datamax': 1.}) - ... viewer.plotMesh() - ... raw_input("Irregular circular mesh. Press to proceed...") - ... except: - ... print "Unable to create a Gist2DViewer" + ... from fipy.viewers.gistViewer.gist2DViewer import Gist2DViewer + ... viewe

Re: example.diffusion.cricle

franck kalala — 2008-08-27T20:53:17

yes Gist is installed, when running this program there a windows intitled c:\windows\system32\cmd.exe which appear and disapear quickly and the follow the message "Unable to create a Gist2DViewer" 2008/8/27 Jonathan Guyer public.gmane.org>

Re: Porting to FiPy-2.0a1. Troubles.

Jonathan Guyer — 2008-08-27T19:57:54

On Aug 27, 2008, at 3:29 PM, I wrote: I forgot a crucial file. Get r2734 or later.

Re: Porting to FiPy-2.0a1. Troubles.

Jonathan Guyer — 2008-08-27T19:29:50

On Aug 25, 2008, at 10:40 AM, I wrote: I've just checked in a fix for this problem. If you update to r2731 and reinstall FiPy, it will hopefully be resolved.

Re: example.diffusion.cricle

Jonathan Guyer — 2008-08-27T14:36:39

On Aug 27, 2008, at 10:08 AM, franck kalala wrote: Have you installed Gist?

Re: example.diffusion.cricle

franck kalala — 2008-08-27T14:08:28

Thanks Jonathan, I have not this error message, but i got this "Unable to create a Gist2DViewer" Regards 2008/8/26 Jonathan Guyer public.gmane.org>

Re: example.diffusion.cricle

Jonathan Guyer — 2008-08-26T20:49:35

On Aug 26, 2008, at 8:34 AM, franck kalala wrote: This error is because FiPy was written for an older version of Gmsh than you have installed. We fixed this problem on trunk a long time ago, but it never got ported back to 1.2. Thanks for reporting the problem. ---- If you are using svn (I don't think you are?) then you can do an svn update You will need to be sure that you have "branches/version-1_2" and not "tags/version-1_2". You can do svn info to check, or you can just do svn switch http://matforge.org/svn/fipy/branches/version-1_2 to be sure. ---- If you are *not* using svn, then you can manually download http://matforge.org/svn/fipy/branches/version-1_2/fipy/meshes/numMesh/gmshImport.py and replace the file of the same name in your fipy/ directory. If your PYTHONPATH points at this working copy of the fipy source, then you're done. If you previously installed fipy, then you'll need to run python setup.py install again.

Re: Solver stagnated error

Daniel Wheeler — 2008-08-26T19:24:32

On Tue, Aug 26, 2008 at 1:58 PM, franck kalala public.gmane.org> wrote: Frank, This is most probably caused by an issue in fipy whereby the variable operators (getGrad, GetFaceGrad etc) are unaware of the boundary conditions assigned to the variable. This is something we are planning to fix in later versions of fipy. There are ways to fix this issue and I'll get back to you in a day or two to help with this. Cheers

Re: Solver stagnated error

franck kalala — 2008-08-26T17:58:57

Re: Surface plot using viewers

franck kalala — 2008-08-26T17:03:17

Thanks it works 2008/8/26 Jonathan Guyer public.gmane.org>

Re: Surface plot using viewers

Jonathan Guyer — 2008-08-26T16:49:25

On Aug 26, 2008, at 11:52 AM, franck kalala wrote: Daniel's code was for trunk/. You'll want to do x = m.getCellCenters()[...,0] y = m.getCellCenters()[...,1] like you've been doing all along.

Re: Surface plot using viewers

franck kalala — 2008-08-26T15:52:54

Thanks for ur code, I try to run it in a shell, but i have this error message, Traceback (most recent call last): File "C:/Python25/surface.py", line 10, in v = CellVariable(mesh=m, value=numerix.exp(-(x**2+y**2))) File "C:\Python25\lib\site-packages\fipy\variables\cellVariable.py", line 74, in __init__ Variable.__init__(self, mesh = mesh, name = name, value = value, unit = unit, array = array) File "C:\Python25\lib\site-packages\fipy\variables\variable.py", line 130, in __init__ self._setValue(value = value, unit = unit, array = array) File "C:\Python25\lib\site-packages\fipy\variables\variable.py", line 512, in _setValue self.value = self._makeValue(value = value, unit = unit, array = array) File "C:\Python25\lib\site-packages\fipy\variables\variable.py", line 542, in _makeValue array[:] = value ValueError: shape mismatch: objects cannot be broadcast to a single shape Regards 2008/8/26 Daniel Wheeler public.gmane.org>

Re: Surface plot using viewers

Daniel Wheeler — 2008-08-26T13:59:18

from fipy import * L = 1. N = 100 m = Grid2D(nx=N, ny=N, dx=L / N, dy=L / N) x, y = m.getCellCenters() v = CellVariable(mesh=m, value=numerix.exp(-(x**2+y**2))) vi = viewers.make(v) vi.plot() raw_input("wait") Hope that helps. Cheers. On Mon, Aug 25, 2008 at 8:25 PM, franck kalala public.gmane.org> wrote:

Re: Surface plot using viewers

Jonathan Guyer — 2008-08-26T13:44:58

On Aug 25, 2008, at 8:25 PM, franck kalala wrote: Create a CellVariable with that value and view it as usual. Alternatively, make CellVariable\s of x and y and then z will automatically be a CellVariable.

Re: example.diffusion.cricle

franck kalala — 2008-08-26T12:34:43

Hello, Thanks now it running, but another error Traceback (most recent call last): File "C:/Python25/hil.py", line 43, in mesh = GmshImporter2D(meshName) File "C:\Python25\lib\site-packages\fipy\meshes\numMesh\gmshImport.py", line 394, in __init__ mesh2D.Mesh2D.__init__(self, **_DataGetter().getData(filename, dimensions = 2, coordDimensions = coordDimensions)) File "C:\Python25\lib\site-packages\fipy\meshes\numMesh\gmshImport.py", line 263, in getData vertexCoords = self._calcVertexCoords(coordDimensions) File "C:\Python25\lib\site-packages\fipy\meshes\numMesh\gmshImport.py", line 286, in _calcVertexCoords numVertices = int(self.inFile.readline()) ValueError: invalid literal for int() with base 10: '2 0 8\n' Regards 2008/8/26 Daniel Stiles public.gmane.org>

Re: example.diffusion.cricle

Daniel Stiles — 2008-08-26T12:07:13

Hello, On this line: 'Plane Surface(11)={10};\n', delete the ending comma, and put in a ] instead, so it looks like: 'Plane Surface(11)={10};\n'] Daniel S On Tue, Aug 26, 2008 at 7:39 AM, franck kalala public.gmane.org>wrote:

example.diffusion.cricle

franck kalala — 2008-08-26T11:39:41

Hi all I am so sorry to disturb you so much, I am very young in python, in fipy, sure I lent many thing from you since I decided to use fipy, I was this example ************************************************************************************************************** cellSize=0.05 radius=1. lines= ['cellSize = '+str(cellSize) + ';\in', 'radius ='+str(radius) + ';\in', 'point(1)={0, 0, 0, cellSize};\n', 'point(2)={-radius, 0, 0, cellSize};\n', 'point(3)={0, radius, 0, cellSize};\n', 'point(4)={radius, 0, 0, cellSize};\n', 'point(5)={0, -radius, 0, cellSize};\n', 'point(6)={2, 1, 3};\n', 'point(7)={3, 1, 4};\n', 'point(8)={4, 1, 5};\n', 'circle(9)={5, 1, 2};\n', 'Line Loop(10)={6, 7, 8, 9};\n', 'Plane Surface(11)={10};\n', import tempfile (f,geomName)=tempfile.mkstemp('.geo') file=open(geoName, 'w') file.writelines(lines) file.close() import os os.close(f) import sys if sys.platform=='win32': meshName=

Surface plot using viewers

franck kalala — 2008-08-26T00:25:52

Hi all How do I plot this surface z=exp(-(x**2+y**2)) using viewers within fipy, for a given 2D mesh. Regards