<br><br><div class="gmail_quote">On Sat, Jul 17, 2010 at 9:48 PM, Charles R Harris <span dir="ltr">&lt;<a href="mailto:charlesr.harris@gmail.com">charlesr.harris@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br><br><div class="gmail_quote"><div><div></div><div class="h5">On Sat, Jul 17, 2010 at 7:38 AM, Ralf Gommers <span dir="ltr">&lt;<a href="mailto:ralf.gommers@googlemail.com" target="_blank">ralf.gommers@googlemail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br><br><div class="gmail_quote"><div><div></div><div>On Sat, Jul 17, 2010 at 8:57 PM, Charles R Harris <span dir="ltr">&lt;<a href="mailto:charlesr.harris@gmail.com" target="_blank">charlesr.harris@gmail.com</a>&gt;</span> wrote:<br>

<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br><br><div class="gmail_quote"><div><div></div><div>On Sat, Jul 17, 2010 at 2:42 AM, Ralf Gommers <span dir="ltr">&lt;<a href="mailto:ralf.gommers@googlemail.com" target="_blank">ralf.gommers@googlemail.com</a>&gt;</span> wrote:<br>


<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
integrate.cumtrapz and numpy.trapz do exactly the same thing and the code is very similar but not identical, as pointed out in <a href="http://projects.scipy.org/scipy/ticket/720" target="_blank">http://projects.scipy.org/scipy/ticket/720</a>. Assuming numpy.trapz is not going anywhere, can we replace the scipy version with it? <br>




<br>For a small bonus, the numpy version is about 10% faster (tested for several array shapes):<br>&gt;&gt;&gt; a = np.arange(1e4).reshape(500, 20)<br>&gt;&gt;&gt; %timeit np.trapz(a, axis=1)<br>10000 loops, best of 3: 182 us per loop<br>




&gt;&gt;&gt; %timeit sp.integrate.cumtrapz(a, axis=1)<br>1000 loops, best of 3: 209 us per loop<br><br></blockquote></div></div><div><br>Replacement seems reasonable to me. We should try to prune duplicate functionality with numpy taking precedence. <br>



<br></div></div></blockquote></div></div><div>Just to double check, that also means deprecating the name cumtrapz right?<br><br></div></div></blockquote></div></div><div><br>Yes, I think so. And maybe a ticket to remove it later although that means having versions in trac for future releases. Maybe we should have a &quot;deprecated&quot; component so we can use such tickets.<br>

<br></div></div></blockquote><div>I was thinking about adding &quot;check all deprecations&quot; as part of the release process (at the beginning). Then there&#39;s less chance to forget it, and less noise in Trac.<br><br>
Ralf<br></div></div><br><br>