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 vote349 views1 answer
Enum case statement in c
C
0 vote417 views1 answer
Ringing a bell using c
C
0 vote361 views1 answer
Block a website on mac
C
0 vote354 views1 answer
Turn a char array into double C
C
0 vote374 views1 answer
Matplotlib legend remove box
C
3 votes407 views1 answer
Error: 'for' loop initial declarations are only allowed in C99 or C11 mode
C
1 vote362 views2 answers
Data types arduino
C
1 vote346 views1 answer
What is system function in c
C
0 vote416 views1 answer
Actionbar content color in android
C
0 vote354 views1 answer
Cyrildewit laravel page view counter package.
C