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

266 votes258 views5 answers
Font-awesome - cdnjs.com - The best FOSS CDN for web
C
1 vote230 views2 answers
Define constant c
C
1 vote203 views1 answer
Grep find and replace
C
0 vote205 views2 answers
Sort linked list C
C
0 vote248 views1 answer
Hostbuilder add environment variables
C
2 votes194 views1 answer
Read from stdin c
C
1 vote214 views1 answer
Lichess puzzles
C
0 vote193 views1 answer
C/c drop mime
C
4 votes205 views3 answers
Multiplication table in c using array
C
0 vote201 views1 answer
Vbl share price
C