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

3 votes229 views2 answers
Word wrap intellij
C
0 vote274 views1 answer
Uri/Beecrowd problem no - 1151 solution in C
C
0 vote285 views1 answer
C sort odd even program
C
1 vote275 views1 answer
Dynamic 2d arr in c
C
0 vote249 views1 answer
C math.h sqrt
C
0 vote273 views1 answer
String insertion into another string
C
0 vote303 views3 answers
Swap two numbers in c
C
1 vote239 views1 answer
Stdio
C
1 vote288 views1 answer
Add border to image android
C
0 vote309 views1 answer
How to remove "Web site created using create-react-app."
C