Many thanks for the patch. The specified commit works beautifully.

Thanks again,
Dave

On Nov 16, 2012, at 12:08 PM, Daniel Wheeler wrote:

> Sorry, that got mangled. See attached file. You only need the changes
> in `fipy/meshes/periodicGrid2D.py`.
> 
> Cheers
> 
> On Fri, Nov 16, 2012 at 11:59 AM, Daniel Wheeler
> <[email protected]> wrote:
>> On Thu, Nov 15, 2012 at 1:17 PM, dely <[email protected]> wrote:
>>> FiPy,
>>> 
>>> I am unable to read pickled dump files from simulations with a 
>>> PeriodicGrid2DLeftRight mesh using FiPy version '3.0-dev5303'. I use the 
>>> fipy.dump function to save the mesh and variables for postprocessing, but I 
>>> receive the following output when I try to unpickle the mesh using 
>>> fipy.dump.read:
>> 
>> Hi Dave,
>> 
>> This is defintely a bug. Thanks for finding it. I filed a ticket for
>> this (see http://matforge.org/fipy/ticket/482) and fixed it. If you
>> want these changes, you can clone the FiPy repository from
>> 
>>  git://code.matforge.org/nist/fipy.git
>> 
>> The commit you'll need is 6e475ec2bf01aae44bc. You could also just
>> patch the class that is broken for you with the following simple
>> changes.
>> 
>>  diff --git a/fipy/meshes/periodicGrid2D.py b/fipy/meshes/periodicGrid2D.py
>>  index bcdbea0..95f8945 100755
>>  --- a/fipy/meshes/periodicGrid2D.py
>>  +++ b/fipy/meshes/periodicGrid2D.py
>>  @@ -44,8 +44,8 @@ from fipy.meshes.grid2D import Grid2D
>>   __all__ = ["PeriodicGrid2D", "PeriodicGrid2DLeftRight",
>> "PeriodicGrid2DTopBottom"]
>> 
>>   class _BasePeriodicGrid2D(Grid2D):
>>  -    def __init__(self, dx = 1., dy = 1., nx = None, ny = None,
>> overlap=2, communicator=parallel):
>>  -        super(_BasePeriodicGrid2D, self).__init__(dx = dx, dy = dy,
>> nx = nx, ny = ny, overlap=overlap, communicator=communicator)
>>  +    def __init__(self, dx = 1., dy = 1., nx = None, ny = None,
>> overlap=2, communicator=parallel, *args, **kwargs):
>>  +        super(_BasePeriodicGrid2D, self).__init__(dx = dx, dy = dy,
>> nx = nx, ny = ny, overlap=overlap, communicator=communicator, *args,
>> **kwargs)
>>           self._nonPeriodicCellVertexIDs = super(_BasePeriodicGrid2D,
>> self)._cellVertexIDs
>>           self._orderedCellVertexIDs_data =
>> super(_BasePeriodicGrid2D, self)._orderedCellVertexIDs
>>           self._nonPeriodicCellFaceIDs =
>> numerix.array(super(_BasePeriodicGrid2D, self).cellFaceIDs)
>> 
>> Good luck!
>> 
>> --
>> Daniel Wheeler
> 
> 
> 
> --
> Daniel Wheeler
> <dave.patch>_______________________________________________
> fipy mailing list
> [email protected]
> http://www.ctcms.nist.gov/fipy
>  [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]


_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
  [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]

Reply via email to