伪代码的高斯消除

Question

我尝试将高斯消除的伪代码 - http://en.wikipedia.org/wiki/Gaussian_elimination#Pseudocode 转换为 C 中的实际实现。我已经做到了，我没有收到任何错误，但我得到的输出与输入相同。有谁知道，我哪里有错误？作为输入，我有 3x4 矩阵 A，其中 num 表示行数（在本例中为 3），num+1 表示列数。

void GE(){
    int i=1;
    int j=1;
    int k,u;

    while(i<=num && j<=num+1){
      int maxi=i;
      for(k=i+1;k<num;k++){
        if(abs(A[k][j])>abs(A[maxi][j])){
            maxi=k;
        }
      }
      if(A[maxi][j]!=0){
        float p;
        p=A[maxi][j];
        A[maxi][j]=A[i][j];
        A[i][j]=p;
        for(u=i+1;u<num;u++){
            A[u][j]=A[u][j]-A[i][j]*A[u][j];
        }
      }//endif
      j++;
      i++;

    }//endwhile

编辑：

FILE *inputMatrix=fopen("file.in","r");
 fscanf(inputMatrix,"%d",&num);
 for(i=0;i<num;++i){
    for(j=0;j<num+1;++j){
        fscanf(inputMatrix,"%f",&A[i][j]);
    }
 }
 fclose(inputMatrix);

file.in:

3
2 1 -1 8
-3 -1 2 -11
-2 1 2 -3

Answer 1

也许你忘了转换 tn C 这一行： “将第 i 行中的每个条目除以 A[i,j]” 类似的东西：

 for (k = j; k <= num+1; ++k) {
        A[i][k] /= A[i][j];
    }

Answer 2

看起来您正在将单元格 A[maxi][j] 与 A[i][j] 交换，但算法要求将整行 maxi 与 i 交换：

swap rows i and maxi, but do not change the value of i
Now A[i,j] will contain the old value of A[maxi,j].

所以你需要一个额外的 for 循环，从 [0, num+1) 开始进行交换：

int x
for(x = 0; x < num+1; x++) {
    p=A[maxi][x];
    A[maxi][x]=A[i][x];
    A[i][x]=p;
}