This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: gsl_linalg_complex_householder_transform
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: Mario Pernici <Mario dot Pernici at mi dot infn dot it>
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Fri, 30 May 2003 15:38:34 +0100
- Subject: Re: gsl_linalg_complex_householder_transform
- References: <Pine.LNX.4.44.0305292229110.7094-100000@ptt.mi.infn.it>
Mario Pernici writes:
> Now my examples for the functions in the section Tridiagonal
> Decomposition of Hermitian Matrices work without modifications of
> the code of gsl. The key point is that gsl_complex
> gsl_linalg_complex_householder_transform (gsl_vector_complex * v)
> returns a complex number tau such that (I - tau^* v v^H) x = -/+
> ||x|| e_1 (see Working Note 72) while I thought the relation was
> (I - tau v v^H) x = -/+ ||x|| e_1 Maybe this should be added in
> the comments to the source code of this function.
Yes, I should have commented it before... I've added some comments
now, and the function gsl_linalg_complex_householder_hv.
Brian