[Maxima] [sage-devel] elliptic_e(0.5, 0.1) differs from Mathematica 7 by about 0.04%.
pip at iszf.irk.ru
Sun Aug 9 23:55:00 CDT 2009
> > 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
redefining ASDF:PERFORM :AROUND (#<STANDARD-CLASS ASDF:LOAD-OP>
#<STANDARD-CLASS ASDF:CL-SOURCE-FILE>) in
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);
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.
More information about the Maxima