This has been excellent in general. A few times it has crashed at line 57:
t_star=fzero(@(t)fixed_point(t,N,I,a2),[0,.1]);
because "??? Error using ==> fzero at 283
The function values at the interval endpoints must differ in sign."
The data doesn't look obviously bad in these cases. A short example vector is [14.0534 13.2851 13.0951 13.1159 14.2221] (this is a shortened version of a length 100 vector that also crashed).
and later
(line 485)
#ifdef _WIN32
unsigned __stdcall transformvolume_gradient_color(double **Args) {
#else
void transformvolume_error_color(double **Args) {
#endif
Those function names should match on both?
After I change this, it compiles but doesn't link. :-(
>> compile_c_files
compiling : bspline_error_2d_double.c
image_interpolation.o: In function `mindex2':
image_interpolation.c:(.text+0x0): multiple definition of `mindex2'
bspline_error_2d_double.o:bspline_error_2d_double.c:(.text+0x0): first defined here