In a worst case scenario, someone can use eigendecomposition (gsl_eigen_hermv) to find [Q,Λ] and from there : Q * inv(Λ) * (Q') can do the trick... Right? ( inv(Λ) should be trivial, as it is diagonal and real)
--- On Thu, 11/17/11, Patrick Alken <al...@colorado.edu> wrote: From: Patrick Alken <al...@colorado.edu> Subject: Re: [Help-gsl] Inverse of complex Hermitian matrices? To: Cc: help-gsl@gnu.org Date: Thursday, November 17, 2011, 7:44 PM Ah oops I was looking at an old git repository On 11/17/2011 06:39 PM, Rhys Ulerich wrote: > On Thu, Nov 17, 2011 at 11:32 AM, Patrick Alken<al...@colorado.edu> wrote: >> It seems that when gsl_linalg_cholesky_invert was added to gsl (v1.12), its >> complex counterpart wasn't. There doesn't seem to be a >> gsl_linalg_complex_cholesky_invert so this is a doc bug. > I see gsl_linalg_complex_cholesky_invert within linalg/choleskyc.c on > trunk. NEWS shows it appeared in 1.15 courtesy of Huan Wu. > > - Rhys _______________________________________________ Help-gsl mailing list Help-gsl@gnu.org https://lists.gnu.org/mailman/listinfo/help-gsl _______________________________________________ Help-gsl mailing list Help-gsl@gnu.org https://lists.gnu.org/mailman/listinfo/help-gsl
