Problem with scipy.integrate.nquad reported on stackoverflow
Could someone with a Windows installation of scipy 0.14.0 run the code at http://stackoverflow.com/questions/27270044/big-integration-error-with-integ... It works for me on Mac OS X. In the comments to the question, Jonathan March reported that the problem occurs with both 32 and 64 bit Windows 7. Does the code work for anyone using Windows? ("Work" means the values reported for 'Intf1corr' and 'Intf2' are the same.) Warren
On Dec 5, 2014, at 5:15 AM, Warren Weckesser <warren.weckesser@gmail.com> wrote:
Could someone with a Windows installation of scipy 0.14.0 run the code at http://stackoverflow.com/questions/27270044/big-integration-error-with-integ...
It works for me on Mac OS X. In the comments to the question, Jonathan March reported that the problem occurs with both 32 and 64 bit Windows 7. Does the code work for anyone using Windows? ("Work" means the values reported for 'Intf1corr' and 'Intf2' are the same.)
Warren _______________________________________________ SciPy-Dev mailing list SciPy-Dev@scipy.org http://mail.scipy.org/mailman/listinfo/scipy-dev
Can someone ask whether or not the old function dblquad gives the same incorrect results? SO won’t let me comment yet, and I don’t think it’s really an answer...
On Fri, Dec 5, 2014 at 11:12 AM, Gmail <charlesnwoods@gmail.com> wrote:
On Dec 5, 2014, at 5:15 AM, Warren Weckesser <warren.weckesser@gmail.com> wrote:
Could someone with a Windows installation of scipy 0.14.0 run the code at http://stackoverflow.com/questions/27270044/big-integration-error-with-integ... ?
It works for me on Mac OS X. In the comments to the question, Jonathan March reported that the problem occurs with both 32 and 64 bit Windows 7. Does the code work for anyone using Windows? ("Work" means the values reported for 'Intf1corr' and 'Intf2' are the same.)
Warren _______________________________________________ SciPy-Dev mailing list SciPy-Dev@scipy.org http://mail.scipy.org/mailman/listinfo/scipy-dev
Can someone ask whether or not the old function dblquad gives the same incorrect results? SO won’t let me comment yet, and I don’t think it’s really an answer...
Good idea. I added a comment with the suggestion to check dblquad. Warren
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@scipy.org http://mail.scipy.org/mailman/listinfo/scipy-dev
On 12/5/14, Gmail <charlesnwoods@gmail.com> wrote:
On Dec 5, 2014, at 5:15 AM, Warren Weckesser <warren.weckesser@gmail.com> wrote:
Could someone with a Windows installation of scipy 0.14.0 run the code at http://stackoverflow.com/questions/27270044/big-integration-error-with-integ...
It works for me on Mac OS X. In the comments to the question, Jonathan March reported that the problem occurs with both 32 and 64 bit Windows 7. Does the code work for anyone using Windows? ("Work" means the values reported for 'Intf1corr' and 'Intf2' are the same.)
Warren _______________________________________________ SciPy-Dev mailing list SciPy-Dev@scipy.org http://mail.scipy.org/mailman/listinfo/scipy-dev
Can someone ask whether or not the old function dblquad gives the same incorrect results? SO won’t let me comment yet, and I don’t think it’s really an answer...
The author of the question reported in a comment that dblquad also gives the wrong result (http://stackoverflow.com/questions/27270044/big-integration-error-with-integ...). Warren
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@scipy.org http://mail.scipy.org/mailman/listinfo/scipy-dev
At least they’re giving consistent results. nquad and dblquad are both just wrappers that call scipy.integrate.quad on itself, so it seems that something about this recursion is interacting badly, maybe with MKL. Maybe the bug can be reproduced with just a call to quad? I’m afraid my knowledge of the .integrate package peters out about here. N
On Dec 8, 2014, at 3:14 PM, Warren Weckesser <warren.weckesser@gmail.com> wrote:
On 12/5/14, Gmail <charlesnwoods@gmail.com> wrote:
On Dec 5, 2014, at 5:15 AM, Warren Weckesser <warren.weckesser@gmail.com> wrote:
Could someone with a Windows installation of scipy 0.14.0 run the code at http://stackoverflow.com/questions/27270044/big-integration-error-with-integ...
It works for me on Mac OS X. In the comments to the question, Jonathan March reported that the problem occurs with both 32 and 64 bit Windows 7. Does the code work for anyone using Windows? ("Work" means the values reported for 'Intf1corr' and 'Intf2' are the same.)
Warren _______________________________________________ SciPy-Dev mailing list SciPy-Dev@scipy.org http://mail.scipy.org/mailman/listinfo/scipy-dev
Can someone ask whether or not the old function dblquad gives the same incorrect results? SO won’t let me comment yet, and I don’t think it’s really an answer...
The author of the question reported in a comment that dblquad also gives the wrong result (http://stackoverflow.com/questions/27270044/big-integration-error-with-integ...).
Warren
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@scipy.org http://mail.scipy.org/mailman/listinfo/scipy-dev
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@scipy.org http://mail.scipy.org/mailman/listinfo/scipy-dev
On Fri, Dec 5, 2014 at 12:15 PM, Warren Weckesser < warren.weckesser@gmail.com> wrote:
Could someone with a Windows installation of scipy 0.14.0 run the code at http://stackoverflow.com/questions/27270044/big-integration-error-with-integ... ?
It works for me on Mac OS X. In the comments to the question, Jonathan March reported that the problem occurs with both 32 and 64 bit Windows 7. Does the code work for anyone using Windows? ("Work" means the values reported for 'Intf1corr' and 'Intf2' are the same.)
I briefly looked at it: * canopy (numpy 1.8.1 and scipy 0.14.0, VS 2008 + ifort + MKL): bug * gholke eggs on python 3.4 w/ numpy 1.9.1 and scipy 0.14.0 (MKL as well): bug * official binaries on 2.7: no bug I suspect either an ifort/MKL bug, or some bad interactions between scipy and this config. I will update once I have time to actually look into the issue. David
Warren _______________________________________________ SciPy-Dev mailing list SciPy-Dev@scipy.org http://mail.scipy.org/mailman/listinfo/scipy-dev
On Fri, Dec 5, 2014 at 11:49 AM, David Cournapeau <cournape@gmail.com> wrote:
On Fri, Dec 5, 2014 at 12:15 PM, Warren Weckesser < warren.weckesser@gmail.com> wrote:
Could someone with a Windows installation of scipy 0.14.0 run the code at http://stackoverflow.com/questions/27270044/big-integration-error-with-integ... ?
It works for me on Mac OS X. In the comments to the question, Jonathan March reported that the problem occurs with both 32 and 64 bit Windows 7. Does the code work for anyone using Windows? ("Work" means the values reported for 'Intf1corr' and 'Intf2' are the same.)
I briefly looked at it:
* canopy (numpy 1.8.1 and scipy 0.14.0, VS 2008 + ifort + MKL): bug * gholke eggs on python 3.4 w/ numpy 1.9.1 and scipy 0.14.0 (MKL as well): bug * official binaries on 2.7: no bug
I suspect either an ifort/MKL bug, or some bad interactions between scipy and this config. I will update once I have time to actually look into the issue.
David
Thanks, David. Warren
Warren _______________________________________________ SciPy-Dev mailing list SciPy-Dev@scipy.org http://mail.scipy.org/mailman/listinfo/scipy-dev
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@scipy.org http://mail.scipy.org/mailman/listinfo/scipy-dev
participants (3)
-
David Cournapeau -
Gmail -
Warren Weckesser