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 vote115 views1 answer
Winautomation block input not working
C
1 vote123 views1 answer
Error: dereferencing pointer to incomplete type
C
2 votes114 views1 answer
Powershell search files for string
C
0 vote120 views1 answer
PATH_MAX
C
0 vote129 views1 answer
Parcel-builder put image in sub folder
C
0 vote117 views2 answers
Int to void* c
C
4 votes112 views1 answer
Typedef c struct
C
0 vote119 views1 answer
Remplir list Java rapidemnet
C
0 vote127 views1 answer
C read csv
C
0 vote114 views1 answer
Turn a char array into double C
C