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 vote368 views1 answer
Router solicitation and advertisement magic is used by
C
0 vote353 views1 answer
Vs code turn off formatter
C
1 vote362 views1 answer
Open a file in from terminal
C
0 vote372 views1 answer
Libdvd-pkg: `apt-get check` failed
C
2 votes356 views1 answer
Terminal size in c
C
0 vote374 views2 answers
Memcpy in c
C
0 vote369 views1 answer
Program optimization
C
10 votes353 views4 answers
Install redis
C
0 vote385 views4 answers
Convert preorder to postorder calculator
C
2 votes327 views1 answer
C string is int
C