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 vote164 views1 answer
Deleting a word with copy fuction c code
C
2 votes119 views2 answers
How to delete virtual hard disk virtualbox
C
0 vote114 views1 answer
C how to check a palindrome string
C
0 vote112 views1 answer
Gemfile.lock`. It is likely that you need to grant write permissions for that path.
C
0 vote102 views1 answer
Gcc options to find out makefiel rules
C
23 votes148 views3 answers
Picasso
C
1 vote109 views1 answer
Which of the following are Cetaceans?
C
1 vote112 views1 answer
Why do you need the return 0 c
C
0 vote108 views1 answer
Pasar a binario recursivo
C
0 vote115 views1 answer
Odd even in c with ternary operator
C