ガウスの消去法

質問

逆行列を求めたいのですが
どうも上手くできません。
下にソースコードを載せますので
どなたかどこが悪いのか
ご教授してくれませんか？
よろしくお願いします。

void MINV(double A[][100],double InvAf2[][100],int G,int R){
int i,j,k;
double temp;//行列値の一時記憶
double B[100][100],C[100][100];//行列の一時記憶
for(k=0;k<R;k++){//k列目を調整する
for(i=0;i<G;i++){//行数を指す
for(j=0;j<R;j++){//列数を指す
temp=A[k][k];
B[i][j]=A[i][j];//InvAf2を計算するのに必要
C[k][j]=InvAf2[k][j];//InvAf2を計算するのに必要(更新した値を使われると困る)
if(i==k){
A[i][j]/=temp;
InvAf2[i][j]/=temp;
}
else{
A[i][j] -= A[k][j] * A[i][k] / temp;
InvAf2[i][j] -= C[k][j] * B[i][k] / temp;
}
}
}
}
}

回答No.1

BLUEPIXY · ベストアンサー

多分、
if(i==k){
A[i][j]/=temp;
InvAf2[i][j]/=temp;
}
は、i==kの時では遅すぎると思う。
あと、B,Cは、別に必要ない。
というか、自動変数（スタックに確保する数）としては、サイズが大きすぎるかもしれません。（それは、呼び出し元のプログラムにも言える）
staticとするか、mallocしてヒープに確保した方がいいと思う。
G,Rは、正則だと思うので、D（次元）だけで良いと思う。
以上の点に注意して書き直してみました。
void MINV(double A[][10],double InvAf2[][10],int D){
int i,j,k;
double temp;
for(k=0;k<D;k++){
temp=A[k][k];
for(j=0;j<D;j++){
A[k][j]/=temp;
InvAf2[k][j]/=temp;
}
for(i=0;i<D;i++){
if(k!=i){
temp=A[i][k];
for(j=0;j<D;j++){
A[i][j]-=A[k][j]*temp;
InvAf2[i][j]-=InvAf2[k][j]*temp;
}
}
}
}
}

回答No.2

BLUEPIXY

＃１です
void MINV(double A[][10],double InvAf2[][10],int D){
は、
void MINV(double A[][100],double InvAf2[][100],int D){
ですね。すみません。

お礼

お礼が遅れて大変申し訳ありませんでした。

おかげさまで、正しく逆行列を求めることが
できました。

今後も、お世話になることがあれば
どうぞよろしくお願いします

ありがとうございました。