Question

这个C代码可能有一百个老式或低效的东西。但主要是，我只是想知道为什么我在日常工作中遇到了一个段错误＆＃34; jacobi_solve＆＃34;。驱动程序代码在这里：

#include <stdio.h>
#include <stdlib.h>


int main() {
  int m=20, n=20;
  int i,j;
  double A[m][n],Anew[m][n];

  for (j=0;j<n;j++){
    for (i=0;i<m;i++){
      A[j][i] = 0.; Anew[j][i] = 0.;
    }
  }

  A[10][10] = 1.;
  printf("%x\n",A); 

  jacobi_solver(A,Anew,m,n);

}

和jacobi_solver在这里：

#include<stdio.h>
#include<math.h>

void jacobi_solver(double **A, double **Anew, int m, int n) {

  double err,tol=1.e-6;
  int iter=0,iter_max=100;
  int i,j;

  err = tol*10.;
  printf("hello\n");
  printf("%f %f %d %d\n",err,tol,iter,iter_max);
  printf("%x\n",A);
  printf("%f\n",A[10][10]);   /* <--- seg fault!  */
  printf("solving ...\n");

  while ( err > tol && iter < iter_max ) {

    for( j = 1; j < n-1; j++) {
      for( i = 1; i < m-1; i++) {
        Anew[j][i] = 0.25 * (A[j][i+1] + A[j][i-1] +
                             A[j-1][i] + A[j+1][i]);
        err = fmax(err, abs(Anew[j][i] - A[j][i]));
      }
    }

    printf("%f\n", err);

    for( j = 1; j < n-1; j++) {
      for( i = 1; i < m-1; i++ ) {
        A[j][i] = Anew[j][i];
      }
    }


    iter++;

  }

}

请注意，如果我拿走print语句，稍后会在主求解循环中出现seg错误。

Answer 1

代码不正确 - 你不能像在C中使用2-D数组一样使用double **（与double *和1-D数组不同，有点令人困惑，但如果你理解内存的工作方式则不然）。我建议你用

中建议的方法之一实现自己的二维数组包装器

How do I work with dynamic multi-dimensional arrays in C?

Answer 2

将求解器函数的界面更改为

void jacobi_solver(int m, int n, double A[n][m], double Anew[n][m]) {
  ...
}

应该只使用2D数组。特别是它不是一个指向矢量的指针表，就像你拥有它一样。必须首先使用数组大小，以便您可以将它们用于数组参数。

这应该适用于符合至少C99的所有编译器。

将数组传递给函数 - 出了什么问题？

2 个答案: