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

学习，那两个涵数什么意思，这是写啥的。

应该还有其他代码吧，程序运行不动的原因应该不在这块，把其他的也帖出来瞧瞧。
或者你传进去(n_mod太大了？

坐等结果

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

这就是开一个二维数组。

你也不去算算1500*1500的Double数据空间有多大

1500*1500*8=18000000≈17.166M

new的是啥？指针存放在哪？

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