一个双精度数组未收到双精度值

Question

我有这段代码（我是从FFT开发出来的），根据Visual Studio，它看起来还不错...但是，当到达该点时，block [a] [b]必须接收该值（即使它不是pow（），而是expo_final，例如，block [a] [b]的值未按照我的期望进行。

来源在这里。

void twiddle_factor(float *d_isub_matrix)
{


    double block[THREAD_SIZE][THREAD_SIZE];
    float *d_osub_matrix = new float[THREAD_SIZE*THREAD_SIZE];
    int a,b,c,d,x,y,z;


    float sum_sines=0.0;
    double expo1,expo2,expo_final;
    float sum_cosines=0.0;
    float sum_sin[THREAD_SIZE][THREAD_SIZE],sum_cos[THREAD_SIZE][THREAD_SIZE];
    float angle=(2*PI)/THREAD_SIZE;

    //put into shared memory the FFT calculation (F(u))

    for(x=0;x<THREAD_SIZE;x++)
    {
        for(y=0;y<THREAD_SIZE;y++)
        {

            for(z=0;z<THREAD_SIZE;z++)
            {
                sum_sines=sum_sines+sin(d_isub_matrix[y*THREAD_SIZE+z]*(angle*(2*z)))+sin(d_isub_matrix[y*THREAD_SIZE+z]*(angle*(2*z+1)));
                sum_cosines=sum_cosines+cos(d_isub_matrix[y*THREAD_SIZE+z]*(angle*(2*z)))+cos(d_isub_matrix[y*THREAD_SIZE+z]*(angle*(2*z+1)));

            }
            sum_sin[x][y]=sum_sines/(2*THREAD_SIZE);
            sum_cos[x][y]=sum_cosines/(2*THREAD_SIZE);

        }
    }

    //getting the value!

    for(a=0;a<THREAD_SIZE;a++)
    {
        for(b=0;b<THREAD_SIZE;a++)
        {
            expo1=exp(sum_sin[a][b]);
            expo2=exp(sum_cos[a][b]);
            expo_final=expo1+expo2;
            block[a][b]=expo_final;
        }

    }

其中定义了ARRAY_SIZE和THREAD_SIZE，它们的值分别为64（ARRAY）和32（THREAD）。 我验证了sum_sines，sum_cosines，sum_sin [] []，sum_cos [] []，angle，expo1，expo2和expo_final的值。 直到最后一个，一切都很好。 它们返回一些负值和正值。 但是我无法获取block [] []的值。 有人知道发生了什么吗？

Answer 1

这是您的问题：

for(b=0;b<THREAD_SIZE;a++)