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 vote178 views1 answer
Cannot get / react router dom
C
0 vote170 views3 answers
Print 2d array in c
C
0 vote166 views1 answer
Fine print serial key
C
0 vote202 views1 answer
Scanf ignore new line
C
0 vote266 views1 answer
Declare type c
C
0 vote169 views1 answer
Reap zombie process in c
C
0 vote172 views1 answer
Solana-test-validator log
C
-1 vote191 views1 answer
C check if file was created
C
0 vote168 views1 answer
Error: argument 1 range [18446744071562067968, 18446744073709551615] exceeds maximum object size 9223372036854775807 [-werror=alloc-size-larger-than=]
C
1 vote178 views2 answers
Slack wordpress
C