Hi Stefan, On 13/10/06, Stefan van der Walt <stefan@sun.ac.za> wrote:
I'd also like to see some tests to verify the behaviour of the new interp2d that is *different* from the old version. It's easier to add it now, while everything is still fresh in our memories.
Will try and do this tomorrow.
Attached is a very simple script (will make it a test later) which uses the test data from netlib->dierckx to check the use fo regrid and surfit. The plot attached shows the results. There could be an error in my dealing with meshgrid (I hate they way it swaps axes), but I think it is right. The surfit part issues a warning on interpolation. If you increase s to say 20 (0 is needed for interpolation) then the warning goes, but then you are smoothing (as can be seen from the corresponding output plot - which was down sampled for file size). I'll make this into a test at some point and also improve the interp2d interface to be more flexible to input array layout - but this will have to be nest week due to time constraints (I need to work...) Hope this demonstrates my point (and that yopu don't find an obvious error in my code...) John
