C++ invert matrix

pockethook picture pockethook · Mar 8, 2012 · Viewed 12.6k times · Source

The following dynamic array contains a non-symmetric n*n matrix (with n <=100):

int **matrix;
matrix = new int*[n];
for (int i = 0; i < n; i++)
    matrix[i] = new int[n];

Is there an extremely easy way to invert it? Ideally I'd only use something from the STL or download a single header file.

Answer

linello picture linello · Mar 8, 2012

Using Eigen.

http://eigen.tuxfamily.org/index.php?title=Main_Page

You can map your array to an Eigen matrix and then perform efficient matrix inversion.

You must only include it.

I add that usually if you have to perform your inversion for linear system solving, it's better to use a matrix decomposition based on the properties of the matrix that you can exploit.

http://eigen.tuxfamily.org/dox/TutorialLinearAlgebra.html