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

0 vote267 views1 answer
Levenshtein c
C
0 vote260 views1 answer
How to reset all values of 2d vector to 0
C
0 vote272 views1 answer
Lsusb command not found
C
0 vote262 views1 answer
Atomic variable c
C
1 vote250 views1 answer
Comnetar toda linha do vs code
C
12 votes290 views2 answers
Owasp
C
0 vote249 views1 answer
Marquee html code with right
C
0 vote254 views7 answers
C++ convert to c online
C
0 vote269 views1 answer
Epita
C
0 vote250 views1 answer
Debian9 remove pack
C