Function box_muller

src/framework/sampling/torch_random.cpp:64–72

    }
    return counter;
}

void box_muller(uint32_t uniform0, uint32_t uniform1, float & normal0, float & normal1) {
    const float radius_input =
        static_cast<float>(uniform0) * kInvTwoPow32 + (kInvTwoPow32 * 0.5F);
    const float angle =
        static_cast<float>(uniform1) * kInvTwoPow32TwoPi + (kInvTwoPow32TwoPi * 0.5F);
    const float radius = std::sqrt(-2.0F * std::log(radius_input));
    normal0 = radius * std::sin(angle);
    normal1 = radius * std::cos(angle);
}