セグメンテーションエラーが消えません

Question

C言語でファイルから行列の値を読み込むときにセグメンテーションエラーが発生しました。

#include <stdio.h>
#include <math.h>
#define N 4
#define M 1
#define nl 10

void lnqgj(double a[][N+M], int n, int m, double epsl, int *isw)
{
  int j;
  int i;
  int k;
  int ip;

  double p;
  double w;
  double a_kk;
  double a_ik;

  for(k = 0 ; k < n ; k++){
    p = 0.0;

    for(i = k;i < n; i++){
      if(p < fabs(a[i][k])){
        p = fabs(a[i][k]);
        ip = i;
      }
    }

    if(p <= epsl){
      *isw = 1;
      printf("\n\t isw = %d\n", *isw);
      return ;
    }

    for(j = k;j < (n + m) ; j++){
      w = a[k][j];
      a[k][j] = a[ip][j];
      a[ip][j] = w;
    }

    a_kk = a[k][k];
    for(j = k;j < (n + m); j++){
      a[k][j] /= a_kk;
    }  

    for (i = 0; i < n ; i++){
      if(i != k){
        a_ik = a[i][k];
        for(j = k; j<(n + m);j++)
          a[i][j] -= (a_ik * a[k][j]);
      }
    }
  }  

  *isw = 0;
  return ;
  }

main(){
  int a;
  int b;
  int s = N;
  int t = M;
  int *isw;
  double min = 1.0e-5;
  double z[s][s+t];
  FILE  *fpin;
  FILE  *fpout;

  if((fpin = fopen("kadai17.dat", "r"))==0){
     fprintf(stderr, "\n\tSorry! I can't open kadai17.dat\n");
     return 0;
  } 
    if((fpout = fopen("kadai17ans.dat", "w"))==0){
     fprintf(stderr, "\n\tSorry! I can't open kadai17ans.dat\n");
     return 0;
  } 

  for( a = 0; a < 4 ;a++){
    printf("\n%d\n",a);
    for( b = 0; b < 5 ;b++){
      fscanf(fpin,"%lf",&(z[a][b]));
      printf("%lf\t",z[a][b]);
      }
    }

  lnqgj(z,s,t,min,isw);

  fprintf(fpout,"isw = %d\n",*isw);
  for( a = 0; a < 4;a++){
    printf("\n\t%lf\n",z[a][4]);
    fprintf(fpout, "\tX%d =%lf\n",a+1 , z[a][4]);
    }

  fclose (fpin);
  fclose (fpout);
  return 0;
}
  
  

試したこと

読み込んでいるファイルの中身はこちらです
3 2 7 1 8
1 5 1 -1 5
4 1 3 -2 7
1 6 4 3 13

printfで読み込まれたデータを表示するのですが、a<4とすると4行目の値が表示されず上のprintfで表示している行数のみが表示されます

0
3.000000 2.000000 7.000000 1.000000 8.000000
1
1.000000 5.000000 1.000000 -1.000000 5.000000
2
4.000000 1.000000 3.000000 -2.000000 7.000000
3
セグメンテーション違反です

Comments

[Orlofsky]

コードは https://teratail.com/help#about-markdown の[コードを入力]を使ってください。

[LouiS0616]

配列zはどのように宣言していますか？

[tkturbo]

配列ｚが3行しかないのでは？

[tkturbo]

z[3][0]に値がセットされていないのでは。

Answer 1

未初期化ポインタ isw による違法メモリアクセスです。
私の手元では次の箇所で落ちました。

void lnqgj(double a[][N+M], int n, int m, double epsl, int *isw)
{
    （省略）
    *isw = 0;      // ここでSegmentation fault
    return;
}

main()関数が、ポインタ変数 isw を初期化しないままlnqgj()関数に渡した、という単純なバグですから、例えば次のようにすれば動きます。

int main()
{
    int a;
    （途中省略）
    FILE  *fpout;
    int foo;       // 変数を追加し

    isw = &foo;    // そこをポイントするよう isw を初期化しておき
    （途中省略）
    lnqgj(z,s,t,min,isw);  // iswをlnqgj()関数に渡す
    （以下省略）

Comments

[edatin]

ありがとうございました。解決しました。一番具体的でわかりやすかったのでベストアンサーにしました。

Answer 2

こんにちは。

ご提示のプログラム部分ではセグメンテーション違反は起きないようです。
それより後のどこかで起きているかも知れません。

printf()の出力はキャッシュされるため、まだキャッシュに溜まっていて画面にでていない可能性があります。printf()は正常に実行され、その後で落ちた場合にいかにも表示された部分までしか実行されていないかのように見え、判断を誤るというミスは割と多いです。

Comments

[edatin]

ご回答ありがとうございます。全体のプログラムを追加したのですが、どこで落ちる可能性があるか教えていただくことはできますでしょうか？

[Chironian]

iswポインタを初期化していないので、あらぬところを指しています。
そこをアクセスしているところで落ちている可能性が高いです。

[edatin]

回答ありがとうございます。返事が遅くなって申し訳ないです。無事解決しました。

Answer 3

  fscanf(fpin,"%lf",&z[a][b]);
 printf("%lf\t",z[a][b]);

双方、%lf ではなくて、%f です

---

とりあえず、

void lnqgj(double a[N][N+M], int n, int m, double epsl, int *isw)

と修正すれば実行できるようになるんじゃないかな

Comments

[edatin]

回答ありがとうございます。実行して見たところ、値がすべて０になり同じ箇所でセグメンテーションエラーが発生してしまいました。

[LouiS0616]

出力はともかく、入力はlfで問題ないように思いますが...
printfにlfを指定するのもC99以降では認められているはずです。

[y_waiwai]

fscanf(fpin,"%f",&(z[a][b]));
としましょう

[y_waiwai]

この関数全体を提示してもらわないと解決できないですね。。

[edatin]

いま、全体を修正したのですが、どうでしょうか。。

[y_waiwai]

回答を追記