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 vote393 views1 answer
Dynamic Memoray alocation For 2D
C
0 vote392 views1 answer
Program optimization
C
0 vote398 views1 answer
Sum average min max in c array
C
0 vote393 views1 answer
C special character display
C
0 vote373 views1 answer
Pygraphviz show
C
0 vote411 views1 answer
Como hacer para que una salida en linux aparezca de poco en poco
C
23 votes461 views6 answers
How to check if a number is even or odd in python
C
0 vote343 views1 answer
11*179*.9*1.35
C
-1 vote413 views1 answer
Matrix of string in c
C
1 vote446 views1 answer
Scranton inhabitants
C