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 vote325 views1 answer
What is wrong with this in c
C
0 vote343 views2 answers
Time include c
C
0 vote347 views1 answer
Rust set toolchain
C
1 vote325 views1 answer
C float remove trailing 0
C
0 vote353 views1 answer
Assets/Scripts/UIManager.cs(65,13): error CS0201: Only assignment, call, increment, decrement, await, and new object expressions can be used as a statement
C
0 vote363 views1 answer
C defined value sum
C
0 vote357 views1 answer
Search array element in c
C
0 vote357 views1 answer
Bp result system
C
0 vote381 views1 answer
Manasa loves maths solution IN C
C
0 vote352 views1 answer
How many archaeologists are there
C