c - 垃圾存储，C 编程

Question

我是新的 C 编程人员。我有一个 781 行的主代码，它失去了控制，因为垃圾值存储在 vector 中。下面显示了主代码的一小部分，其中调用了名为 diff_conv_intermedia1 的子例程。

diff_conv_intermedia1(&factorteta,&N,ID,DIFF,X1_intermedia,Y1_intermedia,X1C_intermedia,Y1C_intermedia,CU1_intermedia,CV1_intermedia,AW1_intermedia,AE1_intermedia,AS1_intermedia,AN1_intermedia,AP1_intermedia,Q1_intermedia,FXI1,FYI1,FI_intermedia1,1,2,1,1);

int q,w;
for(q=1;q<(*factorteta_Ptr)*2+1;q++)
 {
   for(w=1;w<(*N_Ptr)+1;w++)
 {
   printf("%lf\n",AP1_intermedia[q][w]);
 }
 }   

所以子程序如下所示。当我在子例程内部打印结果时，一切正常，但是当我在子例程外部打印结果时，在主代码中，垃圾将作为 AP1_intermedia 存储在 vector 中。我不知道出了什么问题。我对其他子例程重复相同的过程，没有任何错误。

int diff_conv_intermedia1(int *factorteta_Ptr,
                          int *N_Ptr,
                          int ID,
                          double DIFF,
                          double X[(*factorteta_Ptr)*2+1][*N_Ptr+1],
                          double Y[(*factorteta_Ptr)*2+1][*N_Ptr+1],
                          double XC[(*factorteta_Ptr)*2+2][*N_Ptr+2],
                          double YC[(*factorteta_Ptr)*2+2][*N_Ptr+2],
                          double CU[(*factorteta_Ptr)*2+1][*N_Ptr+1],
                          double CV[(*factorteta_Ptr)*2+1][*N_Ptr+1],
                          double AW[(*factorteta_Ptr)*2+1][*N_Ptr+1],
                          double AE[(*factorteta_Ptr)*2+1][*N_Ptr+1],
                          double AS[(*factorteta_Ptr)*2+1][*N_Ptr+1],
                          double AN[(*factorteta_Ptr)*2+1][*N_Ptr+1],
                          double AP[(*factorteta_Ptr)*2+1][*N_Ptr+1],
                          double Q[(*factorteta_Ptr)*2+1][*N_Ptr+1],
                          double FX[(*factorteta_Ptr)*2+1][*N_Ptr+1],
                          double FY[(*factorteta_Ptr)*2+1][*N_Ptr+1],
                          double FI[(*factorteta_Ptr)*2+1][*N_Ptr+1],
                          int WBC,int EBC,int SBC,int NBC)
{
int i,j;
double value,* valuePtr;
double AED, AWD, AND, ASD;
double AEC, AWC, ANC, ASC;
valuePtr = &value;

// Diffusive coefficients

for(i=1;i<(*factorteta_Ptr)*2+1;i++)
{
  for(j=1;j<*N_Ptr+1;j++)
{
  AWD = -DIFF*(Y[i][j-1]-Y[i-1][j-1])/(XC[i][j]-XC[i][j-1]);
  AED = -DIFF*(Y[i][j]-Y[i-1][j])/(XC[i][j+1]-XC[i][j]);
  AND = -DIFF*(X[i][j]-X[i][j-1])/(YC[i+1][j]-YC[i][j]);
  ASD = -DIFF*(X[i-1][j]-X[i-1][j-1])/(YC[i][j]-YC[i-1][j]);

  // Convection term
  if(ID==2)
    {
      max1_or_min2(CU[i][j-1],1,&value);
      AWC=-*valuePtr;
      max1_or_min2(CU[i][j],2,&value);
      AEC=*valuePtr;
      max1_or_min2(CV[i-1][j],1,&value);
      ASC=-*valuePtr;
      max1_or_min2(CV[i][j],2,&value);
      ANC=*valuePtr;
    }
  if(ID==1)
    {
      AWC =-CU[i][j-1]*(1.0-FX[i][j-1]);
      AEC =CU[i][j]*FX[i][j];
      ASC =-CV[i-1][j]*(1.0-FY[i-1][j]);
      ANC =CV[i][j]*FY[i][j];
    }

  // Set Coefficients matrix
  AW[i][j] = AWD+AWC;
  AE[i][j] = AED+AEC;
  AS[i][j] = ASD+ASC;
  AN[i][j] = AND+ANC;
  AP[i][j] = -(AE[i][j]+AW[i][j]+AN[i][j]+AS[i][j]);
  Q[i][j] = 0.0;

}
}
// West Boundary - Inlet B.C
for(i=1;i<(*factorteta_Ptr)*2+1;i++)
 {
  if(WBC==1) Q[i][1] = Q[i][1]-AW[i][1]*FI[i][0];
  if(WBC==2) AP[i][1] = AP[i][1] + AW[i][1];
  AW[i][1] = 0.0;

// East Boundary - (1)Dirichlet (2)ZERP-GRAD Outflow B.C
  if(EBC==1) Q[i][*N_Ptr] = Q[i][*N_Ptr] - AE[i][*N_Ptr]*FI[i][*N_Ptr+1];
  if(EBC==2) AP[i][*N_Ptr] = AP[i][*N_Ptr] + AE[i][*N_Ptr];
  AE[i][*N_Ptr] = 0.0;  
}

// South Boundary - (1)Dirichlet (2)ZERO-GRAD
for(j=1;j<*N_Ptr+1;j++)
{
  if(SBC==1) Q[1][j] = Q[1][j] - AS[1][j]*FI[0][j];
  if(SBC==2) AP[1][j] = AP[1][j] + AS[1][j];
  AS[1][j] = 0.0;

// North Boundary - (1)Dirichlet (2)ZERO-GRAD
  if(NBC==1) Q[(*factorteta_Ptr)*2][j] = Q[(*factorteta_Ptr)*2][j] - AN[(*factorteta_Ptr)*2][j]*FI[(*factorteta_Ptr)*2+1][j];
  if(NBC==2) AP[(*factorteta_Ptr)*2][j] = AP[(*factorteta_Ptr)*2][j] + AN[(*factorteta_Ptr)*2][j];
  AN[(*factorteta_Ptr)*2][j] = 0.0;
}

// Print
int l,k;
for(l=1;l<(*factorteta_Ptr)*2+1;l++)
{
  for(k=1;k<*N_Ptr+1;k++)
{
  printf("%lf %lf %lf %lf\n",AP[l][k],AS[l][k],AN[l][k],FI[l][k]);
}
}

return 0;
}

如果有人愿意，我可以发送所有代码，但有很多扩展。

Answer 1

在你的函数声明中:双 AP[(*factorteta_Ptr)*2+1][*N_Ptr+1]

我不太认为这正在做你认为它正在做的事情。虽然我自己以前没有见过类似的东西，但我相信这是告诉编译器根据其他给定参数为您创建一个可变长度的二维数组。然后，您在函数中填写这些值。但是，因为您不返回该值，也不将其声明为按引用传递，所以当您返回时它会被丢弃，因此工作会丢失，并且 main() 中的数组中有垃圾>。更好的形式是在 main() 中创建此数组，然后通过引用(如 double *AP[][])将其传入，或者在退出时返回此数组，或者修改比这个函数更糟糕的东西，然后将其设为全局函数，以便您可以在任何地方看到它。