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 6d45657a6f4fd99a41aaec63e0c401cc6c3587d5..1b191a1b3f9a5d4696c4d535e94553e6b96f93e6 100644 (file)
--- 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>