This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: cholesky decomposition bug?
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: Pedro Gonnet <pedro at vis dot ethz dot ch>
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Sun, 3 Nov 2002 14:59:28 +0000 (GMT)
- Subject: Re: cholesky decomposition bug?
- References: <1036162079.11311.48.camel@selonia>
Pedro Gonnet writes:
> the c-code that does the decomposition is, abreviated:
>
> A = gsl_matrix_view_array(H,n,n);
> b = gsl_vector_view_array(g,n);
> x = gsl_vector_view_array(x_delta,n);
> gsl_set_error_handler_off();
> if ((res = gsl_linalg_cholesky_decomp(&A.matrix)) != 0) {
> ...
> }
> gsl_linalg_cholesky_solve(&A.matrix,&b.vector,&x.vector);
>
> this stuff works sometimes, so i don't think it's a coding problem on my
> behalf.
Please can you send a complete standalone program which reproduces the
problem, thanks.
Brian