请教各位，如何将程序中计算出的复数数组返回给主函数。
 以下是我的部分代码，请各位高手指出以下我的错误所在。谢谢了。

#include "stdafx.h"
#include <iostream>
#include <string>
#include <complex>
#include <iomanip>
#include <fstream>
#include <conio.h>
#include <math.h>
#include <fftw3.h>
#include <stdio.h>
#define M 128
#define N 128
#define Pi 3.1415926535898 
#define complex_d complex<double>


complex_d BPalgorithm(complex_d G[M][N])
{
complex_d *gg;
complex_d PLm[M][N];
fftw_complex *in, *out;
fftw_plan p;

   
//分配存储空间
gg=new complex_d[M*N];
if(gg==NULL)
{
throw "分配内存失败！";
}
in = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * M);
out = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * M);

// 设置变换计划
p = fftw_plan_dft_1d(M, in, out, FFTW_BACKWARD, FFTW_ESTIMATE);
   
// 设置测试数据
  for(int j=0;j<N;j++)
{
for(int i=0;i<M;i++)
{
  in[i][0]=real(G[i][j]);
  in[i][1]=imag(G[i][j]);
}

fftw_execute(p);  

for(int l=0;l<M;l++)
{
PLm[l][j]=complex_d(out[l][0],out[l][1]);
}

}

gg=PLm; //出现错误（“=”:无法从“std::complex<double> [128][128]”转换“std::complex<double>*”）
return gg;

// 释放内存
delete []gg;
gg = NULL;
  fftw_destory(p);
  fftw_free(in);
fftw_free(out);
}




int main()
{
complex_d G[M][N];
complex_d *g;

 
//分配存储空间
g=new complex_d[M*N];
if(g==NULL)
{
throw "分配内存失败！";

}
  ......

  g=BPalgorithm(G);//这里出现了错误
  ...,...
  }

汗，你又不给所有的代码，你让别人怎么帮你找问题？用眼睛扫描么？？无语了

以数组名做参数，不需要返回，或开辟为全局数组

呵呵，疏忽了。
下面是我的代码。

#include "stdafx.h"
#include <iostream>
#include <string>
#include <complex>
#include <iomanip>
#include <fstream>
#include <conio.h>
#include <math.h>
#include <fftw3.h>
#include <stdio.h>

#define M 128
#define N 128
#define Pi  3.1415926535898
#define complex_d complex<double>
double data[M][2*N];

using namespace std;

void showresult(complex_d G[M][N])
{
    FILE* data = NULL;
    fopen_s(&data, "imagedata.txt", "w");
    if(data == NULL)
    {
        throw "创建数据文件失败。";
    }

    for(int i=0;i<M;i++)
    {
        for(int j=0;j<N;j++)
        {
            fprintf_s(data, "%f\t",10*log10(abs(G[i][j])));
        }
        fprintf_s(data, "\n");
    }
    fclose(data);
}


void dataread( )
{
    using namespace std;
    int a[2];
    double b[N+1];
    double c[1];
    streampos weizhi;
    streampos WZhi;


    ifstream infile("ISAR2D[1].dat",ios::in);

    if(!infile)
    {
        cerr<<"open error!"<<endl;
        exit(1);
    }
    for(int i=0;i<2;i++)
    {
        infile>>a[i];
    }
    //cout<<endl;
    streampos weizhi1 =infile.tellg();                //weizhi1 是 64 64


    infile.seekg(weizhi1);

    for(int i=0;i<N+1;i++)
    {
        infile>>b[i];
    }
    streampos weizhi2 =infile.tellg();               //weizhi2 第二行结束

    infile.seekg(weizhi2);
    infile>>c[0];
    streampos weizhi3=infile.tellg();               //weizhi3 第三行第二个数开始位置
    weizhi=weizhi3-weizhi2;                        //weizhi 求出第一个数所占的空间



    WZhi= weizhi3;
    for(int i=0;i<M;i++)
    {  
        infile.seekg(WZhi);
        for(int j=0;j<2*N;j++)
        {
            infile>>data[i][j];
        }
        streampos weizhi4=infile.tellg(); 
        infile.seekg(weizhi4);
        infile>>c[0];
        streampos weizhi5=infile.tellg(); 
        WZhi=weizhi5;
    }

}

complex_d BPalgorithm(complex_d G[M][N])
{
    complex_d *gg;
    complex_d PLm[M][N];
    fftw_complex *in, *out;
    fftw_plan p;

   
    //分配存储空间
    gg=new complex_d[M*N];
    if(gg==NULL)
    {
        throw "分配内存失败！";
    }
    in  = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * M);
    out = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * M);
   
    // 设置变换计划
    p = fftw_plan_dft_1d(M, in, out, FFTW_BACKWARD, FFTW_ESTIMATE);
   
    // 设置测试数据
    for(int j=0;j<N;j++)
    {
        for(int i=0;i<M;i++)
        {
          in[i][0]=real(G[i][j]);
          in[i][1]=imag(G[i][j]);
        }

        fftw_execute(p);
       
        for(int l=0;l<M;l++)
        {
            complex_d a(out[l][0],out[l][1]);
            PLm[l][j]=a;
        }
        

    }

    gg=PLm;
    return gg;
   
    // 释放内存
    delete []gg;
    gg = NULL;
    fftw_destroy_plan(p);
    fftw_free(in);
    fftw_free(out);
}




int main()
{
    complex_d G[M][N];
    complex_d a;
    complex_d *g;
    complex_d PL[M];
    double  Win[M][N];

    //分配存储空间
    g=new complex_d[M*N];
    if(g==NULL)
    {
        throw "分配内存失败！";
    }
    //读取数据
    dataread( );

    for (int i=0; i<M; i++)
    {
        for (int j=0; j<N; j++)
        {
            complex_d a(data[i][j],data[i][j]);
            G[i][j]=a;
        }
    }
    //hamming窗
    for (int i=0; i<M; i++)
    {
        for (int j=0; j<N; j++)
        {
            Win[i][j]=(0.5-0.5*(cos(2*Pi*j/(8*M-1))))*(0.5-0.5*(cos(2*Pi*i/(8*N-1))));
        }
    }
    //时域加窗
    for (int i=0; i<M; i++)
    {
        for (int j=0; j<N; j++)
        {
            G[i][j]=G[i][j]*Win[i][j];
        }
    }

    g=BPalgorithm(G);
    

    // 显示测试结果
    showresult(G);

}

    complex_d *gg;
    complex_d PLm[M][N];
一个二维数组 一个单指针 怎么能=呢？
complex_d *g;
g=BPalgorithm(G);返回一个complex_d 你怎么能用complex_d *接收呢

我想的是，把这个二维数组complex_d PLm[M][N]的地址赋给*gg，然后返回这个地址。
这样可行否？