Re: [Bug-gsl] problem in cubic solver

[Brian Gough]

At Tue, 28 Apr 2009 17:11:31 +0200,
Lorenzo Moneta wrote:
>   This patch fails the current test for polynomial in gsl, however in  
> my opinion, this is acceptable, because the test condition is too  
> strict. With a slight change of the test coefficient, you will have a  
> similar failure also with the current version.

I've applied the patch and rewritten the angular calculation in terms
of R2 / Q3 so that it is based on the same variables as the condition
(R2 < Q3)--this should avoid the argument of acos being out of range.
Thanks for reporting the problem.

-- 
Brian Gough