Community
cancel
Showing results for 
Search instead for 
Did you mean: 
Shigeo_Kobayashi
Beginner
81 Views

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

0 Kudos
5 Replies
Andrey_N_Intel
Employee
81 Views

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

 

Gennady_F_Intel
Moderator
81 Views

Shigeo,

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

--Gennady



Shigeo_Kobayashi
Beginner
81 Views

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. 

Gennady_F_Intel
Moderator
81 Views

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

Shigeo_Kobayashi
Beginner
81 Views

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

Reply