Re: Test suite?

["Frank Reininghaus" <frank78ac@..............>]

Hi Liam,

thanks for the report. I cannot reproduce this problem. Actually, I wrote
the relevant code (it's just a few lines in gsl-1.11/poly/eval.c and
gsl-1.11/poly/gsl_poly.h) and the test myself and I have no idea what could
go wrong :-( Maybe it would be helpful if you could provide the following
information:

- Which Linux distribution do you use?
- Which compiler do you use (probably gcc?) and which version (output of
'gcc --version' if you're on gcc)?
- What kind of CPU do you have? Is it running in 32- or 64-bit mode?
- Did you change the optimisation options, or do you use the default?
- If you repeat the test, does it always report the same wrong results
(-1.6845078800000004 observed vs ... and -0.30943988
observed vs ...)?

Thanks in advance! If you're subscribed to bug-gsl, it might be best to move
the discussion to that mailing list.

My first idea would be to replace the code in the loop in
gsl_complex_poly_complex_eval () by 'ans = gsl_complex_add (c[i-1],
gsl_complex_mul (x, ans))' and see if the problem persists (although I still
think that the separate calculation of real and imaginary part which I found
to be faster is equivalent). I could also write a program that tests more
polynomials and ask Liam (or somebody else who can reproduce the problem) to
run it to see if I can find a pattern.

I'm sorry for the trouble caused by my code.

Regards
Frank


2008/9/15 Liam Healy <lnp@......................>

> Thanks Frank.
>
> I did make check and it ended with
>
> 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
> ===================
> 1 of 1 tests failed
> ===================
> make[2]: *** [check-TESTS] Error 1
> make[2]: Leaving directory `/home/liam/mathematics/gsl-1.11/poly'
> make[1]: *** [check-am] Error 2
> make[1]: Leaving directory `/home/liam/mathematics/gsl-1.11/poly'
> make: *** [check-recursive] Error 1
>
> It looks like something's wrong either in the test or in the definition.
>