c++ - 为什么此代码不返回空矩阵？

Question

为什么下面粘贴的代码不返回空矩阵，而是返回沿主对角线具有 NaN 的矩阵？ 我哪里做错了？

if(g = 0)
{
for(i=0;i<N;i++)
    {
    for(j=N;j<2*N;j++)
        {
        fnd[i][j]=0;
        }
    }
}

else
{
for(i=0;i<N;i++)
    {
    for(p=0;p<N;p++)
        {
        if(p!=i)
            {
            fnd[i][i+N]+=g*g*3.*fabs(x[i]-x[p])*fabs(x[i+N]-x[p+N]) /fabs(pow(fabs(pow(fabs(x[i]-x[p]),2.)+g*g*pow(fabs(x[i+N]-x[p+N]),2.)+h*h*pow(fabs(x[i+2*N]-x[p+2*N]),2)),2.5) );
            }
        }
    for(j=N;j<2*N;j++)
        {
              //offdiagonal terms<------
        if(i!=j)
            {
            fnd[i][j]=-3.*g*g*fabs(x[i]-x[j-N])* fabs(x[i+N]-x[j])/fabs(pow(fabs(pow(x[i]-x[j-N],2.)+g*g*pow(x[i+N]-x[j],2.)+h*h*pow(x[i+2*N]-x[j+N],2)),2.5) )  ;
            }
          }
    }
}


//-------------------------------------------------------

Answer 1

除了可能的其他内容(这不是完整的代码)，if(g = 0) 应该使用 == 而不是普通的 =.正如所写，g 将始终被分配为零，然后它的计算结果将始终为 false，因此将矩阵初始化为全零的循环不会实际运行。相反，整个 else block 将在 g 为零的情况下运行，我认为这是意外的。