Код обратной матрицы и перемножения матрицна с++

Matrix Matrix::inverted() const
{
    Matrix temp(row, row); // будующая единичная матрица
    Matrix copy(*this); // копия матрицы
    int dim = row;   //размер матрицы
    //только если определитель не равен нулю
    if (det() != 0.0)
    {
        //сформируем единичную матрицу
        for (int i = 0; i < dim; ++i)
         temp.matrix[i][i] = 1.0;
 
        for (int k = 0; k < dim; ++k)
        {
            if (copy.matrix[k][k] == 0.0)
            {
                for (int i = k + 1; i < dim; ++i)
                {
                    if (copy.matrix[i][k] != 0.0)
                    {
                        swap(copy.matrix[k], copy.matrix[i]);
                        swap(temp.matrix[k], temp.matrix[i]);
 
                        break;
                    }
                }
            }
            double div = copy.matrix[k][k];
            for (int j = 0; j < dim; ++j)
            {
                copy.matrix[k][j] /= div;
                temp.matrix[k][j] /= div;
            }
            for (int i = k + 1; i < dim; ++i)
            {
                double multi = copy.matrix[i][k];
                for (int j = 0; j < dim; ++j)
                {
                    copy.matrix[i][j] -= multi * copy.matrix[k][j];
                    temp.matrix[i][j] -= multi * temp.matrix[k][j];
                }
            }
        }
        for (int k = dim - 1; k > 0; --k)
        {
            for (int i = k - 1; i + 1 > 0; --i)
            {
                double multi = copy.matrix[i][k];
 
                for (int j = 0; j < dim; ++j)
                {
                    copy.matrix[i][j] -= multi * copy.matrix[k][j];
                    temp.matrix[i][j] -= multi * temp.matrix[k][j];
                }
            }
        }
    }
    //иначе вернём ноль-матрицу
    else
        temp.resize();
    return temp;
}