Function nc_mul

numexpr3/complex_functions.hpp:395–410

}

static void
nc_mul( npy_intp n, npy_complex64 *a, npy_intp sb1, npy_complex64 *b, npy_intp sb2, npy_complex64 *r)
{
    if( sb1 == sizeof(npy_complex64) && sb2 == sizeof(npy_complex64) ) { // Aligned
        for( npy_intp I = 0; I < n; I++ ) {
            _inline_mul( a[I], b[I], r[I] );
        }
    }
    else {
        sb2 /= sizeof(npy_complex64);
        sb1 /= sizeof(npy_complex64);
        for( npy_intp I = 0; I < n; I++ ) {
            _inline_mul( a[I*sb1], b[I*sb2], r[I] );
        }
    }
}

static void
nc_mul( npy_intp n, npy_complex128 *a, npy_intp sb1, npy_complex128 *b, npy_intp sb2, npy_complex128 *r)