[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
## Re: [Help-gsl] Re: QR Decomposition for a matrix with complex values

**From**: |
Brian Gough |

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

**Date**: |
Wed, 30 Jan 2008 21:52:43 +0000 |

**User-agent**: |
Wanderlust/2.14.0 (Africa) Emacs/22.1 Mule/5.0 (SAKAKI) |

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/