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 vote146 views1 answer
Transform yt video into background overlay
C
1 vote153 views1 answer
":" in struc c
C
4 votes141 views1 answer
Typedef c struct
C
1 vote138 views1 answer
How to know ip address of my pc
C
2 votes98 views2 answers
Swap variables in c
C
0 vote153 views1 answer
Syntax of for loop in c stack over flow
C
1 vote139 views2 answers
Strstr
C
0 vote151 views1 answer
Snprintf with malloc
C
1 vote146 views6 answers
Online linux terminal for c programming
C
0 vote142 views1 answer
Bd number regex
C