
(FPCore (ux uy maxCos) :precision binary32 (let* ((t_0 (+ (- 1.0 ux) (* ux maxCos)))) (* (sin (* (* uy 2.0) PI)) (sqrt (- 1.0 (* t_0 t_0))))))
float code(float ux, float uy, float maxCos) {
float t_0 = (1.0f - ux) + (ux * maxCos);
return sinf(((uy * 2.0f) * ((float) M_PI))) * sqrtf((1.0f - (t_0 * t_0)));
}
function code(ux, uy, maxCos) t_0 = Float32(Float32(Float32(1.0) - ux) + Float32(ux * maxCos)) return Float32(sin(Float32(Float32(uy * Float32(2.0)) * Float32(pi))) * sqrt(Float32(Float32(1.0) - Float32(t_0 * t_0)))) end
function tmp = code(ux, uy, maxCos) t_0 = (single(1.0) - ux) + (ux * maxCos); tmp = sin(((uy * single(2.0)) * single(pi))) * sqrt((single(1.0) - (t_0 * t_0))); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(1 - ux\right) + ux \cdot maxCos\\
\sin \left(\left(uy \cdot 2\right) \cdot \pi\right) \cdot \sqrt{1 - t_0 \cdot t_0}
\end{array}
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 18 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (ux uy maxCos) :precision binary32 (let* ((t_0 (+ (- 1.0 ux) (* ux maxCos)))) (* (sin (* (* uy 2.0) PI)) (sqrt (- 1.0 (* t_0 t_0))))))
float code(float ux, float uy, float maxCos) {
float t_0 = (1.0f - ux) + (ux * maxCos);
return sinf(((uy * 2.0f) * ((float) M_PI))) * sqrtf((1.0f - (t_0 * t_0)));
}
function code(ux, uy, maxCos) t_0 = Float32(Float32(Float32(1.0) - ux) + Float32(ux * maxCos)) return Float32(sin(Float32(Float32(uy * Float32(2.0)) * Float32(pi))) * sqrt(Float32(Float32(1.0) - Float32(t_0 * t_0)))) end
function tmp = code(ux, uy, maxCos) t_0 = (single(1.0) - ux) + (ux * maxCos); tmp = sin(((uy * single(2.0)) * single(pi))) * sqrt((single(1.0) - (t_0 * t_0))); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(1 - ux\right) + ux \cdot maxCos\\
\sin \left(\left(uy \cdot 2\right) \cdot \pi\right) \cdot \sqrt{1 - t_0 \cdot t_0}
\end{array}
\end{array}
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0 (sin (* PI (* 2.0 uy))))
(t_1
(fma
ux
(* 2.0 (- 1.0 maxCos))
(* (+ maxCos -1.0) (* (- 1.0 maxCos) (pow ux 2.0))))))
(cbrt (* (* t_1 (sqrt t_1)) (* t_0 (pow t_0 2.0))))))
float code(float ux, float uy, float maxCos) {
float t_0 = sinf((((float) M_PI) * (2.0f * uy)));
float t_1 = fmaf(ux, (2.0f * (1.0f - maxCos)), ((maxCos + -1.0f) * ((1.0f - maxCos) * powf(ux, 2.0f))));
return cbrtf(((t_1 * sqrtf(t_1)) * (t_0 * powf(t_0, 2.0f))));
}
function code(ux, uy, maxCos) t_0 = sin(Float32(Float32(pi) * Float32(Float32(2.0) * uy))) t_1 = fma(ux, Float32(Float32(2.0) * Float32(Float32(1.0) - maxCos)), Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(Float32(1.0) - maxCos) * (ux ^ Float32(2.0))))) return cbrt(Float32(Float32(t_1 * sqrt(t_1)) * Float32(t_0 * (t_0 ^ Float32(2.0))))) end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(\pi \cdot \left(2 \cdot uy\right)\right)\\
t_1 := \mathsf{fma}\left(ux, 2 \cdot \left(1 - maxCos\right), \left(maxCos + -1\right) \cdot \left(\left(1 - maxCos\right) \cdot {ux}^{2}\right)\right)\\
\sqrt[3]{\left(t_1 \cdot \sqrt{t_1}\right) \cdot \left(t_0 \cdot {t_0}^{2}\right)}
\end{array}
\end{array}
Initial program 55.6%
associate-*l*55.6%
cancel-sign-sub-inv55.6%
+-commutative55.6%
*-commutative55.6%
fma-def55.8%
+-commutative55.8%
associate-+r-55.9%
fma-def55.9%
neg-sub055.9%
associate-+l-55.9%
associate--r-55.9%
metadata-eval55.9%
+-commutative55.9%
*-lft-identity55.9%
*-commutative55.9%
distribute-lft-out--55.9%
Simplified55.9%
Taylor expanded in ux around 0 98.1%
fma-def98.1%
+-commutative98.1%
associate--l+98.2%
sub-neg98.2%
metadata-eval98.2%
distribute-lft-in98.2%
metadata-eval98.2%
+-commutative98.2%
mul-1-neg98.2%
sub-neg98.2%
count-298.2%
*-commutative98.2%
sub-neg98.2%
metadata-eval98.2%
+-commutative98.2%
Simplified98.2%
*-commutative98.2%
add-cbrt-cube98.1%
add-cbrt-cube98.1%
cbrt-unprod98.0%
Applied egg-rr98.2%
Final simplification98.2%
(FPCore (ux uy maxCos)
:precision binary32
(*
(log1p (expm1 (sin (* uy (* 2.0 PI)))))
(sqrt
(fma
ux
(* 2.0 (- 1.0 maxCos))
(* (pow ux 2.0) (* (- 1.0 maxCos) (+ maxCos -1.0)))))))
float code(float ux, float uy, float maxCos) {
return log1pf(expm1f(sinf((uy * (2.0f * ((float) M_PI)))))) * sqrtf(fmaf(ux, (2.0f * (1.0f - maxCos)), (powf(ux, 2.0f) * ((1.0f - maxCos) * (maxCos + -1.0f)))));
}
function code(ux, uy, maxCos) return Float32(log1p(expm1(sin(Float32(uy * Float32(Float32(2.0) * Float32(pi)))))) * sqrt(fma(ux, Float32(Float32(2.0) * Float32(Float32(1.0) - maxCos)), Float32((ux ^ Float32(2.0)) * Float32(Float32(Float32(1.0) - maxCos) * Float32(maxCos + Float32(-1.0))))))) end
\begin{array}{l}
\\
\mathsf{log1p}\left(\mathsf{expm1}\left(\sin \left(uy \cdot \left(2 \cdot \pi\right)\right)\right)\right) \cdot \sqrt{\mathsf{fma}\left(ux, 2 \cdot \left(1 - maxCos\right), {ux}^{2} \cdot \left(\left(1 - maxCos\right) \cdot \left(maxCos + -1\right)\right)\right)}
\end{array}
Initial program 55.6%
associate-*l*55.6%
cancel-sign-sub-inv55.6%
+-commutative55.6%
*-commutative55.6%
fma-def55.8%
+-commutative55.8%
associate-+r-55.9%
fma-def55.9%
neg-sub055.9%
associate-+l-55.9%
associate--r-55.9%
metadata-eval55.9%
+-commutative55.9%
*-lft-identity55.9%
*-commutative55.9%
distribute-lft-out--55.9%
Simplified55.9%
Taylor expanded in ux around 0 98.1%
fma-def98.1%
+-commutative98.1%
associate--l+98.2%
sub-neg98.2%
metadata-eval98.2%
distribute-lft-in98.2%
metadata-eval98.2%
+-commutative98.2%
mul-1-neg98.2%
sub-neg98.2%
count-298.2%
*-commutative98.2%
sub-neg98.2%
metadata-eval98.2%
+-commutative98.2%
Simplified98.2%
log1p-expm1-u98.2%
Applied egg-rr98.2%
Final simplification98.2%
(FPCore (ux uy maxCos)
:precision binary32
(*
(sin (* uy (* 2.0 PI)))
(sqrt
(fma
ux
(* 2.0 (- 1.0 maxCos))
(* (pow ux 2.0) (* (- 1.0 maxCos) (+ maxCos -1.0)))))))
float code(float ux, float uy, float maxCos) {
return sinf((uy * (2.0f * ((float) M_PI)))) * sqrtf(fmaf(ux, (2.0f * (1.0f - maxCos)), (powf(ux, 2.0f) * ((1.0f - maxCos) * (maxCos + -1.0f)))));
}
function code(ux, uy, maxCos) return Float32(sin(Float32(uy * Float32(Float32(2.0) * Float32(pi)))) * sqrt(fma(ux, Float32(Float32(2.0) * Float32(Float32(1.0) - maxCos)), Float32((ux ^ Float32(2.0)) * Float32(Float32(Float32(1.0) - maxCos) * Float32(maxCos + Float32(-1.0))))))) end
\begin{array}{l}
\\
\sin \left(uy \cdot \left(2 \cdot \pi\right)\right) \cdot \sqrt{\mathsf{fma}\left(ux, 2 \cdot \left(1 - maxCos\right), {ux}^{2} \cdot \left(\left(1 - maxCos\right) \cdot \left(maxCos + -1\right)\right)\right)}
\end{array}
Initial program 55.6%
associate-*l*55.6%
cancel-sign-sub-inv55.6%
+-commutative55.6%
*-commutative55.6%
fma-def55.8%
+-commutative55.8%
associate-+r-55.9%
fma-def55.9%
neg-sub055.9%
associate-+l-55.9%
associate--r-55.9%
metadata-eval55.9%
+-commutative55.9%
*-lft-identity55.9%
*-commutative55.9%
distribute-lft-out--55.9%
Simplified55.9%
Taylor expanded in ux around 0 98.1%
fma-def98.1%
+-commutative98.1%
associate--l+98.2%
sub-neg98.2%
metadata-eval98.2%
distribute-lft-in98.2%
metadata-eval98.2%
+-commutative98.2%
mul-1-neg98.2%
sub-neg98.2%
count-298.2%
*-commutative98.2%
sub-neg98.2%
metadata-eval98.2%
+-commutative98.2%
Simplified98.2%
Final simplification98.2%
(FPCore (ux uy maxCos)
:precision binary32
(*
(sin (* 2.0 (* PI uy)))
(sqrt
(+
(* 2.0 (* ux (- 1.0 maxCos)))
(* (pow ux 2.0) (* (- 1.0 maxCos) (+ maxCos -1.0)))))))
float code(float ux, float uy, float maxCos) {
return sinf((2.0f * (((float) M_PI) * uy))) * sqrtf(((2.0f * (ux * (1.0f - maxCos))) + (powf(ux, 2.0f) * ((1.0f - maxCos) * (maxCos + -1.0f)))));
}
function code(ux, uy, maxCos) return Float32(sin(Float32(Float32(2.0) * Float32(Float32(pi) * uy))) * sqrt(Float32(Float32(Float32(2.0) * Float32(ux * Float32(Float32(1.0) - maxCos))) + Float32((ux ^ Float32(2.0)) * Float32(Float32(Float32(1.0) - maxCos) * Float32(maxCos + Float32(-1.0))))))) end
function tmp = code(ux, uy, maxCos) tmp = sin((single(2.0) * (single(pi) * uy))) * sqrt(((single(2.0) * (ux * (single(1.0) - maxCos))) + ((ux ^ single(2.0)) * ((single(1.0) - maxCos) * (maxCos + single(-1.0)))))); end
\begin{array}{l}
\\
\sin \left(2 \cdot \left(\pi \cdot uy\right)\right) \cdot \sqrt{2 \cdot \left(ux \cdot \left(1 - maxCos\right)\right) + {ux}^{2} \cdot \left(\left(1 - maxCos\right) \cdot \left(maxCos + -1\right)\right)}
\end{array}
Initial program 55.6%
associate-*l*55.6%
cancel-sign-sub-inv55.6%
+-commutative55.6%
*-commutative55.6%
fma-def55.8%
+-commutative55.8%
associate-+r-55.9%
fma-def55.9%
neg-sub055.9%
associate-+l-55.9%
associate--r-55.9%
metadata-eval55.9%
+-commutative55.9%
*-lft-identity55.9%
*-commutative55.9%
distribute-lft-out--55.9%
Simplified55.9%
Taylor expanded in ux around 0 98.1%
fma-def98.1%
+-commutative98.1%
associate--l+98.2%
sub-neg98.2%
metadata-eval98.2%
distribute-lft-in98.2%
metadata-eval98.2%
+-commutative98.2%
mul-1-neg98.2%
sub-neg98.2%
count-298.2%
*-commutative98.2%
sub-neg98.2%
metadata-eval98.2%
+-commutative98.2%
Simplified98.2%
Taylor expanded in uy around inf 98.2%
Final simplification98.2%
(FPCore (ux uy maxCos) :precision binary32 (* (sin (* 2.0 (* PI uy))) (sqrt (- (* ux 2.0) (pow ux 2.0)))))
float code(float ux, float uy, float maxCos) {
return sinf((2.0f * (((float) M_PI) * uy))) * sqrtf(((ux * 2.0f) - powf(ux, 2.0f)));
}
function code(ux, uy, maxCos) return Float32(sin(Float32(Float32(2.0) * Float32(Float32(pi) * uy))) * sqrt(Float32(Float32(ux * Float32(2.0)) - (ux ^ Float32(2.0))))) end
function tmp = code(ux, uy, maxCos) tmp = sin((single(2.0) * (single(pi) * uy))) * sqrt(((ux * single(2.0)) - (ux ^ single(2.0)))); end
\begin{array}{l}
\\
\sin \left(2 \cdot \left(\pi \cdot uy\right)\right) \cdot \sqrt{ux \cdot 2 - {ux}^{2}}
\end{array}
Initial program 55.6%
associate-*l*55.6%
cancel-sign-sub-inv55.6%
+-commutative55.6%
*-commutative55.6%
fma-def55.8%
+-commutative55.8%
associate-+r-55.9%
fma-def55.9%
neg-sub055.9%
associate-+l-55.9%
associate--r-55.9%
metadata-eval55.9%
+-commutative55.9%
*-lft-identity55.9%
*-commutative55.9%
distribute-lft-out--55.9%
Simplified55.9%
Taylor expanded in ux around 0 98.1%
fma-def98.1%
+-commutative98.1%
associate--l+98.2%
sub-neg98.2%
metadata-eval98.2%
distribute-lft-in98.2%
metadata-eval98.2%
+-commutative98.2%
mul-1-neg98.2%
sub-neg98.2%
count-298.2%
*-commutative98.2%
sub-neg98.2%
metadata-eval98.2%
+-commutative98.2%
Simplified98.2%
log1p-expm1-u98.2%
Applied egg-rr98.2%
Taylor expanded in maxCos around 0 93.8%
+-commutative93.8%
mul-1-neg93.8%
unsub-neg93.8%
Simplified93.8%
Final simplification93.8%
(FPCore (ux uy maxCos)
:precision binary32
(if (<=
(+ 1.0 (* (+ (- 1.0 ux) (* ux maxCos)) (- (+ ux -1.0) (* ux maxCos))))
0.00022000000171829015)
(* (sin (* uy (* 2.0 PI))) (sqrt (* ux (- (+ 1.0 (- 1.0 maxCos)) maxCos))))
(*
(sin (* 2.0 (* PI uy)))
(sqrt
(+
1.0
(* (+ -1.0 (* ux (- 1.0 maxCos))) (- (+ 1.0 (* ux maxCos)) ux)))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if ((1.0f + (((1.0f - ux) + (ux * maxCos)) * ((ux + -1.0f) - (ux * maxCos)))) <= 0.00022000000171829015f) {
tmp = sinf((uy * (2.0f * ((float) M_PI)))) * sqrtf((ux * ((1.0f + (1.0f - maxCos)) - maxCos)));
} else {
tmp = sinf((2.0f * (((float) M_PI) * uy))) * sqrtf((1.0f + ((-1.0f + (ux * (1.0f - maxCos))) * ((1.0f + (ux * maxCos)) - ux))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (Float32(Float32(1.0) + Float32(Float32(Float32(Float32(1.0) - ux) + Float32(ux * maxCos)) * Float32(Float32(ux + Float32(-1.0)) - Float32(ux * maxCos)))) <= Float32(0.00022000000171829015)) tmp = Float32(sin(Float32(uy * Float32(Float32(2.0) * Float32(pi)))) * sqrt(Float32(ux * Float32(Float32(Float32(1.0) + Float32(Float32(1.0) - maxCos)) - maxCos)))); else tmp = Float32(sin(Float32(Float32(2.0) * Float32(Float32(pi) * uy))) * sqrt(Float32(Float32(1.0) + Float32(Float32(Float32(-1.0) + Float32(ux * Float32(Float32(1.0) - maxCos))) * Float32(Float32(Float32(1.0) + Float32(ux * maxCos)) - ux))))); end return tmp end
function tmp_2 = code(ux, uy, maxCos) tmp = single(0.0); if ((single(1.0) + (((single(1.0) - ux) + (ux * maxCos)) * ((ux + single(-1.0)) - (ux * maxCos)))) <= single(0.00022000000171829015)) tmp = sin((uy * (single(2.0) * single(pi)))) * sqrt((ux * ((single(1.0) + (single(1.0) - maxCos)) - maxCos))); else tmp = sin((single(2.0) * (single(pi) * uy))) * sqrt((single(1.0) + ((single(-1.0) + (ux * (single(1.0) - maxCos))) * ((single(1.0) + (ux * maxCos)) - ux)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;1 + \left(\left(1 - ux\right) + ux \cdot maxCos\right) \cdot \left(\left(ux + -1\right) - ux \cdot maxCos\right) \leq 0.00022000000171829015:\\
\;\;\;\;\sin \left(uy \cdot \left(2 \cdot \pi\right)\right) \cdot \sqrt{ux \cdot \left(\left(1 + \left(1 - maxCos\right)\right) - maxCos\right)}\\
\mathbf{else}:\\
\;\;\;\;\sin \left(2 \cdot \left(\pi \cdot uy\right)\right) \cdot \sqrt{1 + \left(-1 + ux \cdot \left(1 - maxCos\right)\right) \cdot \left(\left(1 + ux \cdot maxCos\right) - ux\right)}\\
\end{array}
\end{array}
if (-.f32 1 (*.f32 (+.f32 (-.f32 1 ux) (*.f32 ux maxCos)) (+.f32 (-.f32 1 ux) (*.f32 ux maxCos)))) < 2.20000002e-4Initial program 32.8%
associate-*l*32.8%
cancel-sign-sub-inv32.8%
+-commutative32.8%
*-commutative32.8%
fma-def32.8%
+-commutative32.8%
associate-+r-33.0%
fma-def33.0%
neg-sub033.0%
associate-+l-33.0%
associate--r-33.0%
metadata-eval33.0%
+-commutative33.0%
*-lft-identity33.0%
*-commutative33.0%
distribute-lft-out--33.0%
Simplified33.0%
Taylor expanded in ux around 0 93.5%
mul-1-neg93.5%
sub-neg93.5%
metadata-eval93.5%
Simplified93.5%
if 2.20000002e-4 < (-.f32 1 (*.f32 (+.f32 (-.f32 1 ux) (*.f32 ux maxCos)) (+.f32 (-.f32 1 ux) (*.f32 ux maxCos)))) Initial program 89.0%
associate-*l*89.0%
cancel-sign-sub-inv89.0%
+-commutative89.0%
*-commutative89.0%
fma-def89.3%
+-commutative89.3%
associate-+r-89.5%
fma-def89.5%
neg-sub089.5%
associate-+l-89.4%
associate--r-89.4%
metadata-eval89.4%
+-commutative89.4%
*-lft-identity89.4%
*-commutative89.4%
distribute-lft-out--89.5%
Simplified89.5%
Taylor expanded in uy around inf 89.1%
Final simplification91.7%
(FPCore (ux uy maxCos) :precision binary32 (if (<= (* 2.0 uy) 0.0032999999821186066) (* 2.0 (* (* PI uy) (sqrt (- (* ux 2.0) (pow ux 2.0))))) (* (sin (* PI (* 2.0 uy))) (sqrt (* ux (- 2.0 (* 2.0 maxCos)))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if ((2.0f * uy) <= 0.0032999999821186066f) {
tmp = 2.0f * ((((float) M_PI) * uy) * sqrtf(((ux * 2.0f) - powf(ux, 2.0f))));
} else {
tmp = sinf((((float) M_PI) * (2.0f * uy))) * sqrtf((ux * (2.0f - (2.0f * maxCos))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (Float32(Float32(2.0) * uy) <= Float32(0.0032999999821186066)) tmp = Float32(Float32(2.0) * Float32(Float32(Float32(pi) * uy) * sqrt(Float32(Float32(ux * Float32(2.0)) - (ux ^ Float32(2.0)))))); else tmp = Float32(sin(Float32(Float32(pi) * Float32(Float32(2.0) * uy))) * sqrt(Float32(ux * Float32(Float32(2.0) - Float32(Float32(2.0) * maxCos))))); end return tmp end
function tmp_2 = code(ux, uy, maxCos) tmp = single(0.0); if ((single(2.0) * uy) <= single(0.0032999999821186066)) tmp = single(2.0) * ((single(pi) * uy) * sqrt(((ux * single(2.0)) - (ux ^ single(2.0))))); else tmp = sin((single(pi) * (single(2.0) * uy))) * sqrt((ux * (single(2.0) - (single(2.0) * maxCos)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;2 \cdot uy \leq 0.0032999999821186066:\\
\;\;\;\;2 \cdot \left(\left(\pi \cdot uy\right) \cdot \sqrt{ux \cdot 2 - {ux}^{2}}\right)\\
\mathbf{else}:\\
\;\;\;\;\sin \left(\pi \cdot \left(2 \cdot uy\right)\right) \cdot \sqrt{ux \cdot \left(2 - 2 \cdot maxCos\right)}\\
\end{array}
\end{array}
if (*.f32 uy 2) < 0.0033Initial program 54.1%
associate-*l*54.1%
cancel-sign-sub-inv54.1%
+-commutative54.1%
*-commutative54.1%
fma-def54.2%
+-commutative54.2%
associate-+r-54.4%
fma-def54.4%
neg-sub054.4%
associate-+l-54.4%
associate--r-54.4%
metadata-eval54.4%
+-commutative54.4%
*-lft-identity54.4%
*-commutative54.4%
distribute-lft-out--54.4%
Simplified54.4%
Taylor expanded in uy around 0 53.6%
Taylor expanded in ux around 0 97.0%
fma-def97.0%
associate--l+97.0%
mul-1-neg97.0%
sub-neg97.0%
metadata-eval97.0%
sub-neg97.0%
metadata-eval97.0%
Simplified97.0%
add-sqr-sqrt96.3%
pow296.3%
*-commutative96.3%
Applied egg-rr96.3%
Taylor expanded in maxCos around 0 93.0%
+-commutative93.0%
mul-1-neg93.0%
unsub-neg93.0%
Simplified93.0%
if 0.0033 < (*.f32 uy 2) Initial program 59.8%
Taylor expanded in ux around 0 73.5%
Final simplification87.7%
(FPCore (ux uy maxCos) :precision binary32 (if (<= ux 0.00011000000085914508) (* (sin (* PI (* 2.0 uy))) (sqrt (* ux (- 2.0 (* 2.0 maxCos))))) (* (sin (* uy (* 2.0 PI))) (sqrt (- 1.0 (* (- 1.0 ux) (- 1.0 ux)))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (ux <= 0.00011000000085914508f) {
tmp = sinf((((float) M_PI) * (2.0f * uy))) * sqrtf((ux * (2.0f - (2.0f * maxCos))));
} else {
tmp = sinf((uy * (2.0f * ((float) M_PI)))) * sqrtf((1.0f - ((1.0f - ux) * (1.0f - ux))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (ux <= Float32(0.00011000000085914508)) tmp = Float32(sin(Float32(Float32(pi) * Float32(Float32(2.0) * uy))) * sqrt(Float32(ux * Float32(Float32(2.0) - Float32(Float32(2.0) * maxCos))))); else tmp = Float32(sin(Float32(uy * Float32(Float32(2.0) * Float32(pi)))) * sqrt(Float32(Float32(1.0) - Float32(Float32(Float32(1.0) - ux) * Float32(Float32(1.0) - ux))))); end return tmp end
function tmp_2 = code(ux, uy, maxCos) tmp = single(0.0); if (ux <= single(0.00011000000085914508)) tmp = sin((single(pi) * (single(2.0) * uy))) * sqrt((ux * (single(2.0) - (single(2.0) * maxCos)))); else tmp = sin((uy * (single(2.0) * single(pi)))) * sqrt((single(1.0) - ((single(1.0) - ux) * (single(1.0) - ux)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;ux \leq 0.00011000000085914508:\\
\;\;\;\;\sin \left(\pi \cdot \left(2 \cdot uy\right)\right) \cdot \sqrt{ux \cdot \left(2 - 2 \cdot maxCos\right)}\\
\mathbf{else}:\\
\;\;\;\;\sin \left(uy \cdot \left(2 \cdot \pi\right)\right) \cdot \sqrt{1 - \left(1 - ux\right) \cdot \left(1 - ux\right)}\\
\end{array}
\end{array}
if ux < 1.10000001e-4Initial program 32.6%
Taylor expanded in ux around 0 93.5%
if 1.10000001e-4 < ux Initial program 88.7%
associate-*l*88.7%
cancel-sign-sub-inv88.7%
+-commutative88.7%
*-commutative88.7%
fma-def89.0%
+-commutative89.0%
associate-+r-89.2%
fma-def89.2%
neg-sub089.2%
associate-+l-89.2%
associate--r-89.2%
metadata-eval89.2%
+-commutative89.2%
*-lft-identity89.2%
*-commutative89.2%
distribute-lft-out--89.2%
Simplified89.2%
Taylor expanded in maxCos around 0 84.2%
Final simplification89.7%
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0 (sin (* uy (* 2.0 PI)))))
(if (<= ux 0.00011000000085914508)
(* t_0 (sqrt (* ux (- (+ 1.0 (- 1.0 maxCos)) maxCos))))
(* t_0 (sqrt (- 1.0 (* (- 1.0 ux) (- 1.0 ux))))))))
float code(float ux, float uy, float maxCos) {
float t_0 = sinf((uy * (2.0f * ((float) M_PI))));
float tmp;
if (ux <= 0.00011000000085914508f) {
tmp = t_0 * sqrtf((ux * ((1.0f + (1.0f - maxCos)) - maxCos)));
} else {
tmp = t_0 * sqrtf((1.0f - ((1.0f - ux) * (1.0f - ux))));
}
return tmp;
}
function code(ux, uy, maxCos) t_0 = sin(Float32(uy * Float32(Float32(2.0) * Float32(pi)))) tmp = Float32(0.0) if (ux <= Float32(0.00011000000085914508)) tmp = Float32(t_0 * sqrt(Float32(ux * Float32(Float32(Float32(1.0) + Float32(Float32(1.0) - maxCos)) - maxCos)))); else tmp = Float32(t_0 * sqrt(Float32(Float32(1.0) - Float32(Float32(Float32(1.0) - ux) * Float32(Float32(1.0) - ux))))); end return tmp end
function tmp_2 = code(ux, uy, maxCos) t_0 = sin((uy * (single(2.0) * single(pi)))); tmp = single(0.0); if (ux <= single(0.00011000000085914508)) tmp = t_0 * sqrt((ux * ((single(1.0) + (single(1.0) - maxCos)) - maxCos))); else tmp = t_0 * sqrt((single(1.0) - ((single(1.0) - ux) * (single(1.0) - ux)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sin \left(uy \cdot \left(2 \cdot \pi\right)\right)\\
\mathbf{if}\;ux \leq 0.00011000000085914508:\\
\;\;\;\;t_0 \cdot \sqrt{ux \cdot \left(\left(1 + \left(1 - maxCos\right)\right) - maxCos\right)}\\
\mathbf{else}:\\
\;\;\;\;t_0 \cdot \sqrt{1 - \left(1 - ux\right) \cdot \left(1 - ux\right)}\\
\end{array}
\end{array}
if ux < 1.10000001e-4Initial program 32.6%
associate-*l*32.6%
cancel-sign-sub-inv32.6%
+-commutative32.6%
*-commutative32.6%
fma-def32.6%
+-commutative32.6%
associate-+r-32.8%
fma-def32.8%
neg-sub032.8%
associate-+l-32.8%
associate--r-32.8%
metadata-eval32.8%
+-commutative32.8%
*-lft-identity32.8%
*-commutative32.8%
distribute-lft-out--32.8%
Simplified32.8%
Taylor expanded in ux around 0 93.6%
mul-1-neg93.6%
sub-neg93.6%
metadata-eval93.6%
Simplified93.6%
if 1.10000001e-4 < ux Initial program 88.7%
associate-*l*88.7%
cancel-sign-sub-inv88.7%
+-commutative88.7%
*-commutative88.7%
fma-def89.0%
+-commutative89.0%
associate-+r-89.2%
fma-def89.2%
neg-sub089.2%
associate-+l-89.2%
associate--r-89.2%
metadata-eval89.2%
+-commutative89.2%
*-lft-identity89.2%
*-commutative89.2%
distribute-lft-out--89.2%
Simplified89.2%
Taylor expanded in maxCos around 0 84.2%
Final simplification89.7%
(FPCore (ux uy maxCos) :precision binary32 (* 2.0 (* uy (* PI (sqrt (- (* ux 2.0) (pow ux 2.0)))))))
float code(float ux, float uy, float maxCos) {
return 2.0f * (uy * (((float) M_PI) * sqrtf(((ux * 2.0f) - powf(ux, 2.0f)))));
}
function code(ux, uy, maxCos) return Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * sqrt(Float32(Float32(ux * Float32(2.0)) - (ux ^ Float32(2.0))))))) end
function tmp = code(ux, uy, maxCos) tmp = single(2.0) * (uy * (single(pi) * sqrt(((ux * single(2.0)) - (ux ^ single(2.0)))))); end
\begin{array}{l}
\\
2 \cdot \left(uy \cdot \left(\pi \cdot \sqrt{ux \cdot 2 - {ux}^{2}}\right)\right)
\end{array}
Initial program 55.6%
associate-*l*55.6%
cancel-sign-sub-inv55.6%
+-commutative55.6%
*-commutative55.6%
fma-def55.8%
+-commutative55.8%
associate-+r-55.9%
fma-def55.9%
neg-sub055.9%
associate-+l-55.9%
associate--r-55.9%
metadata-eval55.9%
+-commutative55.9%
*-lft-identity55.9%
*-commutative55.9%
distribute-lft-out--55.9%
Simplified55.9%
Taylor expanded in uy around 0 46.7%
Taylor expanded in maxCos around 0 45.5%
associate-*l*45.5%
sub-neg45.5%
metadata-eval45.5%
Simplified45.5%
Taylor expanded in ux around 0 77.3%
+-commutative77.3%
mul-1-neg77.3%
unsub-neg77.3%
Simplified77.3%
Final simplification77.3%
(FPCore (ux uy maxCos) :precision binary32 (* 2.0 (* (* PI uy) (sqrt (- (* ux 2.0) (pow ux 2.0))))))
float code(float ux, float uy, float maxCos) {
return 2.0f * ((((float) M_PI) * uy) * sqrtf(((ux * 2.0f) - powf(ux, 2.0f))));
}
function code(ux, uy, maxCos) return Float32(Float32(2.0) * Float32(Float32(Float32(pi) * uy) * sqrt(Float32(Float32(ux * Float32(2.0)) - (ux ^ Float32(2.0)))))) end
function tmp = code(ux, uy, maxCos) tmp = single(2.0) * ((single(pi) * uy) * sqrt(((ux * single(2.0)) - (ux ^ single(2.0))))); end
\begin{array}{l}
\\
2 \cdot \left(\left(\pi \cdot uy\right) \cdot \sqrt{ux \cdot 2 - {ux}^{2}}\right)
\end{array}
Initial program 55.6%
associate-*l*55.6%
cancel-sign-sub-inv55.6%
+-commutative55.6%
*-commutative55.6%
fma-def55.8%
+-commutative55.8%
associate-+r-55.9%
fma-def55.9%
neg-sub055.9%
associate-+l-55.9%
associate--r-55.9%
metadata-eval55.9%
+-commutative55.9%
*-lft-identity55.9%
*-commutative55.9%
distribute-lft-out--55.9%
Simplified55.9%
Taylor expanded in uy around 0 46.7%
Taylor expanded in ux around 0 80.2%
fma-def80.2%
associate--l+80.2%
mul-1-neg80.2%
sub-neg80.2%
metadata-eval80.2%
sub-neg80.2%
metadata-eval80.2%
Simplified80.2%
add-sqr-sqrt79.7%
pow279.7%
*-commutative79.7%
Applied egg-rr79.7%
Taylor expanded in maxCos around 0 77.3%
+-commutative77.3%
mul-1-neg77.3%
unsub-neg77.3%
Simplified77.3%
Final simplification77.3%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= ux 0.00011000000085914508)
(* 2.0 (* (* PI uy) (sqrt (* ux (- (+ 1.0 (- 1.0 maxCos)) maxCos)))))
(*
2.0
(*
(* PI uy)
(sqrt
(+
1.0
(* (+ -1.0 (* ux (- 1.0 maxCos))) (- (+ 1.0 (* ux maxCos)) ux))))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (ux <= 0.00011000000085914508f) {
tmp = 2.0f * ((((float) M_PI) * uy) * sqrtf((ux * ((1.0f + (1.0f - maxCos)) - maxCos))));
} else {
tmp = 2.0f * ((((float) M_PI) * uy) * sqrtf((1.0f + ((-1.0f + (ux * (1.0f - maxCos))) * ((1.0f + (ux * maxCos)) - ux)))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (ux <= Float32(0.00011000000085914508)) tmp = Float32(Float32(2.0) * Float32(Float32(Float32(pi) * uy) * sqrt(Float32(ux * Float32(Float32(Float32(1.0) + Float32(Float32(1.0) - maxCos)) - maxCos))))); else tmp = Float32(Float32(2.0) * Float32(Float32(Float32(pi) * uy) * sqrt(Float32(Float32(1.0) + Float32(Float32(Float32(-1.0) + Float32(ux * Float32(Float32(1.0) - maxCos))) * Float32(Float32(Float32(1.0) + Float32(ux * maxCos)) - ux)))))); end return tmp end
function tmp_2 = code(ux, uy, maxCos) tmp = single(0.0); if (ux <= single(0.00011000000085914508)) tmp = single(2.0) * ((single(pi) * uy) * sqrt((ux * ((single(1.0) + (single(1.0) - maxCos)) - maxCos)))); else tmp = single(2.0) * ((single(pi) * uy) * sqrt((single(1.0) + ((single(-1.0) + (ux * (single(1.0) - maxCos))) * ((single(1.0) + (ux * maxCos)) - ux))))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;ux \leq 0.00011000000085914508:\\
\;\;\;\;2 \cdot \left(\left(\pi \cdot uy\right) \cdot \sqrt{ux \cdot \left(\left(1 + \left(1 - maxCos\right)\right) - maxCos\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left(\pi \cdot uy\right) \cdot \sqrt{1 + \left(-1 + ux \cdot \left(1 - maxCos\right)\right) \cdot \left(\left(1 + ux \cdot maxCos\right) - ux\right)}\right)\\
\end{array}
\end{array}
if ux < 1.10000001e-4Initial program 32.6%
associate-*l*32.6%
cancel-sign-sub-inv32.6%
+-commutative32.6%
*-commutative32.6%
fma-def32.6%
+-commutative32.6%
associate-+r-32.8%
fma-def32.8%
neg-sub032.8%
associate-+l-32.8%
associate--r-32.8%
metadata-eval32.8%
+-commutative32.8%
*-lft-identity32.8%
*-commutative32.8%
distribute-lft-out--32.8%
Simplified32.8%
Taylor expanded in uy around 0 29.7%
Taylor expanded in ux around 0 79.9%
if 1.10000001e-4 < ux Initial program 88.7%
associate-*l*88.7%
cancel-sign-sub-inv88.7%
+-commutative88.7%
*-commutative88.7%
fma-def89.0%
+-commutative89.0%
associate-+r-89.2%
fma-def89.2%
neg-sub089.2%
associate-+l-89.2%
associate--r-89.2%
metadata-eval89.2%
+-commutative89.2%
*-lft-identity89.2%
*-commutative89.2%
distribute-lft-out--89.2%
Simplified89.2%
Taylor expanded in uy around 0 71.2%
Final simplification76.3%
(FPCore (ux uy maxCos) :precision binary32 (if (<= ux 7.000000186963007e-5) (* 2.0 (* (* PI uy) (sqrt (* ux (- (+ 1.0 (- 1.0 maxCos)) maxCos))))) (* 2.0 (* uy (* PI (sqrt (+ 1.0 (- (* ux (- 1.0 ux)) (- 1.0 ux)))))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (ux <= 7.000000186963007e-5f) {
tmp = 2.0f * ((((float) M_PI) * uy) * sqrtf((ux * ((1.0f + (1.0f - maxCos)) - maxCos))));
} else {
tmp = 2.0f * (uy * (((float) M_PI) * sqrtf((1.0f + ((ux * (1.0f - ux)) - (1.0f - ux))))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (ux <= Float32(7.000000186963007e-5)) tmp = Float32(Float32(2.0) * Float32(Float32(Float32(pi) * uy) * sqrt(Float32(ux * Float32(Float32(Float32(1.0) + Float32(Float32(1.0) - maxCos)) - maxCos))))); else tmp = Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * sqrt(Float32(Float32(1.0) + Float32(Float32(ux * Float32(Float32(1.0) - ux)) - Float32(Float32(1.0) - ux))))))); end return tmp end
function tmp_2 = code(ux, uy, maxCos) tmp = single(0.0); if (ux <= single(7.000000186963007e-5)) tmp = single(2.0) * ((single(pi) * uy) * sqrt((ux * ((single(1.0) + (single(1.0) - maxCos)) - maxCos)))); else tmp = single(2.0) * (uy * (single(pi) * sqrt((single(1.0) + ((ux * (single(1.0) - ux)) - (single(1.0) - ux)))))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;ux \leq 7.000000186963007 \cdot 10^{-5}:\\
\;\;\;\;2 \cdot \left(\left(\pi \cdot uy\right) \cdot \sqrt{ux \cdot \left(\left(1 + \left(1 - maxCos\right)\right) - maxCos\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(uy \cdot \left(\pi \cdot \sqrt{1 + \left(ux \cdot \left(1 - ux\right) - \left(1 - ux\right)\right)}\right)\right)\\
\end{array}
\end{array}
if ux < 7.00000019e-5Initial program 31.0%
associate-*l*31.0%
cancel-sign-sub-inv31.0%
+-commutative31.0%
*-commutative31.0%
fma-def31.1%
+-commutative31.1%
associate-+r-31.2%
fma-def31.2%
neg-sub031.2%
associate-+l-31.2%
associate--r-31.2%
metadata-eval31.2%
+-commutative31.2%
*-lft-identity31.2%
*-commutative31.2%
distribute-lft-out--31.2%
Simplified31.2%
Taylor expanded in uy around 0 28.4%
Taylor expanded in ux around 0 80.6%
if 7.00000019e-5 < ux Initial program 87.8%
associate-*l*87.8%
cancel-sign-sub-inv87.8%
+-commutative87.8%
*-commutative87.8%
fma-def88.1%
+-commutative88.1%
associate-+r-88.2%
fma-def88.2%
neg-sub088.2%
associate-+l-88.2%
associate--r-88.2%
metadata-eval88.2%
+-commutative88.2%
*-lft-identity88.2%
*-commutative88.2%
distribute-lft-out--88.2%
Simplified88.2%
Taylor expanded in uy around 0 70.6%
Taylor expanded in maxCos around 0 67.9%
associate-*l*67.9%
sub-neg67.9%
metadata-eval67.9%
Simplified67.9%
*-commutative67.9%
sub-neg67.9%
distribute-lft-in68.6%
*-commutative68.6%
*-un-lft-identity68.6%
Applied egg-rr68.6%
Final simplification75.4%
(FPCore (ux uy maxCos) :precision binary32 (if (<= ux 0.00011000000085914508) (* 2.0 (* uy (* PI (sqrt (* ux 2.0))))) (* 2.0 (* uy (* PI (sqrt (- 1.0 (* (- 1.0 ux) (- 1.0 ux)))))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (ux <= 0.00011000000085914508f) {
tmp = 2.0f * (uy * (((float) M_PI) * sqrtf((ux * 2.0f))));
} else {
tmp = 2.0f * (uy * (((float) M_PI) * sqrtf((1.0f - ((1.0f - ux) * (1.0f - ux))))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (ux <= Float32(0.00011000000085914508)) tmp = Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * sqrt(Float32(ux * Float32(2.0)))))); else tmp = Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * sqrt(Float32(Float32(1.0) - Float32(Float32(Float32(1.0) - ux) * Float32(Float32(1.0) - ux))))))); end return tmp end
function tmp_2 = code(ux, uy, maxCos) tmp = single(0.0); if (ux <= single(0.00011000000085914508)) tmp = single(2.0) * (uy * (single(pi) * sqrt((ux * single(2.0))))); else tmp = single(2.0) * (uy * (single(pi) * sqrt((single(1.0) - ((single(1.0) - ux) * (single(1.0) - ux)))))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;ux \leq 0.00011000000085914508:\\
\;\;\;\;2 \cdot \left(uy \cdot \left(\pi \cdot \sqrt{ux \cdot 2}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(uy \cdot \left(\pi \cdot \sqrt{1 - \left(1 - ux\right) \cdot \left(1 - ux\right)}\right)\right)\\
\end{array}
\end{array}
if ux < 1.10000001e-4Initial program 32.6%
associate-*l*32.6%
cancel-sign-sub-inv32.6%
+-commutative32.6%
*-commutative32.6%
fma-def32.6%
+-commutative32.6%
associate-+r-32.8%
fma-def32.8%
neg-sub032.8%
associate-+l-32.8%
associate--r-32.8%
metadata-eval32.8%
+-commutative32.8%
*-lft-identity32.8%
*-commutative32.8%
distribute-lft-out--32.8%
Simplified32.8%
Taylor expanded in uy around 0 29.7%
Taylor expanded in maxCos around 0 29.7%
associate-*l*29.7%
sub-neg29.7%
metadata-eval29.7%
Simplified29.7%
Taylor expanded in ux around 0 77.7%
if 1.10000001e-4 < ux Initial program 88.7%
associate-*l*88.7%
cancel-sign-sub-inv88.7%
+-commutative88.7%
*-commutative88.7%
fma-def89.0%
+-commutative89.0%
associate-+r-89.2%
fma-def89.2%
neg-sub089.2%
associate-+l-89.2%
associate--r-89.2%
metadata-eval89.2%
+-commutative89.2%
*-lft-identity89.2%
*-commutative89.2%
distribute-lft-out--89.2%
Simplified89.2%
Taylor expanded in uy around 0 71.2%
Taylor expanded in maxCos around 0 68.3%
associate-*l*68.2%
sub-neg68.2%
metadata-eval68.2%
Simplified68.2%
Final simplification73.8%
(FPCore (ux uy maxCos) :precision binary32 (if (<= ux 0.00011000000085914508) (* 2.0 (* uy (* PI (sqrt (* ux 2.0))))) (* 2.0 (* (* PI uy) (sqrt (- 1.0 (* (- 1.0 ux) (- 1.0 ux))))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (ux <= 0.00011000000085914508f) {
tmp = 2.0f * (uy * (((float) M_PI) * sqrtf((ux * 2.0f))));
} else {
tmp = 2.0f * ((((float) M_PI) * uy) * sqrtf((1.0f - ((1.0f - ux) * (1.0f - ux)))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (ux <= Float32(0.00011000000085914508)) tmp = Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * sqrt(Float32(ux * Float32(2.0)))))); else tmp = Float32(Float32(2.0) * Float32(Float32(Float32(pi) * uy) * sqrt(Float32(Float32(1.0) - Float32(Float32(Float32(1.0) - ux) * Float32(Float32(1.0) - ux)))))); end return tmp end
function tmp_2 = code(ux, uy, maxCos) tmp = single(0.0); if (ux <= single(0.00011000000085914508)) tmp = single(2.0) * (uy * (single(pi) * sqrt((ux * single(2.0))))); else tmp = single(2.0) * ((single(pi) * uy) * sqrt((single(1.0) - ((single(1.0) - ux) * (single(1.0) - ux))))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;ux \leq 0.00011000000085914508:\\
\;\;\;\;2 \cdot \left(uy \cdot \left(\pi \cdot \sqrt{ux \cdot 2}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left(\pi \cdot uy\right) \cdot \sqrt{1 - \left(1 - ux\right) \cdot \left(1 - ux\right)}\right)\\
\end{array}
\end{array}
if ux < 1.10000001e-4Initial program 32.6%
associate-*l*32.6%
cancel-sign-sub-inv32.6%
+-commutative32.6%
*-commutative32.6%
fma-def32.6%
+-commutative32.6%
associate-+r-32.8%
fma-def32.8%
neg-sub032.8%
associate-+l-32.8%
associate--r-32.8%
metadata-eval32.8%
+-commutative32.8%
*-lft-identity32.8%
*-commutative32.8%
distribute-lft-out--32.8%
Simplified32.8%
Taylor expanded in uy around 0 29.7%
Taylor expanded in maxCos around 0 29.7%
associate-*l*29.7%
sub-neg29.7%
metadata-eval29.7%
Simplified29.7%
Taylor expanded in ux around 0 77.7%
if 1.10000001e-4 < ux Initial program 88.7%
associate-*l*88.7%
cancel-sign-sub-inv88.7%
+-commutative88.7%
*-commutative88.7%
fma-def89.0%
+-commutative89.0%
associate-+r-89.2%
fma-def89.2%
neg-sub089.2%
associate-+l-89.2%
associate--r-89.2%
metadata-eval89.2%
+-commutative89.2%
*-lft-identity89.2%
*-commutative89.2%
distribute-lft-out--89.2%
Simplified89.2%
Taylor expanded in uy around 0 71.2%
Taylor expanded in maxCos around 0 68.3%
Final simplification73.8%
(FPCore (ux uy maxCos) :precision binary32 (if (<= ux 0.00011000000085914508) (* 2.0 (* (* PI uy) (sqrt (* ux (+ 1.0 (- (- 1.0 maxCos) maxCos)))))) (* 2.0 (* (* PI uy) (sqrt (- 1.0 (* (- 1.0 ux) (- 1.0 ux))))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (ux <= 0.00011000000085914508f) {
tmp = 2.0f * ((((float) M_PI) * uy) * sqrtf((ux * (1.0f + ((1.0f - maxCos) - maxCos)))));
} else {
tmp = 2.0f * ((((float) M_PI) * uy) * sqrtf((1.0f - ((1.0f - ux) * (1.0f - ux)))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (ux <= Float32(0.00011000000085914508)) tmp = Float32(Float32(2.0) * Float32(Float32(Float32(pi) * uy) * sqrt(Float32(ux * Float32(Float32(1.0) + Float32(Float32(Float32(1.0) - maxCos) - maxCos)))))); else tmp = Float32(Float32(2.0) * Float32(Float32(Float32(pi) * uy) * sqrt(Float32(Float32(1.0) - Float32(Float32(Float32(1.0) - ux) * Float32(Float32(1.0) - ux)))))); end return tmp end
function tmp_2 = code(ux, uy, maxCos) tmp = single(0.0); if (ux <= single(0.00011000000085914508)) tmp = single(2.0) * ((single(pi) * uy) * sqrt((ux * (single(1.0) + ((single(1.0) - maxCos) - maxCos))))); else tmp = single(2.0) * ((single(pi) * uy) * sqrt((single(1.0) - ((single(1.0) - ux) * (single(1.0) - ux))))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;ux \leq 0.00011000000085914508:\\
\;\;\;\;2 \cdot \left(\left(\pi \cdot uy\right) \cdot \sqrt{ux \cdot \left(1 + \left(\left(1 - maxCos\right) - maxCos\right)\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left(\pi \cdot uy\right) \cdot \sqrt{1 - \left(1 - ux\right) \cdot \left(1 - ux\right)}\right)\\
\end{array}
\end{array}
if ux < 1.10000001e-4Initial program 32.6%
associate-*l*32.6%
cancel-sign-sub-inv32.6%
+-commutative32.6%
*-commutative32.6%
fma-def32.6%
+-commutative32.6%
associate-+r-32.8%
fma-def32.8%
neg-sub032.8%
associate-+l-32.8%
associate--r-32.8%
metadata-eval32.8%
+-commutative32.8%
*-lft-identity32.8%
*-commutative32.8%
distribute-lft-out--32.8%
Simplified32.8%
Taylor expanded in uy around 0 29.7%
Taylor expanded in ux around 0 79.9%
associate--l+79.8%
mul-1-neg79.8%
sub-neg79.8%
metadata-eval79.8%
Simplified79.8%
if 1.10000001e-4 < ux Initial program 88.7%
associate-*l*88.7%
cancel-sign-sub-inv88.7%
+-commutative88.7%
*-commutative88.7%
fma-def89.0%
+-commutative89.0%
associate-+r-89.2%
fma-def89.2%
neg-sub089.2%
associate-+l-89.2%
associate--r-89.2%
metadata-eval89.2%
+-commutative89.2%
*-lft-identity89.2%
*-commutative89.2%
distribute-lft-out--89.2%
Simplified89.2%
Taylor expanded in uy around 0 71.2%
Taylor expanded in maxCos around 0 68.3%
Final simplification75.1%
(FPCore (ux uy maxCos) :precision binary32 (if (<= ux 0.00011000000085914508) (* 2.0 (* (* PI uy) (sqrt (* ux (- (+ 1.0 (- 1.0 maxCos)) maxCos))))) (* 2.0 (* (* PI uy) (sqrt (- 1.0 (* (- 1.0 ux) (- 1.0 ux))))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (ux <= 0.00011000000085914508f) {
tmp = 2.0f * ((((float) M_PI) * uy) * sqrtf((ux * ((1.0f + (1.0f - maxCos)) - maxCos))));
} else {
tmp = 2.0f * ((((float) M_PI) * uy) * sqrtf((1.0f - ((1.0f - ux) * (1.0f - ux)))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (ux <= Float32(0.00011000000085914508)) tmp = Float32(Float32(2.0) * Float32(Float32(Float32(pi) * uy) * sqrt(Float32(ux * Float32(Float32(Float32(1.0) + Float32(Float32(1.0) - maxCos)) - maxCos))))); else tmp = Float32(Float32(2.0) * Float32(Float32(Float32(pi) * uy) * sqrt(Float32(Float32(1.0) - Float32(Float32(Float32(1.0) - ux) * Float32(Float32(1.0) - ux)))))); end return tmp end
function tmp_2 = code(ux, uy, maxCos) tmp = single(0.0); if (ux <= single(0.00011000000085914508)) tmp = single(2.0) * ((single(pi) * uy) * sqrt((ux * ((single(1.0) + (single(1.0) - maxCos)) - maxCos)))); else tmp = single(2.0) * ((single(pi) * uy) * sqrt((single(1.0) - ((single(1.0) - ux) * (single(1.0) - ux))))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;ux \leq 0.00011000000085914508:\\
\;\;\;\;2 \cdot \left(\left(\pi \cdot uy\right) \cdot \sqrt{ux \cdot \left(\left(1 + \left(1 - maxCos\right)\right) - maxCos\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left(\pi \cdot uy\right) \cdot \sqrt{1 - \left(1 - ux\right) \cdot \left(1 - ux\right)}\right)\\
\end{array}
\end{array}
if ux < 1.10000001e-4Initial program 32.6%
associate-*l*32.6%
cancel-sign-sub-inv32.6%
+-commutative32.6%
*-commutative32.6%
fma-def32.6%
+-commutative32.6%
associate-+r-32.8%
fma-def32.8%
neg-sub032.8%
associate-+l-32.8%
associate--r-32.8%
metadata-eval32.8%
+-commutative32.8%
*-lft-identity32.8%
*-commutative32.8%
distribute-lft-out--32.8%
Simplified32.8%
Taylor expanded in uy around 0 29.7%
Taylor expanded in ux around 0 79.9%
if 1.10000001e-4 < ux Initial program 88.7%
associate-*l*88.7%
cancel-sign-sub-inv88.7%
+-commutative88.7%
*-commutative88.7%
fma-def89.0%
+-commutative89.0%
associate-+r-89.2%
fma-def89.2%
neg-sub089.2%
associate-+l-89.2%
associate--r-89.2%
metadata-eval89.2%
+-commutative89.2%
*-lft-identity89.2%
*-commutative89.2%
distribute-lft-out--89.2%
Simplified89.2%
Taylor expanded in uy around 0 71.2%
Taylor expanded in maxCos around 0 68.3%
Final simplification75.1%
(FPCore (ux uy maxCos) :precision binary32 (* 2.0 (* uy (* PI (sqrt (* ux 2.0))))))
float code(float ux, float uy, float maxCos) {
return 2.0f * (uy * (((float) M_PI) * sqrtf((ux * 2.0f))));
}
function code(ux, uy, maxCos) return Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * sqrt(Float32(ux * Float32(2.0)))))) end
function tmp = code(ux, uy, maxCos) tmp = single(2.0) * (uy * (single(pi) * sqrt((ux * single(2.0))))); end
\begin{array}{l}
\\
2 \cdot \left(uy \cdot \left(\pi \cdot \sqrt{ux \cdot 2}\right)\right)
\end{array}
Initial program 55.6%
associate-*l*55.6%
cancel-sign-sub-inv55.6%
+-commutative55.6%
*-commutative55.6%
fma-def55.8%
+-commutative55.8%
associate-+r-55.9%
fma-def55.9%
neg-sub055.9%
associate-+l-55.9%
associate--r-55.9%
metadata-eval55.9%
+-commutative55.9%
*-lft-identity55.9%
*-commutative55.9%
distribute-lft-out--55.9%
Simplified55.9%
Taylor expanded in uy around 0 46.7%
Taylor expanded in maxCos around 0 45.5%
associate-*l*45.5%
sub-neg45.5%
metadata-eval45.5%
Simplified45.5%
Taylor expanded in ux around 0 64.3%
Final simplification64.3%
herbie shell --seed 2023336
(FPCore (ux uy maxCos)
:name "UniformSampleCone, y"
:precision binary32
:pre (and (and (and (<= 2.328306437e-10 ux) (<= ux 1.0)) (and (<= 2.328306437e-10 uy) (<= uy 1.0))) (and (<= 0.0 maxCos) (<= maxCos 1.0)))
(* (sin (* (* uy 2.0) PI)) (sqrt (- 1.0 (* (+ (- 1.0 ux) (* ux maxCos)) (+ (- 1.0 ux) (* ux maxCos)))))))