Function boxMullerTransform

src/backend/cpu/kernel/random_engine.hpp:259–270

template<typename T>
void boxMullerTransform(data_t<T> *const out1, data_t<T> *const out2,
                        const T r1, const T r2) {
    /*
     * The log of a real value x where 0 < x < 1 is negative.
     */
    using Tc = compute_t<T>;
    Tc r     = sqrt((Tc)(-2.0) * log(static_cast<Tc>(r2)));
    Tc theta = PI_VAL * (static_cast<Tc>(r1));

    *out1 = r * sin(theta);
    *out2 = r * cos(theta);
}

void boxMullerTransform(uint val[4], double *temp) {
    boxMullerTransform<double>(&temp[0], &temp[1], getDoubleNegative11(val, 0),