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 vote227 views1 answer
C gettimeofday example
C
2 votes233 views1 answer
Int8_t is not defined
C
0 vote231 views1 answer
Mysql yyyymm format
C
-1 vote239 views1 answer
Remove last entry from array in Javascript
C
0 vote246 views1 answer
Bd number regex
C
18 votes211 views1 answer
Read files in c
C
0 vote227 views1 answer
Grep C hello world
C
0 vote227 views1 answer
Printf("%3d ",XX);
C
33 votes302 views2 answers
Latex font sizes
C
7 votes180 views2 answers
Install c compiler ubuntu
C