Matrix division works!!

Author: potat-dev

Variants/Lab_12/var_12/main.cpp

@@ -186,4 +186,38 @@ int main() {
 
   det_test.minor_matrix().print();
   det_test.algebraic_additions().print();
+
+  Matrix aaa({
+    {1, 2},
+    {3, 4}
+  });
+  aaa.print();
+
+  aaa.inverse().print();
+
+  Matrix a100({
+    {1,  18, 4,  6,  9 },
+    {11, 17, 13, 14, 8 },
+    {25, 2,  19, 15, 21},
+    {10, 5,  12, 3,  7 },
+    {24, 20, 22, 16, 23}
+  });
+  a100.print();
+
+  cout << a100.determinant() << endl << endl;
+  a100.inverse().print();
+
+  Matrix mult_a({
+    {13, 26},
+    {39, 13}
+  });
+
+  Matrix mult_b({
+    {7, 4},
+    {2, 3}
+  });
+
+  // для проверки
+  // https://ru.wikihow.com/делить-матрицы
+  (mult_a / mult_b).print();
 }

Variants/Lab_12/var_12/matrix.cpp

@@ -176,6 +176,12 @@ Matrix Matrix::algebraic_additions() {
   return temp;
 }
 
+Matrix Matrix::inverse() {
+  Matrix alg_add_T = this->algebraic_additions().transpose();
+  alg_add_T /= this->determinant();
+  return alg_add_T;
+}
+
 // операторы
 
 Matrix& Matrix::operator+=(const Matrix& m) {
@@ -258,6 +264,23 @@ Matrix operator*(Matrix a, const double &b) {
   return a;
 }
 
+Matrix& Matrix::operator/=(const Matrix &b) {
+  *this *= Matrix(b).inverse();
+  return *this;
+  
+  /* поэлементное деление - не гуд
+  for (int x = 0; x < b.arr.size(); x++)
+    for (int y = 0; y < b.arr.size(); y++)
+      arr[y][x] /= b.arr[y][x];
+  return *this;
+  */
+}
+
+Matrix operator/(Matrix a, const Matrix &b) {
+    a /= b;
+    return a;
+}
+
 Matrix& Matrix::operator/=(const double& n) {
   for (int x = 0; x < arr.size(); x++) {
     for (int y = 0; y < arr.size(); y++) {