\sqrt{\frac{u1}{1 - u1}} \cdot \cos \left(6.28318530718 \cdot u2\right)
\left(\sqrt{\frac{u1}{1 - u1}} \cdot \sqrt[3]{{\cos \left(u2 \cdot 6.28318530718\right)}^{2}}\right) \cdot \sqrt[3]{\cos \left(e^{\log 6.28318530718 + \log u2}\right)}
(FPCore (cosTheta_i u1 u2) :precision binary32 (* (sqrt (/ u1 (- 1.0 u1))) (cos (* 6.28318530718 u2))))
(FPCore (cosTheta_i u1 u2) :precision binary32 (* (* (sqrt (/ u1 (- 1.0 u1))) (cbrt (pow (cos (* u2 6.28318530718)) 2.0))) (cbrt (cos (exp (+ (log 6.28318530718) (log u2)))))))
float code(float cosTheta_i, float u1, float u2) {
return sqrtf(u1 / (1.0f - u1)) * cosf(6.28318530718f * u2);
}
float code(float cosTheta_i, float u1, float u2) {
return (sqrtf(u1 / (1.0f - u1)) * cbrtf(powf(cosf(u2 * 6.28318530718f), 2.0f))) * cbrtf(cosf(expf(logf(6.28318530718f) + logf(u2))));
}



Bits error versus cosTheta_i



Bits error versus u1



Bits error versus u2
Results
Initial program 0.3
Applied add-exp-log_binary320.3
Applied add-exp-log_binary320.4
Applied prod-exp_binary320.3
Applied add-cube-cbrt_binary320.5
Applied associate-*r*_binary320.5
Simplified0.4
Applied *-un-lft-identity_binary320.4
Applied cbrt-prod_binary320.4
Applied associate-*l*_binary320.4
Simplified0.3
Final simplification0.3
herbie shell --seed 2022076
(FPCore (cosTheta_i u1 u2)
:name "Trowbridge-Reitz Sample, near normal, slope_x"
:precision binary32
:pre (and (and (and (> cosTheta_i 0.9999) (<= cosTheta_i 1.0)) (and (<= 2.328306437e-10 u1) (<= u1 1.0))) (and (<= 2.328306437e-10 u2) (<= u2 1.0)))
(* (sqrt (/ u1 (- 1.0 u1))) (cos (* 6.28318530718 u2))))