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

1 vote184 views1 answer
Directory folders structure show windows 10 command prompt
C
1 vote203 views1 answer
ERROR: cannot verify nvidia.box.com's certificate, issued by 'CN=GeoTrust RSA CA 2018,OU=www.
C
5 votes212 views2 answers
C strcmp
C
0 vote320 views1 answer
Largest value in u16
C
0 vote248 views1 answer
How to break a loop in c
C
0 vote285 views1 answer
How to arrange a 2d array based on string length in c
C
3 votes224 views2 answers
MongoDb Create
C
0 vote210 views10 answers
Php
C
0 vote193 views1 answer
C input is skipped
C
4 votes222 views2 answers
Grepper vscode
C