| Alternative 1 | |
|---|---|
| Error | 3.2 |
| Cost | 9892 |
\[\begin{array}{l}
t_0 := \cos \left(6.28318530718 \cdot u2\right)\\
\mathbf{if}\;t_0 \leq 0.9999200105667114:\\
\;\;\;\;t_0 \cdot \sqrt{u1}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\frac{u1}{1 - u1}}\\
\end{array}
\]
(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 (* 6.28318530718 u2)) 3.0))))
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((6.28318530718f * u2)), 3.0f));
}
function code(cosTheta_i, u1, u2) return Float32(sqrt(Float32(u1 / Float32(Float32(1.0) - u1))) * cos(Float32(Float32(6.28318530718) * u2))) end
function code(cosTheta_i, u1, u2) return Float32(sqrt(Float32(u1 / Float32(Float32(1.0) - u1))) * cbrt((cos(Float32(Float32(6.28318530718) * u2)) ^ Float32(3.0)))) end
\sqrt{\frac{u1}{1 - u1}} \cdot \cos \left(6.28318530718 \cdot u2\right)
\sqrt{\frac{u1}{1 - u1}} \cdot \sqrt[3]{{\cos \left(6.28318530718 \cdot u2\right)}^{3}}
Results
Initial program 0.3
Applied egg-rr0.3
Final simplification0.3
| Alternative 1 | |
|---|---|
| Error | 3.2 |
| Cost | 9892 |
| Alternative 2 | |
|---|---|
| Error | 2.3 |
| Cost | 6884 |
| Alternative 3 | |
|---|---|
| Error | 0.3 |
| Cost | 6688 |
| Alternative 4 | |
|---|---|
| Error | 6.3 |
| Cost | 3360 |
| Alternative 5 | |
|---|---|
| Error | 11.7 |
| Cost | 3232 |
| Alternative 6 | |
|---|---|
| Error | 29.9 |
| Cost | 32 |

herbie shell --seed 2022297
(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))))