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 vote315 views1 answer
I like to play games
C
1 vote344 views2 answers
Define constant c
C
1 vote314 views1 answer
Operators priority in c
C
0 vote313 views1 answer
Include ‘<stdlib.h>’ or provide a declaration of ‘exit’
C
5 votes246 views1 answer
Printf float format
C
0 vote312 views1 answer
Oracle trunc
C
1 vote287 views1 answer
C null define
C
11 votes329 views6 answers
Tailwindcss install
C
0 vote306 views1 answer
How to zoom in terminal
C
2 votes574 views7 answers
Linked list
C