Gsl matrix invert

gsl_matrix *
invert_a_matrix(gsl_matrix *matrix)
{
    gsl_permutation *p = gsl_permutation_alloc(size);
    int s;

    // Compute the LU decomposition of this matrix
    gsl_linalg_LU_decomp(matrix, p, &s);

    // Compute the  inverse of the LU decomposition
    gsl_matrix *inv = gsl_matrix_alloc(size, size);
    gsl_linalg_LU_invert(matrix, p, inv);

    gsl_permutation_free(p);

    return inv;
}

Source: gist.github.com

Similar codes

10 votes272 views2 answers
Int num = (rand() % (upper - lower + 1)) + lower;
C
4 votes224 views2 answers
Grepper vscode
C
0 vote228 views3 answers
Read from command line c
C
0 vote218 views1 answer
FILE*
C
2 votes236 views1 answer
%.*s in c
C
-1 vote218 views2 answers
Size_t printf
C
0 vote233 views1 answer
Postgres login mac
C
0 vote218 views2 answers
Logical operators in c
C
0 vote212 views1 answer
Unox reclame harmonica tabs
C
2 votes232 views2 answers
How to delete virtual hard disk virtualbox
C