[Help-gsl] Re: Inverse of the bidiagonals from gsl_linalg_cholesky_decomp(...) ?

[John Chludzinski]

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