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 vote409 views3 answers
Console clear c
C
3 votes388 views3 answers
Get the latest field in mongodb collection
C
2 votes362 views1 answer
Vscode arduino [Warning] Output path is not specified. Unable to reuse previously compiled files. Upload could be slow. See README.
C
2 votes391 views5 answers
C array of strings
C
0 vote381 views1 answer
Bootsrap textbox
C
0 vote391 views1 answer
Como somar em C
C
3 votes381 views1 answer
Initialize array in c with 0
C
0 vote376 views1 answer
Function inside struct c
C
1 vote413 views1 answer
And in language c
C
22 votes423 views2 answers
C bind str and int
C