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

1 vote239 views2 answers
Isspace
C
-1 vote248 views2 answers
Printf n characters c
C
0 vote243 views1 answer
Insert image material ui
C
0 vote225 views1 answer
Printf("%d", 10 ? 0 ? 5:1:1:12) what will print
C
0 vote247 views2 answers
Route::resource
C
1 vote235 views2 answers
String input in c
C
0 vote263 views1 answer
Bitwise Operators in C language
C
0 vote234 views1 answer
Pygraphviz show
C
1 vote231 views1 answer
False and true in c
C
0 vote246 views1 answer
Como somar em C
C