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 vote413 views1 answer
Size of Operator in C language
C
14 votes377 views1 answer
Turn a char into an int in c
C
1 vote391 views1 answer
C get time in milliseconds
C
0 vote397 views1 answer
What is a while loop program in C to generate numbers from 0<25?
C
-3 votes388 views5 answers
Return char array in c
C
9 votes411 views9 answers
Selection sort
C
0 vote390 views1 answer
Desktop
C
1 vote387 views1 answer
Get pid c
C
1 vote389 views1 answer
Sdl_renderfillrect
C
0 vote423 views1 answer
How to compile and run a program in c
C