Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

an integral crashes Giac and Sage #31887

Open
mwageringel opened this issue May 31, 2021 · 1 comment
Open

an integral crashes Giac and Sage #31887

mwageringel opened this issue May 31, 2021 · 1 comment

Comments

@mwageringel
Copy link

As reported in this ask-sage question, the following integral crashes with Giac 1.7.0:

    sage: integrand=(x+(1+x)^(1/2))^(1/2)/(x^2+1)/(1+x)^(1/2)
    sage: integrate(integrand, x, algorithm="giac")
    Giac crashed -- automatically restarting.
    sage0*x

When using libgiac, Sage crashes as well:

    sage: from sage.libs.giac import libgiac
    sage: libgiac.integrate(integrand, x)

    [Process exited 139]

With Giac 1.6.0.47 which is in Sage 9.3, one gets some output:

sage: libgiac.integrate(integrand, x)
  ***   Warning: increasing stack size to 2048000.
Bad conditionned root j= 13 value 2.7174396231e+43+1.3659091552e+44*i ratio infinity mindist 5.4348792462e+43
  ***   Warning: increasing stack size to 2048000.
  ***   Warning: increasing stack size to 4096000.
  ***   Warning: increasing stack size to 2048000.
  ***   Warning: increasing stack size to 4096000.
  ***   Warning: increasing stack size to 8192000.
  ***   Warning: increasing stack size to 16384000.
Warning, integration of abs or sign assumes constant sign by intervals (correct if the argument is real):
Check [abs(t_nostep)]
  ***   Warning: increasing stack size to 2048000.
Bad conditionned root j= 13 value 2.7174396231e+43+1.3659091552e+44*i ratio infinity mindist 5.4348792462e+43
  ***   Warning: increasing stack size to 2048000.
  ***   Warning: increasing stack size to 2048000.
  ***   Warning: increasing stack size to 4096000.
  ***   Warning: increasing stack size to 8192000.
  ***   Warning: increasing stack size to 16384000.
Warning, integration of abs or sign assumes constant sign by intervals (correct if the argument is real):
Check [abs(t_nostep)]
---------------------------------------------------------------------------
...
RuntimeError: iquo Error: Bad Argument Type

CC: @sagetrac-parisse

Component: interfaces

Keywords: integrate, giac, integral

Issue created by migration from https://trac.sagemath.org/ticket/31887

@mwageringel mwageringel added this to the sage-9.4 milestone May 31, 2021
@edgarcosta
Copy link
Member

comment:1

Running on OSX, I also get similar issues due to conversion.
Full traceback: https://gist.github.com/edgarcosta/6ba59b0f53a94455f83c2166288cb479
Relevant bit:

TypeError: An error occurred running a Giac command:
INPUT:
sage2:=int(sage0,x):;
OUTPUT:
  ***   Warning: increasing stack size to 2048000.
Unable to convert to real -5.4348792461126750008981111610946592130*E43 Error: Bad Argument Value
Warning, integration of abs or sign assumes constant sign by intervals (correct if the argument is real):
Check [abs(t_nostep)]

Evaluation time: 2.12
"iquo Error: Bad Argument Type"

@mkoeppe mkoeppe modified the milestones: sage-9.4, sage-9.5 Aug 22, 2021
@mkoeppe mkoeppe modified the milestones: sage-9.5, sage-9.6 Dec 18, 2021
@mkoeppe mkoeppe modified the milestones: sage-9.6, sage-9.7 May 3, 2022
@mkoeppe mkoeppe modified the milestones: sage-9.7, sage-9.8 Sep 19, 2022
@mkoeppe mkoeppe removed this from the sage-9.8 milestone Jan 29, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

3 participants