0

前方差分法をコーディングして熱方程式を解こうとしていますが、コードは 13 回目の反復中に常にセグメンテーション フォールトを取得します。そこで、いくつかの printf をループに挿入し、いくつかの変数の値をチェックして、デバッグに行きました。printf に固執すると、ループは以前に障害が発生したポイントを超えて正常に続行されます。ただし、これらのprintfをコメントアウトすると、コードは再びフォールトします...なぜこれが考えられるのか誰にもわかりますか?

/////////////////////////////////////////////////////////////
//
//          Math 448 Assignment 5 Gauss-Seidel
//
/////////////////////////////////////////////////////////////
#include <stdio.h>
#include <iostream>
#include <math.h>

#define I 9
#define J 100
#define M 10


double w[M][J] = {0};
double v[M] = {0};
double K = 0.005;
double H = 0.1;

// J is time, t
// I is distance, x

int main(void){

    // First build up t=0
    int i,j;
    double lambda;
    double t = 0;

        printf("k=%lf\th=%lf\n", K, H);
    lambda = (0.5)*(double)K/(H*H);

    for(i=0;i<=M-1;i++){
        if(i<=5)
        {
            w[i][0] = i*H;
        }else{
            w[i][0] = (double)(1-i*H);
        }
    }

    printf("lam=%lf\n",lambda);

    for(j=1;j<=J;j++){
        v[0] = 0;
        v[I-1] = 0;

        for(i=1;i<=M-1;i++){
            if(t>=0.06){
                            // THESE ARE THE PRINTS THAT CAUSE IT TO WORK
                //printf("i=%d\n1) %lf\n2) %lf\n3) %lf\n",i, w[i-1][j-1], w[i][j-1], w[i+1][j-1]);
                //getchar();
            }
            v[i] = lambda*w[i-1][j-1] + (1-2*lambda)*w[i][j-1] + lambda*w[i+1][j-1];

        }

        t = j*K;
        printf("t=%.4lf\t",t);
        for(i=0;i<M;i++){
            printf("%.4lf\t", v[i]);
        }

        for(i=1;i<M-1;i++){
            w[i][j] = v[i];
        }

        printf("\n");

    }


    /*
            //Psuedocode
        for j=1 to N
            v0 = 0
            vm = 0
            for i=1 to m-1
                vi = LAM*Wi-1 + (1-2LAM)*Wi + LAM*Wi+1
                end for
            t = i*k
            output(v0, v1, v2...)
        end for
    */
    return 0;
}
4

2 に答える 2

1

ここに：

for(j=1;j<=J;j++)

jから0まで反復していますJが、配列には から0までインデックスを付けることができますJ-1。

また、Lol4t0 のコメントに基づいて、ここに:

w[i+1][j-1];

のためにi = M。

于 2013-04-30T19:24:54.313 に答える
0

の配列境界外のスロットにアクセスしようとしていますj == J。インデックスは0...J-1範囲内にあります。

于 2013-04-30T19:24:52.280 に答える