Re: [sage-devel] bug in interface to pari's nffactor

2018-09-19 Thread Bruno Grenet
I've had time to find the origin of the bug: It appears to be a bug in Pari 2.11.0 (cf below). I've submitted a bug report. ? f = x^20 - 15*x^19 + 70989*x^18 - 1646113*x^17 + 3890074283*x^16 - 199035549796*x^15 + 203804256639644*x^14 - 10657741285726487*x^13 + 9779630086245476401*x^12 - 457685

Re: [sage-devel] bug in interface to pari's nffactor

2018-09-14 Thread Dima Pasechnik
One could use 'git bisect' to narrow the problem down to a particular commit. Just in case... On Fri, Sep 14, 2018 at 12:16 PM Bruno Grenet wrote: > > I've also tested on Ubuntu with gcc. Actually, I run new tests: The bug > appears between 8.4beta2 (where it's still ok) and 8.4beta3. I'll try to

Re: [sage-devel] bug in interface to pari's nffactor

2018-09-14 Thread Bruno Grenet
I've also tested on Ubuntu with gcc. Actually, I run new tests: The bug appears between 8.4beta2 (where it's still ok) and 8.4beta3. I'll try to find the the origin of the bug. Bruno Le 14/09/2018 à 10:38, Dima Pasechnik a écrit : It might be a platform dependent bug. I guess John is working

Re: [sage-devel] bug in interface to pari's nffactor

2018-09-14 Thread Dima Pasechnik
It might be a platform dependent bug. I guess John is working on some Ubunty machine, using gcc. Please provide compiler details you're using, too... On Fri, Sep 14, 2018 at 9:25 AM Bruno Grenet wrote: > > Which version of SageMath are you using? I am unable to reproduce the > bug neither on 8

Re: [sage-devel] bug in interface to pari's nffactor

2018-09-14 Thread Bruno Grenet
Which version of SageMath are you using? I am unable to reproduce the bug neither on 8.3beta7 nor on 8.4beta2. Bruno Example with 8.4beta2: sage: x = polygen(QQ) sage: f = x^20 - 15*x^19 + 70989*x^18 - 1646113*x^17 + 3890074283*x^16 - 199035549796*x^15 + 203804256639644*x^14 - 10 : 657741

[sage-devel] bug in interface to pari's nffactor

2018-09-13 Thread John Cremona
I have two polynomials of degree 20 defining the same number field (I obtained the second from the first using pari's polredbest() routine). I am able to use is_isomorphic() to find isomorphisms between them (there are 2) but embeddings() raises an error since roots() does: sage: x = polygen(QQ)