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 vote186 views2 answers
Bash md5sum string
C
0 vote170 views1 answer
FivemStore
C
1 vote158 views1 answer
-> C
C
6 votes152 views1 answer
If statement shorthand c
C
2 votes293 views3 answers
Convert curl to http request with authorization header
C
1 vote111 views1 answer
How to install vscode extension manually
C
0 vote159 views1 answer
Onvert a string into 2d string in c
C
0 vote174 views1 answer
Print 0 1 2 3 4 in c while loop
C
0 vote152 views1 answer
FILE*
C
11 votes162 views6 answers
Tailwindcss install
C