Hello all, I have a set of irregularly gridded data (x,y,z) and I am trying to create an interpolating surface using Thin Plate Splines. I couldn't find any existing Julia routines so I thought I'd just do it my self. Here is my implementation <http://nbviewer.ipython.org/urls/gist.githubusercontent.com/lstagner/89964e63557ddb831e72/raw/49049103742450d80ee28d8876fa1a8b8037f970/tps.ipynb>. As you can see its wrong. I been staring at it for a while now and I am beginning to think I must be hitting some sort of bug or quirk of the language. It either that or I did something wrong. If I get this to work I was thinking about incorporating it into one of the existing interpolation packages.
Can anyone figure out why this is not working? Sources: http://www.geometrictools.com/Documentation/ThinPlateSplines.pdf http://user.engineering.uiowa.edu/~aip/papers/bookstein-89.pdf
