Function boxMullerTransform

src/backend/oneapi/kernel/random_engine_write.hpp:118–131

//
template<typename Tc>
static void boxMullerTransform(cfloat *const cOut, const Tc &r1, const Tc &r2) {
    /*
     * The log of a real value x where 0 < x < 1 is negative.
     */
    Tc r = sycl::sqrt(neg_two<Tc>() * sycl::log(r2));
    Tc s, c;

    // Multiplying by PI instead of 2*PI seems to yeild a better distribution
    // even though the original boxMuller algorithm calls for 2 * PI
    // sincos(two_pi<Tc>() * r1, &s, &c);
    sincospi(r1, &s, &c);
    cOut->real(static_cast<float>(r * s));
    cOut->imag(static_cast<float>(r * c));
}

template<typename Tc>
static void boxMullerTransform(cdouble *const cOut, const Tc &r1,