From: Mikko Rasa <tdb@tdb.fi>
Date: Sun, 2 Jun 2019 13:23:08 +0000 (+0300)
Subject: Fix a memory access error in DynamicMatrix::invert
X-Git-Url: http://git.tdb.fi/?p=libs%2Fmath.git;a=commitdiff_plain;h=459c248db8804b2a51010c89b72d3a0d326267e9

Fix a memory access error in DynamicMatrix::invert

Gauss_jordan returns a reference to the result, which is a local
variable.  It needs to be assigned to *this to retain it.
---

diff --git a/source/linal/dynamicmatrix.h b/source/linal/dynamicmatrix.h
index 6d45657..1b191a1 100644
--- a/source/linal/dynamicmatrix.h
+++ b/source/linal/dynamicmatrix.h
@@ -308,7 +308,8 @@ inline DynamicMatrix<T> &DynamicMatrix<T>::invert()
 	for(unsigned i=0; i<rows_; ++i)
 		r(i, i) = T(1);
 
-	return gauss_jordan(*this, r);
+	gauss_jordan(*this, r);
+	return *this = r;
 }
 
 template<typename T>