template<typename T>
-inline T &invert_matrix(T &m, T &r)
+inline T &gauss_jordan(T &m, T &r)
{
typedef typename T::ElementType V;
using std::abs;
for(unsigned j=i; j-->0; )
r.add_row(i, j, -m.element(j, i));
- return m = r;
+ return r;
}
} // namespace LinAl