Opened 5 years ago

Last modified 13 months ago

#16788 needs_work defect

segfault in numerical_integral() — at Version 7

Reported by: rws Owned by:
Priority: major Milestone: sage-8.4
Component: calculus Keywords: integral, numerics, crash
Cc: kcrisman Merged in:
Authors: Reviewers:
Report Upstream: N/A Work issues:
Branch: Commit:
Dependencies: Stopgaps:

Description (last modified by rws)

sage: numerical_integral(sin(x^2)/(x^2), 1, infinity, max_points=10^10)
---------------------------------------------------------------------------
SignalError                               Traceback (most recent call last)
<ipython-input-7-a164db8ca706> in <module>()
----> 1 numerical_integral(sin(x**Integer(2))/(x**Integer(2)), Integer(1), infinity, max_points=Integer(10)**Integer(10))

/home/ralf/sage/local/lib/python2.7/site-packages/sage/gsl/integration.so in sage.gsl.integration.numerical_integral (build/cythonized/sage/gsl/integration.c:2626)()

/home/ralf/sage/local/lib/python2.7/site-packages/sage/ext/c_lib.so in sage.ext.c_lib.sig_raise_exception (build/cythonized/sage/ext/c_lib.c:1208)()

SignalError: Segmentation fault

Better error messages would be needed here, as well.

The segfault happens on OpenSuSE 12.3 (both 6.2 and 6.3) and SMC (Ubuntu) but not on Mac OS X 10.7

Change History (7)

comment:1 Changed 5 years ago by rws

  • Description modified (diff)

comment:2 Changed 5 years ago by rws

  • Priority changed from major to critical

Segfaults are always critical.

comment:3 Changed 5 years ago by vbraun_spam

  • Milestone changed from sage-6.3 to sage-6.4

comment:4 Changed 5 years ago by kcrisman

  • Cc kcrisman added

What platform is this on? On Mac OS X 10.7 I get (0.285725372660691, 7.519910599701474e-05) for both Sage 5.12 and 6.2.

comment:5 Changed 5 years ago by rws

  • Description modified (diff)

comment:6 Changed 5 years ago by kcrisman

What platform is this on? On Mac OS X 10.7 I get (0.285725372660691, 7.519910599701474e-05) for both Sage 5.12 and 6.2.

This could be a good situation to ask on sage-devel; there are a wide variety of systems there and if everyone tried quickly you might get a lot of feedback. It's too bad the traceback is so short, that's not a lot of info to go on :( as usually I like tracking these things down.

On SMC, which is Ubuntu:

$ sage
┌────────────────────────────────────────────────────────────────────┐
│ Sage Version 6.3, Release Date: 2014-08-10                         │
│ Type "notebook()" for the browser-based notebook interface.        │
│ Type "help()" for help.                                            │
└────────────────────────────────────────────────────────────────────┘
sage: numerical_integral(sin(x^2)/(x^2), 1, infinity, max_points=10^10)
---------------------------------------------------------------------------
SignalError                               Traceback (most recent call last)
<ipython-input-1-a164db8ca706> in <module>()
----> 1 numerical_integral(sin(x**Integer(2))/(x**Integer(2)), Integer(1), infinity, max_points=Integer(10)**Integer(10))
 
/usr/local/sage/sage-6.3.beta6/local/lib/python2.7/site-packages/sage/gsl/integration.so in sage.gsl.integration.numerical_integral (build/cythonized/sage/gsl/integration.c:2626)
()
 
/usr/local/sage/sage-6.3.beta6/local/lib/python2.7/site-packages/sage/ext/c_lib.so in sage.ext.c_lib.sig_raise_exception (build/cythonized/sage/ext/c_lib.c:1208)()
 
SignalError: Segmentation fault
sage: 
Exiting Sage (CPU time 0m0.30s, Wall time 0m7.00s).
~$ uname -a
Linux compute20dc0 3.13.0-32-generic #57-Ubuntu SMP Tue Jul 15 03:51:08 UTC 2014 x86_64 x86_64 x86_64 GNU/Linux

Does yours only happen with Sage 6.3? I don't have that yet locally, though on 6.3.beta5 all is still well.

comment:7 Changed 5 years ago by rws

  • Description modified (diff)
Note: See TracTickets for help on using tickets.