topic dfdInterpolate1D (Data Fitting API) crashes !! in Intel® oneAPI Math Kernel Library

dfdInterpolate1D (Data Fitting API) crashes !!

Shigeo_Kobayashi — Wed, 27 Feb 2013 03:58:35 GMT

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

Hello Shigeo,

Andrey_N_Intel — Wed, 27 Feb 2013 05:57:39 GMT

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,

Gennady_F_Intel — Wed, 27 Feb 2013 06:30:15 GMT

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

Shigeo_Kobayashi — Wed, 27 Feb 2013 08:07:10 GMT

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

Gennady_F_Intel — Wed, 27 Feb 2013 15:28:53 GMT

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

Shigeo_Kobayashi — Thu, 28 Feb 2013 10:11:27 GMT

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!!