Function nc_abs2

numexpr3/complex_functions.hpp:623–637

}
    
static void
nc_abs2( 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] = 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] = x[I*sb1].real*x[I*sb1].real + x[I*sb1].imag*x[I*sb1].imag;
        }
    }
}

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