逆行列を求める関数を作りたい

Question

2次と3次との正方行列の逆行列を求める関数を作っています。2次と3次では、逆行列の求め方が異なるため、一つの関数でどちらも計算できるコードを作ることは難しいと思い、ひとまず、の逆行列を求める関数を作りました。しかし、aの絶対値が正しく認識されないようで、dt2[2][2]={{0,0},{0,0}}と表示されてしまいます。改善すべき点を教えていただきたいです。
やはり、3次は別に定義するしかないですか。

よろしくお願いいたします。

#include <stdio.h>
#include<stdlib.h>

#define ARRAY_SIZE 2

int gyaku(int dt1[2][2], int dt2[2][2]);

int gyaku(int dt1[2][2], int dt2[2][2]) {
	//二次の場合の逆行列
		int a = 0;
		int i, j = 0;
		a = ((dt1[0][0] * dt1[1][1]) - (dt1[0][1] * dt1[1][0]));
		int b = 0;
		b = abs(a);

		if (a = 0) {
			printf("Error!\n");
		}
		else {
			for (i = 0; i < ARRAY_SIZE; i++) {
				for (j = 0; j < ARRAY_SIZE; j++) {
					
					dt2[i][j] = (dt1[i][j] / b);
				}
			}
		}

		return dt2[i][j];
}

int main(void) {
	int dt1[2][2] = { {2,4}, {5,6} };
	int dt2[2][2] = {0,0};

	gyaku(dt1, dt2);

	for (int i = 0; i < ARRAY_SIZE; i++) {
		for (int j = 0; j < ARRAY_SIZE; j++) {
			printf("%2d", dt2[i][j]);
		}
		printf("\n");
	}
}

Answer 1

ちらっとみて
if (a = 0)
は if (a == 0) とか if (b <= delte) # delta は十分小さな値を設定しておく)
の間違いと思われます。

逆行列に関する参考情報

• 行列式の値の求め方、逆行列の作り方の C 言語プログラム

http://thira.plavox.info/blog/2008/06/_c.html

...
逆行列の求め方
何といってもやはり掃き出し法が簡単。2 次から n 次まで対応
...

• 逆行列を求める２通りの方法と例題

https://mathtrain.jp/inversematrix

...
2×2の場合は覚えるのみです
...
一般の n×n の正方行列の逆行列を求める二通りの方法を解説します
...
逆行列の求め方1：掃き出し法
...
逆行列の求め方2：余因子を用いる
...

Comments

[h_proc]

回答ありがとうございます。参考にして作ってみたのですが、nan,infが出てきてしまいます。どこが間違っているのでしょうか。
#include <stdio.h>
#include<stdlib.h>

#define ARRAY_SIZE 3

double gyaku(double dt1[ARRAY_SIZE][ARRAY_SIZE], double dt2[ARRAY_SIZE][ARRAY_SIZE]);

double gyaku(double dt1[ARRAY_SIZE][ARRAY_SIZE], double dt2[ARRAY_SIZE][ARRAY_SIZE]) {
//逆行列の計算
double a = 0;
int i, j = 0;
for (i = 0; i < ARRAY_SIZE; i++) {
for (j = 0; j < ARRAY_SIZE; j++) {
dt2[i][j] = (i == j) ? 1.0 : 0.0;
}
}
for (i = 0; i < ARRAY_SIZE; i++) {
a = 1 / dt1[i][i];
for (j = 0; j < ARRAY_SIZE; j++) {
dt1[i][j] *= a;
dt2[i][j] *= a;
}
for (j = 0; j < ARRAY_SIZE; j++) {
if (i != j) {
a = dt1[j][i];
for (int k = 0; k < ARRAY_SIZE; k++) {
dt1[j][k] -= dt1[i][k] * a;
dt2[j][k] -= dt2[i][k] * a;
}
}
else {
//NR
}
}
}

return dt2[ARRAY_SIZE][ARRAY_SIZE];
}

int main(void) {
double dt1[3][3] = { {1,2,4}, {5,6,8} };
double dt2[3][3] = {0,0,0};

gyaku(dt1, dt2);

for (int i = 0; i < ARRAY_SIZE; i++) {
for (int j = 0; j < ARRAY_SIZE; j++) {
printf("%2f", dt2[i][j]);
}
printf("\n");
}
}

[episteme]

double dt1[3][3] = { {1,2,4}, {5,6,8} }; だと3行目がないよ。マトモなコタエは得られんじゃろ。

おまけ: gyaku の戻り値は何を返したいのかしら? 返すもんがないなら void でええよ。

[h_proc]

つまらないミスをしていました。解決しました。ありがとうございます。