Denis Demidov

Looks like you are using an old (pre C++11) compiler.

I don't have access to Windows/Visual Studio compiler in order to try this, but you should be able to rephrase the problematic line so that your compiler would not complain...

It's a long time since I touched the code, but I think the minimal grid size should be at least 2. The grid will grow as required in order to...

I would try to do `grid = {m0 + 1, n0 + 1}`. The difference could be that they are counting grid cells, and I count grid points.

Each data point only affects 4 grid points around it. So if you have small (5x5x5) initial grid, then each data point contributes to each grid point on this initial...

I am sorry, I don't think I understand the question. Do you want to interpolate vector fields as opposed to scalar (`double`) fields? I think it should be possible to...

There is the `tol` parameter that specifies the target relative approximation error (`||mba(x_i) - z_i|| / ||z||`): https://github.com/ddemidov/mba/blob/420bc7f4d0717a1f66eee05f61122ff2dd28d65a/mba/mba.hpp#L721 The parameter is also exposed to the python wrapper: https://github.com/ddemidov/mba/blob/420bc7f4d0717a1f66eee05f61122ff2dd28d65a/python/pymba.cpp#L30 so you...

I have just tested the `rygg1.dat` example from SINTEF MBA, see results in the following notebook: https://nbviewer.jupyter.org/gist/ddemidov/6fcd134909d674669f6dffc2066bf739 The test point (5,5) has a minor discrepancy, but the data points are...

Note that python interface assumes that all input arrays have C-layout, hence the `.copy()` in cell [2] in the above notebook.

Are large values of u and v still inside your domain (defined with cmin and cmax)? If not then these are outside of interpolation grid and decrease rapidly.