collection of issues found in 1.9-07

Bugs

Modérateur : xcasadmin

Nasser
Messages : 13
Inscription : mer. juil. 11, 2018 11:14 pm

collection of issues found in 1.9-07

Message par Nasser » jeu. juin 16, 2022 4:56 am

Please find current collection of issues found. Some of these might be result of same problem or might be fixed in current source.

Code : Tout sélectionner

>giac
// Using locale /usr/local/share/locale/
// C.UTF-8
// /usr/local/share/locale/
// giac
// UTF-8
// Maximum number of parallel threads 24
Added 0 synonyms
Welcome to giac readline interface, version 1.9.0
(c) 2001,2021 B. Parisse & others
Homepage http://www-fourier.ujf-grenoble.fr/~parisse/giac.html
Released under the GPL license 3.0 or above
See http://www.gnu.org for license details
May contain BSD licensed software parts (lapack, atlas, tinymt)
-------------------------------------------------
Press CTRL and D simultaneously to finish session
Type ?commandname for help
0>> integrate(1/(e*x+d)^(5/2)/(-c*e^2*x^2+c*d^2)^(1/2),x)
Segmentation fault
>

Code : Tout sélectionner

0>> integrate((e*x+d)^(9/2)/(-c*e^2*x^2+c*d^2)^(3/2),x)
Segmentation fault



0>> integrate((e*x+d)^(7/2)/(-c*e^2*x^2+c*d^2)^(3/2),x)
Segmentation fault


0>> integrate(1/(e*x+d)^(5/2)/(-c*e^2*x^2+c*d^2)^(1/2),x)
Segmentation fault



0>> integrate(1/(c*e^2*x^2+2*c*d*e*x+c*d^2)^(3/2),x)
Segmentation fault


0>> integrate(1/(e*x+d)/(c*e^2*x^2+2*c*d*e*x+c*d^2)^(3/2),x)
Unable to build a single algebraic extension for simplifying.
Trying rational simplification only. This might return a wrong answer if simplifying 0/0!
Unable to build a single algebraic extension for simplifying.
Trying rational simplification only. This might return a wrong answer if simplifying 0/0!
Unable to build a single algebraic extension for simplifying.
Trying rational simplification only. This might return a wrong answer if simplifying 0/0!
Unable to build a single algebraic extension for simplifying.
Trying rational simplification only. This might return a wrong answer if simplifying 0/0!
Unable to build a single algebraic extension for simplifying.
Trying rational simplification only. This might return a wrong answer if simplifying 0/0!
Unable to build a single algebraic extension for simplifying.
Trying rational simplification only. This might return a wrong answer if simplifying 0/0!
Unable to build a single algebraic extension for simplifying.
Trying rational simplification only. This might return a wrong answer if simplifying 0/0!
Unable to build a single algebraic extension for simplifying.
Trying rational simplification only. This might return a wrong answer if simplifying 0/0!
Segmentation fault




0>> integrate((e*x+d)^(3/2)*(c*x^2+b*x+a)^p,x)
giac: poly.h:1965: void giac::Tegcdtoabcuv(const giac::tensor<T>&, const giac::tensor<T>&, const giac::tensor<T>&, giac::tensor<T>&, giac::tensor<T>&, giac::tensor<T>&, giac::tensor<T>&) [with T = giac::gen]: Assertion `temp.coord.empty()' failed.
Aborted



0>> integrate((g*x+f)/(e*x+d)^(5/2)/(-c*e^2*x^2-b*e^2*x-b*d*e+c*d^2)^(1/2),x)
Warning, need to choose a branch for the root of a polynomial with parameters. This might be wrong.
The choice was done assuming [b,c,d,exp(1),exp(2)]=[94,26,-27,20,-27]
Warning, need to choose a branch for the root of a polynomial with parameters. This might be wrong.
The choice was done assuming [b,c,d,exp(1),exp(2)]=[45,-9,71,12,-14]
Segmentation fault



0>> integrate(1/x/(c*x^2+b*x+a)^(3/2)/(f*x^2+e*x+d),x)
Segmentation fault



2>> integrate((D*x^3+C*x^2+B*x+A)/(a*x^4+b*x^3+c*x^2+b*x+a),x)
terminate called after throwing an instance of 'std::runtime_error'
  what():  Not invertible Error: Bad Argument Value
terminate called recursively
Aborted


0>> integrate(1/(b*x+a)/(d*x^2+c)^(1/2)/(f*x^2+e)^(1/2),x)
Segmentation fault

0>> integrate((c*x^2+b*x+a)^m*(g*x^3+f*x^2+e*x+d)^n*(-a*d+(a*e*n+b*d*m)*x+(2*a*f*n+b*e*m+b*e*n+2*c*d*m+a*f+b*e+c*d)*x^2+(3*a*g*n+b*f*m+2*b*f*n+2*c*e*m+c*e*n+2*a*g+2*b*f+2*c*e)*x^3+(b*g*m+3*b*g*n+2*c*f*m+2*c*f*n+3*b*g+3*c*f)*x^4+c*g*(4+2*m+3*n)*x^5)/x^2,x)
Simplification assuming d near 0
giac: poly.h:1965: void giac::Tegcdtoabcuv(const giac::tensor<T>&, const giac::tensor<T>&, const giac::tensor<T>&, giac::tensor<T>&, giac::tensor<T>&, giac::tensor<T>&, giac::tensor<T>&) [with T = giac::gen]: Assertion `temp.coord.empty()' failed.
Aborted
>


0>> integrate((d*i*x+c*i)^3*(A+B*log(e*(b*x+a)/(d*x+c)))/(b*g*x+a*g)^3,x)
terminate called after throwing an instance of 'std::runtime_error'
  what():  Not invertible Error: Bad Argument Value
terminate called recursively
Aborted


0>> integrate(x^3*log(c*(a+b/x^2)^p)/(e*x+d),x)
Simplification assuming a near 0
Simplification assuming d near 0
terminate called after throwing an instance of 'std::runtime_error'
  what():  Not invertible Error: Bad Argument Value
Aborted
>

0>> integrate(x^3*log(c*(a+b/x^3)^p)/(e*x+d),x)
Simplification assuming a near 0
Simplification assuming d near 0
terminate called after throwing an instance of 'std::runtime_error'
  what():  Not invertible Error: Bad Argument Value
terminate called recursively
Aborted

parisse
Messages : 5734
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: collection of issues found in 1.9-07

Message par parisse » jeu. juin 16, 2022 11:27 am

Hi,

Could you please just send the inputs that raise a bug? (Each line ending with a ;)
Something like

Code : Tout sélectionner

integrate((e*x+d)^(9/2)/(-c*e^2*x^2+c*d^2)^(3/2),x);
integrate(1/(c*e^2*x^2+2*c*d*e*x+c*d^2)^(3/2),x);
This is a lot easier for me, because I can copy/paste the whole input at once in a file and run icas on the file, check if the bugs are already fixed, or run the file from gdb for those that are not yet fixed.
Thanks!

lukamar
Messages : 331
Inscription : ven. juin 30, 2017 9:55 am
Localisation : Zagreb, Croatia

Re: collection of issues found in 1.9-07

Message par lukamar » sam. juin 18, 2022 1:22 pm

It seems that using e as a symbol in the expression e*x+d is the problem. In giac, it is the base of natural logarithm (a predefined constant). Nevertheless, using it shouldn't crash icas.

parisse
Messages : 5734
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: collection of issues found in 1.9-07

Message par parisse » sam. juin 18, 2022 6:44 pm

parse_e(1) will prevent parsing e to exp(1). From C++ parse_e(true,&ct);

parisse
Messages : 5734
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: collection of issues found in 1.9-07

Message par parisse » dim. juin 19, 2022 6:07 am

https://www-fourier.univ-grenoble-alpes ... .0.tar.bz2 was updated with fixes. Some integrals still hang, but they should not segfault anymore...

Répondre