Yet another Maple bug (asympt, vv 11/10) discovered at Cyber Tester's CAS exact challenge ascent
Vladimir Bondarenko
[ ... ]
"Actually, the challenges are a tiny fraction of what we hope
could become the definitive, world's largest benchmark test
for new powerful computer algebra systems to be created within
the next 20-30 years."
[ ... ]
We see yet another boon of immediate nature brought by the
Cyber Tester's exact challenge top climbing. (An in-passing
remark for the advanced readers, don't get disappointed too
early, right now we are still at the piedmont... thanks for your
patience.)
While trying to beat our challenges, we the customers now and
then encounter new bugs thus hopefully helping the computer
algebra systems manufacturers to improve their systems.
On Apr 9, 10:40 am, "G. A. Edgar" <e...@math.ohio-state.edu.invalid>
writes:
http://groups.google.com/group/sci.math.symbolic/msg/e4ff1eb81c8abdf7?hl=en&
GAE> restart;
GAE> (t/(1-t))^a/t:
GAE> I1 := Int(%,t=0..1/2):
GAE> with(student):
GAE> changevar(u=t/(1-t),I1,u):
GAE> I2 := simplify(%):
GAE> value(%):
GAE> convert(%,Psi);
GAE> asympt(%,a);
GAE> 1/2*Psi(1/2+1/2*a)-1/2*Psi(1/2*a)
GAE> 1/2/a+1/4/a^2-1/8/a^4+O(1/a^6)
Maple 11 and Maple 10 return
Pi*csc(Pi*a)-1/(a-1)*(-a+1)/a-1/(a-1)*(-a
+1)*(1/2*Psi(1/2-1/2*a)-1/2*Psi(-1/2*a))
Error, (in asympt) unable to compute series
Maple crisis, alive & kickin'!
http://maple.bug-list.org/maple-crisis.php
Best wishes,
Vladimir Bondarenko
VM and GEMM architect
Co-founder, CEO, Mathematical Director
http://www.cybertester.com/ Cyber Tester, LLC
http://maple.bug-list.org/ Maple Bugs Encyclopaedia
http://www.CAS-testing.org/ CAS Testing
Axel Vogt
Hm ... may be a problem with the student package?
For me this works:
restart; interface(version);
(t/(1-t))^a/t:
I1 := Int(%,t=0..1/2):
student[changevar](u=t/(1-t),I1,u):
I2 := simplify(%):
value(%):
#% assuming 0 < a;
v:=convert(%,Psi);
asympt(%,a);
Classic Worksheet Interface, Maple 10.06, Windows, Oct 2 2006 Build ID 255401
v := 1/2 Psi(1/2 + a/2) - 1/2 Psi(a/2)
1 1 1 1
--- + ---- - ---- + O(----)
2 a 2 4 6
4 a 8 a a
Vladimir Bondarenko
AV> For me this works
Gee!
Congratulations!
Together, we discovered an INTERMITTENT bug in Maple 11/Maple 10.
We again encountered a non-deterministic behavior in Maple kernel
I pointed out in the Maple crisis link ... but this bug was absent
in earlier versions...
The Maple crisis is expanding...
http://maple.bug-list.org/maple-crisis.php
So:
..................................................................................
> restart; interface(version);(t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
`Classic Worksheet Interface, Maple 10.06, Windows, Oct 2 2006 Build
ID 255401`
Error, (in asympt) unable to compute series
> restart; interface(version);(t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
`Classic Worksheet Interface, Maple 10.06, Windows, Oct 2 2006 Build
ID 255401`
1/2/a+1/4/a^2-1/8/a^4+O(1/a^6)
> restart; interface(version);(t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
`Classic Worksheet Interface, Maple 10.06, Windows, Oct 2 2006 Build
ID 255401`
Error, (in asympt) unable to compute series
..................................................................................
> restart; interface(version);(t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
`Classic Worksheet Interface, Maple 11.00, Windows, Feb 16 2007 Build
ID 277223`
Error, (in asympt) unable to compute series
> restart; interface(version);(t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
`Classic Worksheet Interface, Maple 11.00, Windows, Feb 16 2007 Build
ID 277223`
Error, (in asympt) unable to compute series
> restart; interface(version);(t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
`Classic Worksheet Interface, Maple 11.00, Windows, Feb 16 2007 Build
ID 277223`
Error, (in asympt) unable to compute series
> restart; interface(version);(t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
`Classic Worksheet Interface, Maple 11.00, Windows, Feb 16 2007 Build
ID 277223`
Error, (in asympt) unable to compute series
> restart; interface(version);(t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
`Classic Worksheet Interface, Maple 11.00, Windows, Feb 16 2007 Build
ID 277223`
Error, (in asympt) unable to compute series
> restart; interface(version);(t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
`Classic Worksheet Interface, Maple 11.00, Windows, Feb 16 2007 Build
ID 277223`
1/2/a+1/4/a^2-1/8/a^4+O(1/a^6)
..................................................................................
On Apr 9, 12:52 pm, Axel Vogt <&...@axelvogt.de> wrote:
> Vladimir Bondarenko wrote:
> >http://groups.google.com/group/sci.math.symbolic/msg/aaa2ae83b7f110da...
>
> > [ ... ]
>
> > "Actually, the challenges are a tiny fraction of what we hope
> > could become the definitive, world's largest benchmark test
> > for new powerful computer algebra systems to be created within
> > the next 20-30 years."
>
> > [ ... ]
>
> > We see yet another boon of immediate nature brought by the
> > Cyber Tester's exact challenge top climbing. (An in-passing
> > remark for the advanced readers, don't get disappointed too
> > early, right now we are still at the piedmont... thanks for your
> > patience.)
>
> > While trying to beat our challenges, we the customers now and
> > then encounter new bugs thus hopefully helping the computer
> > algebra systems manufacturers to improve their systems.
>
> > On Apr 9, 10:40 am, "G. A. Edgar" <e...@math.ohio-state.edu.invalid>
> > writes:
>
> >http://groups.google.com/group/sci.math.symbolic/msg/e4ff1eb81c8abdf7...
> 4 a 8 a a- Hide quoted text -
>
> - Show quoted text -
Axel Vogt
Vladimir Bondarenko
AV> well ... for me it works.
Please inject the following input 100 or better 1000 times:
restart; (t/(1-t))^a/t:
I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
You always (100/1000 times) see
1/2/a+1/4/a^2-1/8/a^4+O(1/a^6)
and never
Error, (in asympt) unable to compute series
?
If this is the case, what is your CPU and operating system?
Could anyone else save AV reproduce this experiment and share
the results?
Axel Vogt
> On Apr 9, 1:23 pm, Axel Vogt <&...@axelvogt.de> wrote:
>
> AV> well ... for me it works.
>
> Please inject the following input 100 or better 1000 times:
>
> restart; (t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
>
> You always (100/1000 times) see
>
> 1/2/a+1/4/a^2-1/8/a^4+O(1/a^6)
>
> and never
>
> Error, (in asympt) unable to compute series
>
> ?
Seems I missed your point, sorry ... so explicitly: yes, it
is switching in the results and occasionally refuses to do
the series.
Eric Gisse
wrote:
Buh?
That's...interesting.
>
> restart; (t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
1 1 1 /1 \
--- + ---- - ---- + O|--|
2 a 2 4 | 6|
4 a 8 a \a /
>
> restart; (t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
1 1 1 /1 \
--- + ---- - ---- + O|--|
2 a 2 4 | 6|
4 a 8 a \a /
>
> restart; (t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
1 1 1 /1 \
--- + ---- - ---- + O|--|
2 a 2 4 | 6|
4 a 8 a \a /
>
> restart; (t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
Error, (in asympt) unable to compute series
>
> restart; (t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
Error, (in asympt) unable to compute series
>
> restart; (t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
Error, (in asympt) unable to compute series
>
> restart; (t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
Error, (in asympt) unable to compute series
>
> restart; (t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
Error, (in asympt) unable to compute series
>
> restart; (t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
Error, (in asympt) unable to compute series
>
> restart; (t/(1-t))^a/t:
> I1 := Int(%,t=0..1/2):student[changevar](u=t/(1-t),I1,u):
> I2 := simplify(%):value(%):v:=convert(%,Psi): asympt(%,a);
Error, (in asympt) unable to compute series