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 vote129 views1 answer
What is c
C
1 vote133 views1 answer
C bit access union
C
0 vote140 views1 answer
How to use arry
C
18 votes132 views9 answers
While loop
C
5 votes129 views1 answer
Type change in c
C
-3 votes154 views2 answers
How to return a char array from a function in C
C
1 vote170 views3 answers
Ruby inject
C
3 votes94 views1 answer
Scala list
C
0 vote125 views1 answer
::template
C
0 vote135 views1 answer
Fgets em c
C