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

2 votes139 views2 answers
End task cmd
C
0 vote210 views1 answer
Largest value in u32
C
1 vote209 views1 answer
TypeError: Incorrect type of self (must be 'Feature2D' or its derivative)
C
0 vote195 views1 answer
C program boilerplate code
C
45 votes344 views3 answers
String concatenation in c
C
2 votes196 views5 answers
C array of strings
C
0 vote193 views1 answer
Pendu langage c
C
0 vote196 views1 answer
Sorted characters in a string in c
C
4 votes204 views2 answers
C get time
C
5 votes200 views2 answers
Latex noindent
C