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 vote204 views2 answers
Pointeur de pointeur en language c
C
17 votes201 views2 answers
What is chocolatey
C
0 vote244 views5 answers
Online c compiler with mpi
C
1 vote207 views1 answer
Markdown empty line
C
0 vote200 views1 answer
MongoDb Read
C
2 votes205 views1 answer
Can we write a program without main in c
C
0 vote199 views1 answer
What is char in programming
C
1 vote211 views1 answer
Access denied creating xampp-control.ini
C
9 votes222 views1 answer
C how to get an integer from user input
C
0 vote208 views2 answers
Print integer to stdout using write or putchar?
C