On Mon, Oct 24, 2016 at 12:20 AM, Zhekai Deng
<zhekaideng2...@u.northwestern.edu> wrote:
>
> This is not a issue if mesh system is 1D, and boundary cell is just 1. Or
> maybe on the rectangular mesh system. However, does any one has any idea to
> how to transfer it properly ?
Hi Zhekai; I haven't read through everything but I would guess that
you need to interpolate from one mesh to another. To do that you can
use one of Scipy's inerpolation functions such as
https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.interpolate.griddata.html
We also built that capability into FiPy, but it's probably better to
use Scipy's interpolation at this point.
https://github.com/usnistgov/fipy/blob/develop/fipy/variables/cellVariable.py#L167
If the mesh points are aligned, it is still a good idea to use the
interpolation approach as it removes the issue of worrying about mesh
ordering.
--
Daniel Wheeler
_______________________________________________
fipy mailing list
fipy@nist.gov
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
