Opened 6 years ago
Last modified 2 months ago
#17606 new defect
Calculation involving nintegral works with Sage 6.3 but not with Sage 6.4
Reported by: | alyasoff | Owned by: | |
---|---|---|---|
Priority: | major | Milestone: | sage-6.4 |
Component: | calculus | Keywords: | nintegral, integral |
Cc: | Merged in: | ||
Authors: | Andrew Lyasoff | Reviewers: | |
Report Upstream: | N/A | Work issues: | |
Branch: | Commit: | ||
Dependencies: | Stopgaps: |
Description (last modified by )
The following code:
var('u') assume(u>0) fu=1/8*(I*sqrt(2)*(erf(1/2*sqrt(2)*(I*pi + 2*u + 1)) - erf(1/2*sqrt(2)*(-I*pi + 2*u + 1)))*e^(u + 1/2) + sqrt(2)*(-I*erf(1/2*sqrt(2)*(I*pi - 2*u + 1)) + I*erf(1/2*sqrt(2)*(-I*pi - 2*u + 1)))*e^(u + 1/2) - (sqrt(2)*e^(1/8) - sqrt(2)*e^(2*u + 1/8))*(erf(1/4*sqrt(2)*(2*I*pi + 4*u + 1)) + erf(1/4*sqrt(2)*(-2*I*pi + 4*u + 1))) - (sqrt(2)*e^(1/8) - sqrt(2)*e^(2*u + 1/8))*(erf(1/4*sqrt(2)*(2*I*pi - 4*u + 1)) + erf(1/4*sqrt(2)*(-2*I*pi - 4*u + 1))))*e^(-u)/sqrt(pi) fu.nintegral(u,0,2)
returns:
(1.4275778409741315, 3.656408002774788e-14, 63, 0)
with Sage 6.3, but with Sage 6.4 returns:
Traceback (most recent call last): File "<stdin>", line 1, in <module> File "_sage_input_4.py", line 10, in <module> exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("ZnUubmludGVncmFsKHUsMCwyKQ=="),globals())+"\\n"); execfile(os.path.abspath("___code___.py")) File "", line 1, in <module> File "/tmp/tmpbwhja8/___code___.py", line 3, in <module> exec compile(u'fu.nintegral(u,_sage_const_0 ,_sage_const_2 ) File "", line 1, in <module> File "sage/symbolic/expression.pyx", line 10244, in sage.symbolic.expression.Expression.nintegral (build/cythonized/sage/symbolic/expression.cpp:51090) File "/usr/lib/sagemath/local/lib/python2.7/site-packages/sage/calculus/calculus.py", line 792, in nintegral raise TypeError(err) TypeError: ECL says: Error executing code in Maxima:
Attachments (2)
Change History (7)
Changed 6 years ago by
comment:1 Changed 6 years ago by
- Description modified (diff)
comment:2 Changed 6 years ago by
- Component changed from numerical to calculus
comment:3 Changed 15 months ago by
comment:4 Changed 2 months ago by
- Keywords integral added
comment:5 Changed 2 months ago by
The lines give me an error with 9.2b12 (on MacOS 10.15.2). Full error message is attached.
RuntimeError: ECL says: Error executing code in Maxima: quad_qags: Cannot numerically evaluate errexp1 at 1.0 During handling of the above exception, another exception occurred: ... TypeError: ECL says: Error executing code in Maxima: quad_qags: Cannot numerically evaluate errexp1 at 1.0. During handling of the above exception, another exception occurred: ... TypeError: ECL says: Error executing code in Maxima: quad_qags: Cannot numerically evaluate errexp1 at 1.0
No error with 9.1 (on my laptop and on CoCalc
).
Note: See
TracTickets for help on using
tickets.
works fine in 8.9.b7