[SciPy-Dev] Problem with scipy.integrate.nquad reported on stackoverflow

[Gmail]

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 at gmail.com> wrote:
> 
> On 12/5/14, Gmail <charlesnwoods at gmail.com> wrote:
>> 
>>> On Dec 5, 2014, at 5:15 AM, Warren Weckesser <warren.weckesser at 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-integrate-ndquad?
>>> 
>>> 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 at 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-integrate-nquad).
> 
> Warren
> 
> 
>> _______________________________________________
>> SciPy-Dev mailing list
>> SciPy-Dev at scipy.org
>> http://mail.scipy.org/mailman/listinfo/scipy-dev
>> 
> _______________________________________________
> SciPy-Dev mailing list
> SciPy-Dev at scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-dev