c言語 数値計算 ガウスの消去法 意図した数値が出ない

Question

ガウスの消去法を用いて Ax = y のxを求めようとしています。

A = {{1 -2 3}, {2 1 0}, {1 2 -1}}
y = {1 5 5}です。

答えはx = {1 3 2}となります。
しかし私が書いたコードだと x = {2 3 2} となってしまいます。

どこがおかしいのか分からない状態です。
初学者であるため、コードが読みにくいかもしれません。
お手数おかけしますが、よろしくお願いいたします．

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include<string.h>
#include<unistd.h>

#define N 3 // N×Nの行列

void partial_pivoding(int, double arr[][N], double vec[]); 

int main(){
    double a[][N] = {{1., -2., 3.}, {2., 1., 0.}, {1., 2., -1.}};
    double y[] = {1., 5., 5.}, x[N];
    int k, i, j;
    double c, sum = 0.;

    for(k=0; k<N; k++){
        partial_pivoding(k, a, y);　//部分ピボット選択

        for(i=k+1; i<N; i++){　//前進代入
            c = *(*(a+i)+k) / *(*(a+k)+k);

            for(j=k; j<N; j++) *(*(a+i)+j) -= c * *(*(a+k)+j);
            
            *(y+i) -= c * *(y+k);
        }
    }

    *(x+N-1) = *(y+N-1) / *(*(a+N-1)+N-1);
    
//ここまでの動作が正常なのは確認済

    for(i=N-2; i>=0; i--){　//後退代入　おそらく、ここがおかしい
        
        for(j=i+1; j<N; j++) sum += *(*(a+i)+j) * *(x+j);

        *(x+i) = (*(y+i) - sum) / *(*(a+i)+i);
    }

    printf("result\n");
    for(i=0; i<N; i++) printf("%lf\n", *(x+i));　//x表示
}

void partial_pivoding(int j, double arr[][N], double vec[]){
    int i, k;
    double t_vec[N], t_scal;

    for(i=j+1; i<N; i++){
        if(*(*(arr+j)+j)<*(*(arr+i)+j)){

            for(k=j; k<N; k++){
                *(t_vec+k) = *(*(arr+i)+k);
                *(*(arr+i)+k) = *(*(arr+j)+k);
                *(*(arr+j)+k) = *(t_vec+k);
            }

            t_scal = *(vec+j);
            *(vec+j) = *(vec+i);
            *(vec+i) = t_scal;
        }
    }
    
}

Answer 1

C言語のコードを組むなら、デバッグ環境を揃えましょう
コードの任意の場所で実行を止めて、変数のナカミを見ることができます
また、１行づつ実行させて、動作を確認できます

そうすれば、自分が組んだコードがどう動いてるのかわからないってことはなくなるかと思います

Comments

[drilllaser]

アドバイスありがとうございます。
おすすめの環境などはありますか？
わたしは今VScodeで行っています。

[y_waiwai]

Windowsなら、VisualStudioが最強ですね
あとEclipseなんかも。

Answer 2

すみません。後退代入の中でのsumを初期化してませんでした。