#10682: sum fails with lower bound != 0 or 1 (upgrade maxima to 5.26)
---------------------------------------------------------+------------------
Reporter: fmaltey | Owner:
burcin
Type: defect | Status:
needs_info
Priority: critical | Milestone:
sage-5.0
Component: symbolics | Keywords:
maxima 5.26.0 binomial sum
Work_issues: domain issue, error handling in nintegral | Upstream: N/A
Reviewer: Jean-Pierre Flori, | Author:
Dima Pasechnik
Merged: | Dependencies:
---------------------------------------------------------+------------------
Comment(by nbruin):
OK, this poking around in maxima's source is really quite unproductive,
but I think I have located the source of the change in behaviour. With
maxima 5.26, observe:
{{{
sage: x.nintegrate(x,0,1,1e-14)
(0.0, 0.0, 0, 6)
sage: from sage.libs.ecl import *
sage: ecl_eval("(setf f2cl-lib::*stop-signals-error-p* T)")
<ECL: T>
sage: x.nintegrate(x,0,1,1e-14)
ValueError: Maxima (via quadpack) cannot compute the integral to that
precision
}}}
whereas with maxima 5.23.2 we have the opposite:
{{{
sage: x.nintegrate(x,0,1,1e-14)
ValueError: Maxima (via quadpack) cannot compute the integral to that
precision
sage: from sage.libs.ecl import *
sage: ecl_eval("(setf f2cl-lib::*stop-signals-error-p* T)")
<ECL: T>
sage: x.nintegrate(x,0,1,1e-14)
(0.0, 0.0, 0, 6)
}}}
I think
[http://maxima.git.sourceforge.net/git/gitweb.cgi?p=maxima/maxima;a=commitdiff;h=9067349233cc217786983b07c23a2b78cf7d27bc
this commit] is responsible]:
{{{
(defun stop (&optional arg)
(when arg
(format cl::*error-output* "~A~%" arg))
- (unless *stop-signals-error-p*
+ (when *stop-signals-error-p*
(cerror "Continue anyway" "STOP reached")))
}}}
In later commits, the following is added:
{{{
+;;; Revision 1.116 2011/02/20 20:51:04 rtoy
+;;; Oops. STOP should signal an error if *STOP-SIGNALS-ERROR-P* is
+;;; non-NIL.
}}}
so I think the maxima 5.26 behaviour is as intended. I think you can give
Dima's patch a positive review now. Everything is explained and any
improvements should be handled on future tickets.
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/10682#comment:65>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
