strided_expression.c

static int __mini_mangle___pyx_array_expression5strided(Py_ssize_t const *const CYTHON_RESTRICT __mini_mangle___pyx_shape, double *const CYTHON_RESTRICT __pyx_op1_data, Py_ssize_t const *const CYTHON_RESTRICT __pyx_op1_strides, double const *const CYTHON_RESTRICT __pyx_op2_data, Py_ssize_t const *const CYTHON_RESTRICT __pyx_op2_strides, double const *const CYTHON_RESTRICT __pyx_op3_data, Py_ssize_t const *const CYTHON_RESTRICT __pyx_op3_strides) {
    Py_ssize_t __mini_mangle_temp0;
    #ifdef _OPENMP
    #pragma omp parallel for if(1)
    #endif
    for (__mini_mangle_temp0 = 0; __mini_mangle_temp0 < (__mini_mangle___pyx_shape[0]); __mini_mangle_temp0++) {
        double *CYTHON_RESTRICT __mini_mangle_temp1;
        double const *CYTHON_RESTRICT __mini_mangle_temp2;
        double const *CYTHON_RESTRICT __mini_mangle_temp3;
        Py_ssize_t __mini_mangle_temp4;
        __mini_mangle_temp1 = ((double *) (((char *) __pyx_op1_data) + (__mini_mangle_temp0 * (__pyx_op1_strides[0]))));
        __mini_mangle_temp2 = ((double *) (((char *) __pyx_op2_data) + (__mini_mangle_temp0 * (__pyx_op2_strides[0]))));
        __mini_mangle_temp3 = ((double *) (((char *) __pyx_op3_data) + (__mini_mangle_temp0 * (__pyx_op3_strides[0]))));
        for (__mini_mangle_temp4 = 0; __mini_mangle_temp4 < (__mini_mangle___pyx_shape[1]); __mini_mangle_temp4++) {
            (__mini_mangle_temp1[__mini_mangle_temp4]) = ((__mini_mangle_temp2[__mini_mangle_temp4]) + (*((double const *CYTHON_RESTRICT) (((char *) __mini_mangle_temp3) + (__mini_mangle_temp4 * (__pyx_op3_strides[1]))))));
        }
    }
    return 0;
}