Re: [sage-devel] linbox 64-bit charpoly

2016-10-03 Thread parisse
The probabilistic early termination does not take much time here, the charpoly stabilizes at about 85% of the primes required to reach the Hadamard bound. Testing with a few random matrices, I often get stabilization at about 80% (+/-10%), in this situation I think it's best to wait a little mo

Re: [sage-devel] linbox 64-bit charpoly

2016-09-28 Thread Clement Pernet
Hi, Update: I think I finally found the bug that led some rare computations to hang forever: givaro's random iterator was seed from the 6 digits of the current time microseconds, and could, with proba 10^-6 be seeded with 0, and the congurential generator would then always output 0, causing the

Re: [sage-devel] linbox 64-bit charpoly

2016-09-27 Thread parisse
Le mercredi 28 septembre 2016 03:13:11 UTC+2, Jonathan Bober a écrit : > > > Ah, yes, I'm wrong again, as the multimodular in Flint is pretty new. I > didn't look at what Sage has until now (flint 2.5.2, which looks likes it > uses a fairly simple O(n^4) algorithm). I had previously looked at t

Re: [sage-devel] linbox 64-bit charpoly

2016-09-27 Thread Jonathan Bober
On Tue, Sep 27, 2016 at 8:34 PM, 'Bill Hart' via sage-devel < sage-devel@googlegroups.com> wrote: > > > On Tuesday, 27 September 2016 20:53:28 UTC+2, Jonathan Bober wrote: >> >> On Tue, Sep 27, 2016 at 7:18 PM, 'Bill Hart' via sage-devel < >> sage-...@googlegroups.com> wrote: >> >>> I'm pretty sur

Re: [sage-devel] linbox 64-bit charpoly

2016-09-27 Thread Jonathan Bober
On Tue, Sep 27, 2016 at 8:02 PM, William Stein wrote: > On Tue, Sep 27, 2016 at 11:53 AM, Jonathan Bober > wrote: > > On Tue, Sep 27, 2016 at 7:18 PM, 'Bill Hart' via sage-devel > > wrote: > >> > >> I'm pretty sure the charpoly routine in Flint is much more recent that 2 > >> years. Are you ref

Re: [sage-devel] linbox 64-bit charpoly

2016-09-27 Thread 'Bill Hart' via sage-devel
On Tuesday, 27 September 2016 20:53:28 UTC+2, Jonathan Bober wrote: > > On Tue, Sep 27, 2016 at 7:18 PM, 'Bill Hart' via sage-devel < > sage-...@googlegroups.com > wrote: > >> I'm pretty sure the charpoly routine in Flint is much more recent that 2 >> years. Are you referring to a Sage implement

Re: [sage-devel] linbox 64-bit charpoly

2016-09-27 Thread William Stein
On Tue, Sep 27, 2016 at 11:53 AM, Jonathan Bober wrote: > On Tue, Sep 27, 2016 at 7:18 PM, 'Bill Hart' via sage-devel > wrote: >> >> I'm pretty sure the charpoly routine in Flint is much more recent that 2 >> years. Are you referring to a Sage implementation on top of Flint arithmetic >> or somet

Re: [sage-devel] linbox 64-bit charpoly

2016-09-27 Thread Jonathan Bober
On Tue, Sep 27, 2016 at 7:18 PM, 'Bill Hart' via sage-devel < sage-devel@googlegroups.com> wrote: > I'm pretty sure the charpoly routine in Flint is much more recent that 2 > years. Are you referring to a Sage implementation on top of Flint > arithmetic or something? > It is just a problem with S

Re: [sage-devel] linbox 64-bit charpoly

2016-09-27 Thread 'Bill Hart' via sage-devel
I'm pretty sure the charpoly routine in Flint is much more recent that 2 years. Are you referring to a Sage implementation on top of Flint arithmetic or something? The only timing that I can find right at the moment had us about 5x faster than Sage. It's not in a released version of Flint thoug

Re: [sage-devel] linbox 64-bit charpoly

2016-09-26 Thread Jonathan Bober
On Tue, Sep 27, 2016 at 4:18 AM, William Stein wrote: > On Mon, Sep 26, 2016 at 6:55 PM, Jonathan Bober wrote: > > On Mon, Sep 26, 2016 at 11:52 PM, William Stein > wrote: > >> > >> On Mon, Sep 26, 2016 at 3:27 PM, Jonathan Bober > wrote: > >> > In the matrix_integer_dense charpoly() function,

Re: [sage-devel] linbox 64-bit charpoly

2016-09-26 Thread William Stein
On Mon, Sep 26, 2016 at 6:55 PM, Jonathan Bober wrote: > On Mon, Sep 26, 2016 at 11:52 PM, William Stein wrote: >> >> On Mon, Sep 26, 2016 at 3:27 PM, Jonathan Bober wrote: >> > In the matrix_integer_dense charpoly() function, there is a note in the >> > docstring which says "Linbox charpoly dis

Re: [sage-devel] linbox 64-bit charpoly

2016-09-26 Thread Jonathan Bober
On Mon, Sep 26, 2016 at 11:52 PM, William Stein wrote: > On Mon, Sep 26, 2016 at 3:27 PM, Jonathan Bober wrote: > > In the matrix_integer_dense charpoly() function, there is a note in the > > docstring which says "Linbox charpoly disabled on 64-bit machines, since > it > > hangs in many cases."

Re: [sage-devel] linbox 64-bit charpoly

2016-09-26 Thread William Stein
On Mon, Sep 26, 2016 at 3:27 PM, Jonathan Bober wrote: > In the matrix_integer_dense charpoly() function, there is a note in the > docstring which says "Linbox charpoly disabled on 64-bit machines, since it > hangs in many cases." > > As far as I can tell, that is not true, in the sense that (1) I

[sage-devel] linbox 64-bit charpoly

2016-09-26 Thread Jonathan Bober
In the matrix_integer_dense charpoly() function, there is a note in the docstring which says "Linbox charpoly disabled on 64-bit machines, since it hangs in many cases." As far as I can tell, that is not true, in the sense that (1) I have 64-bit machines, and Linbox charpoly is not disabled, (2) c