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 vote178 views1 answer
Print 0 1 2 3 4 in c while loop
C
0 vote151 views1 answer
How to ban websites on mac
C
0 vote172 views1 answer
<img src="https://www.google.com/url?sa=i&url=https%3A%2F%2Ftwitter.com%2Fimgatennis%3Flang%3Dar-x-fm&psig=AOvVaw1QH9um9Q1ktPXKUufXt8Av&ust=1636491440205000&source=images&cd=vfe&ved=0CAsQjRxqFwoTCKi9v-rTifQCFQAAAAAdAAA
C
5 votes181 views2 answers
Convert string to float c
C
0 vote158 views1 answer
Object-fit not working
C
0 vote201 views1 answer
KneesDev
C
0 vote160 views1 answer
C read a whole string from a file
C
0 vote172 views1 answer
Algorithm that flips sentences and numbers
C
1 vote177 views1 answer
Material-ui Roboto font
C
0 vote166 views1 answer
Mark rober
C