On 11/14/14, 1:23 PM, serge Guelton wrote:
On Fri, Nov 14, 2014 at 12:02:15PM -0800, Mehdi AMINI wrote:Hi, from pythran/pythonic/patch/README.rst """ The implementation of std::complex is very slow, due to the complex limited range (see `-fcx-limited-range`) feature. Numpy does not implement it, so we have to conform to numpy's version. """Let's focus on the multiplication.
Are you kidding, I explicitly linked the Python division implementation!
So, Python/Numpy behavior and C++ behavior are different with respect to complex multiplication. What we did is provide a std::complex<T> that match Python/Numpy behavior. We cannot ask nt2 to follow numpy's behvaior, so... we hijacked std::complex :-/
You have to ask NT2 to be able to supply your own complex implementation, hijacking std::complex is a no-no! (non standard I assume)
Mehdi