Karim Belabas on Sun, 02 Jun 2013 15:14:01 +0200


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

Re: weird error in thueinit


* Max Alekseyev [2013-06-02 13:31]:
> In my lengthy computations involving many thueinit() / thue() calls, I
> suddenly get the following error:
> 
>   *** thueinit: incorrect priority in nfroots: variable x >= x
> 
> in the call "thueinit(p)" with apperently innocent p = -29*x^4 +
> 184*x^3 - 294*x^2 + 184*x - 29.
> Right after this, calling thueinit(p) causes the same error, while
> explicit calling thueinit(-29*x^4 + 184*x^3 - 294*x^2 + 184*x - 29)
> works fine.

This is apparently generated when calling 
  thueinit(t_POL with automorphisms)
  -> bnfinit
  -> nfroots (to compute automorphisms in order to speed up bnfinit)

The "incorrect priority" is either 
1) a programming mistake on your side. Try to 

  - replace thueinit(p) by thueinit(simplify(p))

  - add write("some_logfile", getrand()) statements at various places
    Maybe resetting the random seed to this know value will help
    reproduce the problem.

2) the result of internal data corruption somewhere, which will be very
hard to track down unless you can supply a script that generates the
problem from scratch. Can you experiment and try to generate such a script ?

Cheers,

    K.B.
--
Karim Belabas, IMB (UMR 5251)  Tel: (+33) (0)5 40 00 26 17
Universite Bordeaux 1          Fax: (+33) (0)5 40 00 69 50
351, cours de la Liberation    http://www.math.u-bordeaux1.fr/~kbelabas/
F-33405 Talence (France)       http://pari.math.u-bordeaux1.fr/  [PARI/GP]
`