
(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
(*
(sqrt
(fma
(* ux ux)
(* (+ maxCos -1.0) (- 1.0 maxCos))
(* ux (fma maxCos -2.0 2.0))))
(sin (* (* 2.0 uy) PI))))
float code(float ux, float uy, float maxCos) {
return sqrtf(fmaf((ux * ux), ((maxCos + -1.0f) * (1.0f - maxCos)), (ux * fmaf(maxCos, -2.0f, 2.0f)))) * sinf(((2.0f * uy) * ((float) M_PI)));
}
function code(ux, uy, maxCos) return Float32(sqrt(fma(Float32(ux * ux), Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), Float32(ux * fma(maxCos, Float32(-2.0), Float32(2.0))))) * sin(Float32(Float32(Float32(2.0) * uy) * Float32(pi)))) end
\begin{array}{l}
\\
\sqrt{\mathsf{fma}\left(ux \cdot ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), ux \cdot \mathsf{fma}\left(maxCos, -2, 2\right)\right)} \cdot \sin \left(\left(2 \cdot uy\right) \cdot \pi\right)
\end{array}
Initial program 57.5%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites98.2%
lift-*.f32N/A
*-commutativeN/A
lower-*.f3298.2
Applied rewrites98.3%
Final simplification98.3%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= (* 2.0 uy) 0.02199999988079071)
(*
(sqrt
(fma
(* ux ux)
(* (+ maxCos -1.0) (- 1.0 maxCos))
(* ux (fma maxCos -2.0 2.0))))
(* uy (fma (* -1.3333333333333333 (* uy uy)) (* PI (* PI PI)) (* 2.0 PI))))
(* (sin (* (* 2.0 uy) PI)) (sqrt (* (* ux ux) (+ -1.0 (/ 2.0 ux)))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if ((2.0f * uy) <= 0.02199999988079071f) {
tmp = sqrtf(fmaf((ux * ux), ((maxCos + -1.0f) * (1.0f - maxCos)), (ux * fmaf(maxCos, -2.0f, 2.0f)))) * (uy * fmaf((-1.3333333333333333f * (uy * uy)), (((float) M_PI) * (((float) M_PI) * ((float) M_PI))), (2.0f * ((float) M_PI))));
} else {
tmp = sinf(((2.0f * uy) * ((float) M_PI))) * sqrtf(((ux * ux) * (-1.0f + (2.0f / ux))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (Float32(Float32(2.0) * uy) <= Float32(0.02199999988079071)) tmp = Float32(sqrt(fma(Float32(ux * ux), Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), Float32(ux * fma(maxCos, Float32(-2.0), Float32(2.0))))) * Float32(uy * fma(Float32(Float32(-1.3333333333333333) * Float32(uy * uy)), Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))), Float32(Float32(2.0) * Float32(pi))))); else tmp = Float32(sin(Float32(Float32(Float32(2.0) * uy) * 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}
\mathbf{if}\;2 \cdot uy \leq 0.02199999988079071:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(ux \cdot ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), ux \cdot \mathsf{fma}\left(maxCos, -2, 2\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)\\
\mathbf{else}:\\
\;\;\;\;\sin \left(\left(2 \cdot uy\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.0219999999Initial program 60.5%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites98.3%
lift-*.f32N/A
*-commutativeN/A
lower-*.f3298.3
Applied rewrites98.4%
Taylor expanded in uy around 0
lower-*.f32N/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
cube-multN/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f3298.5
Applied rewrites98.5%
if 0.0219999999 < (*.f32 uy #s(literal 2 binary32)) Initial program 47.2%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites97.6%
Taylor expanded in maxCos around 0
Applied rewrites92.3%
Final simplification97.1%
(FPCore (ux uy maxCos)
:precision binary32
(*
(sin (* (* 2.0 uy) PI))
(sqrt
(*
ux
(+
(- 1.0 maxCos)
(fma ux (fma maxCos (- 1.0 maxCos) (+ maxCos -1.0)) (- 1.0 maxCos)))))))
float code(float ux, float uy, float maxCos) {
return sinf(((2.0f * uy) * ((float) M_PI))) * sqrtf((ux * ((1.0f - maxCos) + fmaf(ux, fmaf(maxCos, (1.0f - maxCos), (maxCos + -1.0f)), (1.0f - maxCos)))));
}
function code(ux, uy, maxCos) return Float32(sin(Float32(Float32(Float32(2.0) * uy) * Float32(pi))) * sqrt(Float32(ux * Float32(Float32(Float32(1.0) - maxCos) + fma(ux, fma(maxCos, Float32(Float32(1.0) - maxCos), Float32(maxCos + Float32(-1.0))), Float32(Float32(1.0) - maxCos)))))) end
\begin{array}{l}
\\
\sin \left(\left(2 \cdot uy\right) \cdot \pi\right) \cdot \sqrt{ux \cdot \left(\left(1 - maxCos\right) + \mathsf{fma}\left(ux, \mathsf{fma}\left(maxCos, 1 - maxCos, maxCos + -1\right), 1 - maxCos\right)\right)}
\end{array}
Initial program 57.5%
lift--.f32N/A
sub-negN/A
+-commutativeN/A
lift-*.f32N/A
lift-+.f32N/A
distribute-rgt-inN/A
distribute-neg-inN/A
associate-+l+N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
lower-fma.f32N/A
Applied rewrites57.5%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
associate-*r*N/A
lower-fma.f32N/A
mul-1-negN/A
lower-neg.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f32N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f32N/A
lower-/.f32N/A
unpow2N/A
lower-*.f32N/A
lower-/.f3256.9
Applied rewrites56.9%
Taylor expanded in ux around 0
lower-*.f32N/A
associate-+r+N/A
lower-+.f32N/A
mul-1-negN/A
sub-negN/A
lower--.f32N/A
+-commutativeN/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
distribute-neg-inN/A
mul-1-negN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f32N/A
Applied rewrites98.2%
Final simplification98.2%
(FPCore (ux uy maxCos) :precision binary32 (* (sin (* (* 2.0 uy) 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(((2.0f * uy) * ((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(Float32(2.0) * uy) * 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(2 \cdot uy\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 57.5%
Taylor expanded in ux around 0
lower-*.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
lower-fma.f32N/A
Applied rewrites98.2%
Final simplification98.2%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= (* 2.0 uy) 0.02199999988079071)
(*
(sqrt
(fma
(* ux ux)
(* (+ maxCos -1.0) (- 1.0 maxCos))
(* ux (fma maxCos -2.0 2.0))))
(* uy (fma (* -1.3333333333333333 (* uy uy)) (* PI (* PI PI)) (* 2.0 PI))))
(* (sin (* (* 2.0 uy) PI)) (sqrt (fma ux 2.0 (* ux (- ux)))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if ((2.0f * uy) <= 0.02199999988079071f) {
tmp = sqrtf(fmaf((ux * ux), ((maxCos + -1.0f) * (1.0f - maxCos)), (ux * fmaf(maxCos, -2.0f, 2.0f)))) * (uy * fmaf((-1.3333333333333333f * (uy * uy)), (((float) M_PI) * (((float) M_PI) * ((float) M_PI))), (2.0f * ((float) M_PI))));
} else {
tmp = sinf(((2.0f * uy) * ((float) M_PI))) * sqrtf(fmaf(ux, 2.0f, (ux * -ux)));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (Float32(Float32(2.0) * uy) <= Float32(0.02199999988079071)) tmp = Float32(sqrt(fma(Float32(ux * ux), Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), Float32(ux * fma(maxCos, Float32(-2.0), Float32(2.0))))) * Float32(uy * fma(Float32(Float32(-1.3333333333333333) * Float32(uy * uy)), Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))), Float32(Float32(2.0) * Float32(pi))))); else tmp = Float32(sin(Float32(Float32(Float32(2.0) * uy) * Float32(pi))) * sqrt(fma(ux, Float32(2.0), Float32(ux * Float32(-ux))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;2 \cdot uy \leq 0.02199999988079071:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(ux \cdot ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), ux \cdot \mathsf{fma}\left(maxCos, -2, 2\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)\\
\mathbf{else}:\\
\;\;\;\;\sin \left(\left(2 \cdot uy\right) \cdot \pi\right) \cdot \sqrt{\mathsf{fma}\left(ux, 2, ux \cdot \left(-ux\right)\right)}\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.0219999999Initial program 60.5%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites98.3%
lift-*.f32N/A
*-commutativeN/A
lower-*.f3298.3
Applied rewrites98.4%
Taylor expanded in uy around 0
lower-*.f32N/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
cube-multN/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f3298.5
Applied rewrites98.5%
if 0.0219999999 < (*.f32 uy #s(literal 2 binary32)) Initial program 47.2%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites97.6%
lift-*.f32N/A
*-commutativeN/A
lower-*.f3297.6
Applied rewrites97.8%
Taylor expanded in maxCos around 0
Applied rewrites92.3%
Final simplification97.1%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= (* 2.0 uy) 0.02199999988079071)
(*
(sqrt
(fma
(* ux ux)
(* (+ maxCos -1.0) (- 1.0 maxCos))
(* ux (fma maxCos -2.0 2.0))))
(* uy (fma (* -1.3333333333333333 (* uy uy)) (* PI (* PI PI)) (* 2.0 PI))))
(* (sin (* (* 2.0 uy) PI)) (sqrt (* ux (- 2.0 ux))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if ((2.0f * uy) <= 0.02199999988079071f) {
tmp = sqrtf(fmaf((ux * ux), ((maxCos + -1.0f) * (1.0f - maxCos)), (ux * fmaf(maxCos, -2.0f, 2.0f)))) * (uy * fmaf((-1.3333333333333333f * (uy * uy)), (((float) M_PI) * (((float) M_PI) * ((float) M_PI))), (2.0f * ((float) M_PI))));
} else {
tmp = sinf(((2.0f * uy) * ((float) M_PI))) * sqrtf((ux * (2.0f - ux)));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (Float32(Float32(2.0) * uy) <= Float32(0.02199999988079071)) tmp = Float32(sqrt(fma(Float32(ux * ux), Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), Float32(ux * fma(maxCos, Float32(-2.0), Float32(2.0))))) * Float32(uy * fma(Float32(Float32(-1.3333333333333333) * Float32(uy * uy)), Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))), Float32(Float32(2.0) * Float32(pi))))); else tmp = Float32(sin(Float32(Float32(Float32(2.0) * uy) * Float32(pi))) * sqrt(Float32(ux * Float32(Float32(2.0) - ux)))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;2 \cdot uy \leq 0.02199999988079071:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(ux \cdot ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), ux \cdot \mathsf{fma}\left(maxCos, -2, 2\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)\\
\mathbf{else}:\\
\;\;\;\;\sin \left(\left(2 \cdot uy\right) \cdot \pi\right) \cdot \sqrt{ux \cdot \left(2 - ux\right)}\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.0219999999Initial program 60.5%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites98.3%
lift-*.f32N/A
*-commutativeN/A
lower-*.f3298.3
Applied rewrites98.4%
Taylor expanded in uy around 0
lower-*.f32N/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
cube-multN/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f3298.5
Applied rewrites98.5%
if 0.0219999999 < (*.f32 uy #s(literal 2 binary32)) Initial program 47.2%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites97.6%
lift-*.f32N/A
*-commutativeN/A
lower-*.f3297.6
Applied rewrites97.8%
Taylor expanded in maxCos around 0
Applied rewrites92.3%
Applied rewrites92.2%
Final simplification97.1%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= (* 2.0 uy) 0.041999999433755875)
(*
(sqrt
(fma
(* ux ux)
(* (+ maxCos -1.0) (- 1.0 maxCos))
(* ux (fma maxCos -2.0 2.0))))
(* uy (fma (* -1.3333333333333333 (* uy uy)) (* PI (* PI PI)) (* 2.0 PI))))
(* (sin (* (* 2.0 uy) PI)) (sqrt (* ux 2.0)))))
float code(float ux, float uy, float maxCos) {
float tmp;
if ((2.0f * uy) <= 0.041999999433755875f) {
tmp = sqrtf(fmaf((ux * ux), ((maxCos + -1.0f) * (1.0f - maxCos)), (ux * fmaf(maxCos, -2.0f, 2.0f)))) * (uy * fmaf((-1.3333333333333333f * (uy * uy)), (((float) M_PI) * (((float) M_PI) * ((float) M_PI))), (2.0f * ((float) M_PI))));
} else {
tmp = sinf(((2.0f * uy) * ((float) M_PI))) * sqrtf((ux * 2.0f));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (Float32(Float32(2.0) * uy) <= Float32(0.041999999433755875)) tmp = Float32(sqrt(fma(Float32(ux * ux), Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), Float32(ux * fma(maxCos, Float32(-2.0), Float32(2.0))))) * Float32(uy * fma(Float32(Float32(-1.3333333333333333) * Float32(uy * uy)), Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))), Float32(Float32(2.0) * Float32(pi))))); else tmp = Float32(sin(Float32(Float32(Float32(2.0) * uy) * Float32(pi))) * sqrt(Float32(ux * Float32(2.0)))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;2 \cdot uy \leq 0.041999999433755875:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(ux \cdot ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), ux \cdot \mathsf{fma}\left(maxCos, -2, 2\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)\\
\mathbf{else}:\\
\;\;\;\;\sin \left(\left(2 \cdot uy\right) \cdot \pi\right) \cdot \sqrt{ux \cdot 2}\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.0419999994Initial program 60.2%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites98.3%
lift-*.f32N/A
*-commutativeN/A
lower-*.f3298.3
Applied rewrites98.4%
Taylor expanded in uy around 0
lower-*.f32N/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
cube-multN/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f3298.1
Applied rewrites98.1%
if 0.0419999994 < (*.f32 uy #s(literal 2 binary32)) Initial program 46.0%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites97.5%
lift-*.f32N/A
*-commutativeN/A
lower-*.f3297.5
Applied rewrites97.7%
Taylor expanded in maxCos around 0
Applied rewrites91.5%
Taylor expanded in ux around 0
Applied rewrites79.8%
Final simplification94.6%
(FPCore (ux uy maxCos)
:precision binary32
(*
(sqrt
(fma
(* ux ux)
(* (+ maxCos -1.0) (- 1.0 maxCos))
(* ux (fma maxCos -2.0 2.0))))
(* uy (fma (* -1.3333333333333333 (* uy uy)) (* PI (* PI PI)) (* 2.0 PI)))))
float code(float ux, float uy, float maxCos) {
return sqrtf(fmaf((ux * ux), ((maxCos + -1.0f) * (1.0f - maxCos)), (ux * fmaf(maxCos, -2.0f, 2.0f)))) * (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(Float32(ux * ux), Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), Float32(ux * fma(maxCos, Float32(-2.0), Float32(2.0))))) * 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(ux \cdot ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), ux \cdot \mathsf{fma}\left(maxCos, -2, 2\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 57.5%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites98.2%
lift-*.f32N/A
*-commutativeN/A
lower-*.f3298.2
Applied rewrites98.3%
Taylor expanded in uy around 0
lower-*.f32N/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
cube-multN/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f3289.0
Applied rewrites89.0%
Final simplification89.0%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= (* 2.0 uy) 0.0014250000240281224)
(*
(sqrt
(fma
(* ux ux)
(* (+ maxCos -1.0) (- 1.0 maxCos))
(* ux (fma maxCos -2.0 2.0))))
(* 2.0 (* uy PI)))
(*
(* uy (fma (* -1.3333333333333333 (* uy uy)) (* PI (* PI PI)) (* 2.0 PI)))
(sqrt (* ux (fma -2.0 maxCos 2.0))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if ((2.0f * uy) <= 0.0014250000240281224f) {
tmp = sqrtf(fmaf((ux * ux), ((maxCos + -1.0f) * (1.0f - maxCos)), (ux * fmaf(maxCos, -2.0f, 2.0f)))) * (2.0f * (uy * ((float) M_PI)));
} else {
tmp = (uy * fmaf((-1.3333333333333333f * (uy * uy)), (((float) M_PI) * (((float) M_PI) * ((float) M_PI))), (2.0f * ((float) M_PI)))) * sqrtf((ux * fmaf(-2.0f, maxCos, 2.0f)));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (Float32(Float32(2.0) * uy) <= Float32(0.0014250000240281224)) tmp = Float32(sqrt(fma(Float32(ux * ux), Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), Float32(ux * fma(maxCos, Float32(-2.0), Float32(2.0))))) * Float32(Float32(2.0) * Float32(uy * Float32(pi)))); else tmp = Float32(Float32(uy * fma(Float32(Float32(-1.3333333333333333) * Float32(uy * uy)), Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))), Float32(Float32(2.0) * Float32(pi)))) * sqrt(Float32(ux * fma(Float32(-2.0), maxCos, Float32(2.0))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;2 \cdot uy \leq 0.0014250000240281224:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(ux \cdot ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), ux \cdot \mathsf{fma}\left(maxCos, -2, 2\right)\right)} \cdot \left(2 \cdot \left(uy \cdot \pi\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\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) \cdot \sqrt{ux \cdot \mathsf{fma}\left(-2, maxCos, 2\right)}\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.00142500002Initial program 60.7%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites98.3%
lift-*.f32N/A
*-commutativeN/A
lower-*.f3298.3
Applied rewrites98.4%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3297.5
Applied rewrites97.5%
if 0.00142500002 < (*.f32 uy #s(literal 2 binary32)) Initial program 51.1%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites97.9%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3249.3
Applied rewrites49.3%
Taylor expanded in ux around 0
cancel-sign-sub-invN/A
metadata-evalN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3247.4
Applied rewrites47.4%
Taylor expanded in uy around 0
lower-*.f32N/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
cube-multN/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f3261.9
Applied rewrites61.9%
Final simplification85.4%
(FPCore (ux uy maxCos)
:precision binary32
(if (<=
(+ 1.0 (* (+ (- 1.0 ux) (* ux maxCos)) (- (+ ux -1.0) (* ux maxCos))))
0.00019999999494757503)
(* (* 2.0 (* uy PI)) (sqrt (* ux (fma -2.0 maxCos 2.0))))
(* 2.0 (* (* uy PI) (sqrt (fma (- 1.0 ux) (+ ux -1.0) 1.0))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if ((1.0f + (((1.0f - ux) + (ux * maxCos)) * ((ux + -1.0f) - (ux * maxCos)))) <= 0.00019999999494757503f) {
tmp = (2.0f * (uy * ((float) M_PI))) * sqrtf((ux * fmaf(-2.0f, maxCos, 2.0f)));
} else {
tmp = 2.0f * ((uy * ((float) M_PI)) * sqrtf(fmaf((1.0f - ux), (ux + -1.0f), 1.0f)));
}
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.00019999999494757503)) tmp = Float32(Float32(Float32(2.0) * Float32(uy * Float32(pi))) * sqrt(Float32(ux * fma(Float32(-2.0), maxCos, Float32(2.0))))); else tmp = Float32(Float32(2.0) * Float32(Float32(uy * Float32(pi)) * sqrt(fma(Float32(Float32(1.0) - ux), Float32(ux + Float32(-1.0)), Float32(1.0))))); end return 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.00019999999494757503:\\
\;\;\;\;\left(2 \cdot \left(uy \cdot \pi\right)\right) \cdot \sqrt{ux \cdot \mathsf{fma}\left(-2, maxCos, 2\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left(uy \cdot \pi\right) \cdot \sqrt{\mathsf{fma}\left(1 - ux, ux + -1, 1\right)}\right)\\
\end{array}
\end{array}
if (-.f32 #s(literal 1 binary32) (*.f32 (+.f32 (-.f32 #s(literal 1 binary32) ux) (*.f32 ux maxCos)) (+.f32 (-.f32 #s(literal 1 binary32) ux) (*.f32 ux maxCos)))) < 1.99999995e-4Initial program 36.3%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites98.1%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3278.0
Applied rewrites78.0%
Taylor expanded in ux around 0
cancel-sign-sub-invN/A
metadata-evalN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3274.2
Applied rewrites74.2%
if 1.99999995e-4 < (-.f32 #s(literal 1 binary32) (*.f32 (+.f32 (-.f32 #s(literal 1 binary32) ux) (*.f32 ux maxCos)) (+.f32 (-.f32 #s(literal 1 binary32) ux) (*.f32 ux maxCos)))) Initial program 88.8%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-sqrt.f32N/A
sub-negN/A
+-commutativeN/A
unpow2N/A
distribute-rgt-neg-inN/A
lower-fma.f32N/A
Applied rewrites79.2%
Taylor expanded in maxCos around 0
Applied rewrites76.1%
Final simplification75.0%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= maxCos 1.4999999621068127e-5)
(*
(* uy (fma (* -1.3333333333333333 (* uy uy)) (* PI (* PI PI)) (* 2.0 PI)))
(sqrt (fma ux 2.0 (* ux (- ux)))))
(*
(* 2.0 (* uy PI))
(sqrt
(*
(* ux ux)
(fma (+ maxCos -1.0) (- 1.0 maxCos) (/ (fma maxCos -2.0 2.0) ux)))))))
float code(float ux, float uy, float maxCos) {
float tmp;
if (maxCos <= 1.4999999621068127e-5f) {
tmp = (uy * fmaf((-1.3333333333333333f * (uy * uy)), (((float) M_PI) * (((float) M_PI) * ((float) M_PI))), (2.0f * ((float) M_PI)))) * sqrtf(fmaf(ux, 2.0f, (ux * -ux)));
} else {
tmp = (2.0f * (uy * ((float) M_PI))) * sqrtf(((ux * ux) * fmaf((maxCos + -1.0f), (1.0f - maxCos), (fmaf(maxCos, -2.0f, 2.0f) / ux))));
}
return tmp;
}
function code(ux, uy, maxCos) tmp = Float32(0.0) if (maxCos <= Float32(1.4999999621068127e-5)) tmp = Float32(Float32(uy * fma(Float32(Float32(-1.3333333333333333) * Float32(uy * uy)), Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))), Float32(Float32(2.0) * Float32(pi)))) * sqrt(fma(ux, Float32(2.0), Float32(ux * Float32(-ux))))); else tmp = Float32(Float32(Float32(2.0) * Float32(uy * Float32(pi))) * sqrt(Float32(Float32(ux * ux) * fma(Float32(maxCos + Float32(-1.0)), Float32(Float32(1.0) - maxCos), Float32(fma(maxCos, Float32(-2.0), Float32(2.0)) / ux))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;maxCos \leq 1.4999999621068127 \cdot 10^{-5}:\\
\;\;\;\;\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) \cdot \sqrt{\mathsf{fma}\left(ux, 2, ux \cdot \left(-ux\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \left(uy \cdot \pi\right)\right) \cdot \sqrt{\left(ux \cdot ux\right) \cdot \mathsf{fma}\left(maxCos + -1, 1 - maxCos, \frac{\mathsf{fma}\left(maxCos, -2, 2\right)}{ux}\right)}\\
\end{array}
\end{array}
if maxCos < 1.49999996e-5Initial program 57.9%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites98.2%
lift-*.f32N/A
*-commutativeN/A
lower-*.f3298.2
Applied rewrites98.3%
Taylor expanded in maxCos around 0
Applied rewrites97.5%
Taylor expanded in uy around 0
lower-*.f32N/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
cube-multN/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f3288.9
Applied rewrites88.9%
if 1.49999996e-5 < maxCos Initial program 54.2%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites97.9%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3280.5
Applied rewrites80.5%
Final simplification87.9%
(FPCore (ux uy maxCos)
:precision binary32
(*
(sqrt
(fma
(* ux ux)
(* (+ maxCos -1.0) (- 1.0 maxCos))
(* ux (fma maxCos -2.0 2.0))))
(* 2.0 (* uy PI))))
float code(float ux, float uy, float maxCos) {
return sqrtf(fmaf((ux * ux), ((maxCos + -1.0f) * (1.0f - maxCos)), (ux * fmaf(maxCos, -2.0f, 2.0f)))) * (2.0f * (uy * ((float) M_PI)));
}
function code(ux, uy, maxCos) return Float32(sqrt(fma(Float32(ux * ux), Float32(Float32(maxCos + Float32(-1.0)) * Float32(Float32(1.0) - maxCos)), Float32(ux * fma(maxCos, Float32(-2.0), Float32(2.0))))) * Float32(Float32(2.0) * Float32(uy * Float32(pi)))) end
\begin{array}{l}
\\
\sqrt{\mathsf{fma}\left(ux \cdot ux, \left(maxCos + -1\right) \cdot \left(1 - maxCos\right), ux \cdot \mathsf{fma}\left(maxCos, -2, 2\right)\right)} \cdot \left(2 \cdot \left(uy \cdot \pi\right)\right)
\end{array}
Initial program 57.5%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites98.2%
lift-*.f32N/A
*-commutativeN/A
lower-*.f3298.2
Applied rewrites98.3%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3281.1
Applied rewrites81.1%
Final simplification81.1%
(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 57.5%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites98.2%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3281.1
Applied rewrites81.1%
Taylor expanded in ux around 0
cancel-sign-sub-invN/A
metadata-evalN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3264.8
Applied rewrites64.8%
Taylor expanded in ux around 0
lower-*.f32N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+l+N/A
metadata-evalN/A
cancel-sign-sub-invN/A
Applied rewrites81.0%
Final simplification81.0%
(FPCore (ux uy maxCos) :precision binary32 (* (sqrt (fma ux 2.0 (* ux (- ux)))) (* 2.0 (* uy PI))))
float code(float ux, float uy, float maxCos) {
return sqrtf(fmaf(ux, 2.0f, (ux * -ux))) * (2.0f * (uy * ((float) M_PI)));
}
function code(ux, uy, maxCos) return Float32(sqrt(fma(ux, Float32(2.0), Float32(ux * Float32(-ux)))) * Float32(Float32(2.0) * Float32(uy * Float32(pi)))) end
\begin{array}{l}
\\
\sqrt{\mathsf{fma}\left(ux, 2, ux \cdot \left(-ux\right)\right)} \cdot \left(2 \cdot \left(uy \cdot \pi\right)\right)
\end{array}
Initial program 57.5%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites98.2%
lift-*.f32N/A
*-commutativeN/A
lower-*.f3298.2
Applied rewrites98.3%
Taylor expanded in maxCos around 0
Applied rewrites92.6%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3277.2
Applied rewrites77.2%
(FPCore (ux uy maxCos) :precision binary32 (* (* 2.0 (* uy PI)) (sqrt (* ux (fma -2.0 maxCos 2.0)))))
float code(float ux, float uy, float maxCos) {
return (2.0f * (uy * ((float) M_PI))) * sqrtf((ux * fmaf(-2.0f, maxCos, 2.0f)));
}
function code(ux, uy, maxCos) return Float32(Float32(Float32(2.0) * Float32(uy * Float32(pi))) * sqrt(Float32(ux * fma(Float32(-2.0), maxCos, Float32(2.0))))) end
\begin{array}{l}
\\
\left(2 \cdot \left(uy \cdot \pi\right)\right) \cdot \sqrt{ux \cdot \mathsf{fma}\left(-2, maxCos, 2\right)}
\end{array}
Initial program 57.5%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.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
Applied rewrites98.2%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3281.1
Applied rewrites81.1%
Taylor expanded in ux around 0
cancel-sign-sub-invN/A
metadata-evalN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3264.8
Applied rewrites64.8%
(FPCore (ux uy maxCos) :precision binary32 (* 2.0 (* (* ux maxCos) (* uy (* PI (sqrt -1.0))))))
float code(float ux, float uy, float maxCos) {
return 2.0f * ((ux * maxCos) * (uy * (((float) M_PI) * sqrtf(-1.0f))));
}
function code(ux, uy, maxCos) return Float32(Float32(2.0) * Float32(Float32(ux * maxCos) * Float32(uy * Float32(Float32(pi) * sqrt(Float32(-1.0)))))) end
function tmp = code(ux, uy, maxCos) tmp = single(2.0) * ((ux * maxCos) * (uy * (single(pi) * sqrt(single(-1.0))))); end
\begin{array}{l}
\\
2 \cdot \left(\left(ux \cdot maxCos\right) \cdot \left(uy \cdot \left(\pi \cdot \sqrt{-1}\right)\right)\right)
\end{array}
Initial program 57.5%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-sqrt.f32N/A
sub-negN/A
+-commutativeN/A
unpow2N/A
distribute-rgt-neg-inN/A
lower-fma.f32N/A
Applied rewrites51.7%
Taylor expanded in maxCos around inf
Applied rewrites-0.0%
Final simplification-0.0%
herbie shell --seed 2024221
(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)))))))