ajdawson/gridfill

Fill missing values in a grid by solving Poisson's equation via an iterative relaxation scheme.

gridfill's Issues

Fill fails for dtype=np.float32

Fortran code is expecting double precision, so the easiest thing to do is probably cast the input grids to np.float64 before passing to the fortran layer.

resmax and numiter are not the maximum for multiple grids

poisson_fill_grids returns the number of iterations numiter and the maximum residual resmax from the last grid operated on. It should return residual and iteration count for all grids.

In xarray

By the way @ajdawson your code is very useful - thanks for sharing. I think it could be a great addition to xarray. They currently have an ffill method but it only operates along one single dimension as far as I can tell. See example below for quick comparison. @shoyer ?

fill not found in new version

Hi, just noticed that this library is updated to a modernized version.

If I install the new version using python setup.py install, I can not find the fill function if I import from gridfill import fill

However this works if I use the old version.(c0ae1f9). Do you know why this happens?