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 vote306 views1 answer
C static
C
0 vote330 views1 answer
Sdl_rect
C
0 vote326 views1 answer
Warning: function returns address of local variable [-Wreturn-local-addr]
C
27 votes354 views6 answers
Hello world c
C
0 vote293 views1 answer
Entity framework core discard changes
C
0 vote326 views1 answer
Libreoffice reference cell in different sheet with sheet name with space
C
0 vote327 views1 answer
Code to reverse the words in a sentnce
C
0 vote318 views2 answers
Memcpy in c
C
23 votes373 views6 answers
How to check if a number is even or odd in python
C
2 votes328 views1 answer
Int8_t is not defined
C