On Tue, Jun 10, 2014 at 07:13:03AM +0200, Julien Puydt wrote: > Hi, > > Peter Bruin bissected the problem down to a commit (8cfb9c6b, > incorrect rounding in mulur), see comment #91 on the sage bug.
Hello Julien, First, I need to clarify that PARI never promised to return identical results on 32 and 64 bit platforms, and I do not see anything wrong with commit 8cfb9c6b. Secondly, on the face of it, the issue is catastrophic cancellation in computel.gp, in the loop for(y=1,ncoeff,Gvec[d+y]=1.*sum(j=1,d,subst(recG[j],n,y)*Gvec[d+y-j])); Cheers, -- Bill. <[email protected]> Imagine a large red swirl here. -- To UNSUBSCRIBE, email to [email protected] with a subject of "unsubscribe". Trouble? Contact [email protected]
