[sage-devel] elliptic_e(0.5, 0.1) differs from Mathematica 7 by about 0.04%.



> > Dr. David Kirkby wrote:
> > > I've built Sage on a sun4u SPARC machine, using Maxima 5.19.0 and ECL
> > > 9.8.1
> > >
> > > Sage fails about 15 tests, one of which is below. I computed the result
> > > in Mathematica too, and find the result on the SPARC does differ from
> > > what Sage expects, and the result Sage expects is much closer to the
> > > result from Mathematica.
> > >
> > >      sage: elliptic_eu (0.5, 0.1)
> > > Expected:
> > >      0.496054551287
> > > Got:
> > >      0.495848403419
> >
> > sorry, I posted the wrong result - there is this failure, but the one I
> > meant to post was:
> >
> >      sage: elliptic_e(0.5, 0.1)
> > Expected:
> >       0.498011394499
> > Got:
> >       0.497801100392
> >
> > Mathematica gives: 0.49801139449883153311546104061744810584963105068054
>
> [va at colombo py-q]$ rmaxima
> Maxima 5.18.1 http://maxima.sourceforge.net
> Using Lisp SBCL 1.0.29
> Distributed under the GNU Public License. See the file COPYING.
> Dedicated to the memory of William Schelter.
> The function bug_report() provides bug reporting information.
> (%i1) elliptic_e(0.5, 0.1);
> (%o1)                          .4980113944988315
>
Sorry last was for the previous version of Maxima

[va at colombo maxima-5.19.0]$ rlwrap ./maxima-local
STYLE-WARNING:
   redefining ASDF:PERFORM :AROUND (#<STANDARD-CLASS ASDF:LOAD-OP>
                                    #<STANDARD-CLASS ASDF:CL-SOURCE-FILE>) in 
DEFMETHOD
Maxima 5.19.0 http://maxima.sourceforge.net
Using Lisp SBCL 1.0.29
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) elliptic_e(0.5, 0.1);
(%o1)                          .4980113944988315

I wonder how sage people resist to use the high perfomance lisps like cmucl or 
sbcl (which has windows version as well) to compile maxima. 

best luck
Valery