1) How did I store the lower triangular matrix into a SIZExSIZE matrix:
for (int i = 0; i < SIZE; i++)
{
for (int j = 0; j <= i; j++)
{
*(L->data + (i*SIZE + j)) = *(M->data + (i*SIZE + j));
}
}2) How did I compute the inverse if that (lower triangular) matrix: gsl_linalg_LU_decomp(L, p, &s); gsl_linalg_LU_invert(L, p, Linv); ---John _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
