Why this code crashes ?
I confirmed the crash under following environment:
Microsoft Windows 7 (64bits)
Intel(R) Core(TM) i7-2600S CPU 2.80GHz、2801 Mhz 8Gb memories.
Intel(R) Math Kernel Library Version 11.0.2 Product Build 20130123 for 32-bit applications.
I saw no crash on 'Intel(R) Math Kernel Library Version 11.0.1 Product Build 20121016 for 32-bit applications' (but not exactly the same code bellow).
void TestSpline()
{
int NY = 1;
double x[] = {1.0,2.0};
double y[] = {0.0,2.0};
int NX = sizeof(x)/sizeof(x[0]);
int order = 2;
double boundaryVals[3];
int ndorder = 1;
int dorder[] = { 1 };
int type = 0;
int bc_type = 0;
double xv = 0.0;
double yv = 0.0;
double *coef = (double*)mkl_malloc(8*10240,16);
DFTaskPtr task; int status;
boundaryVals[0] = 0; boundaryVals[1] = 0; boundaryVals[2] = 3.333;
status = dfdNewTask1D( &task,NX,x, DF_NO_HINT, NY, y, DF_MATRIX_STORAGE_ROWS );
status = dfdEditPPSpline1D(task, order, type, bc_type , boundaryVals, DF_NO_IC, NULL, coef, DF_NO_HINT);
status = dfdConstruct1D( task, DF_PP_SPLINE, DF_METHOD_STD );
for(int i=0;i<=6;++i) {
xv = i; // Crashes when xv==x[max indeces] !!!!!
status = dfdInterpolate1D( task, DF_INTERP, DF_METHOD_PP, 1, &xv, DF_NON_UNIFORM_PARTITION, ndorder,dorder, NULL,
&yv,DF_MATRIX_STORAGE_COLS, NULL );
printf("status = %d x=%g y=%g\n",status,xv,yv);
}
}
