<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><br class=""><div><blockquote type="cite" class=""><div class="">Am 21.03.2015 um 05:13 schrieb Travis Oliphant <<a href="mailto:travis@continuum.io" class="">travis@continuum.io</a>>:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class="">I have not tested this, but I suspect that the MATLAB routine is using the erf implementation from the Intel Math Kernel Libraries (MKL).    <div class=""><br class=""></div><div class="">There is a function in MKL called vdErf that takes a vector of doubles and is likely tuned to the hardware.    This could be linked to NumPy with similar speed benefits. <div class=""><div class=""><div class=""><br class=""></div><div class="">-Travis</div><div class=""><br class=""></div></div></div></div></div></div></blockquote></div><div class=""><br class=""></div>The uvml module <a href="https://github.com/geggo/uvml" class="">https://github.com/geggo/uvml</a> exposes the fast MKL/VML erf implementation to numpy. Unfortunately, no binaries available.<div class=""><br class=""></div><div class="">Gregor</div></body></html>