splie2（用于计算函数在节点处关于y的二阶偏导数）.cpp

用C++编写的双三次样条插值算法

void splie2(double x1a[], double x2a[], double ya[][11], int m, int n, double y2a[][11])
{
    double ytmp[101], y2tmp[101];
	int j,k;
    for (j = 1; j<=m; j++)
	{
        for (k = 1; k<=n; k++)
            ytmp[k] = ya[j][k];
        spline(x2a, ytmp, n, 1e+30, 1e+30, y2tmp);
        for (k = 1; k<=n; k++)
            y2a[j][k] = y2tmp[k];
    }
}