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 vote335 views1 answer
Worst case algorithm
C
0 vote336 views1 answer
Yt-project annotate_scale
C
0 vote346 views1 answer
Ecs Invalid 'cpu' setting for task
C
6 votes346 views2 answers
C substring
C
0 vote340 views1 answer
Table de hachage en c
C
0 vote342 views1 answer
Bullseye lxc network problem
C
0 vote361 views1 answer
Postgres login mac
C
0 vote342 views1 answer
How to create calculator with switch in c
C
0 vote321 views1 answer
What is char in programming
C
0 vote347 views1 answer
Docker-compose file break line
C