First possible successful mpi sampling code.
Moved Lua code out of sampling.cpp, the call is made from mara.cpp now. sampling.lua test performs correctly on 10,000 points x 16 mpi procs on laptop.
Modififed subgrids.cpp SubgridAtPosition not to iterate over non-existent dimensions.
Generalized sampling to work in 1,2, or 3 dimensions. Mildly tested in 3, nothing else yet.
Fixed a bug in mara.cpp init_prim, index offset incorrect. Pausing to pull the bug fix into master.
Have bilinear interpolation with sampling working in 2d, not tested in 3d or with MPI. Fixed a potential bug where MPI non-blocking sends would have referenced out-of-date data. Exposed bounary.ApplyBoundaries to Lua.
Moved bilinear interp into a function, trilinear interp function is a shell, thinking about how to write a test for it.
Filled out 3d trilinear interp function, still no test written for it.
Brought in bug fix to lua_vis.cpp from master.
Needed cstring in sampling.cpp.
Made live colormap update work, pressing 'c' reloads the present texture.
Merge branch 'master' into lua_glfw
Merged in updates from master branch.
Sampling test is setup for runs on a cluster.
Fixed one bug in the 3d sampling call: (Nd == 1) should have been (Nd == 3).
Added beginning of lines rendering.
Sampling points one-at-a-time from Lua is quite slow. Wrote test for speed of equivalent C function.
Forgot to initialize prim array in test/sampling.lua on last commit.
Caught a bug in the 3d sampling, trilinear_interp has a y index that should have been z. Also added the bulk-sampling routine and a test which (works) in serial, untested in parallel. By (works), I mean there is still a bug in the sampling or interpolation which crops up near boundaries.
Actually, that interpolation thing was not a bug. It happens because outflow boundaries violate the linear dependence assumed by the test.