[sage-support] Re: Why this numerical integral bombs?
On Apr 26, 1:36 am, mabshoff michael.absh...@mathematik.uni- dortmund.de wrote: I do not see how that rules out clisp 2.47 - or am I just not getting what you are driving at :) I meant if I can recreate the problem from 3.4 source as well as 3.4.1 source then doesn't that mean I can recreate the problem in both version of clisp? Hence clisp is innocent? cs --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---
[sage-support] Re: Why this numerical integral bombs?
On Apr 26, 11:05 pm, Chris Seberino cseber...@gmail.com wrote: On Apr 26, 1:36 am, mabshoff michael.absh...@mathematik.uni- dortmund.de wrote: I do not see how that rules out clisp 2.47 - or am I just not getting what you are driving at :) I meant if I can recreate the problem from 3.4 source as well as 3.4.1 source then doesn't that mean I can recreate the problem in both version of clisp? Hence clisp is innocent? No, my guess is that either clisp release triggers the same problem. Any recent clisp source release is incredibly easy to miscompile and we are building it with -O0 with any compiler to maximize our chances to make it work. But no clisp release since 2.41 has ever even build on the build farm or SkyNet, much less various OSX releases or Solaris (at least not with anything beyond gcc 4.x). That is the reason we are switching and numerous other projects have prevented me from doing the clisp - ecl switch, but finally its time has come since that is the only way to support Solaris properly. cs Cheers, Michael --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---
[sage-support] Re: Why this numerical integral bombs?
mabshoff wrote: On Apr 25, 12:23 am, Paul Zimmermann paul.zimmerm...@loria.fr wrote: Michael, Hi Paul, And Sage 3.4: | Sage Version 3.4, Release Date: 2009-03-11 | | Type notebook() for the GUI, and license() for information.| -- sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 So I cannot reproduce this. What platform are you on and are you using a binary, build from source, etc? the example that failed has 'exp' instead of 'cos'. However with Sage 3.4 it works for me: -- | Sage Version 3.4, Release Date: 2009-03-11 | | Type notebook() for the GUI, and license() for information.| -- sage: numerical_integral(sin(pi*exp(x/2)),0,2)[0] -0.43734547482524966 Paul Zimmermann I ran this on 3.4, 3.4.1 and a couple other rc releases in between: sage: for i in (0..200): : numerical_integral(sin(pi*exp(x/2)),0,2)[0] : and I could not get a failure. Note that the most likely culprit here is clisp IMHO. Just today I found out that clisp 2.47 doesn't build on the latest Mandriva release in 32 and 64 bit flavors. Oh well ... We should have an Ubuntu 9.04 test image in the not too distant future, so maybe we will hit it there. I just upgraded to Jaunty over the weekend, and then compiled 3.4.1 from scratch. Both expressions work for me: $ cat /etc/lsb-release DISTRIB_ID=Ubuntu DISTRIB_RELEASE=9.04 DISTRIB_CODENAME=jaunty DISTRIB_DESCRIPTION=Ubuntu 9.04 $ sage -- | Sage Version 3.4.1, Release Date: 2009-04-21 | | Type notebook() for the GUI, and license() for information.| -- sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*exp(x/2)),0,2)[0] -0.4373454748252496 sage: numerical_integral(sin(pi*exp(x/2)),0,2)[0] -0.4373454748252496 sage: numerical_integral(sin(pi*exp(x/2)),0,2)[0] -0.4373454748252496 sage: -Jason --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---
[sage-support] Re: Why this numerical integral bombs?
On Apr 27, 2:33 am, Jason Grout jason-s...@creativetrax.com wrote: I just upgraded to Jaunty over the weekend, and then compiled 3.4.1 from scratch. Both expressions work for me: Hmm. Does Sage leave a log build regarding the compilation? Perhaps my build had some errors yours doesn't have I should look for? cs --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---
[sage-support] Re: Why this numerical integral bombs?
Chris Seberino wrote: On Apr 27, 2:33 am, Jason Grout jason-s...@creativetrax.com wrote: I just upgraded to Jaunty over the weekend, and then compiled 3.4.1 from scratch. Both expressions work for me: Hmm. Does Sage leave a log build regarding the compilation? Perhaps my build had some errors yours doesn't have I should look for? Yes, in install.log in your root sage directory, if I recall correctly. Do you see any errors in there (it's a *huge* text file, though). Jason --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---
[sage-support] Re: Why this numerical integral bombs?
On Apr 25, 10:42 pm, Chris Seberino cseber...@gmail.com wrote: On Apr 25, 7:02 pm, mabshoff michael.absh...@mathematik.uni- dortmund.de wrote: What should I do until Sage 4.0? Is there a way to tweak my 3.4.1 to use clisp 2.47? Do you mean clisp 2.46 since 3.4.1 ships with clisp 2.47 while Sage 3.4 used clisp 2.46? Since *I* was able to reproduce the problem on 3.4 on Ubuntu Jaunty haven't I ruled out that clisp 2.47 is not the problem? * Above you asked about Sage 3.4.1 which ships clisp-2.47. * Sage 3.4 ships clisp-2.46 which you mention you tried on Ubuntu 9.04 in the message right above. I do not see how that rules out clisp 2.47 - or am I just not getting what you are driving at :) Is there something else I can try? I am not sure. cs Cheers, Michael --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---
[sage-support] Re: Why this numerical integral bombs?
On Apr 24, 7:09 pm, mabshoff michael.absh...@mathematik.uni- dortmund.de wrote: So I cannot reproduce this. What platform are you on and are you using a binary, build from source, etc? I'm on an Ubuntu 9.04 ThinkPad T41 laptop. I built my Sage install from 3.4 source myself as there was no binary yet for this latest Ubuntu. cs --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---
[sage-support] Re: Why this numerical integral bombs?
Michael, And Sage 3.4: | Sage Version 3.4, Release Date: 2009-03-11 | | Type notebook() for the GUI, and license() for information.| -- sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 So I cannot reproduce this. What platform are you on and are you using a binary, build from source, etc? the example that failed has 'exp' instead of 'cos'. However with Sage 3.4 it works for me: -- | Sage Version 3.4, Release Date: 2009-03-11 | | Type notebook() for the GUI, and license() for information.| -- sage: numerical_integral(sin(pi*exp(x/2)),0,2)[0] -0.43734547482524966 Paul Zimmermann --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---
[sage-support] Re: Why this numerical integral bombs?
On Apr 25, 12:23 am, Paul Zimmermann paul.zimmerm...@loria.fr wrote: Michael, Hi Paul, And Sage 3.4: | Sage Version 3.4, Release Date: 2009-03-11 | | Type notebook() for the GUI, and license() for information. | -- sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 So I cannot reproduce this. What platform are you on and are you using a binary, build from source, etc? the example that failed has 'exp' instead of 'cos'. However with Sage 3.4 it works for me: -- | Sage Version 3.4, Release Date: 2009-03-11 | | Type notebook() for the GUI, and license() for information. | -- sage: numerical_integral(sin(pi*exp(x/2)),0,2)[0] -0.43734547482524966 Paul Zimmermann I ran this on 3.4, 3.4.1 and a couple other rc releases in between: sage: for i in (0..200): : numerical_integral(sin(pi*exp(x/2)),0,2)[0] : and I could not get a failure. Note that the most likely culprit here is clisp IMHO. Just today I found out that clisp 2.47 doesn't build on the latest Mandriva release in 32 and 64 bit flavors. Oh well ... We should have an Ubuntu 9.04 test image in the not too distant future, so maybe we will hit it there. Cheers, Michael --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---
[sage-support] Re: Why this numerical integral bombs?
On Apr 25, 2:10 am, mabshoff michael.absh...@mathematik.uni- dortmund.de wrote: I ran this on 3.4, 3.4.1 and a couple other rc releases in between: sage: for i in (0..200): : numerical_integral(sin(pi*exp(x/2)),0,2)[0] : and I could not get a failure. Hmm. That is weird. I just compiled 3.4.1 from source and got the same error on the newer version too. We should have an Ubuntu 9.04 test image in the not too distant future, so maybe we will hit it there. That seems to be the only delta between you guys and me so chances are that's the problem. cs --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---
[sage-support] Re: Why this numerical integral bombs?
On Apr 25, 2:17 pm, Chris Seberino cseber...@gmail.com wrote: Hi Chris, On Apr 25, 2:10 am, mabshoff michael.absh...@mathematik.uni- SNIP Hmm. That is weird. I just compiled 3.4.1 from source and got the same error on the newer version too. Ok. SNIP That seems to be the only delta between you guys and me so chances are that's the problem. Not really, we switched from clisp 2.46 to 2.47 and if you look at https://groups.google.com/group/sage-devel/t/753a8d2f3c052183 for some people who have had issues for a long, long time we could never get to the bottom of all the sudden they completely disappear. Clisp seems to be working or not working depending on compiler release, target CPU, etc, so I am not surprised you are seeing oddities no one can reproduce. We are in progress of switching to ecl as the basis for Maxima which will hopefully serve as a more stable basis than clisp. Given that we want 3.4.2 out early next week the ecl switch will probably miss the boat on that one, but since it is essential to switch for Sage 4.0 to get Solaris support working it will be in 4.0. cs Cheers, Michael --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---
[sage-support] Re: Why this numerical integral bombs?
On Apr 25, 5:36 pm, mabshoff michael.absh...@mathematik.uni- dortmund.de wrote: Not really, we switched from clisp 2.46 to 2.47 and if you look at https://groups.google.com/group/sage-devel/t/753a8d2f3c052183 for some people who have had issues for a long, long time we could never get to the bottom of all the sudden they completely disappear. What should I do until Sage 4.0? Is there a way to tweak my 3.4.1 to use clisp 2.47? cs --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---
[sage-support] Re: Why this numerical integral bombs?
On Apr 25, 6:38 pm, Chris Seberino cseber...@gmail.com wrote: On Apr 25, 5:36 pm, mabshoff michael.absh...@mathematik.uni- Hi Chris, dortmund.de wrote: Not really, we switched from clisp 2.46 to 2.47 and if you look at https://groups.google.com/group/sage-devel/t/753a8d2f3c052183 for some people who have had issues for a long, long time we could never get to the bottom of all the sudden they completely disappear. What should I do until Sage 4.0? Is there a way to tweak my 3.4.1 to use clisp 2.47? Do you mean clisp 2.46 since 3.4.1 ships with clisp 2.47 while Sage 3.4 used clisp 2.46? You can force the use of clisp-2.46.spkg by running ./sage -f /patch/ to/clisp-2.46.spkg, but that will break the pure clisp mode AFAIK. It will also not use the readline-off mode we hacked into clisp 2.47, but if you never saw any Maxima hangs it should not matter to you too. You must also force a rebuild of the Maxima.spkg (via ./sage -f ../maxima- x.y.z.spkg after rebuilding clisp to make it use that version of clisp you just build without trouble). Note that the above change will likely be overwritten by updated Maxima and ecl.spkgs in Sage 4.0, but not 3.4.2. Maxima itself has a disable-readline option, but it seems broken with Maxima 5.16.3+clisp-2.[46|47], I don't know what the status of 5.18.1 which was just released is. cs Cheers, Michael --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---
[sage-support] Re: Why this numerical integral bombs?
On Apr 25, 7:02 pm, mabshoff michael.absh...@mathematik.uni- dortmund.de wrote: What should I do until Sage 4.0? Is there a way to tweak my 3.4.1 to use clisp 2.47? Do you mean clisp 2.46 since 3.4.1 ships with clisp 2.47 while Sage 3.4 used clisp 2.46? Since *I* was able to reproduce the problem on 3.4 on Ubuntu Jaunty haven't I ruled out that clisp 2.47 is not the problem? Is there something else I can try? cs --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---
[sage-support] Re: Why this numerical integral bombs?
On Apr 24, 7:03 pm, Chris Seberino cseber...@gmail.com wrote: Hi, Notice first and second are very similar except for the cos() and exp (). Why first one ok but second bombs? sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 I guess you are using Sage 3.4? Works for me in my current merge tree: -- | Sage Version 3.4.2.alpha0, Release Date: 2009-04-24| | Type notebook() for the GUI, and license() for information.| -- sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: 3.4.1: -- | Sage Version 3.4.1, Release Date: 2009-04-21 | | Type notebook() for the GUI, and license() for information.| -- sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: And Sage 3.4: | Sage Version 3.4, Release Date: 2009-03-11 | | Type notebook() for the GUI, and license() for information.| -- sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: numerical_integral(sin(pi*cos(x/2)),0,2)[0] 0.85397903781471396 sage: So I cannot reproduce this. What platform are you on and are you using a binary, build from source, etc? Cheers, Michael --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~--~~~~--~~--~--~---