Function nc_abs

numexpr3/complex_functions.hpp:591–605

static void
nc_abs( npy_intp n, npy_complex64 *x, npy_intp sb1, npy_float32 *r )
{
    if( sb1 == sizeof(npy_complex64) ) { // Aligned
        for( npy_intp I = 0; I < n; I++ ) {
             r[I] = sqrtf( x[I].real*x[I].real + x[I].imag*x[I].imag);
        }
    }
    else {
        sb1 /= sizeof(npy_complex64);
        for( npy_intp I = 0; I < n; I++ ) {
             r[I] = sqrtf( x[I*sb1].real*x[I*sb1].real + x[I*sb1].imag*x[I*sb1].imag);
        }
    }
}

static void
nc_abs( npy_intp n, npy_complex128 *x, npy_intp sb1, npy_float64 *r )