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/ _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
