
(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 16 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
(*
(sin (cbrt (* (pow (* uy 2.0) 3.0) (pow PI 3.0))))
(sqrt
(fma
ux
(+ 1.0 (- (- 1.0 maxCos) maxCos))
(* (- 1.0 maxCos) (* (+ maxCos -1.0) (* ux ux)))))))
float code(float ux, float uy, float maxCos) {
return sinf(cbrtf((powf((uy * 2.0f), 3.0f) * powf(((float) M_PI), 3.0f)))) * sqrtf(fmaf(ux, (1.0f + ((1.0f - maxCos) - maxCos)), ((1.0f - maxCos) * ((maxCos + -1.0f) * (ux * ux)))));
}
function code(ux, uy, maxCos) return Float32(sin(cbrt(Float32((Float32(uy * Float32(2.0)) ^ Float32(3.0)) * (Float32(pi) ^ Float32(3.0))))) * sqrt(fma(ux, Float32(Float32(1.0) + Float32(Float32(Float32(1.0) - maxCos) - maxCos)), Float32(Float32(Float32(1.0) - maxCos) * Float32(Float32(maxCos + Float32(-1.0)) * Float32(ux * ux)))))) end
\begin{array}{l}
\\
\sin \left(\sqrt[3]{{\left(uy \cdot 2\right)}^{3} \cdot {\pi}^{3}}\right) \cdot \sqrt{\mathsf{fma}\left(ux, 1 + \left(\left(1 - maxCos\right) - maxCos\right), \left(1 - maxCos\right) \cdot \left(\left(maxCos + -1\right) \cdot \left(ux \cdot ux\right)\right)\right)}
\end{array}
Initial program 55.7%
associate-*l*55.7%
sub-neg55.7%
+-commutative55.7%
distribute-rgt-neg-in55.7%
fma-def55.9%
+-commutative55.9%
associate-+r-55.8%
fma-def55.9%
neg-sub055.9%
+-commutative55.9%
associate-+r-55.7%
associate--r-55.7%
neg-sub055.7%
+-commutative55.7%
sub-neg55.7%
fma-def55.7%
Simplified55.7%
Taylor expanded in ux around 0 98.2%
+-commutative98.2%
fma-def98.2%
associate--l+98.3%
mul-1-neg98.3%
sub-neg98.3%
metadata-eval98.3%
distribute-neg-in98.3%
metadata-eval98.3%
+-commutative98.3%
sub-neg98.3%
sub-neg98.3%
metadata-eval98.3%
*-commutative98.3%
unpow298.3%
Simplified98.3%
associate-*r*98.3%
add-cbrt-cube98.3%
add-cbrt-cube98.3%
cbrt-unprod98.3%
pow398.4%
pow398.4%
Applied egg-rr98.4%
Final simplification98.4%
(FPCore (ux uy maxCos)
:precision binary32
(*
(sin (* uy (* 2.0 PI)))
(pow
(fma
ux
(+ 1.0 (- 1.0 (+ maxCos maxCos)))
(* (- 1.0 maxCos) (* (+ maxCos -1.0) (* ux ux))))
0.5)))
float code(float ux, float uy, float maxCos) {
return sinf((uy * (2.0f * ((float) M_PI)))) * powf(fmaf(ux, (1.0f + (1.0f - (maxCos + maxCos))), ((1.0f - maxCos) * ((maxCos + -1.0f) * (ux * ux)))), 0.5f);
}
function code(ux, uy, maxCos) return Float32(sin(Float32(uy * Float32(Float32(2.0) * Float32(pi)))) * (fma(ux, Float32(Float32(1.0) + Float32(Float32(1.0) - Float32(maxCos + maxCos))), Float32(Float32(Float32(1.0) - maxCos) * Float32(Float32(maxCos + Float32(-1.0)) * Float32(ux * ux)))) ^ Float32(0.5))) end
\begin{array}{l}
\\
\sin \left(uy \cdot \left(2 \cdot \pi\right)\right) \cdot {\left(\mathsf{fma}\left(ux, 1 + \left(1 - \left(maxCos + maxCos\right)\right), \left(1 - maxCos\right) \cdot \left(\left(maxCos + -1\right) \cdot \left(ux \cdot ux\right)\right)\right)\right)}^{0.5}
\end{array}
Initial program 55.7%
associate-*l*55.7%
sub-neg55.7%
+-commutative55.7%
distribute-rgt-neg-in55.7%
fma-def55.9%
+-commutative55.9%
associate-+r-55.8%
fma-def55.9%
neg-sub055.9%
+-commutative55.9%
associate-+r-55.7%
associate--r-55.7%
neg-sub055.7%
+-commutative55.7%
sub-neg55.7%
fma-def55.7%
Simplified55.7%
Taylor expanded in ux around 0 98.2%
+-commutative98.2%
fma-def98.2%
associate--l+98.3%
mul-1-neg98.3%
sub-neg98.3%
metadata-eval98.3%
distribute-neg-in98.3%
metadata-eval98.3%
+-commutative98.3%
sub-neg98.3%
sub-neg98.3%
metadata-eval98.3%
*-commutative98.3%
unpow298.3%
Simplified98.3%
pow1/298.3%
associate--l-98.3%
associate-*r*98.3%
Applied egg-rr98.3%
Final simplification98.3%
(FPCore (ux uy maxCos)
:precision binary32
(*
(sqrt
(fma
(- 2.0 (+ maxCos maxCos))
ux
(* (- 1.0 maxCos) (* ux (* ux (+ maxCos -1.0))))))
(sin (* (* uy 2.0) PI))))
float code(float ux, float uy, float maxCos) {
return sqrtf(fmaf((2.0f - (maxCos + maxCos)), ux, ((1.0f - maxCos) * (ux * (ux * (maxCos + -1.0f)))))) * sinf(((uy * 2.0f) * ((float) M_PI)));
}
function code(ux, uy, maxCos) return Float32(sqrt(fma(Float32(Float32(2.0) - Float32(maxCos + maxCos)), ux, Float32(Float32(Float32(1.0) - maxCos) * Float32(ux * Float32(ux * Float32(maxCos + Float32(-1.0))))))) * sin(Float32(Float32(uy * Float32(2.0)) * Float32(pi)))) end
\begin{array}{l}
\\
\sqrt{\mathsf{fma}\left(2 - \left(maxCos + maxCos\right), ux, \left(1 - maxCos\right) \cdot \left(ux \cdot \left(ux \cdot \left(maxCos + -1\right)\right)\right)\right)} \cdot \sin \left(\left(uy \cdot 2\right) \cdot \pi\right)
\end{array}
Initial program 55.7%
associate-*l*55.7%
sub-neg55.7%
+-commutative55.7%
distribute-rgt-neg-in55.7%
fma-def55.9%
+-commutative55.9%
associate-+r-55.8%
fma-def55.9%
neg-sub055.9%
+-commutative55.9%
associate-+r-55.7%
associate--r-55.7%
neg-sub055.7%
+-commutative55.7%
sub-neg55.7%
fma-def55.7%
Simplified55.7%
Taylor expanded in ux around 0 98.2%
+-commutative98.2%
fma-def98.2%
associate--l+98.3%
mul-1-neg98.3%
sub-neg98.3%
metadata-eval98.3%
distribute-neg-in98.3%
metadata-eval98.3%
+-commutative98.3%
sub-neg98.3%
sub-neg98.3%
metadata-eval98.3%
*-commutative98.3%
unpow298.3%
Simplified98.3%
log1p-expm1-u98.3%
Applied egg-rr98.3%
Taylor expanded in uy around inf 98.3%
fma-def98.3%
count-298.3%
sub-neg98.3%
metadata-eval98.3%
*-commutative98.3%
unpow298.3%
associate-*r*98.3%
associate-*r*98.3%
Simplified98.3%
Final simplification98.3%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= maxCos 6.199999916134402e-5)
(* (sin (* (* uy 2.0) PI)) (sqrt (fma 2.0 ux (* ux (- ux)))))
(*
2.0
(*
(sqrt
(fma
(- 2.0 (+ maxCos maxCos))
ux
(* (- 1.0 maxCos) (* ux (* ux (+ maxCos -1.0))))))
(* uy PI)))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (maxCos <= 6.199999916134402e-5f) {
tmp = sinf(((uy * 2.0f) * ((float) M_PI))) * sqrtf(fmaf(2.0f, ux, (ux * -ux)));
} else {
tmp = 2.0f * (sqrtf(fmaf((2.0f - (maxCos + maxCos)), ux, ((1.0f - maxCos) * (ux * (ux * (maxCos + -1.0f)))))) * (uy * ((float) M_PI)));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (maxCos <= Float32(6.199999916134402e-5)) tmp = Float32(sin(Float32(Float32(uy * Float32(2.0)) * Float32(pi))) * sqrt(fma(Float32(2.0), ux, Float32(ux * Float32(-ux))))); else tmp = Float32(Float32(2.0) * Float32(sqrt(fma(Float32(Float32(2.0) - Float32(maxCos + maxCos)), ux, Float32(Float32(Float32(1.0) - maxCos) * Float32(ux * Float32(ux * Float32(maxCos + Float32(-1.0))))))) * Float32(uy * Float32(pi)))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;maxCos \leq 6.199999916134402 \cdot 10^{-5}:\\
\;\;\;\;\sin \left(\left(uy \cdot 2\right) \cdot \pi\right) \cdot \sqrt{\mathsf{fma}\left(2, ux, ux \cdot \left(-ux\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\sqrt{\mathsf{fma}\left(2 - \left(maxCos + maxCos\right), ux, \left(1 - maxCos\right) \cdot \left(ux \cdot \left(ux \cdot \left(maxCos + -1\right)\right)\right)\right)} \cdot \left(uy \cdot \pi\right)\right)\\
\end{array}
\end{array}
if maxCos < 6.19999992e-5Initial program 56.0%
associate-*l*56.0%
sub-neg56.0%
+-commutative56.0%
distribute-rgt-neg-in56.0%
fma-def56.1%
+-commutative56.1%
associate-+r-56.0%
fma-def56.0%
neg-sub056.0%
+-commutative56.0%
associate-+r-56.0%
associate--r-56.0%
neg-sub056.0%
+-commutative56.0%
sub-neg56.0%
fma-def56.0%
Simplified56.0%
Taylor expanded in ux around 0 98.3%
+-commutative98.3%
fma-def98.4%
associate--l+98.4%
mul-1-neg98.4%
sub-neg98.4%
metadata-eval98.4%
distribute-neg-in98.4%
metadata-eval98.4%
+-commutative98.4%
sub-neg98.4%
sub-neg98.4%
metadata-eval98.4%
*-commutative98.4%
unpow298.4%
Simplified98.4%
Taylor expanded in maxCos around 0 97.7%
associate-*r*97.7%
+-commutative97.7%
fma-def97.7%
mul-1-neg97.7%
unpow297.7%
distribute-rgt-neg-in97.7%
Simplified97.7%
if 6.19999992e-5 < maxCos Initial program 54.1%
associate-*l*54.1%
sub-neg54.1%
+-commutative54.1%
distribute-rgt-neg-in54.1%
fma-def54.3%
+-commutative54.3%
associate-+r-53.9%
fma-def54.5%
neg-sub054.5%
+-commutative54.5%
associate-+r-53.0%
associate--r-53.0%
neg-sub053.0%
+-commutative53.0%
sub-neg53.0%
fma-def53.0%
Simplified53.0%
Taylor expanded in ux around 0 97.3%
+-commutative97.3%
fma-def97.3%
associate--l+98.1%
mul-1-neg98.1%
sub-neg98.1%
metadata-eval98.1%
distribute-neg-in98.1%
metadata-eval98.1%
+-commutative98.1%
sub-neg98.1%
sub-neg98.1%
metadata-eval98.1%
*-commutative98.1%
unpow298.1%
Simplified98.1%
log1p-expm1-u98.0%
Applied egg-rr98.0%
Taylor expanded in uy around 0 79.4%
fma-def79.5%
count-279.5%
sub-neg79.5%
metadata-eval79.5%
*-commutative79.5%
unpow279.5%
associate-*r*79.5%
Simplified79.5%
Final simplification95.6%
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0 (+ (- 1.0 ux) (* ux maxCos))))
(if (<= t_0 0.999779999256134)
(*
(sin (* (* uy 2.0) PI))
(sqrt (+ 1.0 (* t_0 (- (+ ux -1.0) (* ux maxCos))))))
(* (sin (* uy (* 2.0 PI))) (sqrt (* ux (- 2.0 (* 2.0 maxCos))))))))
float code(float ux, float uy, float maxCos) {
float t_0 = (1.0f - ux) + (ux * maxCos);
float tmp;
if (t_0 <= 0.999779999256134f) {
tmp = sinf(((uy * 2.0f) * ((float) M_PI))) * sqrtf((1.0f + (t_0 * ((ux + -1.0f) - (ux * maxCos)))));
} else {
tmp = sinf((uy * (2.0f * ((float) M_PI)))) * sqrtf((ux * (2.0f - (2.0f * maxCos))));
}
return tmp;
}
function code(ux, uy, maxCos) t_0 = Float32(Float32(Float32(1.0) - ux) + Float32(ux * maxCos)) tmp = Float32(0.0) if (t_0 <= Float32(0.999779999256134)) tmp = Float32(sin(Float32(Float32(uy * Float32(2.0)) * Float32(pi))) * sqrt(Float32(Float32(1.0) + Float32(t_0 * Float32(Float32(ux + Float32(-1.0)) - Float32(ux * maxCos)))))); else tmp = Float32(sin(Float32(uy * Float32(Float32(2.0) * Float32(pi)))) * sqrt(Float32(ux * Float32(Float32(2.0) - Float32(Float32(2.0) * maxCos))))); end return tmp end
function tmp_2 = code(ux, uy, maxCos) t_0 = (single(1.0) - ux) + (ux * maxCos); tmp = single(0.0); if (t_0 <= single(0.999779999256134)) tmp = sin(((uy * single(2.0)) * single(pi))) * sqrt((single(1.0) + (t_0 * ((ux + single(-1.0)) - (ux * maxCos))))); else tmp = sin((uy * (single(2.0) * single(pi)))) * sqrt((ux * (single(2.0) - (single(2.0) * maxCos)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(1 - ux\right) + ux \cdot maxCos\\
\mathbf{if}\;t_0 \leq 0.999779999256134:\\
\;\;\;\;\sin \left(\left(uy \cdot 2\right) \cdot \pi\right) \cdot \sqrt{1 + t_0 \cdot \left(\left(ux + -1\right) - ux \cdot maxCos\right)}\\
\mathbf{else}:\\
\;\;\;\;\sin \left(uy \cdot \left(2 \cdot \pi\right)\right) \cdot \sqrt{ux \cdot \left(2 - 2 \cdot maxCos\right)}\\
\end{array}
\end{array}
if (+.f32 (-.f32 1 ux) (*.f32 ux maxCos)) < 0.999779999Initial program 90.6%
if 0.999779999 < (+.f32 (-.f32 1 ux) (*.f32 ux maxCos)) Initial program 34.5%
associate-*l*34.5%
+-commutative34.5%
associate-+r-34.5%
fma-def34.5%
+-commutative34.5%
associate-+r-34.4%
fma-def34.4%
Simplified34.4%
Taylor expanded in ux around 0 92.2%
Final simplification91.6%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= uy 0.0009500000160187483)
(*
2.0
(*
(sqrt
(fma
(- 2.0 (+ maxCos maxCos))
ux
(* (- 1.0 maxCos) (* ux (* ux (+ maxCos -1.0))))))
(* uy PI)))
(* (sin (* uy (* 2.0 PI))) (sqrt (* ux (- 2.0 (* 2.0 maxCos)))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (uy <= 0.0009500000160187483f) {
tmp = 2.0f * (sqrtf(fmaf((2.0f - (maxCos + maxCos)), ux, ((1.0f - maxCos) * (ux * (ux * (maxCos + -1.0f)))))) * (uy * ((float) M_PI)));
} else {
tmp = sinf((uy * (2.0f * ((float) M_PI)))) * sqrtf((ux * (2.0f - (2.0f * maxCos))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (uy <= Float32(0.0009500000160187483)) tmp = Float32(Float32(2.0) * Float32(sqrt(fma(Float32(Float32(2.0) - Float32(maxCos + maxCos)), ux, Float32(Float32(Float32(1.0) - maxCos) * Float32(ux * Float32(ux * Float32(maxCos + Float32(-1.0))))))) * Float32(uy * Float32(pi)))); else tmp = Float32(sin(Float32(uy * Float32(Float32(2.0) * Float32(pi)))) * sqrt(Float32(ux * Float32(Float32(2.0) - Float32(Float32(2.0) * maxCos))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;uy \leq 0.0009500000160187483:\\
\;\;\;\;2 \cdot \left(\sqrt{\mathsf{fma}\left(2 - \left(maxCos + maxCos\right), ux, \left(1 - maxCos\right) \cdot \left(ux \cdot \left(ux \cdot \left(maxCos + -1\right)\right)\right)\right)} \cdot \left(uy \cdot \pi\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin \left(uy \cdot \left(2 \cdot \pi\right)\right) \cdot \sqrt{ux \cdot \left(2 - 2 \cdot maxCos\right)}\\
\end{array}
\end{array}
if uy < 9.50000016e-4Initial program 59.6%
associate-*l*59.6%
sub-neg59.6%
+-commutative59.6%
distribute-rgt-neg-in59.6%
fma-def59.6%
+-commutative59.6%
associate-+r-59.4%
fma-def59.5%
neg-sub059.5%
+-commutative59.5%
associate-+r-59.4%
associate--r-59.4%
neg-sub059.4%
+-commutative59.4%
sub-neg59.4%
fma-def59.4%
Simplified59.4%
Taylor expanded in ux around 0 98.5%
+-commutative98.5%
fma-def98.6%
associate--l+98.7%
mul-1-neg98.7%
sub-neg98.7%
metadata-eval98.7%
distribute-neg-in98.7%
metadata-eval98.7%
+-commutative98.7%
sub-neg98.7%
sub-neg98.7%
metadata-eval98.7%
*-commutative98.7%
unpow298.7%
Simplified98.7%
log1p-expm1-u98.7%
Applied egg-rr98.7%
Taylor expanded in uy around 0 97.1%
fma-def97.2%
count-297.2%
sub-neg97.2%
metadata-eval97.2%
*-commutative97.2%
unpow297.2%
associate-*r*97.2%
Simplified97.2%
if 9.50000016e-4 < uy Initial program 49.2%
associate-*l*49.2%
+-commutative49.2%
associate-+r-49.3%
fma-def49.3%
+-commutative49.3%
associate-+r-49.0%
fma-def49.0%
Simplified49.0%
Taylor expanded in ux around 0 80.7%
Final simplification91.1%
(FPCore (ux uy maxCos) :precision binary32 (if (<= ux 0.00039999998989515007) (* (sin (* uy (* 2.0 PI))) (sqrt (* ux (- 2.0 (* 2.0 maxCos))))) (* (sqrt (+ 1.0 (* (- 1.0 ux) (+ ux -1.0)))) (sin (* 2.0 (* uy PI))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (ux <= 0.00039999998989515007f) {
tmp = sinf((uy * (2.0f * ((float) M_PI)))) * sqrtf((ux * (2.0f - (2.0f * maxCos))));
} else {
tmp = sqrtf((1.0f + ((1.0f - ux) * (ux + -1.0f)))) * sinf((2.0f * (uy * ((float) M_PI))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (ux <= Float32(0.00039999998989515007)) tmp = Float32(sin(Float32(uy * Float32(Float32(2.0) * Float32(pi)))) * sqrt(Float32(ux * Float32(Float32(2.0) - Float32(Float32(2.0) * maxCos))))); else tmp = Float32(sqrt(Float32(Float32(1.0) + Float32(Float32(Float32(1.0) - ux) * Float32(ux + Float32(-1.0))))) * sin(Float32(Float32(2.0) * Float32(uy * Float32(pi))))); end return tmp end
function tmp_2 = code(ux, uy, maxCos) tmp = single(0.0); if (ux <= single(0.00039999998989515007)) tmp = sin((uy * (single(2.0) * single(pi)))) * sqrt((ux * (single(2.0) - (single(2.0) * maxCos)))); else tmp = sqrt((single(1.0) + ((single(1.0) - ux) * (ux + single(-1.0))))) * sin((single(2.0) * (uy * single(pi)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;ux \leq 0.00039999998989515007:\\
\;\;\;\;\sin \left(uy \cdot \left(2 \cdot \pi\right)\right) \cdot \sqrt{ux \cdot \left(2 - 2 \cdot maxCos\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{1 + \left(1 - ux\right) \cdot \left(ux + -1\right)} \cdot \sin \left(2 \cdot \left(uy \cdot \pi\right)\right)\\
\end{array}
\end{array}
if ux < 3.9999999e-4Initial program 35.4%
associate-*l*35.4%
+-commutative35.4%
associate-+r-35.4%
fma-def35.4%
+-commutative35.4%
associate-+r-35.2%
fma-def35.2%
Simplified35.2%
Taylor expanded in ux around 0 91.7%
if 3.9999999e-4 < ux Initial program 90.9%
associate-*l*90.9%
sub-neg90.9%
+-commutative90.9%
distribute-rgt-neg-in90.9%
fma-def91.4%
+-commutative91.4%
associate-+r-91.2%
fma-def91.3%
neg-sub091.3%
+-commutative91.3%
associate-+r-91.1%
associate--r-91.1%
neg-sub091.1%
+-commutative91.1%
sub-neg91.1%
fma-def91.1%
Simplified91.1%
Taylor expanded in maxCos around 0 86.0%
Final simplification89.6%
(FPCore (ux uy maxCos) :precision binary32 (if (<= uy 0.0005099999834783375) (* 2.0 (* (sqrt (fma 2.0 ux (* ux (- ux)))) (* uy PI))) (* (sin (* uy (* 2.0 PI))) (sqrt (* ux (- 2.0 (* 2.0 maxCos)))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (uy <= 0.0005099999834783375f) {
tmp = 2.0f * (sqrtf(fmaf(2.0f, ux, (ux * -ux))) * (uy * ((float) M_PI)));
} else {
tmp = sinf((uy * (2.0f * ((float) M_PI)))) * sqrtf((ux * (2.0f - (2.0f * maxCos))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (uy <= Float32(0.0005099999834783375)) tmp = Float32(Float32(2.0) * Float32(sqrt(fma(Float32(2.0), ux, Float32(ux * Float32(-ux)))) * Float32(uy * Float32(pi)))); else tmp = Float32(sin(Float32(uy * Float32(Float32(2.0) * Float32(pi)))) * sqrt(Float32(ux * Float32(Float32(2.0) - Float32(Float32(2.0) * maxCos))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;uy \leq 0.0005099999834783375:\\
\;\;\;\;2 \cdot \left(\sqrt{\mathsf{fma}\left(2, ux, ux \cdot \left(-ux\right)\right)} \cdot \left(uy \cdot \pi\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin \left(uy \cdot \left(2 \cdot \pi\right)\right) \cdot \sqrt{ux \cdot \left(2 - 2 \cdot maxCos\right)}\\
\end{array}
\end{array}
if uy < 5.09999983e-4Initial program 60.2%
associate-*l*60.2%
sub-neg60.2%
+-commutative60.2%
distribute-rgt-neg-in60.2%
fma-def60.3%
+-commutative60.3%
associate-+r-60.1%
fma-def60.2%
neg-sub060.2%
+-commutative60.2%
associate-+r-60.1%
associate--r-60.1%
neg-sub060.1%
+-commutative60.1%
sub-neg60.1%
fma-def60.1%
Simplified60.1%
Taylor expanded in uy around 0 59.7%
Taylor expanded in ux around 0 62.9%
associate--l+62.8%
*-commutative62.8%
sub-neg62.8%
metadata-eval62.8%
unpow262.8%
fma-neg62.8%
associate--l+62.9%
mul-1-neg62.9%
sub-neg62.9%
metadata-eval62.9%
metadata-eval62.9%
Simplified62.9%
Taylor expanded in maxCos around 0 91.5%
+-commutative91.5%
fma-def91.5%
mul-1-neg91.5%
unpow291.5%
distribute-rgt-neg-in91.5%
Simplified91.5%
if 5.09999983e-4 < uy Initial program 48.8%
associate-*l*48.8%
+-commutative48.8%
associate-+r-48.9%
fma-def48.9%
+-commutative48.9%
associate-+r-48.6%
fma-def48.6%
Simplified48.6%
Taylor expanded in ux around 0 80.8%
Final simplification87.3%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= ux 0.00013000000035390258)
(* 2.0 (* PI (* uy (sqrt (* ux (fma maxCos -2.0 2.0))))))
(*
2.0
(*
(* uy PI)
(sqrt
(+
1.0
(* (+ ux (- -1.0 (* ux maxCos))) (- (+ 1.0 (* ux maxCos)) ux))))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (ux <= 0.00013000000035390258f) {
tmp = 2.0f * (((float) M_PI) * (uy * sqrtf((ux * fmaf(maxCos, -2.0f, 2.0f)))));
} else {
tmp = 2.0f * ((uy * ((float) M_PI)) * sqrtf((1.0f + ((ux + (-1.0f - (ux * maxCos))) * ((1.0f + (ux * maxCos)) - ux)))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (ux <= Float32(0.00013000000035390258)) tmp = Float32(Float32(2.0) * Float32(Float32(pi) * Float32(uy * sqrt(Float32(ux * fma(maxCos, Float32(-2.0), Float32(2.0))))))); else tmp = Float32(Float32(2.0) * Float32(Float32(uy * Float32(pi)) * sqrt(Float32(Float32(1.0) + Float32(Float32(ux + Float32(Float32(-1.0) - Float32(ux * maxCos))) * Float32(Float32(Float32(1.0) + Float32(ux * maxCos)) - ux)))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;ux \leq 0.00013000000035390258:\\
\;\;\;\;2 \cdot \left(\pi \cdot \left(uy \cdot \sqrt{ux \cdot \mathsf{fma}\left(maxCos, -2, 2\right)}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left(uy \cdot \pi\right) \cdot \sqrt{1 + \left(ux + \left(-1 - ux \cdot maxCos\right)\right) \cdot \left(\left(1 + ux \cdot maxCos\right) - ux\right)}\right)\\
\end{array}
\end{array}
if ux < 1.3e-4Initial program 33.3%
associate-*l*33.3%
sub-neg33.3%
+-commutative33.3%
distribute-rgt-neg-in33.3%
fma-def33.3%
+-commutative33.3%
associate-+r-33.4%
fma-def33.4%
neg-sub033.4%
+-commutative33.4%
associate-+r-33.2%
associate--r-33.2%
neg-sub033.2%
+-commutative33.2%
sub-neg33.2%
fma-def33.2%
Simplified33.2%
Taylor expanded in uy around 0 30.4%
Taylor expanded in ux around 0 33.9%
associate--l+33.9%
*-commutative33.9%
sub-neg33.9%
metadata-eval33.9%
unpow233.9%
fma-neg33.9%
associate--l+33.9%
mul-1-neg33.9%
sub-neg33.9%
metadata-eval33.9%
metadata-eval33.9%
Simplified33.9%
Taylor expanded in ux around 0 72.7%
Taylor expanded in uy around 0 72.7%
*-commutative72.7%
cancel-sign-sub-inv72.7%
metadata-eval72.7%
associate-*r*72.7%
*-commutative72.7%
*-commutative72.7%
metadata-eval72.7%
cancel-sign-sub-inv72.7%
*-commutative72.7%
cancel-sign-sub-inv72.7%
metadata-eval72.7%
+-commutative72.7%
*-commutative72.7%
fma-def72.7%
Simplified72.7%
if 1.3e-4 < ux Initial program 89.6%
associate-*l*89.6%
sub-neg89.6%
+-commutative89.6%
distribute-rgt-neg-in89.6%
fma-def90.0%
+-commutative90.0%
associate-+r-89.7%
fma-def89.8%
neg-sub089.8%
+-commutative89.8%
associate-+r-89.6%
associate--r-89.6%
neg-sub089.6%
+-commutative89.6%
sub-neg89.6%
fma-def89.6%
Simplified89.6%
Taylor expanded in uy around 0 75.7%
Final simplification73.9%
(FPCore (ux uy maxCos) :precision binary32 (* 2.0 (* uy (* PI (sqrt (fma 2.0 ux (* ux (- ux))))))))
float code(float ux, float uy, float maxCos) {
return 2.0f * (uy * (((float) M_PI) * sqrtf(fmaf(2.0f, ux, (ux * -ux)))));
}
function code(ux, uy, maxCos) return Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * sqrt(fma(Float32(2.0), ux, Float32(ux * Float32(-ux))))))) end
\begin{array}{l}
\\
2 \cdot \left(uy \cdot \left(\pi \cdot \sqrt{\mathsf{fma}\left(2, ux, ux \cdot \left(-ux\right)\right)}\right)\right)
\end{array}
Initial program 55.7%
associate-*l*55.7%
sub-neg55.7%
+-commutative55.7%
distribute-rgt-neg-in55.7%
fma-def55.9%
+-commutative55.9%
associate-+r-55.8%
fma-def55.9%
neg-sub055.9%
+-commutative55.9%
associate-+r-55.7%
associate--r-55.7%
neg-sub055.7%
+-commutative55.7%
sub-neg55.7%
fma-def55.7%
Simplified55.7%
Taylor expanded in uy around 0 48.4%
Taylor expanded in ux around 0 51.1%
associate--l+51.0%
*-commutative51.0%
sub-neg51.0%
metadata-eval51.0%
unpow251.0%
fma-neg51.0%
associate--l+51.0%
mul-1-neg51.0%
sub-neg51.0%
metadata-eval51.0%
metadata-eval51.0%
Simplified51.0%
Taylor expanded in maxCos around 0 73.7%
associate-*l*73.6%
+-commutative73.6%
fma-def73.6%
mul-1-neg73.6%
unpow273.6%
distribute-rgt-neg-in73.6%
Simplified73.6%
Final simplification73.6%
(FPCore (ux uy maxCos) :precision binary32 (* 2.0 (* (sqrt (fma 2.0 ux (* ux (- ux)))) (* uy PI))))
float code(float ux, float uy, float maxCos) {
return 2.0f * (sqrtf(fmaf(2.0f, ux, (ux * -ux))) * (uy * ((float) M_PI)));
}
function code(ux, uy, maxCos) return Float32(Float32(2.0) * Float32(sqrt(fma(Float32(2.0), ux, Float32(ux * Float32(-ux)))) * Float32(uy * Float32(pi)))) end
\begin{array}{l}
\\
2 \cdot \left(\sqrt{\mathsf{fma}\left(2, ux, ux \cdot \left(-ux\right)\right)} \cdot \left(uy \cdot \pi\right)\right)
\end{array}
Initial program 55.7%
associate-*l*55.7%
sub-neg55.7%
+-commutative55.7%
distribute-rgt-neg-in55.7%
fma-def55.9%
+-commutative55.9%
associate-+r-55.8%
fma-def55.9%
neg-sub055.9%
+-commutative55.9%
associate-+r-55.7%
associate--r-55.7%
neg-sub055.7%
+-commutative55.7%
sub-neg55.7%
fma-def55.7%
Simplified55.7%
Taylor expanded in uy around 0 48.4%
Taylor expanded in ux around 0 51.1%
associate--l+51.0%
*-commutative51.0%
sub-neg51.0%
metadata-eval51.0%
unpow251.0%
fma-neg51.0%
associate--l+51.0%
mul-1-neg51.0%
sub-neg51.0%
metadata-eval51.0%
metadata-eval51.0%
Simplified51.0%
Taylor expanded in maxCos around 0 73.7%
+-commutative73.7%
fma-def73.7%
mul-1-neg73.7%
unpow273.7%
distribute-rgt-neg-in73.7%
Simplified73.7%
Final simplification73.7%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= ux 0.00013000000035390258)
(* 2.0 (* (* uy PI) (sqrt (+ (* -2.0 (* ux maxCos)) (* 2.0 ux)))))
(*
2.0
(*
(* uy PI)
(sqrt
(+
1.0
(* (+ ux (- -1.0 (* ux maxCos))) (- (+ 1.0 (* ux maxCos)) ux))))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (ux <= 0.00013000000035390258f) {
tmp = 2.0f * ((uy * ((float) M_PI)) * sqrtf(((-2.0f * (ux * maxCos)) + (2.0f * ux))));
} else {
tmp = 2.0f * ((uy * ((float) M_PI)) * sqrtf((1.0f + ((ux + (-1.0f - (ux * maxCos))) * ((1.0f + (ux * maxCos)) - ux)))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (ux <= Float32(0.00013000000035390258)) tmp = Float32(Float32(2.0) * Float32(Float32(uy * Float32(pi)) * sqrt(Float32(Float32(Float32(-2.0) * Float32(ux * maxCos)) + Float32(Float32(2.0) * ux))))); else tmp = Float32(Float32(2.0) * Float32(Float32(uy * Float32(pi)) * sqrt(Float32(Float32(1.0) + Float32(Float32(ux + Float32(Float32(-1.0) - Float32(ux * 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.00013000000035390258)) tmp = single(2.0) * ((uy * single(pi)) * sqrt(((single(-2.0) * (ux * maxCos)) + (single(2.0) * ux)))); else tmp = single(2.0) * ((uy * single(pi)) * sqrt((single(1.0) + ((ux + (single(-1.0) - (ux * maxCos))) * ((single(1.0) + (ux * maxCos)) - ux))))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;ux \leq 0.00013000000035390258:\\
\;\;\;\;2 \cdot \left(\left(uy \cdot \pi\right) \cdot \sqrt{-2 \cdot \left(ux \cdot maxCos\right) + 2 \cdot ux}\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left(uy \cdot \pi\right) \cdot \sqrt{1 + \left(ux + \left(-1 - ux \cdot maxCos\right)\right) \cdot \left(\left(1 + ux \cdot maxCos\right) - ux\right)}\right)\\
\end{array}
\end{array}
if ux < 1.3e-4Initial program 33.3%
associate-*l*33.3%
sub-neg33.3%
+-commutative33.3%
distribute-rgt-neg-in33.3%
fma-def33.3%
+-commutative33.3%
associate-+r-33.4%
fma-def33.4%
neg-sub033.4%
+-commutative33.4%
associate-+r-33.2%
associate--r-33.2%
neg-sub033.2%
+-commutative33.2%
sub-neg33.2%
fma-def33.2%
Simplified33.2%
Taylor expanded in uy around 0 30.4%
Taylor expanded in ux around 0 33.9%
associate--l+33.9%
*-commutative33.9%
sub-neg33.9%
metadata-eval33.9%
unpow233.9%
fma-neg33.9%
associate--l+33.9%
mul-1-neg33.9%
sub-neg33.9%
metadata-eval33.9%
metadata-eval33.9%
Simplified33.9%
Taylor expanded in ux around 0 72.7%
Taylor expanded in maxCos around 0 72.7%
if 1.3e-4 < ux Initial program 89.6%
associate-*l*89.6%
sub-neg89.6%
+-commutative89.6%
distribute-rgt-neg-in89.6%
fma-def90.0%
+-commutative90.0%
associate-+r-89.7%
fma-def89.8%
neg-sub089.8%
+-commutative89.8%
associate-+r-89.6%
associate--r-89.6%
neg-sub089.6%
+-commutative89.6%
sub-neg89.6%
fma-def89.6%
Simplified89.6%
Taylor expanded in uy around 0 75.7%
Final simplification73.9%
(FPCore (ux uy maxCos) :precision binary32 (if (<= ux 0.00014000000373926014) (* 2.0 (* (* uy PI) (sqrt (+ (* -2.0 (* ux maxCos)) (* 2.0 ux))))) (* 2.0 (* (* uy PI) (sqrt (+ 1.0 (* (- 1.0 ux) (+ ux -1.0))))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (ux <= 0.00014000000373926014f) {
tmp = 2.0f * ((uy * ((float) M_PI)) * sqrtf(((-2.0f * (ux * maxCos)) + (2.0f * ux))));
} else {
tmp = 2.0f * ((uy * ((float) M_PI)) * sqrtf((1.0f + ((1.0f - ux) * (ux + -1.0f)))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (ux <= Float32(0.00014000000373926014)) tmp = Float32(Float32(2.0) * Float32(Float32(uy * Float32(pi)) * sqrt(Float32(Float32(Float32(-2.0) * Float32(ux * maxCos)) + Float32(Float32(2.0) * ux))))); else tmp = Float32(Float32(2.0) * Float32(Float32(uy * Float32(pi)) * sqrt(Float32(Float32(1.0) + Float32(Float32(Float32(1.0) - ux) * Float32(ux + Float32(-1.0))))))); end return tmp end
function tmp_2 = code(ux, uy, maxCos) tmp = single(0.0); if (ux <= single(0.00014000000373926014)) tmp = single(2.0) * ((uy * single(pi)) * sqrt(((single(-2.0) * (ux * maxCos)) + (single(2.0) * ux)))); else tmp = single(2.0) * ((uy * single(pi)) * sqrt((single(1.0) + ((single(1.0) - ux) * (ux + single(-1.0)))))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;ux \leq 0.00014000000373926014:\\
\;\;\;\;2 \cdot \left(\left(uy \cdot \pi\right) \cdot \sqrt{-2 \cdot \left(ux \cdot maxCos\right) + 2 \cdot ux}\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left(uy \cdot \pi\right) \cdot \sqrt{1 + \left(1 - ux\right) \cdot \left(ux + -1\right)}\right)\\
\end{array}
\end{array}
if ux < 1.40000004e-4Initial program 33.5%
associate-*l*33.5%
sub-neg33.5%
+-commutative33.5%
distribute-rgt-neg-in33.5%
fma-def33.5%
+-commutative33.5%
associate-+r-33.5%
fma-def33.5%
neg-sub033.5%
+-commutative33.5%
associate-+r-33.4%
associate--r-33.4%
neg-sub033.4%
+-commutative33.4%
sub-neg33.4%
fma-def33.4%
Simplified33.4%
Taylor expanded in uy around 0 30.4%
Taylor expanded in ux around 0 33.9%
associate--l+33.9%
*-commutative33.9%
sub-neg33.9%
metadata-eval33.9%
unpow233.9%
fma-neg33.9%
associate--l+33.9%
mul-1-neg33.9%
sub-neg33.9%
metadata-eval33.9%
metadata-eval33.9%
Simplified33.9%
Taylor expanded in ux around 0 72.4%
Taylor expanded in maxCos around 0 72.4%
if 1.40000004e-4 < ux Initial program 89.8%
associate-*l*89.8%
sub-neg89.8%
+-commutative89.8%
distribute-rgt-neg-in89.8%
fma-def90.2%
+-commutative90.2%
associate-+r-90.0%
fma-def90.1%
neg-sub090.1%
+-commutative90.1%
associate-+r-89.9%
associate--r-89.9%
neg-sub089.9%
+-commutative89.9%
sub-neg89.9%
fma-def89.9%
Simplified89.9%
Taylor expanded in uy around 0 76.1%
Taylor expanded in maxCos around 0 72.8%
Final simplification72.5%
(FPCore (ux uy maxCos) :precision binary32 (* 2.0 (* (* uy PI) (sqrt (+ (* -2.0 (* ux maxCos)) (* 2.0 ux))))))
float code(float ux, float uy, float maxCos) {
return 2.0f * ((uy * ((float) M_PI)) * sqrtf(((-2.0f * (ux * maxCos)) + (2.0f * ux))));
}
function code(ux, uy, maxCos) return Float32(Float32(2.0) * Float32(Float32(uy * Float32(pi)) * sqrt(Float32(Float32(Float32(-2.0) * Float32(ux * maxCos)) + Float32(Float32(2.0) * ux))))) end
function tmp = code(ux, uy, maxCos) tmp = single(2.0) * ((uy * single(pi)) * sqrt(((single(-2.0) * (ux * maxCos)) + (single(2.0) * ux)))); end
\begin{array}{l}
\\
2 \cdot \left(\left(uy \cdot \pi\right) \cdot \sqrt{-2 \cdot \left(ux \cdot maxCos\right) + 2 \cdot ux}\right)
\end{array}
Initial program 55.7%
associate-*l*55.7%
sub-neg55.7%
+-commutative55.7%
distribute-rgt-neg-in55.7%
fma-def55.9%
+-commutative55.9%
associate-+r-55.8%
fma-def55.9%
neg-sub055.9%
+-commutative55.9%
associate-+r-55.7%
associate--r-55.7%
neg-sub055.7%
+-commutative55.7%
sub-neg55.7%
fma-def55.7%
Simplified55.7%
Taylor expanded in uy around 0 48.4%
Taylor expanded in ux around 0 51.1%
associate--l+51.0%
*-commutative51.0%
sub-neg51.0%
metadata-eval51.0%
unpow251.0%
fma-neg51.0%
associate--l+51.0%
mul-1-neg51.0%
sub-neg51.0%
metadata-eval51.0%
metadata-eval51.0%
Simplified51.0%
Taylor expanded in ux around 0 61.5%
Taylor expanded in maxCos around 0 61.5%
Final simplification61.5%
(FPCore (ux uy maxCos) :precision binary32 (* 2.0 (* (* uy PI) (sqrt (* ux (- 2.0 (* 2.0 maxCos)))))))
float code(float ux, float uy, float maxCos) {
return 2.0f * ((uy * ((float) M_PI)) * sqrtf((ux * (2.0f - (2.0f * maxCos)))));
}
function code(ux, uy, maxCos) return Float32(Float32(2.0) * Float32(Float32(uy * Float32(pi)) * sqrt(Float32(ux * Float32(Float32(2.0) - Float32(Float32(2.0) * maxCos)))))) end
function tmp = code(ux, uy, maxCos) tmp = single(2.0) * ((uy * single(pi)) * sqrt((ux * (single(2.0) - (single(2.0) * maxCos))))); end
\begin{array}{l}
\\
2 \cdot \left(\left(uy \cdot \pi\right) \cdot \sqrt{ux \cdot \left(2 - 2 \cdot maxCos\right)}\right)
\end{array}
Initial program 55.7%
associate-*l*55.7%
sub-neg55.7%
+-commutative55.7%
distribute-rgt-neg-in55.7%
fma-def55.9%
+-commutative55.9%
associate-+r-55.8%
fma-def55.9%
neg-sub055.9%
+-commutative55.9%
associate-+r-55.7%
associate--r-55.7%
neg-sub055.7%
+-commutative55.7%
sub-neg55.7%
fma-def55.7%
Simplified55.7%
Taylor expanded in uy around 0 48.4%
Taylor expanded in ux around 0 51.1%
associate--l+51.0%
*-commutative51.0%
sub-neg51.0%
metadata-eval51.0%
unpow251.0%
fma-neg51.0%
associate--l+51.0%
mul-1-neg51.0%
sub-neg51.0%
metadata-eval51.0%
metadata-eval51.0%
Simplified51.0%
Taylor expanded in ux around 0 61.5%
Final simplification61.5%
(FPCore (ux uy maxCos) :precision binary32 (* 2.0 (* (* uy PI) (sqrt (* 2.0 ux)))))
float code(float ux, float uy, float maxCos) {
return 2.0f * ((uy * ((float) M_PI)) * sqrtf((2.0f * ux)));
}
function code(ux, uy, maxCos) return Float32(Float32(2.0) * Float32(Float32(uy * Float32(pi)) * sqrt(Float32(Float32(2.0) * ux)))) end
function tmp = code(ux, uy, maxCos) tmp = single(2.0) * ((uy * single(pi)) * sqrt((single(2.0) * ux))); end
\begin{array}{l}
\\
2 \cdot \left(\left(uy \cdot \pi\right) \cdot \sqrt{2 \cdot ux}\right)
\end{array}
Initial program 55.7%
associate-*l*55.7%
sub-neg55.7%
+-commutative55.7%
distribute-rgt-neg-in55.7%
fma-def55.9%
+-commutative55.9%
associate-+r-55.8%
fma-def55.9%
neg-sub055.9%
+-commutative55.9%
associate-+r-55.7%
associate--r-55.7%
neg-sub055.7%
+-commutative55.7%
sub-neg55.7%
fma-def55.7%
Simplified55.7%
Taylor expanded in uy around 0 48.4%
Taylor expanded in ux around 0 51.1%
associate--l+51.0%
*-commutative51.0%
sub-neg51.0%
metadata-eval51.0%
unpow251.0%
fma-neg51.0%
associate--l+51.0%
mul-1-neg51.0%
sub-neg51.0%
metadata-eval51.0%
metadata-eval51.0%
Simplified51.0%
Taylor expanded in ux around 0 61.5%
Taylor expanded in maxCos around 0 58.9%
Final simplification58.9%
herbie shell --seed 2023188
(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)))))))