<div dir="ltr"><div><div>Hi Bernard,<br><br></div>Looks like you're on to something - two other people have raised this discrepancy before: <a href="https://github.com/numpy/numpy/issues/2588">https://github.com/numpy/numpy/issues/2588</a>. Unfortunately, when it comes to resolving the discrepancy one of the previous comments takes the opposite view. Namely, that the docstring is correct and the code is wrong.<br>
<br></div><div>Do different domains use different conventions here? Are there some references to back up one stance or another?<br><br></div><div>But all else being equal, I'm guessing there'll be far more appetite for updating the documentation than the code.<br>
</div><div><br></div>Regards,<br>Richard Hattersley<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On 7 October 2013 22:09, Bernhard Spinnler <span dir="ltr"><<a href="mailto:Bernhard.Spinnler@gmx.net" target="_blank">Bernhard.Spinnler@gmx.net</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">The numpy.correlate documentation says:<br>
<br>
        correlate(a, v) = z[k] = sum_n a[n] * conj(v[n+k])<br>
<br>
In [1]: a = [1, 2]<br>
<br>
In [2]: v = [2, 1j]<br>
<br>
In [3]: z = correlate(a, v, 'full')<br>
<br>
In [4]: z<br>
Out[4]: array([ 0.-1.j,  2.-2.j,  4.+0.j])<br>
<br>
However, according to the documentation, z should be<br>
<br>
        z[-1] = a[1] * conj(v[0]) = 4.+0.j<br>
        z[0]  = a[0] * conj(v[0]) + a[1] * conj(v[1]) = 2.-2.j<br>
        z[1] = a[0] * conj(v[1]) = 0.-1.j<br>
<br>
which is the time reversed version of what correlate() calculates.<br>
<br>
IMHO, the correlate() code is correct. The correct formula in the docs (which is also the correlation formula in standard text books) should be<br>
<br>
        z[k] = sum_n a[n+k] * conj(v[n])<br>
<br>
Cheers,<br>
        Bernhard<br>
_______________________________________________<br>
NumPy-Discussion mailing list<br>
<a href="mailto:NumPy-Discussion@scipy.org">NumPy-Discussion@scipy.org</a><br>
<a href="http://mail.scipy.org/mailman/listinfo/numpy-discussion" target="_blank">http://mail.scipy.org/mailman/listinfo/numpy-discussion</a><br>
</blockquote></div><br></div>