Function _inline_div

numexpr3/complex_functions.hpp:437–444

The naive complex division function here is not very accurate and sometimes fails the autotest assert. Probably there's something better... https://arxiv.org/pdf/1210.4539.pdf or https://arxiv.org/pdf/1608.07596.pdf Most complex division algorithms seem to have a branch however, which would take a bit of work to avoid.

// Most complex division algorithms seem to have a branch however, which
// would take a bit of work to avoid.
static inline void
_inline_div( npy_complex64 a, npy_complex64 b, npy_complex64 &r)
{
    npy_float32 d = 1.0f/(b.real*b.real + b.imag*b.imag);
    r.real = (a.real*b.real + a.imag*b.imag)*d;
    r.imag = (a.imag*b.real - a.real*b.imag)*d;
 
}
static inline void
_inline_div( npy_complex128 a, npy_complex128 b, npy_complex128 &r)
{