On Thu, Sep 18, 2008 at 3:36 AM, Frank Reininghaus
<[EMAIL PROTECTED]> wrote:
> Hi Liam,
>
> thanks for your reply!
>
> 2008/9/18 Liam Healy <[EMAIL PROTECTED]>
>
>> I modified the function in poly/eval.c almost like this, but using z
>> instead of x as that looks more correct and it is what is in the
>> commented-out portion
>>
>
> It would be great if you could check this (adding #include
> <gsl/gsl_complex_math.h> as Steven suggested) and tell us if that fixes the
> problem. You might have to make the same change in the file
> gsl-1.11/poly/gsl_poly.h for the inlined version of the function.
>
> Thanks a lot!
> Frank
I made all these changes and reran make, make check and got the same answer
make[2]: Entering directory `/home/liam/mathematics/gsl-1.11/poly'
FAIL: y.real, gsl_complex_poly_complex_eval ({-2.31 + 0.44i, 4.21 -
3.19i, 0.93 + 1.04i, -0.42 + 0.68i}, 0.49 + 0.95i)
(-1.6845078800000004 observed vs 1.82462012000000007 expected) [9]
FAIL: y.imag, gsl_complex_poly_complex_eval ({-2.31 + 0.44i, 4.21 -
3.19i, 0.93 + 1.04i, -0.42 + 0.68i}, 0.49 + 0.95i) (-0.30943988
observed vs 2.30389411999999982 expected) [10]
FAIL: test
Liam
_______________________________________________
Bug-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/bug-gsl
