Re: [sage-devel] Modular forms slowdown

2013-01-08 Thread Jeroen Demeyer
On 2013-01-07 20:50, Jeroen Demeyer wrote: > On 2013-01-07 12:24, Jeroen Demeyer wrote: >> I am running the command >> >> sage: time B = >> ModularFormsRing(Gamma0(3)).modular_forms_of_weight(80).cuspidal_submodule().q_expansion_basis(prec=30) >> >> sage-5.4: 17.94 s >> sage-5.5: 56.28 s > > I can

Re: [sage-devel] Modular forms slowdown

2013-01-08 Thread Jeroen Demeyer
On 2013-01-08 09:11, Jeroen Demeyer wrote: > On 2013-01-07 23:43, Volker Braun wrote: >> But the command is mostly linear algebra over QQ. Its possible that the >> old GCC does something very stupid, but I doubt that you can speed up >> integer computations by a factor of 5x just by optimizing code

Re: [sage-devel] Modular forms slowdown

2013-01-08 Thread Jeroen Demeyer
On 2013-01-07 23:43, Volker Braun wrote: > But the command is mostly linear algebra over QQ. Its possible that the > old GCC does something very stupid, but I doubt that you can speed up > integer computations by a factor of 5x just by optimizing code paths. So > it pretty much must be the blas ins

Re: [sage-devel] Modular forms slowdown

2013-01-07 Thread Volker Braun
But the command is mostly linear algebra over QQ. Its possible that the old GCC does something very stupid, but I doubt that you can speed up integer computations by a factor of 5x just by optimizing code paths. So it pretty much must be the blas inside linbox. On Monday, January 7, 2013 9:1

Re: [sage-devel] Modular forms slowdown

2013-01-07 Thread Jeroen Demeyer
On 2013-01-07 21:59, Volker Braun wrote: > Its pretty clear that the newer GCC correctly vectorizes ATLAS for the > CPU whereas the old one did not. I doubt it's ATLAS. I didn't find any particular difference in the ATLAS install logs. It's curious to see that ATLAS somehow brute-forces CPU flags,

Re: [sage-devel] Modular forms slowdown

2013-01-07 Thread Volker Braun
Its pretty clear that the newer GCC correctly vectorizes ATLAS for the CPU whereas the old one did not. On Monday, January 7, 2013 8:50:14 PM UTC, Snark wrote: > > >> > ModularFormsRing(Gamma0(3)).modular_forms_of_weight(80).cuspidal_submodule().q_expansion_basis(prec=30) > > > I'm still impre

Re: [sage-devel] Modular forms slowdown

2013-01-07 Thread Julien Puydt
Le 07/01/2013 20:50, Jeroen Demeyer a écrit : On 2013-01-07 12:24, Jeroen Demeyer wrote: I am running the command sage: time B = ModularFormsRing(Gamma0(3)).modular_forms_of_weight(80).cuspidal_submodule().q_expansion_basis(prec=30) sage-5.4: 17.94 s sage-5.5: 56.28 s I can confirm that buil

Re: [sage-devel] Modular forms slowdown

2013-01-07 Thread Jeroen Demeyer
On 2013-01-07 12:24, Jeroen Demeyer wrote: > I am running the command > > sage: time B = > ModularFormsRing(Gamma0(3)).modular_forms_of_weight(80).cuspidal_submodule().q_expansion_basis(prec=30) > > sage-5.4: 17.94 s > sage-5.5: 56.28 s I can confirm that building Sage with GCC-4.2.4 gives the s

Re: [sage-devel] Modular forms slowdown

2013-01-07 Thread Jeroen Demeyer
On 2013-01-07 15:27, John Cremona wrote: > On 7 January 2013 14:27, Jeroen Demeyer wrote: >> On 2013-01-07 12:24, Jeroen Demeyer wrote: >>> I am running the command >>> >>> sage: time B = >>> ModularFormsRing(Gamma0(3)).modular_forms_of_weight(80).cuspidal_submodule().q_expansion_basis(prec=30) >>

Re: [sage-devel] Modular forms slowdown

2013-01-07 Thread John Cremona
On 7 January 2013 14:27, Jeroen Demeyer wrote: > On 2013-01-07 12:24, Jeroen Demeyer wrote: >> I am running the command >> >> sage: time B = >> ModularFormsRing(Gamma0(3)).modular_forms_of_weight(80).cuspidal_submodule().q_expansion_basis(prec=30) >> >> sage-5.4: 17.94 s >> sage-5.5: 56.28 s > > F

Re: [sage-devel] Modular forms slowdown

2013-01-07 Thread Jeroen Demeyer
On 2013-01-07 12:24, Jeroen Demeyer wrote: > I am running the command > > sage: time B = > ModularFormsRing(Gamma0(3)).modular_forms_of_weight(80).cuspidal_submodule().q_expansion_basis(prec=30) > > sage-5.4: 17.94 s > sage-5.5: 56.28 s False alarm, the cause seems to be differing GCC versions o

Re: [sage-devel] Modular forms slowdown

2013-01-07 Thread Jeroen Demeyer
On 2013-01-07 15:14, David Loeffler wrote: > That's very strange! The time-critical part of that computation is > almost entirely in doing linear algebra with the coefficients of power > series -- it bypasses most of the modular symbols code, which would be > very inefficient in this example since

Re: [sage-devel] Modular forms slowdown

2013-01-07 Thread David Loeffler
That's very strange! The time-critical part of that computation is almost entirely in doing linear algebra with the coefficients of power series -- it bypasses most of the modular symbols code, which would be very inefficient in this example since the level is small and the weight large -- and oddl

Re: [sage-devel] Modular forms slowdown

2013-01-07 Thread Jeroen Demeyer
Actually, the slowdown is in 5.4 -> 5.4.1, which narrows it down quite a lot: http://boxen.math.washington.edu/home/release/sage-5.4.1/tickets.html On 2013-01-07 12:24, Jeroen Demeyer wrote: > I am running the command > > sage: time B = > ModularFormsRing(Gamma0(3)).modular_forms_of_weight(80).c

[sage-devel] Modular forms slowdown

2013-01-07 Thread Jeroen Demeyer
I am running the command sage: time B = ModularFormsRing(Gamma0(3)).modular_forms_of_weight(80).cuspidal_submodule().q_expansion_basis(prec=30) sage-5.4: 17.94 s sage-5.5: 56.28 s -- You received this message because you are subscribed to the Google Groups "sage-devel" group. To post to this g