dfdInterpolate1D (Data Fitting API) crashes !!

dfdInterpolate1D (Data Fitting API) crashes !!

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);
    }
}

6 posts / 0 new
Last post
For more complete information about compiler optimizations, see our Optimization Notice.

Hello Shigeo,

We confirm this is the bug in Intel(R) MKL11.0.2 which is caused by code optimizations of the interpolation function for a scalar case. There is a recent discussion of the same issue here: http://software.intel.com/en-us/forums/topic/369956. We expect fix for the bug in the next version of the library.

Can you please clarify if this bug blocks development of your application?

Thanks, Andrey

 

Shigeo,

This an known issue. Please see the disscussion here:  http://software.intel.com/en-us/forums/topic/369956

--Gennady

Hi Gennady,
Thank you for your post. Yes yes, it must be the same(which I could not find through topic titles...)!
If fixed version come soon,then I can wait for it or I will revert to update 1. 

yes, we are planning to add the fix into the next update. We will let you know then the fix would be available.

Hi Andrey and Gennady,
I can wait for new version of course!  What I want to know is 'Am I wrong ?' and you already answered me, thanks!!

Leave a Comment

Please sign in to add a comment. Not a member? Join today