[Numpy-discussion] Problem with correlate
Charles R Harris
charlesr.harris@gmail....
Sun May 31 14:47:55 CDT 2009
On Sun, May 31, 2009 at 11:54 AM, rob steed <rjsteed@talk21.com> wrote:
>
> Hi,
> After my previous email, I have opened a ticket #1117 (correlate not order
> dependent)
>
> I have found that the correlate function is defined in multiarraymodule.c
> and
> that inputs are being swapped using the following code
>
> n1 = ap1->dimensions[0];
> n2 = ap2->dimensions[0];
> if (n1 < n2) {
> ret = ap1;
> ap1 = ap2;
> ap2 = ret;
> ret = NULL;
> i = n1;
> n1 = n2;
> n2 = i;
> }
>
> I do not know the code well enough to see whether this could just be
> removed (I don't know c either).
> Maybe the algorithmn requires the inputs to be length ordered? I will try
> to work it out.
>
If the correlation algorithm doesn't use an fft and is done explicitly, then
the maximum overlap for any shift is the length of the shortest input.
Swapping the arrays makes that logic easier to implement, but it isn't
necessary.
Chuck
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/numpy-discussion/attachments/20090531/2653accd/attachment.html
More information about the Numpy-discussion
mailing list