On Thu, Jun 23, 2022 at 09:00:40PM -0700, 'Nasser M. Abbasi' via FriCAS -
computer algebra system wrote:
> Thanks for version 1.3.8.
>
> I've only did partial test (less than 50%) but regression report generated
> so far shows there are some integrals no longer evaluate in 1.3.8 but did
> in 1.3.7. This could be due to bug fix. But thought to let you know just in
> case.
>
> example 1
> ========
> using 1.3.7
> integrate((F^((-a*x+1)^(1/2)/(a*x+1)^(1/2)))^n/(-a^2*x^2+1),x);
> gives
> ((-1)*Ei((n*log(F)*((-1)*a*x+1)^(1/2)*(a*x+1)^(1/2))/(a*x+1)))/a
> But does not evaluate (returns the input) in 1.3.8
>
> example 2
> ========
> using 1.3.7
> integrate(F^(3*(-a*x+1)^(1/2)/(a*x+1)^(1/2))/(-a^2*x^2+1),x)
> gives
> ((-1)*Ei((3*log(F)*((-1)*a*x+1)^(1/2)*(a*x+1)^(1/2))/(a*x+1)))/a
> But does not evaluate (returns the input) in 1.3.8
>
>
> example 3
> ========
> using 1.3.7
> integrate(F^(2*(-a*x+1)^(1/2)/(a*x+1)^(1/2))/(-a^2*x^2+1),x)
> gives
> ((-1)*Ei((2*log(F)*((-1)*a*x+1)^(1/2)*(a*x+1)^(1/2))/(a*x+1)))/a
> But does not evaluate (returns the input) in 1.3.8
>
> example 4
> ========
> using 1.3.7
> integrate(F^(2*(-a*x+1)^(1/2)/(a*x+1)^(1/2))/(-a^2*x^2+1),x)
> gives
> ((-1)*Ei((2*log(F)*((-1)*a*x+1)^(1/2)*(a*x+1)^(1/2))/(a*x+1)))/a
> But does not evaluate (returns the input) in 1.3.8
>
Thanks for the report. In 1.3.7 integrals above were done
by pattern matching integrator. In 1.3.8 pattern matching
integrator is removed, algorithmic integrater is supposed
to handle them, but fails.
> example 5
> ========
> using 1.3.7
> integrate(sec(d*x+c)*(a+b*sin(d*x+c))^(5/2),x)
> gives this long output
<snip>
> But does not evaluate in 1.3.8
This looks like an issue with new code to generate elliptic
functions, I will investigate.
> The following all timeout now when they did not on 1.3.7. I need to re-run
> these again to make sure the timeout is due to version number change and
> not anything else
>
> integrate((j*x^4+i*x^3+h*x^2+g*x+f)/(c*x^2+b*x+a)^(5/2),x,)
> integrate((j*x^4+i*x^3+h*x^2+g*x+f)/(-c*x^2+b*x+a)^(5/2),x)
> integrate(x*(e*x^3+d)/(c*x^6+b*x^3+a),x)
> integrate((e*x^3+d)/(c*x^6+b*x^3+a),x)
> integrate(x^2*(e*x^4+d)/(c*x^8+b*x^4+a),x)
> integrate((e*x^4+d)/(c*x^8+b*x^4+a),x)
> integrate((e*x)^m*(a+b*log(c*log(d*x)^p)),x)
The first six are long running, 137.66 sec, 135.69, 213.04 sec, 110.06 sec,
254.77 and 81.66 sec respectively on my machine. The last
one needs 0.22 sec on my machine.
It is possible that 1.3.8 is a bit slower on them than 1.3.7
and there is substantial variation from run to run, so for the
first six it is possible that 1.3.8 did not do them within limit.
Or some timing oddity, expecially for the last one.
--
Waldek Hebisch
--
You received this message because you are subscribed to the Google Groups
"FriCAS - computer algebra system" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/fricas-devel/20220624081839.GA18167%40fricas.math.uni.wroc.pl.
