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 vote411 views1 answer
Unable to locate package dos2unix
C
0 vote397 views1 answer
Gitlab ci heroku
C
0 vote451 views1 answer
Iterating through a linked list
C
10 votes444 views2 answers
Malloc c include
C
-1 vote387 views1 answer
Free array in c
C
0 vote416 views1 answer
Pasar a binario recursivo
C
0 vote410 views1 answer
Two way communication between child and parent processes in C using pipes
C
1 vote401 views1 answer
Install kubernetes kubectl on mac
C
1 vote387 views1 answer
C iterate string
C
0 vote405 views1 answer
Homeworld 3
C