C言語です

質問

#include<stdio.h>
int main(void){
int matA[4][4], matB[4][4], matC[4][4];
int i, j, k;
/* 行列A の読み込み */
for(i = 0; i < 4; ++i){
for(j = 0; j < 4; ++j){
scanf("%d",&matA[i][j]);
}
}
/* 行列B の読み込み */
for(i = 0; i < 4; ++i){
for(j = 0; j < 4; ++j){
scanf("%d",&matB[i][j]);
}
}
/*行列の積の計算*/
for(i = 0; i < 4; ++i){
for(j = 0; j < 4; ++j){
matC[i][j]=0; /*初期化が必要*/
for(k = 0; k < 4; ++k){
matC[i][j]+=matA[i][k]*matB[k][j];
}
}
}
/*結果の表示*/
for(i = 0; i < 4; ++i){
for(j = 0; j < 4; ++j){
printf("%d\t",matC[i][j]);
}
printf("\n");
}
return 0;
}

このプログラムについてですが、
/*行列の積の計算*/
for(i = 0; i < 4; ++i){
for(j = 0; j < 4; ++j){
matC[i][j]=0; /*初期化が必要*/
for(k = 0; k < 4; ++k){
matC[i][j]+=matA[i][k]*matB[k][j];
}
}
}

この部分の初期化がfor(k = 0; k < 4; ++k){の部分の後ろだったら
実行結果が異なります。
なぜ、
for(k = 0; k < 4; ++k){
matC[i][j]=0; /*初期化が必要*/

この順番だといけないんですか？

回答No.1

hirotn · ベストアンサー

+= 演算子を使っているので、積算しています（計算結果をmatCに足し続ける）

しかし、for(k = 0; k < 4; ++k){のループの度にmatCを初期化してしまうので、最後の掛け算しか残っていないことになる。

と読めます。

回答No.4

KEIS050162

計算の過程が知りたいとき、部分的にデバッグ用に、printf() 文とgetch()　などを組み合せて追加し、処理の過程などを確認してみると良いでしょう。デバッグが終わったら、コメントにしてしまえばOKです。

for(k = 0; k < 4; ++k){
　　　matC[i][j]+=matA[i][k]*matB[k][j];

　　　printf( "matc{ %d ][ %d ] = %d\n". i, i, matc[i][i] );
　　　getch();

デバッグが終わったら、コメントにする。

//　　　printf( "matc{ %d ][ %d ] = %d\n". i, i, matc[i][i] );
//　　　getch();

色んな場所に上記の様なprintf（）を入れ込んで試してみてください。

回答No.3

Wr5

ステップ実行のやりかたとか調べるか、pritnf()などで変数の内容出力してみるとかしないんですか？
質問乱立しまくりのようですが……。

for(i = 0; i < 4; ++i){
　for(j = 0; j < 4; ++j){
　　matC[i][j]=0; /*初期化が必要*/
　　for(k = 0; k < 4; ++k){
　　　matC[i][j]+=matA[i][k]*matB[k][j];
　　}
　}
}
を
for(i = 0; i < 4; ++i){
　for(j = 0; j < 4; ++j){
　　for(k = 0; k < 4; ++k){
　　　matC[i][j]=0; /*初期化が必要*/
　　　matC[i][j]+=matA[i][k]*matB[k][j];
　　}
　}
}
にしたら、結果が異なるのはな～んでだろ～？

って質問でいいんですかね？
for(k = 0; k < 4; ++k)のループで
matC[i][j]=0;により0が入ります。
次にmatC[i][j]+=matA[i][k]*matB[k][j];の演算結果を代入します。
次のループに入って、
matC[i][j]=0;により0が入ります。
# つまり「以前の演算結果はご破算にしろ」と言っています。
それを繰り返して、結果として有効になるのは最後の演算結果だけです。
「そうしろ」というコードなのですから、CPUは言われたとおりに忠実に実行しただけです。

回答No.2

hirotn

matC[i][j]+=matA[i][k]*matB[k][j];を、
matC[i][j]=matA[i][k]*matB[k][j];　として試行してみてください。実行結果が異なる場合とおなじになりませんか？