X-Git-Url: http://git.tdb.fi/?a=blobdiff_plain;f=source%2Flinal%2Fsquarematrix.h;h=b8531aafcb8359bcda6ea35628a9f561f77b677b;hb=fb1e198;hp=0a97e61172b6489c58e91296421e946cfddfb9d8;hpb=5797fe0a5296952cb8f8643fdc6cabddee19a554;p=libs%2Fmath.git diff --git a/source/linal/squarematrix.h b/source/linal/squarematrix.h index 0a97e61..b8531aa 100644 --- a/source/linal/squarematrix.h +++ b/source/linal/squarematrix.h @@ -47,14 +47,16 @@ template SquareMatrix &SquareMatrix::invert() { SquareMatrix r = identity(); - return invert_matrix(*this, r); + gauss_jordan(*this, r); + return *this = r; } template inline SquareMatrix invert(const SquareMatrix &m) { - SquareMatrix r = m; - return r.invert(); + SquareMatrix temp = m; + SquareMatrix r = SquareMatrix::identity(); + return gauss_jordan(temp, r); } } // namespace LinAl