Olivier Ramare on Sun, 07 Oct 2001 19:54:21 +0200


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Un bug avec polsturm


Bonjour,

  un petit probleme avec polsturm ?
  Voici la session en entier :

(* le polynome gg est en fait: X*intformal((4*X*(1-X))^14) + 1  *)
---------------------------------------------------------------------------------------------------------------------------------
Appele avec : /usr/local/bin/gp -s 10000000 -p 500000 -emacs

                       GP/PARI CALCULATOR Version 2.1.0 (released)
                     i686 running linux (ix86 kernel) 32-bit version
                     (readline v4.1 enabled, extended help available)

                            Copyright (C) 2000 The PARI Group

PARI/GP is free software, covered by the GNU General Public License, and
comes WITHOUT 
ANY WARRANTY WHATSOEVER.

Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.

   realprecision = 28 significant digits
   seriesprecision = 16 significant terms
   format = g0.28

parisize = 10000000, primelimit = 500000
? gg=268435456/29*X^30 - 134217728*X^29 + 24427626496/27*X^28 -
3758096384*X^27 + 268703891456/25*X^26 - 67175972864/3*X^25 +
806111674368/23*X^24 - 41875931136*X^23 + 38386270208*X^22 -
134351945728/5*X^21 + 268703891456/19*X^20 - 48855252992/9*X^19 +
24427626496/17*X^18 - 234881024*X^17 + 268435456/15*X^16 + 1
%1 = 268435456/29*X^30 - 134217728*X^29 + 24427626496/27*X^28 -
3758096384*X^27 + 268703891456/25*X^26 - 67175972864/3*X^25 +
806111674368/23*X^24 - 41875931136*X^23 + 38386270208*X^22 -
134351945728/5*X^21 + 268703891456/19*X^20 - 48855252992/9*X^19 +
24427626496/17*X^18 - 234881024*X^17 + 268435456/15*X^16 + 1
? polsturm(gg,0,0.5)
  ***   bug in GP (Segmentation Fault), please report
? 
---------------------------------------------------------------------------------------------------------------------------------
  Bon courage donc !
  Amities,
           Olivier