N次正方行列AとN次元ベクトルxの積y=Axを計算するプログラム

Question

c言語でN次正方行列AとN次元ベクトルxの積y=Axを計算するプログラムを作っているのですが、
最初のscanfで数値を入力したらsegmentation faultと表示されてしまいます。何が原因なのでしょうか。

#include<stdio.h>
#include<stdlib.h>
void PrintVector(int *y, int n){
  int i;
  for(i=0;i<n;i++){
    printf("%6d",*(y+i));
  }
}
void PrintMatrix(int *a, int n){
  int i,j;
  int *ai;
  ai= a;
  for( i = 0 ; i < n ; i++){
    for( j = 0 ; j < n ; j++){
      printf("%6d", *(ai+j));
    }
    ai+=n;
    printf("\n");
  }
  printf("\n");
}
int main(){
  int i, j,num;
  int *ai;
int *a, *x, *y;  
printf("N・・");
  scanf("%d",&num);
  
  ai= (int *) malloc(num*sizeof(int));
  ai=a;
  for(i=0;i<num;i++){
    for(j=0;j<num;j++){
      printf("A[%d][%d]=",i,j);
      scanf("%d",ai+j);
     }
   ai++;
  }
 
 for(i=0;i<num;i++){
      printf("X[%d]=",i);
      scanf("%d",x+i);
    }
  printf("A = \n"); PrintMatrix(a,num);
  printf("x = \n"); PrintVector(x,num);
  for(i=0;i<num;i++){
    *(y+1)=0;
    for(j=0;j<num;j++){
      *(y+i)+=*(ai+j)**(x+j);
    }
    ai+num;
  }
  free(ai);
 
  return 0;
}

Answer 1

コンパイラの警告レベルを最大にして警告内容を確認してみてはどうでしょうか。おかしなところがみつかるかと思います。

Comments

[YOshim]

自分の環境だと最初のscanfではエラーは出ませんでした。2回目のscanfで入力しても反応が芳しくなかったですが、
質問のエラー発生は2回目のscanfではないでしょうか。

Answer 2

3点修正したら正常と思われる表示が得られました。

main関数内にai = aとあるが、aは実態が伴っていないのでおかしい。
ai++と操作しているのでaをaiの先頭番地として記憶させたいのであれば
a = aiとするのが正しい。

xとyがnullポインタになっている。
x,yのサイズを動的に変えたいなら、ai同様にx,yそれぞれにmallocでメモリを割り当てるのがよい。
動的に変える必要がないなら、int x[100],y[100]でもよい。

main関数の最初のfor文中にai++とあるが、これはai+=numの誤りだと思われます。

Comments

[ksyiwk]

上記の3点を修正したのですが、行列表示ができません。

[YOshim]

修正箇所が違うのではないでしょうか。
main関数内だけですよ。

あとエラー内容が分かればより具体的な指摘ができると思います。

[ksyiwk]

すみませんPrintMatrix内もa = aiとしていました

Answer 3

こんにちわ。

ai= (int *) malloc(num*sizeof(int));

上記のコードで確保される整数型配列のサイズはNであり、N*Nではありません。
なので、その先のコードで領域外参照が発生しているのではないでしょうか。