im2=allocate_f_2d(n_mod,n_mod);
double **allocate_f_2d(int r,int c)
/*给二维浮点数组申请储存空间*/
{
    double **a=new double*[r];
    for(int i=0;i<r;i++)
        a[i]=new double[c];
    return a;
}

是的，和你说的一样。后来发现是其他调用函数时，实参形参没有对应，导致实际开内存比这还要大。

这就是开一个二维数组。

new的是double类型，存放在a[]；还有double类型指针，放在指针a里面。