int main(){
double value;
ae_int_t m = 2;
ae_int_t n = 3;
real_1d_array x = "[0,1,2]";
real_1d_array y = "[0,1]";
real_2d_array f = "[[1,2],[1,2,3]]";
spline2dinterpolant s;
spline2dbuildbicubic(x,y,f,m,n,s);
value = spline2dcalc(s,.5,.5);
printf("value=%f\n",value);
return 0;
}

int main(int argc, char *argv[])
{
    QCoreApplication proga(argc, argv);

alglib::ae_int_t  max_x, max_y;
alglib::spline2dinterpolant c;

    max_x=15;
    max_y=15;


XX.setlength(max_x);
YY.setlength(max_y);
FF.setlength(max_y, max_x);

for (int i=0; i<max_x-1; i++)
{
    XX[i]=1.0*i;
}

for (int j=0; j<max_y-1; j++)
{
    YY[j]=1.0+j;
}

for (int i=0; i<max_x-1; i++)
for (int j=0; j<max_y-1; j++)
{
    FF[j][i]=YY[j]*XX[i];
}
/**********************************************/
spline2dbuildbicubic(XX, YY, FF, max_x, max_y, &c);
/**********************************************/
return proga.exec();
}

void interpn(double *outgrid, double *gridx, double *gridy, double *gridxy, double *agrid1, double *agrid2, int ngrid, int ngrid2) {
alglib::real_1d_array gridxx;
gridxx.setcontent(ngrid,gridx);

alglib::real_1d_array gridyy;
gridyy.setcontent(ngrid2,gridy);

alglib::real_1d_array ingrid;
ingrid.setcontent(ngrid*ngrid2,gridxy);

alglib::spline1dinterpolant splone;

// build spline
alglib::spline2dbuildbicubicv(gridxx,ngrid,gridyy,ngrid2,ingrid,1,splone);

for (int ii=0; ii<ngrid; ii++) {
    for (int jj=0; jj<ngrid2; jj++) {
        for (int kk=0; kk<nstates_as; kk++){
            for (int ll=0; ll<nstates_is; ll++){
                outgrid[ii*ngrid2+jj+(kk+ll)*(ngrid*ngrid2)] = alglib::spline2dcalc(splone, agrid1[ii*ngrid2+jj], agrid2[ii*ngrid2+jj]);
           }
        }
    }
}

}
.
.
.
interpn(&k2prime_bu[0][0][0][0],kkgrid,kmgrid,tempmat,&kprime[0][0][0][0],&kmprime[0][0][0][0],ngrid,ngrid2);

alglib::spline2dbuildbicubicv(gridxx,ngrid,gridyy,ngrid2,ingrid,1,splone);