Re: [Help-gsl] Re: QR Decomposition for a matrix with complex values

[Brian Gough]

At Tue, 29 Jan 2008 11:46:17 +0000,
Jonny Taylor wrote:
> 
> OK, it turns out the problem was my confusion over exactly how the  
> complex householder transformations work. With that sorted, I can get  
> it to work.
> 
> The comment in householdercomplex.c:
> "To compute H'v, pass the conjugate of tau as the first argument to  
> gsl_linalg_householder_hm() rather than tau itself"
> was the bit I hadn't appreciated. Shouldn't this be mentioned  
> somewhere in the manual? It's only by reading the source that you  
> discover this pretty crucial point.

Thanks. I'll put that in the manual.

Would you be able to write some tests (similar to the existing ones in
linalg/test.c) so your complex QR can be included?

-- 
Brian Gough

GNU Scientific Library -
http://www.gnu.org/software/gsl/