
(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 23 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 (* (* uy 2.0) PI))
(sqrt
(fma
(sqrt (* ux (* ux (fma maxCos -2.0 2.0))))
(sqrt (fma maxCos -2.0 2.0))
(* (fma ux maxCos (- ux)) (fma maxCos (- ux) ux))))))
float code(float ux, float uy, float maxCos) {
return sinf(((uy * 2.0f) * ((float) M_PI))) * sqrtf(fmaf(sqrtf((ux * (ux * fmaf(maxCos, -2.0f, 2.0f)))), sqrtf(fmaf(maxCos, -2.0f, 2.0f)), (fmaf(ux, maxCos, -ux) * fmaf(maxCos, -ux, ux))));
}
function code(ux, uy, maxCos) return Float32(sin(Float32(Float32(uy * Float32(2.0)) * Float32(pi))) * sqrt(fma(sqrt(Float32(ux * Float32(ux * fma(maxCos, Float32(-2.0), Float32(2.0))))), sqrt(fma(maxCos, Float32(-2.0), Float32(2.0))), Float32(fma(ux, maxCos, Float32(-ux)) * fma(maxCos, Float32(-ux), ux))))) end
\begin{array}{l}
\\
\sin \left(\left(uy \cdot 2\right) \cdot \pi\right) \cdot \sqrt{\mathsf{fma}\left(\sqrt{ux \cdot \left(ux \cdot \mathsf{fma}\left(maxCos, -2, 2\right)\right)}, \sqrt{\mathsf{fma}\left(maxCos, -2, 2\right)}, \mathsf{fma}\left(ux, maxCos, -ux\right) \cdot \mathsf{fma}\left(maxCos, -ux, ux\right)\right)}
\end{array}
Initial program 56.6%
Taylor expanded in ux around inf
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
sub-negN/A
+-commutativeN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-inN/A
metadata-evalN/A
associate-+l+N/A
Simplified98.3%
Applied egg-rr98.4%
Final simplification98.4%
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0
(sqrt
(*
ux
(fma
ux
(* (+ maxCos -1.0) (- 1.0 maxCos))
(fma maxCos -2.0 2.0))))))
(if (<= (* uy 2.0) 0.02879999950528145)
(fma
(* t_0 -1.3333333333333333)
(* uy (* uy (* uy (* PI (* PI PI)))))
(* uy (* t_0 (* 2.0 PI))))
(* (sin (* (* uy 2.0) PI)) (sqrt (* (* ux ux) (+ -1.0 (/ 2.0 ux))))))))
float code(float ux, float uy, float maxCos) {
float t_0 = sqrtf((ux * fmaf(ux, ((maxCos + -1.0f) * (1.0f - maxCos)), fmaf(maxCos, -2.0f, 2.0f))));
float tmp;
if ((uy * 2.0f) <= 0.02879999950528145f) {
tmp = fmaf((t_0 * -1.3333333333333333f), (uy * (uy * (uy * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))))), (uy * (t_0 * (2.0f * ((float) M_PI)))));
} else {
tmp = sinf(((uy * 2.0f) * ((float) M_PI))) * sqrtf(((ux * ux) * (-1.0f + (2.0f / ux))));
}
return tmp;
}
function code(ux, uy, maxCos) t_0 = sqrt(Float32(ux * fma(ux, Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), fma(maxCos, Float32(-2.0), Float32(2.0))))) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(0.02879999950528145)) tmp = fma(Float32(t_0 * Float32(-1.3333333333333333)), Float32(uy * Float32(uy * Float32(uy * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))))), Float32(uy * Float32(t_0 * Float32(Float32(2.0) * Float32(pi))))); else tmp = Float32(sin(Float32(Float32(uy * Float32(2.0)) * Float32(pi))) * sqrt(Float32(Float32(ux * ux) * Float32(Float32(-1.0) + Float32(Float32(2.0) / ux))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{ux \cdot \mathsf{fma}\left(ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), \mathsf{fma}\left(maxCos, -2, 2\right)\right)}\\
\mathbf{if}\;uy \cdot 2 \leq 0.02879999950528145:\\
\;\;\;\;\mathsf{fma}\left(t\_0 \cdot -1.3333333333333333, uy \cdot \left(uy \cdot \left(uy \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right)\right), uy \cdot \left(t\_0 \cdot \left(2 \cdot \pi\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin \left(\left(uy \cdot 2\right) \cdot \pi\right) \cdot \sqrt{\left(ux \cdot ux\right) \cdot \left(-1 + \frac{2}{ux}\right)}\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.0287999995Initial program 56.3%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
cancel-sign-sub-invN/A
+-commutativeN/A
metadata-evalN/A
associate-+l+N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified98.4%
Taylor expanded in uy around 0
Simplified98.5%
+-commutativeN/A
distribute-lft-inN/A
Applied egg-rr98.6%
if 0.0287999995 < (*.f32 uy #s(literal 2 binary32)) Initial program 58.3%
Taylor expanded in maxCos around 0
sub-negN/A
+-commutativeN/A
unpow2N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
accelerator-lowering-fma.f32N/A
--lowering--.f32N/A
mul-1-negN/A
neg-lowering-neg.f32N/A
--lowering--.f3255.4
Simplified55.4%
Taylor expanded in ux around inf
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f32N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f3291.0
Simplified91.0%
Final simplification97.3%
(FPCore (ux uy maxCos)
:precision binary32
(*
(sin (* (* uy 2.0) PI))
(sqrt
(fma
(fma ux (* (+ maxCos -1.0) (- 1.0 maxCos)) (* maxCos -2.0))
ux
(* 2.0 ux)))))
float code(float ux, float uy, float maxCos) {
return sinf(((uy * 2.0f) * ((float) M_PI))) * sqrtf(fmaf(fmaf(ux, ((maxCos + -1.0f) * (1.0f - maxCos)), (maxCos * -2.0f)), ux, (2.0f * ux)));
}
function code(ux, uy, maxCos) return Float32(sin(Float32(Float32(uy * Float32(2.0)) * Float32(pi))) * sqrt(fma(fma(ux, Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), Float32(maxCos * Float32(-2.0))), ux, Float32(Float32(2.0) * ux)))) end
\begin{array}{l}
\\
\sin \left(\left(uy \cdot 2\right) \cdot \pi\right) \cdot \sqrt{\mathsf{fma}\left(\mathsf{fma}\left(ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), maxCos \cdot -2\right), ux, 2 \cdot ux\right)}
\end{array}
Initial program 56.6%
Taylor expanded in ux around inf
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
sub-negN/A
+-commutativeN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-inN/A
metadata-evalN/A
associate-+l+N/A
Simplified98.3%
distribute-lft-inN/A
associate-*r*N/A
*-commutativeN/A
div-invN/A
associate-*l*N/A
inv-powN/A
pow2N/A
pow-prod-upN/A
metadata-evalN/A
unpow1N/A
*-commutativeN/A
distribute-lft-inN/A
associate-+r+N/A
distribute-rgt-inN/A
accelerator-lowering-fma.f32N/A
Applied egg-rr98.4%
(FPCore (ux uy maxCos) :precision binary32 (* (sin (* (* uy 2.0) PI)) (sqrt (* ux (fma ux (* (+ maxCos -1.0) (- 1.0 maxCos)) (fma maxCos -2.0 2.0))))))
float code(float ux, float uy, float maxCos) {
return sinf(((uy * 2.0f) * ((float) M_PI))) * sqrtf((ux * fmaf(ux, ((maxCos + -1.0f) * (1.0f - maxCos)), fmaf(maxCos, -2.0f, 2.0f))));
}
function code(ux, uy, maxCos) return Float32(sin(Float32(Float32(uy * Float32(2.0)) * Float32(pi))) * sqrt(Float32(ux * fma(ux, Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), fma(maxCos, Float32(-2.0), Float32(2.0)))))) end
\begin{array}{l}
\\
\sin \left(\left(uy \cdot 2\right) \cdot \pi\right) \cdot \sqrt{ux \cdot \mathsf{fma}\left(ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), \mathsf{fma}\left(maxCos, -2, 2\right)\right)}
\end{array}
Initial program 56.6%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
cancel-sign-sub-invN/A
+-commutativeN/A
metadata-evalN/A
associate-+l+N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified98.3%
(FPCore (ux uy maxCos) :precision binary32 (* (sin (* (* uy 2.0) PI)) (sqrt (fma ux (- 2.0 ux) (* (fma 2.0 ux -2.0) (* ux maxCos))))))
float code(float ux, float uy, float maxCos) {
return sinf(((uy * 2.0f) * ((float) M_PI))) * sqrtf(fmaf(ux, (2.0f - ux), (fmaf(2.0f, ux, -2.0f) * (ux * maxCos))));
}
function code(ux, uy, maxCos) return Float32(sin(Float32(Float32(uy * Float32(2.0)) * Float32(pi))) * sqrt(fma(ux, Float32(Float32(2.0) - ux), Float32(fma(Float32(2.0), ux, Float32(-2.0)) * Float32(ux * maxCos))))) end
\begin{array}{l}
\\
\sin \left(\left(uy \cdot 2\right) \cdot \pi\right) \cdot \sqrt{\mathsf{fma}\left(ux, 2 - ux, \mathsf{fma}\left(2, ux, -2\right) \cdot \left(ux \cdot maxCos\right)\right)}
\end{array}
Initial program 56.6%
Taylor expanded in ux around inf
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
sub-negN/A
+-commutativeN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-inN/A
metadata-evalN/A
associate-+l+N/A
Simplified98.3%
distribute-lft-inN/A
associate-*r*N/A
*-commutativeN/A
div-invN/A
associate-*l*N/A
inv-powN/A
pow2N/A
pow-prod-upN/A
metadata-evalN/A
unpow1N/A
*-commutativeN/A
distribute-lft-inN/A
associate-+r+N/A
distribute-rgt-inN/A
accelerator-lowering-fma.f32N/A
Applied egg-rr98.4%
Taylor expanded in maxCos around 0
associate-+r+N/A
+-commutativeN/A
mul-1-negN/A
unpow2N/A
distribute-lft-neg-inN/A
mul-1-negN/A
distribute-rgt-inN/A
accelerator-lowering-fma.f32N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
sub-negN/A
metadata-evalN/A
accelerator-lowering-fma.f3297.8
Simplified97.8%
Final simplification97.8%
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0
(sqrt
(*
ux
(fma
ux
(* (+ maxCos -1.0) (- 1.0 maxCos))
(fma maxCos -2.0 2.0))))))
(if (<= (* uy 2.0) 0.02879999950528145)
(fma
(* t_0 -1.3333333333333333)
(* uy (* uy (* uy (* PI (* PI PI)))))
(* uy (* t_0 (* 2.0 PI))))
(* (sin (* (* uy 2.0) PI)) (sqrt (* ux (- 2.0 ux)))))))
float code(float ux, float uy, float maxCos) {
float t_0 = sqrtf((ux * fmaf(ux, ((maxCos + -1.0f) * (1.0f - maxCos)), fmaf(maxCos, -2.0f, 2.0f))));
float tmp;
if ((uy * 2.0f) <= 0.02879999950528145f) {
tmp = fmaf((t_0 * -1.3333333333333333f), (uy * (uy * (uy * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))))), (uy * (t_0 * (2.0f * ((float) M_PI)))));
} else {
tmp = sinf(((uy * 2.0f) * ((float) M_PI))) * sqrtf((ux * (2.0f - ux)));
}
return tmp;
}
function code(ux, uy, maxCos) t_0 = sqrt(Float32(ux * fma(ux, Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), fma(maxCos, Float32(-2.0), Float32(2.0))))) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(0.02879999950528145)) tmp = fma(Float32(t_0 * Float32(-1.3333333333333333)), Float32(uy * Float32(uy * Float32(uy * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))))), Float32(uy * Float32(t_0 * Float32(Float32(2.0) * Float32(pi))))); else tmp = Float32(sin(Float32(Float32(uy * Float32(2.0)) * Float32(pi))) * sqrt(Float32(ux * Float32(Float32(2.0) - ux)))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{ux \cdot \mathsf{fma}\left(ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), \mathsf{fma}\left(maxCos, -2, 2\right)\right)}\\
\mathbf{if}\;uy \cdot 2 \leq 0.02879999950528145:\\
\;\;\;\;\mathsf{fma}\left(t\_0 \cdot -1.3333333333333333, uy \cdot \left(uy \cdot \left(uy \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right)\right), uy \cdot \left(t\_0 \cdot \left(2 \cdot \pi\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin \left(\left(uy \cdot 2\right) \cdot \pi\right) \cdot \sqrt{ux \cdot \left(2 - ux\right)}\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.0287999995Initial program 56.3%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
cancel-sign-sub-invN/A
+-commutativeN/A
metadata-evalN/A
associate-+l+N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified98.4%
Taylor expanded in uy around 0
Simplified98.5%
+-commutativeN/A
distribute-lft-inN/A
Applied egg-rr98.6%
if 0.0287999995 < (*.f32 uy #s(literal 2 binary32)) Initial program 58.3%
Taylor expanded in maxCos around 0
sub-negN/A
+-commutativeN/A
unpow2N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
accelerator-lowering-fma.f32N/A
--lowering--.f32N/A
mul-1-negN/A
neg-lowering-neg.f32N/A
--lowering--.f3255.4
Simplified55.4%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f3290.8
Simplified90.8%
Final simplification97.3%
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0
(sqrt
(*
ux
(fma
ux
(* (+ maxCos -1.0) (- 1.0 maxCos))
(fma maxCos -2.0 2.0))))))
(if (<= (* uy 2.0) 0.05999999865889549)
(fma
(* t_0 -1.3333333333333333)
(* uy (* uy (* uy (* PI (* PI PI)))))
(* uy (* t_0 (* 2.0 PI))))
(* (sin (* (* uy 2.0) PI)) (sqrt (* 2.0 ux))))))
float code(float ux, float uy, float maxCos) {
float t_0 = sqrtf((ux * fmaf(ux, ((maxCos + -1.0f) * (1.0f - maxCos)), fmaf(maxCos, -2.0f, 2.0f))));
float tmp;
if ((uy * 2.0f) <= 0.05999999865889549f) {
tmp = fmaf((t_0 * -1.3333333333333333f), (uy * (uy * (uy * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))))), (uy * (t_0 * (2.0f * ((float) M_PI)))));
} else {
tmp = sinf(((uy * 2.0f) * ((float) M_PI))) * sqrtf((2.0f * ux));
}
return tmp;
}
function code(ux, uy, maxCos) t_0 = sqrt(Float32(ux * fma(ux, Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), fma(maxCos, Float32(-2.0), Float32(2.0))))) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(0.05999999865889549)) tmp = fma(Float32(t_0 * Float32(-1.3333333333333333)), Float32(uy * Float32(uy * Float32(uy * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))))), Float32(uy * Float32(t_0 * Float32(Float32(2.0) * Float32(pi))))); else tmp = Float32(sin(Float32(Float32(uy * Float32(2.0)) * Float32(pi))) * sqrt(Float32(Float32(2.0) * ux))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{ux \cdot \mathsf{fma}\left(ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), \mathsf{fma}\left(maxCos, -2, 2\right)\right)}\\
\mathbf{if}\;uy \cdot 2 \leq 0.05999999865889549:\\
\;\;\;\;\mathsf{fma}\left(t\_0 \cdot -1.3333333333333333, uy \cdot \left(uy \cdot \left(uy \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right)\right), uy \cdot \left(t\_0 \cdot \left(2 \cdot \pi\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin \left(\left(uy \cdot 2\right) \cdot \pi\right) \cdot \sqrt{2 \cdot ux}\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.0599999987Initial program 56.8%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
cancel-sign-sub-invN/A
+-commutativeN/A
metadata-evalN/A
associate-+l+N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified98.5%
Taylor expanded in uy around 0
Simplified98.3%
+-commutativeN/A
distribute-lft-inN/A
Applied egg-rr98.4%
if 0.0599999987 < (*.f32 uy #s(literal 2 binary32)) Initial program 55.5%
Taylor expanded in maxCos around 0
sub-negN/A
+-commutativeN/A
unpow2N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
accelerator-lowering-fma.f32N/A
--lowering--.f32N/A
mul-1-negN/A
neg-lowering-neg.f32N/A
--lowering--.f3252.2
Simplified52.2%
Taylor expanded in ux around 0
*-lowering-*.f3272.4
Simplified72.4%
Final simplification94.3%
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0
(sqrt
(*
ux
(fma
ux
(* (+ maxCos -1.0) (- 1.0 maxCos))
(fma maxCos -2.0 2.0))))))
(fma
(* t_0 (* -1.3333333333333333 (* uy (* uy (* PI (* PI PI))))))
uy
(* uy (* t_0 (* 2.0 PI))))))
float code(float ux, float uy, float maxCos) {
float t_0 = sqrtf((ux * fmaf(ux, ((maxCos + -1.0f) * (1.0f - maxCos)), fmaf(maxCos, -2.0f, 2.0f))));
return fmaf((t_0 * (-1.3333333333333333f * (uy * (uy * (((float) M_PI) * (((float) M_PI) * ((float) M_PI))))))), uy, (uy * (t_0 * (2.0f * ((float) M_PI)))));
}
function code(ux, uy, maxCos) t_0 = sqrt(Float32(ux * fma(ux, Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), fma(maxCos, Float32(-2.0), Float32(2.0))))) return fma(Float32(t_0 * Float32(Float32(-1.3333333333333333) * Float32(uy * Float32(uy * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))))))), uy, Float32(uy * Float32(t_0 * Float32(Float32(2.0) * Float32(pi))))) end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{ux \cdot \mathsf{fma}\left(ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), \mathsf{fma}\left(maxCos, -2, 2\right)\right)}\\
\mathsf{fma}\left(t\_0 \cdot \left(-1.3333333333333333 \cdot \left(uy \cdot \left(uy \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right)\right)\right), uy, uy \cdot \left(t\_0 \cdot \left(2 \cdot \pi\right)\right)\right)
\end{array}
\end{array}
Initial program 56.6%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
cancel-sign-sub-invN/A
+-commutativeN/A
metadata-evalN/A
associate-+l+N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified98.3%
Taylor expanded in uy around 0
Simplified88.8%
+-commutativeN/A
distribute-lft-inN/A
*-commutativeN/A
Applied egg-rr88.9%
Final simplification88.9%
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0
(sqrt
(*
ux
(fma
ux
(* (+ maxCos -1.0) (- 1.0 maxCos))
(fma maxCos -2.0 2.0))))))
(fma
(* t_0 -1.3333333333333333)
(* uy (* uy (* uy (* PI (* PI PI)))))
(* uy (* t_0 (* 2.0 PI))))))
float code(float ux, float uy, float maxCos) {
float t_0 = sqrtf((ux * fmaf(ux, ((maxCos + -1.0f) * (1.0f - maxCos)), fmaf(maxCos, -2.0f, 2.0f))));
return fmaf((t_0 * -1.3333333333333333f), (uy * (uy * (uy * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))))), (uy * (t_0 * (2.0f * ((float) M_PI)))));
}
function code(ux, uy, maxCos) t_0 = sqrt(Float32(ux * fma(ux, Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), fma(maxCos, Float32(-2.0), Float32(2.0))))) return fma(Float32(t_0 * Float32(-1.3333333333333333)), Float32(uy * Float32(uy * Float32(uy * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))))), Float32(uy * Float32(t_0 * Float32(Float32(2.0) * Float32(pi))))) end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{ux \cdot \mathsf{fma}\left(ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), \mathsf{fma}\left(maxCos, -2, 2\right)\right)}\\
\mathsf{fma}\left(t\_0 \cdot -1.3333333333333333, uy \cdot \left(uy \cdot \left(uy \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right)\right), uy \cdot \left(t\_0 \cdot \left(2 \cdot \pi\right)\right)\right)
\end{array}
\end{array}
Initial program 56.6%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
cancel-sign-sub-invN/A
+-commutativeN/A
metadata-evalN/A
associate-+l+N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified98.3%
Taylor expanded in uy around 0
Simplified88.8%
+-commutativeN/A
distribute-lft-inN/A
Applied egg-rr88.9%
Final simplification88.9%
(FPCore (ux uy maxCos)
:precision binary32
(*
uy
(fma
2.0
(*
PI
(sqrt
(fma
(fma ux maxCos (- ux))
(* ux (- 1.0 maxCos))
(* ux (fma maxCos -2.0 2.0)))))
(*
-1.3333333333333333
(*
(sqrt
(* ux (fma ux (* (+ maxCos -1.0) (- 1.0 maxCos)) (fma -2.0 maxCos 2.0))))
(* (* PI (* PI PI)) (* uy uy)))))))
float code(float ux, float uy, float maxCos) {
return uy * fmaf(2.0f, (((float) M_PI) * sqrtf(fmaf(fmaf(ux, maxCos, -ux), (ux * (1.0f - maxCos)), (ux * fmaf(maxCos, -2.0f, 2.0f))))), (-1.3333333333333333f * (sqrtf((ux * fmaf(ux, ((maxCos + -1.0f) * (1.0f - maxCos)), fmaf(-2.0f, maxCos, 2.0f)))) * ((((float) M_PI) * (((float) M_PI) * ((float) M_PI))) * (uy * uy)))));
}
function code(ux, uy, maxCos) return Float32(uy * fma(Float32(2.0), Float32(Float32(pi) * sqrt(fma(fma(ux, maxCos, Float32(-ux)), Float32(ux * Float32(Float32(1.0) - maxCos)), Float32(ux * fma(maxCos, Float32(-2.0), Float32(2.0)))))), Float32(Float32(-1.3333333333333333) * Float32(sqrt(Float32(ux * fma(ux, Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), fma(Float32(-2.0), maxCos, Float32(2.0))))) * Float32(Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))) * Float32(uy * uy)))))) end
\begin{array}{l}
\\
uy \cdot \mathsf{fma}\left(2, \pi \cdot \sqrt{\mathsf{fma}\left(\mathsf{fma}\left(ux, maxCos, -ux\right), ux \cdot \left(1 - maxCos\right), ux \cdot \mathsf{fma}\left(maxCos, -2, 2\right)\right)}, -1.3333333333333333 \cdot \left(\sqrt{ux \cdot \mathsf{fma}\left(ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), \mathsf{fma}\left(-2, maxCos, 2\right)\right)} \cdot \left(\left(\pi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(uy \cdot uy\right)\right)\right)\right)
\end{array}
Initial program 56.6%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
cancel-sign-sub-invN/A
+-commutativeN/A
metadata-evalN/A
associate-+l+N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified98.3%
Taylor expanded in uy around 0
Simplified88.8%
distribute-rgt-inN/A
associate-*r*N/A
distribute-rgt-inN/A
*-commutativeN/A
neg-mul-1N/A
associate-*l*N/A
*-commutativeN/A
+-commutativeN/A
*-commutativeN/A
+-commutativeN/A
rem-square-sqrtN/A
sqrt-prodN/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
sqrt-unprodN/A
Applied egg-rr88.9%
Final simplification88.9%
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0
(sqrt
(*
ux
(fma
ux
(* (+ maxCos -1.0) (- 1.0 maxCos))
(fma maxCos -2.0 2.0))))))
(*
uy
(fma
(* 2.0 t_0)
PI
(* t_0 (* -1.3333333333333333 (* uy (* uy (* PI (* PI PI))))))))))
float code(float ux, float uy, float maxCos) {
float t_0 = sqrtf((ux * fmaf(ux, ((maxCos + -1.0f) * (1.0f - maxCos)), fmaf(maxCos, -2.0f, 2.0f))));
return uy * fmaf((2.0f * t_0), ((float) M_PI), (t_0 * (-1.3333333333333333f * (uy * (uy * (((float) M_PI) * (((float) M_PI) * ((float) M_PI))))))));
}
function code(ux, uy, maxCos) t_0 = sqrt(Float32(ux * fma(ux, Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), fma(maxCos, Float32(-2.0), Float32(2.0))))) return Float32(uy * fma(Float32(Float32(2.0) * t_0), Float32(pi), Float32(t_0 * Float32(Float32(-1.3333333333333333) * Float32(uy * Float32(uy * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))))))))) end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{ux \cdot \mathsf{fma}\left(ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), \mathsf{fma}\left(maxCos, -2, 2\right)\right)}\\
uy \cdot \mathsf{fma}\left(2 \cdot t\_0, \pi, t\_0 \cdot \left(-1.3333333333333333 \cdot \left(uy \cdot \left(uy \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right)\right)\right)\right)
\end{array}
\end{array}
Initial program 56.6%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
cancel-sign-sub-invN/A
+-commutativeN/A
metadata-evalN/A
associate-+l+N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified98.3%
Taylor expanded in uy around 0
Simplified88.8%
associate-*r*N/A
accelerator-lowering-fma.f32N/A
Applied egg-rr88.8%
Final simplification88.8%
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0
(sqrt
(*
ux
(fma
ux
(* (+ maxCos -1.0) (- 1.0 maxCos))
(fma -2.0 maxCos 2.0))))))
(*
uy
(fma
2.0
(* PI t_0)
(* -1.3333333333333333 (* t_0 (* (* PI (* PI PI)) (* uy uy))))))))
float code(float ux, float uy, float maxCos) {
float t_0 = sqrtf((ux * fmaf(ux, ((maxCos + -1.0f) * (1.0f - maxCos)), fmaf(-2.0f, maxCos, 2.0f))));
return uy * fmaf(2.0f, (((float) M_PI) * t_0), (-1.3333333333333333f * (t_0 * ((((float) M_PI) * (((float) M_PI) * ((float) M_PI))) * (uy * uy)))));
}
function code(ux, uy, maxCos) t_0 = sqrt(Float32(ux * fma(ux, Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), fma(Float32(-2.0), maxCos, Float32(2.0))))) return Float32(uy * fma(Float32(2.0), Float32(Float32(pi) * t_0), Float32(Float32(-1.3333333333333333) * Float32(t_0 * Float32(Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))) * Float32(uy * uy)))))) end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{ux \cdot \mathsf{fma}\left(ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), \mathsf{fma}\left(-2, maxCos, 2\right)\right)}\\
uy \cdot \mathsf{fma}\left(2, \pi \cdot t\_0, -1.3333333333333333 \cdot \left(t\_0 \cdot \left(\left(\pi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(uy \cdot uy\right)\right)\right)\right)
\end{array}
\end{array}
Initial program 56.6%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
cancel-sign-sub-invN/A
+-commutativeN/A
metadata-evalN/A
associate-+l+N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified98.3%
Taylor expanded in uy around 0
Simplified88.8%
Final simplification88.8%
(FPCore (ux uy maxCos)
:precision binary32
(*
(sqrt
(fma
(sqrt (* ux (* ux (fma maxCos -2.0 2.0))))
(sqrt (fma maxCos -2.0 2.0))
(* (fma ux maxCos (- ux)) (fma maxCos (- ux) ux))))
(* uy (fma (* -1.3333333333333333 (* uy uy)) (* PI (* PI PI)) (* 2.0 PI)))))
float code(float ux, float uy, float maxCos) {
return sqrtf(fmaf(sqrtf((ux * (ux * fmaf(maxCos, -2.0f, 2.0f)))), sqrtf(fmaf(maxCos, -2.0f, 2.0f)), (fmaf(ux, maxCos, -ux) * fmaf(maxCos, -ux, ux)))) * (uy * fmaf((-1.3333333333333333f * (uy * uy)), (((float) M_PI) * (((float) M_PI) * ((float) M_PI))), (2.0f * ((float) M_PI))));
}
function code(ux, uy, maxCos) return Float32(sqrt(fma(sqrt(Float32(ux * Float32(ux * fma(maxCos, Float32(-2.0), Float32(2.0))))), sqrt(fma(maxCos, Float32(-2.0), Float32(2.0))), Float32(fma(ux, maxCos, Float32(-ux)) * fma(maxCos, Float32(-ux), ux)))) * Float32(uy * fma(Float32(Float32(-1.3333333333333333) * Float32(uy * uy)), Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))), Float32(Float32(2.0) * Float32(pi))))) end
\begin{array}{l}
\\
\sqrt{\mathsf{fma}\left(\sqrt{ux \cdot \left(ux \cdot \mathsf{fma}\left(maxCos, -2, 2\right)\right)}, \sqrt{\mathsf{fma}\left(maxCos, -2, 2\right)}, \mathsf{fma}\left(ux, maxCos, -ux\right) \cdot \mathsf{fma}\left(maxCos, -ux, ux\right)\right)} \cdot \left(uy \cdot \mathsf{fma}\left(-1.3333333333333333 \cdot \left(uy \cdot uy\right), \pi \cdot \left(\pi \cdot \pi\right), 2 \cdot \pi\right)\right)
\end{array}
Initial program 56.6%
Taylor expanded in ux around inf
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
sub-negN/A
+-commutativeN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-inN/A
metadata-evalN/A
associate-+l+N/A
Simplified98.3%
Applied egg-rr98.4%
Taylor expanded in uy around 0
*-lowering-*.f32N/A
associate-*r*N/A
accelerator-lowering-fma.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
cube-multN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3288.8
Simplified88.8%
Final simplification88.8%
(FPCore (ux uy maxCos)
:precision binary32
(*
(sqrt
(fma
(fma ux (* (+ maxCos -1.0) (- 1.0 maxCos)) (* maxCos -2.0))
ux
(* 2.0 ux)))
(* uy (fma (* -1.3333333333333333 (* uy uy)) (* PI (* PI PI)) (* 2.0 PI)))))
float code(float ux, float uy, float maxCos) {
return sqrtf(fmaf(fmaf(ux, ((maxCos + -1.0f) * (1.0f - maxCos)), (maxCos * -2.0f)), ux, (2.0f * ux))) * (uy * fmaf((-1.3333333333333333f * (uy * uy)), (((float) M_PI) * (((float) M_PI) * ((float) M_PI))), (2.0f * ((float) M_PI))));
}
function code(ux, uy, maxCos) return Float32(sqrt(fma(fma(ux, Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), Float32(maxCos * Float32(-2.0))), ux, Float32(Float32(2.0) * ux))) * Float32(uy * fma(Float32(Float32(-1.3333333333333333) * Float32(uy * uy)), Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))), Float32(Float32(2.0) * Float32(pi))))) end
\begin{array}{l}
\\
\sqrt{\mathsf{fma}\left(\mathsf{fma}\left(ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), maxCos \cdot -2\right), ux, 2 \cdot ux\right)} \cdot \left(uy \cdot \mathsf{fma}\left(-1.3333333333333333 \cdot \left(uy \cdot uy\right), \pi \cdot \left(\pi \cdot \pi\right), 2 \cdot \pi\right)\right)
\end{array}
Initial program 56.6%
Taylor expanded in ux around inf
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
sub-negN/A
+-commutativeN/A
+-commutativeN/A
distribute-neg-inN/A
distribute-lft-neg-inN/A
metadata-evalN/A
associate-+l+N/A
Simplified98.3%
distribute-lft-inN/A
associate-*r*N/A
*-commutativeN/A
div-invN/A
associate-*l*N/A
inv-powN/A
pow2N/A
pow-prod-upN/A
metadata-evalN/A
unpow1N/A
*-commutativeN/A
distribute-lft-inN/A
associate-+r+N/A
distribute-rgt-inN/A
accelerator-lowering-fma.f32N/A
Applied egg-rr98.4%
Taylor expanded in uy around 0
*-lowering-*.f32N/A
associate-*r*N/A
accelerator-lowering-fma.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
cube-multN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3288.8
Simplified88.8%
Final simplification88.8%
(FPCore (ux uy maxCos) :precision binary32 (* (sqrt (* ux (fma ux (* (+ maxCos -1.0) (- 1.0 maxCos)) (fma maxCos -2.0 2.0)))) (* uy (fma -1.3333333333333333 (* (* PI (* PI PI)) (* uy uy)) (* 2.0 PI)))))
float code(float ux, float uy, float maxCos) {
return sqrtf((ux * fmaf(ux, ((maxCos + -1.0f) * (1.0f - maxCos)), fmaf(maxCos, -2.0f, 2.0f)))) * (uy * fmaf(-1.3333333333333333f, ((((float) M_PI) * (((float) M_PI) * ((float) M_PI))) * (uy * uy)), (2.0f * ((float) M_PI))));
}
function code(ux, uy, maxCos) return Float32(sqrt(Float32(ux * fma(ux, Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), fma(maxCos, Float32(-2.0), Float32(2.0))))) * Float32(uy * fma(Float32(-1.3333333333333333), Float32(Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))) * Float32(uy * uy)), Float32(Float32(2.0) * Float32(pi))))) end
\begin{array}{l}
\\
\sqrt{ux \cdot \mathsf{fma}\left(ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), \mathsf{fma}\left(maxCos, -2, 2\right)\right)} \cdot \left(uy \cdot \mathsf{fma}\left(-1.3333333333333333, \left(\pi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(uy \cdot uy\right), 2 \cdot \pi\right)\right)
\end{array}
Initial program 56.6%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
cancel-sign-sub-invN/A
+-commutativeN/A
metadata-evalN/A
associate-+l+N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified98.3%
Taylor expanded in uy around 0
*-lowering-*.f32N/A
accelerator-lowering-fma.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
cube-multN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3288.7
Simplified88.7%
Final simplification88.7%
(FPCore (ux uy maxCos)
:precision binary32
(*
uy
(*
2.0
(*
PI
(sqrt
(*
ux
(+
(fma -2.0 maxCos 2.0)
(* (fma maxCos (- ux) ux) (+ maxCos -1.0)))))))))
float code(float ux, float uy, float maxCos) {
return uy * (2.0f * (((float) M_PI) * sqrtf((ux * (fmaf(-2.0f, maxCos, 2.0f) + (fmaf(maxCos, -ux, ux) * (maxCos + -1.0f)))))));
}
function code(ux, uy, maxCos) return Float32(uy * Float32(Float32(2.0) * Float32(Float32(pi) * sqrt(Float32(ux * Float32(fma(Float32(-2.0), maxCos, Float32(2.0)) + Float32(fma(maxCos, Float32(-ux), ux) * Float32(maxCos + Float32(-1.0))))))))) end
\begin{array}{l}
\\
uy \cdot \left(2 \cdot \left(\pi \cdot \sqrt{ux \cdot \left(\mathsf{fma}\left(-2, maxCos, 2\right) + \mathsf{fma}\left(maxCos, -ux, ux\right) \cdot \left(maxCos + -1\right)\right)}\right)\right)
\end{array}
Initial program 56.6%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
cancel-sign-sub-invN/A
+-commutativeN/A
metadata-evalN/A
associate-+l+N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified98.3%
Taylor expanded in uy around 0
Simplified88.8%
Taylor expanded in uy around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
Simplified82.8%
Final simplification82.8%
(FPCore (ux uy maxCos) :precision binary32 (* (sqrt (* ux (fma ux (* (+ maxCos -1.0) (- 1.0 maxCos)) (fma maxCos -2.0 2.0)))) (* 2.0 (* uy PI))))
float code(float ux, float uy, float maxCos) {
return sqrtf((ux * fmaf(ux, ((maxCos + -1.0f) * (1.0f - maxCos)), fmaf(maxCos, -2.0f, 2.0f)))) * (2.0f * (uy * ((float) M_PI)));
}
function code(ux, uy, maxCos) return Float32(sqrt(Float32(ux * fma(ux, Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), fma(maxCos, Float32(-2.0), Float32(2.0))))) * Float32(Float32(2.0) * Float32(uy * Float32(pi)))) end
\begin{array}{l}
\\
\sqrt{ux \cdot \mathsf{fma}\left(ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), \mathsf{fma}\left(maxCos, -2, 2\right)\right)} \cdot \left(2 \cdot \left(uy \cdot \pi\right)\right)
\end{array}
Initial program 56.6%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
cancel-sign-sub-invN/A
+-commutativeN/A
metadata-evalN/A
associate-+l+N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified98.3%
Taylor expanded in uy around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3282.7
Simplified82.7%
Final simplification82.7%
(FPCore (ux uy maxCos) :precision binary32 (* 2.0 (* (* uy PI) (sqrt (fma (fma ux maxCos (- 1.0 ux)) (fma maxCos (- ux) ux) ux)))))
float code(float ux, float uy, float maxCos) {
return 2.0f * ((uy * ((float) M_PI)) * sqrtf(fmaf(fmaf(ux, maxCos, (1.0f - ux)), fmaf(maxCos, -ux, ux), ux)));
}
function code(ux, uy, maxCos) return Float32(Float32(2.0) * Float32(Float32(uy * Float32(pi)) * sqrt(fma(fma(ux, maxCos, Float32(Float32(1.0) - ux)), fma(maxCos, Float32(-ux), ux), ux)))) end
\begin{array}{l}
\\
2 \cdot \left(\left(uy \cdot \pi\right) \cdot \sqrt{\mathsf{fma}\left(\mathsf{fma}\left(ux, maxCos, 1 - ux\right), \mathsf{fma}\left(maxCos, -ux, ux\right), ux\right)}\right)
\end{array}
Initial program 56.6%
Taylor expanded in uy around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
sqrt-lowering-sqrt.f32N/A
sub-negN/A
+-commutativeN/A
unpow2N/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified49.9%
distribute-lft-inN/A
associate-+l+N/A
accelerator-lowering-fma.f32N/A
accelerator-lowering-fma.f32N/A
--lowering--.f32N/A
distribute-rgt-inN/A
neg-mul-1N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
accelerator-lowering-fma.f32N/A
neg-lowering-neg.f32N/A
+-lowering-+.f32N/A
Applied egg-rr56.3%
Taylor expanded in maxCos around 0
Simplified78.1%
(FPCore (ux uy maxCos) :precision binary32 (* uy (* 2.0 (* PI (sqrt (* ux (- 2.0 ux)))))))
float code(float ux, float uy, float maxCos) {
return uy * (2.0f * (((float) M_PI) * sqrtf((ux * (2.0f - ux)))));
}
function code(ux, uy, maxCos) return Float32(uy * Float32(Float32(2.0) * Float32(Float32(pi) * sqrt(Float32(ux * Float32(Float32(2.0) - ux)))))) end
function tmp = code(ux, uy, maxCos) tmp = uy * (single(2.0) * (single(pi) * sqrt((ux * (single(2.0) - ux))))); end
\begin{array}{l}
\\
uy \cdot \left(2 \cdot \left(\pi \cdot \sqrt{ux \cdot \left(2 - ux\right)}\right)\right)
\end{array}
Initial program 56.6%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
cancel-sign-sub-invN/A
+-commutativeN/A
metadata-evalN/A
associate-+l+N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified98.3%
Taylor expanded in uy around 0
Simplified88.8%
Taylor expanded in maxCos around 0
*-lowering-*.f32N/A
distribute-rgt-inN/A
mul-1-negN/A
distribute-lft-neg-inN/A
unpow2N/A
mul-1-negN/A
+-commutativeN/A
sqrt-lowering-sqrt.f32N/A
+-commutativeN/A
mul-1-negN/A
unpow2N/A
distribute-lft-neg-inN/A
mul-1-negN/A
distribute-rgt-inN/A
*-lowering-*.f32N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
PI-lowering-PI.f3282.9
Simplified82.9%
Taylor expanded in uy around 0
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
sub-negN/A
mul-1-negN/A
sqrt-lowering-sqrt.f32N/A
*-lowering-*.f32N/A
mul-1-negN/A
sub-negN/A
--lowering--.f3277.8
Simplified77.8%
(FPCore (ux uy maxCos) :precision binary32 (* 2.0 (* (sqrt (* ux (- 2.0 ux))) (* uy PI))))
float code(float ux, float uy, float maxCos) {
return 2.0f * (sqrtf((ux * (2.0f - ux))) * (uy * ((float) M_PI)));
}
function code(ux, uy, maxCos) return Float32(Float32(2.0) * Float32(sqrt(Float32(ux * Float32(Float32(2.0) - ux))) * Float32(uy * Float32(pi)))) end
function tmp = code(ux, uy, maxCos) tmp = single(2.0) * (sqrt((ux * (single(2.0) - ux))) * (uy * single(pi))); end
\begin{array}{l}
\\
2 \cdot \left(\sqrt{ux \cdot \left(2 - ux\right)} \cdot \left(uy \cdot \pi\right)\right)
\end{array}
Initial program 56.6%
Taylor expanded in uy around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
sqrt-lowering-sqrt.f32N/A
sub-negN/A
+-commutativeN/A
unpow2N/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified49.9%
distribute-lft-inN/A
associate-+l+N/A
accelerator-lowering-fma.f32N/A
accelerator-lowering-fma.f32N/A
--lowering--.f32N/A
distribute-rgt-inN/A
neg-mul-1N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
accelerator-lowering-fma.f32N/A
neg-lowering-neg.f32N/A
+-lowering-+.f32N/A
Applied egg-rr56.3%
Taylor expanded in maxCos around 0
+-commutativeN/A
accelerator-lowering-fma.f32N/A
--lowering--.f3277.7
Simplified77.7%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
mul-1-negN/A
sub-negN/A
--lowering--.f3277.7
Simplified77.7%
Final simplification77.7%
(FPCore (ux uy maxCos) :precision binary32 (* 2.0 (* (sqrt (* 2.0 ux)) (* uy PI))))
float code(float ux, float uy, float maxCos) {
return 2.0f * (sqrtf((2.0f * ux)) * (uy * ((float) M_PI)));
}
function code(ux, uy, maxCos) return Float32(Float32(2.0) * Float32(sqrt(Float32(Float32(2.0) * ux)) * Float32(uy * Float32(pi)))) end
function tmp = code(ux, uy, maxCos) tmp = single(2.0) * (sqrt((single(2.0) * ux)) * (uy * single(pi))); end
\begin{array}{l}
\\
2 \cdot \left(\sqrt{2 \cdot ux} \cdot \left(uy \cdot \pi\right)\right)
\end{array}
Initial program 56.6%
Taylor expanded in uy around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
sqrt-lowering-sqrt.f32N/A
sub-negN/A
+-commutativeN/A
unpow2N/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified49.9%
distribute-lft-inN/A
associate-+l+N/A
accelerator-lowering-fma.f32N/A
accelerator-lowering-fma.f32N/A
--lowering--.f32N/A
distribute-rgt-inN/A
neg-mul-1N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
accelerator-lowering-fma.f32N/A
neg-lowering-neg.f32N/A
+-lowering-+.f32N/A
Applied egg-rr56.3%
Taylor expanded in maxCos around 0
+-commutativeN/A
accelerator-lowering-fma.f32N/A
--lowering--.f3277.7
Simplified77.7%
Taylor expanded in ux around 0
*-lowering-*.f3264.2
Simplified64.2%
Final simplification64.2%
(FPCore (ux uy maxCos) :precision binary32 (* 2.0 (* (* uy PI) (* ux maxCos))))
float code(float ux, float uy, float maxCos) {
return 2.0f * ((uy * ((float) M_PI)) * (ux * maxCos));
}
function code(ux, uy, maxCos) return Float32(Float32(2.0) * Float32(Float32(uy * Float32(pi)) * Float32(ux * maxCos))) end
function tmp = code(ux, uy, maxCos) tmp = single(2.0) * ((uy * single(pi)) * (ux * maxCos)); end
\begin{array}{l}
\\
2 \cdot \left(\left(uy \cdot \pi\right) \cdot \left(ux \cdot maxCos\right)\right)
\end{array}
Initial program 56.6%
Taylor expanded in uy around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
sqrt-lowering-sqrt.f32N/A
sub-negN/A
+-commutativeN/A
unpow2N/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified49.9%
distribute-lft-inN/A
associate-+l+N/A
accelerator-lowering-fma.f32N/A
accelerator-lowering-fma.f32N/A
--lowering--.f32N/A
distribute-rgt-inN/A
neg-mul-1N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
accelerator-lowering-fma.f32N/A
neg-lowering-neg.f32N/A
+-lowering-+.f32N/A
Applied egg-rr56.3%
Applied egg-rr16.4%
Taylor expanded in maxCos around inf
associate-*r*N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3210.0
Simplified10.0%
Final simplification10.0%
(FPCore (ux uy maxCos) :precision binary32 (* -2.0 (* (* uy PI) (* ux maxCos))))
float code(float ux, float uy, float maxCos) {
return -2.0f * ((uy * ((float) M_PI)) * (ux * maxCos));
}
function code(ux, uy, maxCos) return Float32(Float32(-2.0) * Float32(Float32(uy * Float32(pi)) * Float32(ux * maxCos))) end
function tmp = code(ux, uy, maxCos) tmp = single(-2.0) * ((uy * single(pi)) * (ux * maxCos)); end
\begin{array}{l}
\\
-2 \cdot \left(\left(uy \cdot \pi\right) \cdot \left(ux \cdot maxCos\right)\right)
\end{array}
Initial program 56.6%
Taylor expanded in uy around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
sqrt-lowering-sqrt.f32N/A
sub-negN/A
+-commutativeN/A
unpow2N/A
distribute-rgt-neg-inN/A
accelerator-lowering-fma.f32N/A
Simplified49.9%
distribute-lft-inN/A
associate-+l+N/A
accelerator-lowering-fma.f32N/A
accelerator-lowering-fma.f32N/A
--lowering--.f32N/A
distribute-rgt-inN/A
neg-mul-1N/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
accelerator-lowering-fma.f32N/A
neg-lowering-neg.f32N/A
+-lowering-+.f32N/A
Applied egg-rr56.3%
Applied egg-rr16.4%
Taylor expanded in maxCos around -inf
*-lowering-*.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f326.0
Simplified6.0%
Final simplification6.0%
herbie shell --seed 2024204
(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)))))))