This is the code which accompanies the paper "Divergence-free interpolation of vector fields from point values - exact div B=0 in numerical simulations". The code has been designed to be simple and minimal, so that the reader can appreciate that the radial basis function based method is conceptually no more compicated than Taylor-series based finite-difference approaches. A consequence of these design choices is that the demonstration code here is utterly unsuited to solving any real astrophysical problems.
07/06/2011 - Lauri Perkkio has brough to my attention a typo in equation 7, where the expression for phi_12 is missing a factor of xy. You can see the correct expression in line 87 of the file rbfvmod.py in the original sample code.
Back to my main website at www.colinmcnally.ca.