File kernel/src/utils/SiconosAlgebra/InvertMatrix.hpp

Go to the source code of this file

The following code inverts the matrix input using LU-decomposition with backsubstitution of unit vectors.

Reference: Numerical Recipies in C, 2nd ed., by Press, Teukolsky, Vetterling & Flannery.

you can solve Ax=b using three lines of ublas code:

permutation_matrix<> piv; lu_factorize(A, piv); lu_substitute(A, piv, x);


template<class T, class U, class V>
bool InvertMatrix(const boost::numeric::ublas::matrix<T, U, V> &input, boost::numeric::ublas::matrix<T, U, V> &inverse)

Matrix inversion routine.

Uses lu_factorize and lu_substitute in uBLAS to invert a matrix

  • input: source matrix

  • output: inverted matrix.