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 vote174 views1 answer
Router solicitation and advertisement magic is used by
C
0 vote184 views1 answer
Strlen stop before '\0' or at
C
0 vote149 views1 answer
Calling of a void in c
C
0 vote175 views1 answer
Sdl_rect
C
5 votes154 views1 answer
Array loop in c
C
1 vote160 views1 answer
Difference between malloc and calloc and realloc
C
0 vote163 views1 answer
Xor in c
C
0 vote185 views1 answer
Disable gnu++11 option
C
2 votes170 views1 answer
Numeri in ordine crescente C
C
0 vote149 views3 answers
Script in c
C