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 vote141 views1 answer
Winautomation block input not working
C
0 vote169 views1 answer
4 byte alignment c code
C
4 votes131 views2 answers
Maximo comun divisor
C
0 vote144 views1 answer
Ubuntu get to local disk
C
39 votes321 views4 answers
Celsius to fahrenheit formula
C
0 vote131 views1 answer
Len of str vbs
C
5 votes140 views1 answer
Https://www.codegrepper.com/search.php?q=convert%20string%20to%20float%20c
C
0 vote142 views1 answer
Java Node<E>
C
0 vote142 views5 answers
[email protected]
C
22 votes150 views6 answers
ModuleNotFoundError: No module named 'cv2'
C