
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* (* (- 1.0 ux) maxCos) ux))
(t_1 (sqrt (- 1.0 (* t_0 t_0))))
(t_2 (* (* uy 2.0) PI)))
(+ (+ (* (* (cos t_2) t_1) xi) (* (* (sin t_2) t_1) yi)) (* t_0 zi))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ((1.0f - ux) * maxCos) * ux;
float t_1 = sqrtf((1.0f - (t_0 * t_0)));
float t_2 = (uy * 2.0f) * ((float) M_PI);
return (((cosf(t_2) * t_1) * xi) + ((sinf(t_2) * t_1) * yi)) + (t_0 * zi);
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(Float32(Float32(1.0) - ux) * maxCos) * ux) t_1 = sqrt(Float32(Float32(1.0) - Float32(t_0 * t_0))) t_2 = Float32(Float32(uy * Float32(2.0)) * Float32(pi)) return Float32(Float32(Float32(Float32(cos(t_2) * t_1) * xi) + Float32(Float32(sin(t_2) * t_1) * yi)) + Float32(t_0 * zi)) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = ((single(1.0) - ux) * maxCos) * ux; t_1 = sqrt((single(1.0) - (t_0 * t_0))); t_2 = (uy * single(2.0)) * single(pi); tmp = (((cos(t_2) * t_1) * xi) + ((sin(t_2) * t_1) * yi)) + (t_0 * zi); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\left(1 - ux\right) \cdot maxCos\right) \cdot ux\\
t_1 := \sqrt{1 - t\_0 \cdot t\_0}\\
t_2 := \left(uy \cdot 2\right) \cdot \pi\\
\left(\left(\cos t\_2 \cdot t\_1\right) \cdot xi + \left(\sin t\_2 \cdot t\_1\right) \cdot yi\right) + t\_0 \cdot zi
\end{array}
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 25 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* (* (- 1.0 ux) maxCos) ux))
(t_1 (sqrt (- 1.0 (* t_0 t_0))))
(t_2 (* (* uy 2.0) PI)))
(+ (+ (* (* (cos t_2) t_1) xi) (* (* (sin t_2) t_1) yi)) (* t_0 zi))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ((1.0f - ux) * maxCos) * ux;
float t_1 = sqrtf((1.0f - (t_0 * t_0)));
float t_2 = (uy * 2.0f) * ((float) M_PI);
return (((cosf(t_2) * t_1) * xi) + ((sinf(t_2) * t_1) * yi)) + (t_0 * zi);
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(Float32(Float32(1.0) - ux) * maxCos) * ux) t_1 = sqrt(Float32(Float32(1.0) - Float32(t_0 * t_0))) t_2 = Float32(Float32(uy * Float32(2.0)) * Float32(pi)) return Float32(Float32(Float32(Float32(cos(t_2) * t_1) * xi) + Float32(Float32(sin(t_2) * t_1) * yi)) + Float32(t_0 * zi)) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = ((single(1.0) - ux) * maxCos) * ux; t_1 = sqrt((single(1.0) - (t_0 * t_0))); t_2 = (uy * single(2.0)) * single(pi); tmp = (((cos(t_2) * t_1) * xi) + ((sin(t_2) * t_1) * yi)) + (t_0 * zi); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\left(1 - ux\right) \cdot maxCos\right) \cdot ux\\
t_1 := \sqrt{1 - t\_0 \cdot t\_0}\\
t_2 := \left(uy \cdot 2\right) \cdot \pi\\
\left(\left(\cos t\_2 \cdot t\_1\right) \cdot xi + \left(\sin t\_2 \cdot t\_1\right) \cdot yi\right) + t\_0 \cdot zi
\end{array}
\end{array}
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* ux (* (- 1.0 ux) maxCos)))
(t_1 (sqrt (+ (* t_0 (* ux (* maxCos (+ ux -1.0)))) 1.0)))
(t_2 (* (* uy 2.0) PI)))
(+ (+ (* (* (cos t_2) t_1) xi) (* (* t_1 (sin t_2)) yi)) (* t_0 zi))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ux * ((1.0f - ux) * maxCos);
float t_1 = sqrtf(((t_0 * (ux * (maxCos * (ux + -1.0f)))) + 1.0f));
float t_2 = (uy * 2.0f) * ((float) M_PI);
return (((cosf(t_2) * t_1) * xi) + ((t_1 * sinf(t_2)) * yi)) + (t_0 * zi);
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(ux * Float32(Float32(Float32(1.0) - ux) * maxCos)) t_1 = sqrt(Float32(Float32(t_0 * Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0))))) + Float32(1.0))) t_2 = Float32(Float32(uy * Float32(2.0)) * Float32(pi)) return Float32(Float32(Float32(Float32(cos(t_2) * t_1) * xi) + Float32(Float32(t_1 * sin(t_2)) * yi)) + Float32(t_0 * zi)) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = ux * ((single(1.0) - ux) * maxCos); t_1 = sqrt(((t_0 * (ux * (maxCos * (ux + single(-1.0))))) + single(1.0))); t_2 = (uy * single(2.0)) * single(pi); tmp = (((cos(t_2) * t_1) * xi) + ((t_1 * sin(t_2)) * yi)) + (t_0 * zi); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := ux \cdot \left(\left(1 - ux\right) \cdot maxCos\right)\\
t_1 := \sqrt{t\_0 \cdot \left(ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\right) + 1}\\
t_2 := \left(uy \cdot 2\right) \cdot \pi\\
\left(\left(\cos t\_2 \cdot t\_1\right) \cdot xi + \left(t\_1 \cdot \sin t\_2\right) \cdot yi\right) + t\_0 \cdot zi
\end{array}
\end{array}
Initial program 98.7%
Final simplification98.7%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* ux (* (- 1.0 ux) maxCos))))
(+
(* t_0 zi)
(+
(*
(*
(cos (* (* uy 2.0) PI))
(sqrt (+ (* t_0 (* ux (* maxCos (+ ux -1.0)))) 1.0)))
xi)
(* yi (sin (* 2.0 (* uy PI))))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ux * ((1.0f - ux) * maxCos);
return (t_0 * zi) + (((cosf(((uy * 2.0f) * ((float) M_PI))) * sqrtf(((t_0 * (ux * (maxCos * (ux + -1.0f)))) + 1.0f))) * xi) + (yi * sinf((2.0f * (uy * ((float) M_PI))))));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(ux * Float32(Float32(Float32(1.0) - ux) * maxCos)) return Float32(Float32(t_0 * zi) + Float32(Float32(Float32(cos(Float32(Float32(uy * Float32(2.0)) * Float32(pi))) * sqrt(Float32(Float32(t_0 * Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0))))) + Float32(1.0)))) * xi) + Float32(yi * sin(Float32(Float32(2.0) * Float32(uy * Float32(pi))))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = ux * ((single(1.0) - ux) * maxCos); tmp = (t_0 * zi) + (((cos(((uy * single(2.0)) * single(pi))) * sqrt(((t_0 * (ux * (maxCos * (ux + single(-1.0))))) + single(1.0)))) * xi) + (yi * sin((single(2.0) * (uy * single(pi)))))); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := ux \cdot \left(\left(1 - ux\right) \cdot maxCos\right)\\
t\_0 \cdot zi + \left(\left(\cos \left(\left(uy \cdot 2\right) \cdot \pi\right) \cdot \sqrt{t\_0 \cdot \left(ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\right) + 1}\right) \cdot xi + yi \cdot \sin \left(2 \cdot \left(uy \cdot \pi\right)\right)\right)
\end{array}
\end{array}
Initial program 98.7%
Taylor expanded in ux around 0
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3298.6
Applied rewrites98.6%
Final simplification98.6%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* 2.0 (* uy PI)))
(t_1
(sqrt
(fma
(* (* maxCos maxCos) (* ux ux))
(* (- 1.0 ux) (+ ux -1.0))
1.0))))
(if (<= (* uy 2.0) 0.009499999694526196)
(fma
(* zi (* (- 1.0 ux) maxCos))
ux
(fma
uy
(fma
uy
(*
t_1
(fma
-1.3333333333333333
(* (* uy yi) (* PI (* PI PI)))
(* -2.0 (* xi (* PI PI)))))
(* t_1 (* 2.0 (* PI yi))))
(* xi t_1)))
(fma maxCos (* ux zi) (fma xi (cos t_0) (* yi (sin t_0)))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = 2.0f * (uy * ((float) M_PI));
float t_1 = sqrtf(fmaf(((maxCos * maxCos) * (ux * ux)), ((1.0f - ux) * (ux + -1.0f)), 1.0f));
float tmp;
if ((uy * 2.0f) <= 0.009499999694526196f) {
tmp = fmaf((zi * ((1.0f - ux) * maxCos)), ux, fmaf(uy, fmaf(uy, (t_1 * fmaf(-1.3333333333333333f, ((uy * yi) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))), (-2.0f * (xi * (((float) M_PI) * ((float) M_PI)))))), (t_1 * (2.0f * (((float) M_PI) * yi)))), (xi * t_1)));
} else {
tmp = fmaf(maxCos, (ux * zi), fmaf(xi, cosf(t_0), (yi * sinf(t_0))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(2.0) * Float32(uy * Float32(pi))) t_1 = sqrt(fma(Float32(Float32(maxCos * maxCos) * Float32(ux * ux)), Float32(Float32(Float32(1.0) - ux) * Float32(ux + Float32(-1.0))), Float32(1.0))) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(0.009499999694526196)) tmp = fma(Float32(zi * Float32(Float32(Float32(1.0) - ux) * maxCos)), ux, fma(uy, fma(uy, Float32(t_1 * fma(Float32(-1.3333333333333333), Float32(Float32(uy * yi) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))), Float32(Float32(-2.0) * Float32(xi * Float32(Float32(pi) * Float32(pi)))))), Float32(t_1 * Float32(Float32(2.0) * Float32(Float32(pi) * yi)))), Float32(xi * t_1))); else tmp = fma(maxCos, Float32(ux * zi), fma(xi, cos(t_0), Float32(yi * sin(t_0)))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 2 \cdot \left(uy \cdot \pi\right)\\
t_1 := \sqrt{\mathsf{fma}\left(\left(maxCos \cdot maxCos\right) \cdot \left(ux \cdot ux\right), \left(1 - ux\right) \cdot \left(ux + -1\right), 1\right)}\\
\mathbf{if}\;uy \cdot 2 \leq 0.009499999694526196:\\
\;\;\;\;\mathsf{fma}\left(zi \cdot \left(\left(1 - ux\right) \cdot maxCos\right), ux, \mathsf{fma}\left(uy, \mathsf{fma}\left(uy, t\_1 \cdot \mathsf{fma}\left(-1.3333333333333333, \left(uy \cdot yi\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), -2 \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right)\right), t\_1 \cdot \left(2 \cdot \left(\pi \cdot yi\right)\right)\right), xi \cdot t\_1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(maxCos, ux \cdot zi, \mathsf{fma}\left(xi, \cos t\_0, yi \cdot \sin t\_0\right)\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.00949999969Initial program 99.2%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3299.2
Applied rewrites99.2%
Taylor expanded in uy around 0
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites94.5%
Applied rewrites94.7%
Taylor expanded in uy around 0
Applied rewrites99.4%
if 0.00949999969 < (*.f32 uy #s(literal 2 binary32)) Initial program 97.2%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3296.7
Applied rewrites96.7%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3295.3
Applied rewrites95.3%
Final simplification98.4%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* 2.0 (* uy PI)))
(t_1
(sqrt
(fma
(* (* maxCos maxCos) (* ux ux))
(* (- 1.0 ux) (+ ux -1.0))
1.0))))
(if (<= (* uy 2.0) 0.009499999694526196)
(fma
(* zi (* (- 1.0 ux) maxCos))
ux
(fma
uy
(fma
uy
(*
t_1
(fma
-1.3333333333333333
(* (* uy yi) (* PI (* PI PI)))
(* -2.0 (* xi (* PI PI)))))
(* t_1 (* 2.0 (* PI yi))))
(* xi t_1)))
(fma xi (cos t_0) (fma yi (sin t_0) (* zi (* ux maxCos)))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = 2.0f * (uy * ((float) M_PI));
float t_1 = sqrtf(fmaf(((maxCos * maxCos) * (ux * ux)), ((1.0f - ux) * (ux + -1.0f)), 1.0f));
float tmp;
if ((uy * 2.0f) <= 0.009499999694526196f) {
tmp = fmaf((zi * ((1.0f - ux) * maxCos)), ux, fmaf(uy, fmaf(uy, (t_1 * fmaf(-1.3333333333333333f, ((uy * yi) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))), (-2.0f * (xi * (((float) M_PI) * ((float) M_PI)))))), (t_1 * (2.0f * (((float) M_PI) * yi)))), (xi * t_1)));
} else {
tmp = fmaf(xi, cosf(t_0), fmaf(yi, sinf(t_0), (zi * (ux * maxCos))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(2.0) * Float32(uy * Float32(pi))) t_1 = sqrt(fma(Float32(Float32(maxCos * maxCos) * Float32(ux * ux)), Float32(Float32(Float32(1.0) - ux) * Float32(ux + Float32(-1.0))), Float32(1.0))) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(0.009499999694526196)) tmp = fma(Float32(zi * Float32(Float32(Float32(1.0) - ux) * maxCos)), ux, fma(uy, fma(uy, Float32(t_1 * fma(Float32(-1.3333333333333333), Float32(Float32(uy * yi) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))), Float32(Float32(-2.0) * Float32(xi * Float32(Float32(pi) * Float32(pi)))))), Float32(t_1 * Float32(Float32(2.0) * Float32(Float32(pi) * yi)))), Float32(xi * t_1))); else tmp = fma(xi, cos(t_0), fma(yi, sin(t_0), Float32(zi * Float32(ux * maxCos)))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 2 \cdot \left(uy \cdot \pi\right)\\
t_1 := \sqrt{\mathsf{fma}\left(\left(maxCos \cdot maxCos\right) \cdot \left(ux \cdot ux\right), \left(1 - ux\right) \cdot \left(ux + -1\right), 1\right)}\\
\mathbf{if}\;uy \cdot 2 \leq 0.009499999694526196:\\
\;\;\;\;\mathsf{fma}\left(zi \cdot \left(\left(1 - ux\right) \cdot maxCos\right), ux, \mathsf{fma}\left(uy, \mathsf{fma}\left(uy, t\_1 \cdot \mathsf{fma}\left(-1.3333333333333333, \left(uy \cdot yi\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), -2 \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right)\right), t\_1 \cdot \left(2 \cdot \left(\pi \cdot yi\right)\right)\right), xi \cdot t\_1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(xi, \cos t\_0, \mathsf{fma}\left(yi, \sin t\_0, zi \cdot \left(ux \cdot maxCos\right)\right)\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.00949999969Initial program 99.2%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3299.2
Applied rewrites99.2%
Taylor expanded in uy around 0
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites94.5%
Applied rewrites94.7%
Taylor expanded in uy around 0
Applied rewrites99.4%
if 0.00949999969 < (*.f32 uy #s(literal 2 binary32)) Initial program 97.2%
Taylor expanded in ux around 0
+-commutativeN/A
associate-+l+N/A
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-fma.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f3295.2
Applied rewrites95.2%
Final simplification98.4%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (let* ((t_0 (* 2.0 (* uy PI)))) (fma xi (cos t_0) (fma yi (sin t_0) (* maxCos (* (- 1.0 ux) (* ux zi)))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = 2.0f * (uy * ((float) M_PI));
return fmaf(xi, cosf(t_0), fmaf(yi, sinf(t_0), (maxCos * ((1.0f - ux) * (ux * zi)))));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(2.0) * Float32(uy * Float32(pi))) return fma(xi, cos(t_0), fma(yi, sin(t_0), Float32(maxCos * Float32(Float32(Float32(1.0) - ux) * Float32(ux * zi))))) end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 2 \cdot \left(uy \cdot \pi\right)\\
\mathsf{fma}\left(xi, \cos t\_0, \mathsf{fma}\left(yi, \sin t\_0, maxCos \cdot \left(\left(1 - ux\right) \cdot \left(ux \cdot zi\right)\right)\right)\right)
\end{array}
\end{array}
Initial program 98.7%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-fma.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
Applied rewrites98.5%
Final simplification98.5%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* 2.0 (* uy PI)))
(t_1
(sqrt
(fma
(* (* maxCos maxCos) (* ux ux))
(* (- 1.0 ux) (+ ux -1.0))
1.0))))
(if (<= (* uy 2.0) 0.04500000178813934)
(fma
(* zi (* (- 1.0 ux) maxCos))
ux
(fma
uy
(fma
uy
(*
t_1
(fma
-1.3333333333333333
(* (* uy yi) (* PI (* PI PI)))
(* -2.0 (* xi (* PI PI)))))
(* t_1 (* 2.0 (* PI yi))))
(* xi t_1)))
(fma xi (cos t_0) (* yi (sin t_0))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = 2.0f * (uy * ((float) M_PI));
float t_1 = sqrtf(fmaf(((maxCos * maxCos) * (ux * ux)), ((1.0f - ux) * (ux + -1.0f)), 1.0f));
float tmp;
if ((uy * 2.0f) <= 0.04500000178813934f) {
tmp = fmaf((zi * ((1.0f - ux) * maxCos)), ux, fmaf(uy, fmaf(uy, (t_1 * fmaf(-1.3333333333333333f, ((uy * yi) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))), (-2.0f * (xi * (((float) M_PI) * ((float) M_PI)))))), (t_1 * (2.0f * (((float) M_PI) * yi)))), (xi * t_1)));
} else {
tmp = fmaf(xi, cosf(t_0), (yi * sinf(t_0)));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(2.0) * Float32(uy * Float32(pi))) t_1 = sqrt(fma(Float32(Float32(maxCos * maxCos) * Float32(ux * ux)), Float32(Float32(Float32(1.0) - ux) * Float32(ux + Float32(-1.0))), Float32(1.0))) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(0.04500000178813934)) tmp = fma(Float32(zi * Float32(Float32(Float32(1.0) - ux) * maxCos)), ux, fma(uy, fma(uy, Float32(t_1 * fma(Float32(-1.3333333333333333), Float32(Float32(uy * yi) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))), Float32(Float32(-2.0) * Float32(xi * Float32(Float32(pi) * Float32(pi)))))), Float32(t_1 * Float32(Float32(2.0) * Float32(Float32(pi) * yi)))), Float32(xi * t_1))); else tmp = fma(xi, cos(t_0), Float32(yi * sin(t_0))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 2 \cdot \left(uy \cdot \pi\right)\\
t_1 := \sqrt{\mathsf{fma}\left(\left(maxCos \cdot maxCos\right) \cdot \left(ux \cdot ux\right), \left(1 - ux\right) \cdot \left(ux + -1\right), 1\right)}\\
\mathbf{if}\;uy \cdot 2 \leq 0.04500000178813934:\\
\;\;\;\;\mathsf{fma}\left(zi \cdot \left(\left(1 - ux\right) \cdot maxCos\right), ux, \mathsf{fma}\left(uy, \mathsf{fma}\left(uy, t\_1 \cdot \mathsf{fma}\left(-1.3333333333333333, \left(uy \cdot yi\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), -2 \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right)\right), t\_1 \cdot \left(2 \cdot \left(\pi \cdot yi\right)\right)\right), xi \cdot t\_1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(xi, \cos t\_0, yi \cdot \sin t\_0\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.0450000018Initial program 99.2%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3299.2
Applied rewrites99.2%
Taylor expanded in uy around 0
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites92.5%
Applied rewrites92.7%
Taylor expanded in uy around 0
Applied rewrites99.0%
if 0.0450000018 < (*.f32 uy #s(literal 2 binary32)) Initial program 96.8%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3291.9
Applied rewrites91.9%
Final simplification97.6%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(+
(* (* ux (* (- 1.0 ux) maxCos)) zi)
(*
yi
(*
(sqrt (fma (* maxCos maxCos) (* (* ux ux) (* (- 1.0 ux) (+ ux -1.0))) 1.0))
(fma
(fma (* -2.0 (* uy uy)) (* PI PI) 1.0)
(/ xi yi)
(sin (* 2.0 (* uy PI))))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return ((ux * ((1.0f - ux) * maxCos)) * zi) + (yi * (sqrtf(fmaf((maxCos * maxCos), ((ux * ux) * ((1.0f - ux) * (ux + -1.0f))), 1.0f)) * fmaf(fmaf((-2.0f * (uy * uy)), (((float) M_PI) * ((float) M_PI)), 1.0f), (xi / yi), sinf((2.0f * (uy * ((float) M_PI)))))));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(Float32(ux * Float32(Float32(Float32(1.0) - ux) * maxCos)) * zi) + Float32(yi * Float32(sqrt(fma(Float32(maxCos * maxCos), Float32(Float32(ux * ux) * Float32(Float32(Float32(1.0) - ux) * Float32(ux + Float32(-1.0)))), Float32(1.0))) * fma(fma(Float32(Float32(-2.0) * Float32(uy * uy)), Float32(Float32(pi) * Float32(pi)), Float32(1.0)), Float32(xi / yi), sin(Float32(Float32(2.0) * Float32(uy * Float32(pi)))))))) end
\begin{array}{l}
\\
\left(ux \cdot \left(\left(1 - ux\right) \cdot maxCos\right)\right) \cdot zi + yi \cdot \left(\sqrt{\mathsf{fma}\left(maxCos \cdot maxCos, \left(ux \cdot ux\right) \cdot \left(\left(1 - ux\right) \cdot \left(ux + -1\right)\right), 1\right)} \cdot \mathsf{fma}\left(\mathsf{fma}\left(-2 \cdot \left(uy \cdot uy\right), \pi \cdot \pi, 1\right), \frac{xi}{yi}, \sin \left(2 \cdot \left(uy \cdot \pi\right)\right)\right)\right)
\end{array}
Initial program 98.7%
Taylor expanded in yi around inf
lower-*.f32N/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites98.6%
Taylor expanded in uy around 0
Applied rewrites92.6%
Final simplification92.6%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(*
zi
(fma
maxCos
(* ux (- 1.0 ux))
(/
(*
(sqrt
(fma (* maxCos maxCos) (* (* ux ux) (* (- 1.0 ux) (+ ux -1.0))) 1.0))
(fma
xi
(fma (* -2.0 (* uy uy)) (* PI PI) 1.0)
(* yi (sin (* 2.0 (* uy PI))))))
zi))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return zi * fmaf(maxCos, (ux * (1.0f - ux)), ((sqrtf(fmaf((maxCos * maxCos), ((ux * ux) * ((1.0f - ux) * (ux + -1.0f))), 1.0f)) * fmaf(xi, fmaf((-2.0f * (uy * uy)), (((float) M_PI) * ((float) M_PI)), 1.0f), (yi * sinf((2.0f * (uy * ((float) M_PI))))))) / zi));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(zi * fma(maxCos, Float32(ux * Float32(Float32(1.0) - ux)), Float32(Float32(sqrt(fma(Float32(maxCos * maxCos), Float32(Float32(ux * ux) * Float32(Float32(Float32(1.0) - ux) * Float32(ux + Float32(-1.0)))), Float32(1.0))) * fma(xi, fma(Float32(Float32(-2.0) * Float32(uy * uy)), Float32(Float32(pi) * Float32(pi)), Float32(1.0)), Float32(yi * sin(Float32(Float32(2.0) * Float32(uy * Float32(pi))))))) / zi))) end
\begin{array}{l}
\\
zi \cdot \mathsf{fma}\left(maxCos, ux \cdot \left(1 - ux\right), \frac{\sqrt{\mathsf{fma}\left(maxCos \cdot maxCos, \left(ux \cdot ux\right) \cdot \left(\left(1 - ux\right) \cdot \left(ux + -1\right)\right), 1\right)} \cdot \mathsf{fma}\left(xi, \mathsf{fma}\left(-2 \cdot \left(uy \cdot uy\right), \pi \cdot \pi, 1\right), yi \cdot \sin \left(2 \cdot \left(uy \cdot \pi\right)\right)\right)}{zi}\right)
\end{array}
Initial program 98.7%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3298.6
Applied rewrites98.6%
Taylor expanded in zi around -inf
Applied rewrites98.6%
Taylor expanded in uy around 0
Applied rewrites92.7%
Final simplification92.7%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* PI (* PI PI))))
(if (<= (* uy 2.0) 0.04500000178813934)
(*
zi
(fma
maxCos
(* ux (- 1.0 ux))
(/
(*
(sqrt
(fma (* maxCos maxCos) (* (* ux ux) (* (- 1.0 ux) (+ ux -1.0))) 1.0))
(fma
uy
(fma
uy
(fma
-1.3333333333333333
(* (* uy yi) t_0)
(* (* PI PI) (* xi -2.0)))
(* 2.0 (* PI yi)))
xi))
zi)))
(fma
xi
(cos (* 2.0 (* uy PI)))
(* yi (* uy (fma (* -1.3333333333333333 (* uy uy)) t_0 (* 2.0 PI))))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ((float) M_PI) * (((float) M_PI) * ((float) M_PI));
float tmp;
if ((uy * 2.0f) <= 0.04500000178813934f) {
tmp = zi * fmaf(maxCos, (ux * (1.0f - ux)), ((sqrtf(fmaf((maxCos * maxCos), ((ux * ux) * ((1.0f - ux) * (ux + -1.0f))), 1.0f)) * fmaf(uy, fmaf(uy, fmaf(-1.3333333333333333f, ((uy * yi) * t_0), ((((float) M_PI) * ((float) M_PI)) * (xi * -2.0f))), (2.0f * (((float) M_PI) * yi))), xi)) / zi));
} else {
tmp = fmaf(xi, cosf((2.0f * (uy * ((float) M_PI)))), (yi * (uy * fmaf((-1.3333333333333333f * (uy * uy)), t_0, (2.0f * ((float) M_PI))))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(0.04500000178813934)) tmp = Float32(zi * fma(maxCos, Float32(ux * Float32(Float32(1.0) - ux)), Float32(Float32(sqrt(fma(Float32(maxCos * maxCos), Float32(Float32(ux * ux) * Float32(Float32(Float32(1.0) - ux) * Float32(ux + Float32(-1.0)))), Float32(1.0))) * fma(uy, fma(uy, fma(Float32(-1.3333333333333333), Float32(Float32(uy * yi) * t_0), Float32(Float32(Float32(pi) * Float32(pi)) * Float32(xi * Float32(-2.0)))), Float32(Float32(2.0) * Float32(Float32(pi) * yi))), xi)) / zi))); else tmp = fma(xi, cos(Float32(Float32(2.0) * Float32(uy * Float32(pi)))), Float32(yi * Float32(uy * fma(Float32(Float32(-1.3333333333333333) * Float32(uy * uy)), t_0, Float32(Float32(2.0) * Float32(pi)))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \pi \cdot \left(\pi \cdot \pi\right)\\
\mathbf{if}\;uy \cdot 2 \leq 0.04500000178813934:\\
\;\;\;\;zi \cdot \mathsf{fma}\left(maxCos, ux \cdot \left(1 - ux\right), \frac{\sqrt{\mathsf{fma}\left(maxCos \cdot maxCos, \left(ux \cdot ux\right) \cdot \left(\left(1 - ux\right) \cdot \left(ux + -1\right)\right), 1\right)} \cdot \mathsf{fma}\left(uy, \mathsf{fma}\left(uy, \mathsf{fma}\left(-1.3333333333333333, \left(uy \cdot yi\right) \cdot t\_0, \left(\pi \cdot \pi\right) \cdot \left(xi \cdot -2\right)\right), 2 \cdot \left(\pi \cdot yi\right)\right), xi\right)}{zi}\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(xi, \cos \left(2 \cdot \left(uy \cdot \pi\right)\right), yi \cdot \left(uy \cdot \mathsf{fma}\left(-1.3333333333333333 \cdot \left(uy \cdot uy\right), t\_0, 2 \cdot \pi\right)\right)\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.0450000018Initial program 99.2%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3299.2
Applied rewrites99.2%
Taylor expanded in zi around -inf
Applied rewrites99.0%
Taylor expanded in uy around 0
Applied rewrites98.7%
if 0.0450000018 < (*.f32 uy #s(literal 2 binary32)) Initial program 96.8%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3296.1
Applied rewrites96.1%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3291.9
Applied rewrites91.9%
Taylor expanded in uy around 0
Applied rewrites64.8%
Final simplification92.1%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(if (<= (* uy 2.0) 0.04600000008940697)
(*
zi
(fma
maxCos
(* ux (- 1.0 ux))
(/
(*
(sqrt
(fma (* maxCos maxCos) (* (* ux ux) (* (- 1.0 ux) (+ ux -1.0))) 1.0))
(fma
uy
(fma
uy
(fma
-1.3333333333333333
(* (* uy yi) (* PI (* PI PI)))
(* (* PI PI) (* xi -2.0)))
(* 2.0 (* PI yi)))
xi))
zi)))
(fma
xi
(fma (* -2.0 (* uy uy)) (* PI PI) 1.0)
(* yi (sin (* 2.0 (* uy PI)))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float tmp;
if ((uy * 2.0f) <= 0.04600000008940697f) {
tmp = zi * fmaf(maxCos, (ux * (1.0f - ux)), ((sqrtf(fmaf((maxCos * maxCos), ((ux * ux) * ((1.0f - ux) * (ux + -1.0f))), 1.0f)) * fmaf(uy, fmaf(uy, fmaf(-1.3333333333333333f, ((uy * yi) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))), ((((float) M_PI) * ((float) M_PI)) * (xi * -2.0f))), (2.0f * (((float) M_PI) * yi))), xi)) / zi));
} else {
tmp = fmaf(xi, fmaf((-2.0f * (uy * uy)), (((float) M_PI) * ((float) M_PI)), 1.0f), (yi * sinf((2.0f * (uy * ((float) M_PI))))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(0.04600000008940697)) tmp = Float32(zi * fma(maxCos, Float32(ux * Float32(Float32(1.0) - ux)), Float32(Float32(sqrt(fma(Float32(maxCos * maxCos), Float32(Float32(ux * ux) * Float32(Float32(Float32(1.0) - ux) * Float32(ux + Float32(-1.0)))), Float32(1.0))) * fma(uy, fma(uy, fma(Float32(-1.3333333333333333), Float32(Float32(uy * yi) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))), Float32(Float32(Float32(pi) * Float32(pi)) * Float32(xi * Float32(-2.0)))), Float32(Float32(2.0) * Float32(Float32(pi) * yi))), xi)) / zi))); else tmp = fma(xi, fma(Float32(Float32(-2.0) * Float32(uy * uy)), Float32(Float32(pi) * Float32(pi)), Float32(1.0)), Float32(yi * sin(Float32(Float32(2.0) * Float32(uy * Float32(pi)))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;uy \cdot 2 \leq 0.04600000008940697:\\
\;\;\;\;zi \cdot \mathsf{fma}\left(maxCos, ux \cdot \left(1 - ux\right), \frac{\sqrt{\mathsf{fma}\left(maxCos \cdot maxCos, \left(ux \cdot ux\right) \cdot \left(\left(1 - ux\right) \cdot \left(ux + -1\right)\right), 1\right)} \cdot \mathsf{fma}\left(uy, \mathsf{fma}\left(uy, \mathsf{fma}\left(-1.3333333333333333, \left(uy \cdot yi\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), \left(\pi \cdot \pi\right) \cdot \left(xi \cdot -2\right)\right), 2 \cdot \left(\pi \cdot yi\right)\right), xi\right)}{zi}\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(xi, \mathsf{fma}\left(-2 \cdot \left(uy \cdot uy\right), \pi \cdot \pi, 1\right), yi \cdot \sin \left(2 \cdot \left(uy \cdot \pi\right)\right)\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.0460000001Initial program 99.2%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3299.2
Applied rewrites99.2%
Taylor expanded in zi around -inf
Applied rewrites99.0%
Taylor expanded in uy around 0
Applied rewrites98.6%
if 0.0460000001 < (*.f32 uy #s(literal 2 binary32)) Initial program 96.8%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3296.0
Applied rewrites96.0%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3291.8
Applied rewrites91.8%
Taylor expanded in uy around 0
Applied rewrites63.1%
Final simplification91.8%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* (- 1.0 ux) maxCos)))
(if (<= (* uy 2.0) 0.00011999999696854502)
(fma
(* zi t_0)
ux
(*
(sqrt (fma t_0 (* ux (* ux (* maxCos (+ ux -1.0)))) 1.0))
(fma uy (* 2.0 (* PI yi)) xi)))
(fma
xi
(fma (* -2.0 (* uy uy)) (* PI PI) 1.0)
(* yi (sin (* 2.0 (* uy PI))))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = (1.0f - ux) * maxCos;
float tmp;
if ((uy * 2.0f) <= 0.00011999999696854502f) {
tmp = fmaf((zi * t_0), ux, (sqrtf(fmaf(t_0, (ux * (ux * (maxCos * (ux + -1.0f)))), 1.0f)) * fmaf(uy, (2.0f * (((float) M_PI) * yi)), xi)));
} else {
tmp = fmaf(xi, fmaf((-2.0f * (uy * uy)), (((float) M_PI) * ((float) M_PI)), 1.0f), (yi * sinf((2.0f * (uy * ((float) M_PI))))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(Float32(1.0) - ux) * maxCos) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(0.00011999999696854502)) tmp = fma(Float32(zi * t_0), ux, Float32(sqrt(fma(t_0, Float32(ux * Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0))))), Float32(1.0))) * fma(uy, Float32(Float32(2.0) * Float32(Float32(pi) * yi)), xi))); else tmp = fma(xi, fma(Float32(Float32(-2.0) * Float32(uy * uy)), Float32(Float32(pi) * Float32(pi)), Float32(1.0)), Float32(yi * sin(Float32(Float32(2.0) * Float32(uy * Float32(pi)))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(1 - ux\right) \cdot maxCos\\
\mathbf{if}\;uy \cdot 2 \leq 0.00011999999696854502:\\
\;\;\;\;\mathsf{fma}\left(zi \cdot t\_0, ux, \sqrt{\mathsf{fma}\left(t\_0, ux \cdot \left(ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\right), 1\right)} \cdot \mathsf{fma}\left(uy, 2 \cdot \left(\pi \cdot yi\right), xi\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(xi, \mathsf{fma}\left(-2 \cdot \left(uy \cdot uy\right), \pi \cdot \pi, 1\right), yi \cdot \sin \left(2 \cdot \left(uy \cdot \pi\right)\right)\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 1.19999997e-4Initial program 99.3%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3299.3
Applied rewrites99.3%
Taylor expanded in uy around 0
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites99.3%
Applied rewrites99.5%
if 1.19999997e-4 < (*.f32 uy #s(literal 2 binary32)) Initial program 97.9%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3297.7
Applied rewrites97.7%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3293.6
Applied rewrites93.6%
Taylor expanded in uy around 0
Applied rewrites80.7%
Final simplification91.2%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(if (<= yi 1.0000000168623835e-16)
(*
zi
(fma
maxCos
(* ux (- 1.0 ux))
(/
(*
(sqrt
(fma (* maxCos maxCos) (* (* ux ux) (* (- 1.0 ux) (+ ux -1.0))) 1.0))
(fma uy (fma 2.0 (* PI yi) (* (* xi (* PI PI)) (* uy -2.0))) xi))
zi)))
(fma xi 1.0 (* yi (sin (* 2.0 (* uy PI)))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float tmp;
if (yi <= 1.0000000168623835e-16f) {
tmp = zi * fmaf(maxCos, (ux * (1.0f - ux)), ((sqrtf(fmaf((maxCos * maxCos), ((ux * ux) * ((1.0f - ux) * (ux + -1.0f))), 1.0f)) * fmaf(uy, fmaf(2.0f, (((float) M_PI) * yi), ((xi * (((float) M_PI) * ((float) M_PI))) * (uy * -2.0f))), xi)) / zi));
} else {
tmp = fmaf(xi, 1.0f, (yi * sinf((2.0f * (uy * ((float) M_PI))))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) tmp = Float32(0.0) if (yi <= Float32(1.0000000168623835e-16)) tmp = Float32(zi * fma(maxCos, Float32(ux * Float32(Float32(1.0) - ux)), Float32(Float32(sqrt(fma(Float32(maxCos * maxCos), Float32(Float32(ux * ux) * Float32(Float32(Float32(1.0) - ux) * Float32(ux + Float32(-1.0)))), Float32(1.0))) * fma(uy, fma(Float32(2.0), Float32(Float32(pi) * yi), Float32(Float32(xi * Float32(Float32(pi) * Float32(pi))) * Float32(uy * Float32(-2.0)))), xi)) / zi))); else tmp = fma(xi, Float32(1.0), Float32(yi * sin(Float32(Float32(2.0) * Float32(uy * Float32(pi)))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;yi \leq 1.0000000168623835 \cdot 10^{-16}:\\
\;\;\;\;zi \cdot \mathsf{fma}\left(maxCos, ux \cdot \left(1 - ux\right), \frac{\sqrt{\mathsf{fma}\left(maxCos \cdot maxCos, \left(ux \cdot ux\right) \cdot \left(\left(1 - ux\right) \cdot \left(ux + -1\right)\right), 1\right)} \cdot \mathsf{fma}\left(uy, \mathsf{fma}\left(2, \pi \cdot yi, \left(xi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(uy \cdot -2\right)\right), xi\right)}{zi}\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(xi, 1, yi \cdot \sin \left(2 \cdot \left(uy \cdot \pi\right)\right)\right)\\
\end{array}
\end{array}
if yi < 1.00000002e-16Initial program 98.6%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3298.5
Applied rewrites98.5%
Taylor expanded in zi around -inf
Applied rewrites98.6%
Taylor expanded in uy around 0
Applied rewrites87.5%
if 1.00000002e-16 < yi Initial program 98.9%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3298.8
Applied rewrites98.8%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3295.8
Applied rewrites95.8%
Taylor expanded in uy around 0
Applied rewrites91.5%
Final simplification88.4%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(if (<= (* uy 2.0) 0.20000000298023224)
(*
zi
(fma
maxCos
(* ux (- 1.0 ux))
(/
(*
(sqrt
(fma (* maxCos maxCos) (* (* ux ux) (* (- 1.0 ux) (+ ux -1.0))) 1.0))
(fma uy (fma 2.0 (* PI yi) (* (* xi (* PI PI)) (* uy -2.0))) xi))
zi)))
(* (sin (* (* uy 2.0) PI)) yi)))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float tmp;
if ((uy * 2.0f) <= 0.20000000298023224f) {
tmp = zi * fmaf(maxCos, (ux * (1.0f - ux)), ((sqrtf(fmaf((maxCos * maxCos), ((ux * ux) * ((1.0f - ux) * (ux + -1.0f))), 1.0f)) * fmaf(uy, fmaf(2.0f, (((float) M_PI) * yi), ((xi * (((float) M_PI) * ((float) M_PI))) * (uy * -2.0f))), xi)) / zi));
} else {
tmp = sinf(((uy * 2.0f) * ((float) M_PI))) * yi;
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(0.20000000298023224)) tmp = Float32(zi * fma(maxCos, Float32(ux * Float32(Float32(1.0) - ux)), Float32(Float32(sqrt(fma(Float32(maxCos * maxCos), Float32(Float32(ux * ux) * Float32(Float32(Float32(1.0) - ux) * Float32(ux + Float32(-1.0)))), Float32(1.0))) * fma(uy, fma(Float32(2.0), Float32(Float32(pi) * yi), Float32(Float32(xi * Float32(Float32(pi) * Float32(pi))) * Float32(uy * Float32(-2.0)))), xi)) / zi))); else tmp = Float32(sin(Float32(Float32(uy * Float32(2.0)) * Float32(pi))) * yi); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;uy \cdot 2 \leq 0.20000000298023224:\\
\;\;\;\;zi \cdot \mathsf{fma}\left(maxCos, ux \cdot \left(1 - ux\right), \frac{\sqrt{\mathsf{fma}\left(maxCos \cdot maxCos, \left(ux \cdot ux\right) \cdot \left(\left(1 - ux\right) \cdot \left(ux + -1\right)\right), 1\right)} \cdot \mathsf{fma}\left(uy, \mathsf{fma}\left(2, \pi \cdot yi, \left(xi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(uy \cdot -2\right)\right), xi\right)}{zi}\right)\\
\mathbf{else}:\\
\;\;\;\;\sin \left(\left(uy \cdot 2\right) \cdot \pi\right) \cdot yi\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.200000003Initial program 99.1%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3299.1
Applied rewrites99.1%
Taylor expanded in zi around -inf
Applied rewrites99.0%
Taylor expanded in uy around 0
Applied rewrites92.4%
if 0.200000003 < (*.f32 uy #s(literal 2 binary32)) Initial program 95.3%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3294.0
Applied rewrites94.0%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3290.6
Applied rewrites90.6%
Taylor expanded in xi around 0
Applied rewrites52.4%
Final simplification88.0%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* (- 1.0 ux) maxCos)))
(if (<= (* uy 2.0) 7.000000186963007e-5)
(fma
(* zi t_0)
ux
(*
(sqrt (fma t_0 (* ux (* ux (* maxCos (+ ux -1.0)))) 1.0))
(fma uy (* 2.0 (* PI yi)) xi)))
(fma
uy
(fma
2.0
(* PI yi)
(*
uy
(fma
-1.3333333333333333
(* (* uy yi) (* PI (* PI PI)))
(* -2.0 (* xi (* PI PI))))))
xi))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = (1.0f - ux) * maxCos;
float tmp;
if ((uy * 2.0f) <= 7.000000186963007e-5f) {
tmp = fmaf((zi * t_0), ux, (sqrtf(fmaf(t_0, (ux * (ux * (maxCos * (ux + -1.0f)))), 1.0f)) * fmaf(uy, (2.0f * (((float) M_PI) * yi)), xi)));
} else {
tmp = fmaf(uy, fmaf(2.0f, (((float) M_PI) * yi), (uy * fmaf(-1.3333333333333333f, ((uy * yi) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))), (-2.0f * (xi * (((float) M_PI) * ((float) M_PI))))))), xi);
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(Float32(1.0) - ux) * maxCos) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(7.000000186963007e-5)) tmp = fma(Float32(zi * t_0), ux, Float32(sqrt(fma(t_0, Float32(ux * Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0))))), Float32(1.0))) * fma(uy, Float32(Float32(2.0) * Float32(Float32(pi) * yi)), xi))); else tmp = fma(uy, fma(Float32(2.0), Float32(Float32(pi) * yi), Float32(uy * fma(Float32(-1.3333333333333333), Float32(Float32(uy * yi) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))), Float32(Float32(-2.0) * Float32(xi * Float32(Float32(pi) * Float32(pi))))))), xi); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(1 - ux\right) \cdot maxCos\\
\mathbf{if}\;uy \cdot 2 \leq 7.000000186963007 \cdot 10^{-5}:\\
\;\;\;\;\mathsf{fma}\left(zi \cdot t\_0, ux, \sqrt{\mathsf{fma}\left(t\_0, ux \cdot \left(ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\right), 1\right)} \cdot \mathsf{fma}\left(uy, 2 \cdot \left(\pi \cdot yi\right), xi\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(uy, \mathsf{fma}\left(2, \pi \cdot yi, uy \cdot \mathsf{fma}\left(-1.3333333333333333, \left(uy \cdot yi\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), -2 \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right)\right)\right), xi\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 7.00000019e-5Initial program 99.4%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3299.4
Applied rewrites99.4%
Taylor expanded in uy around 0
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites99.3%
Applied rewrites99.5%
if 7.00000019e-5 < (*.f32 uy #s(literal 2 binary32)) Initial program 98.0%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3297.7
Applied rewrites97.7%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3293.8
Applied rewrites93.8%
Taylor expanded in uy around 0
Applied rewrites55.6%
Taylor expanded in uy around 0
Applied rewrites72.6%
Final simplification87.1%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* (- 1.0 ux) maxCos)))
(if (<= (* uy 2.0) 7.000000186963007e-5)
(fma
(* ux t_0)
zi
(*
(sqrt (fma t_0 (* ux (* ux (* maxCos (+ ux -1.0)))) 1.0))
(fma uy (* 2.0 (* PI yi)) xi)))
(fma
uy
(fma
2.0
(* PI yi)
(*
uy
(fma
-1.3333333333333333
(* (* uy yi) (* PI (* PI PI)))
(* -2.0 (* xi (* PI PI))))))
xi))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = (1.0f - ux) * maxCos;
float tmp;
if ((uy * 2.0f) <= 7.000000186963007e-5f) {
tmp = fmaf((ux * t_0), zi, (sqrtf(fmaf(t_0, (ux * (ux * (maxCos * (ux + -1.0f)))), 1.0f)) * fmaf(uy, (2.0f * (((float) M_PI) * yi)), xi)));
} else {
tmp = fmaf(uy, fmaf(2.0f, (((float) M_PI) * yi), (uy * fmaf(-1.3333333333333333f, ((uy * yi) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))), (-2.0f * (xi * (((float) M_PI) * ((float) M_PI))))))), xi);
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(Float32(1.0) - ux) * maxCos) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(7.000000186963007e-5)) tmp = fma(Float32(ux * t_0), zi, Float32(sqrt(fma(t_0, Float32(ux * Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0))))), Float32(1.0))) * fma(uy, Float32(Float32(2.0) * Float32(Float32(pi) * yi)), xi))); else tmp = fma(uy, fma(Float32(2.0), Float32(Float32(pi) * yi), Float32(uy * fma(Float32(-1.3333333333333333), Float32(Float32(uy * yi) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))), Float32(Float32(-2.0) * Float32(xi * Float32(Float32(pi) * Float32(pi))))))), xi); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(1 - ux\right) \cdot maxCos\\
\mathbf{if}\;uy \cdot 2 \leq 7.000000186963007 \cdot 10^{-5}:\\
\;\;\;\;\mathsf{fma}\left(ux \cdot t\_0, zi, \sqrt{\mathsf{fma}\left(t\_0, ux \cdot \left(ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\right), 1\right)} \cdot \mathsf{fma}\left(uy, 2 \cdot \left(\pi \cdot yi\right), xi\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(uy, \mathsf{fma}\left(2, \pi \cdot yi, uy \cdot \mathsf{fma}\left(-1.3333333333333333, \left(uy \cdot yi\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), -2 \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right)\right)\right), xi\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 7.00000019e-5Initial program 99.4%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3299.4
Applied rewrites99.4%
Taylor expanded in uy around 0
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites99.3%
Applied rewrites99.4%
if 7.00000019e-5 < (*.f32 uy #s(literal 2 binary32)) Initial program 98.0%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3297.7
Applied rewrites97.7%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3293.8
Applied rewrites93.8%
Taylor expanded in uy around 0
Applied rewrites55.6%
Taylor expanded in uy around 0
Applied rewrites72.6%
Final simplification87.1%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(if (<= (* uy 2.0) 7.000000186963007e-5)
(fma
(sqrt (fma (* maxCos maxCos) (* (* ux ux) (* (- 1.0 ux) (+ ux -1.0))) 1.0))
(fma 2.0 (* uy (* PI yi)) xi)
(* maxCos (* (- 1.0 ux) (* ux zi))))
(fma
uy
(fma
2.0
(* PI yi)
(*
uy
(fma
-1.3333333333333333
(* (* uy yi) (* PI (* PI PI)))
(* -2.0 (* xi (* PI PI))))))
xi)))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float tmp;
if ((uy * 2.0f) <= 7.000000186963007e-5f) {
tmp = fmaf(sqrtf(fmaf((maxCos * maxCos), ((ux * ux) * ((1.0f - ux) * (ux + -1.0f))), 1.0f)), fmaf(2.0f, (uy * (((float) M_PI) * yi)), xi), (maxCos * ((1.0f - ux) * (ux * zi))));
} else {
tmp = fmaf(uy, fmaf(2.0f, (((float) M_PI) * yi), (uy * fmaf(-1.3333333333333333f, ((uy * yi) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))), (-2.0f * (xi * (((float) M_PI) * ((float) M_PI))))))), xi);
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(7.000000186963007e-5)) tmp = fma(sqrt(fma(Float32(maxCos * maxCos), Float32(Float32(ux * ux) * Float32(Float32(Float32(1.0) - ux) * Float32(ux + Float32(-1.0)))), Float32(1.0))), fma(Float32(2.0), Float32(uy * Float32(Float32(pi) * yi)), xi), Float32(maxCos * Float32(Float32(Float32(1.0) - ux) * Float32(ux * zi)))); else tmp = fma(uy, fma(Float32(2.0), Float32(Float32(pi) * yi), Float32(uy * fma(Float32(-1.3333333333333333), Float32(Float32(uy * yi) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))), Float32(Float32(-2.0) * Float32(xi * Float32(Float32(pi) * Float32(pi))))))), xi); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;uy \cdot 2 \leq 7.000000186963007 \cdot 10^{-5}:\\
\;\;\;\;\mathsf{fma}\left(\sqrt{\mathsf{fma}\left(maxCos \cdot maxCos, \left(ux \cdot ux\right) \cdot \left(\left(1 - ux\right) \cdot \left(ux + -1\right)\right), 1\right)}, \mathsf{fma}\left(2, uy \cdot \left(\pi \cdot yi\right), xi\right), maxCos \cdot \left(\left(1 - ux\right) \cdot \left(ux \cdot zi\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(uy, \mathsf{fma}\left(2, \pi \cdot yi, uy \cdot \mathsf{fma}\left(-1.3333333333333333, \left(uy \cdot yi\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), -2 \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right)\right)\right), xi\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 7.00000019e-5Initial program 99.4%
Taylor expanded in uy around 0
+-commutativeN/A
associate-+r+N/A
associate-*r*N/A
distribute-rgt-outN/A
Applied rewrites99.2%
if 7.00000019e-5 < (*.f32 uy #s(literal 2 binary32)) Initial program 98.0%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3297.7
Applied rewrites97.7%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3293.8
Applied rewrites93.8%
Taylor expanded in uy around 0
Applied rewrites55.6%
Taylor expanded in uy around 0
Applied rewrites72.6%
Final simplification87.0%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* (- 1.0 ux) maxCos)))
(if (<= (* uy 2.0) 7.000000186963007e-5)
(fma
(* zi t_0)
ux
(*
(fma uy (* 2.0 (* PI yi)) xi)
(sqrt (fma t_0 (* (- ux) (* ux maxCos)) 1.0))))
(fma
uy
(fma
2.0
(* PI yi)
(*
uy
(fma
-1.3333333333333333
(* (* uy yi) (* PI (* PI PI)))
(* -2.0 (* xi (* PI PI))))))
xi))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = (1.0f - ux) * maxCos;
float tmp;
if ((uy * 2.0f) <= 7.000000186963007e-5f) {
tmp = fmaf((zi * t_0), ux, (fmaf(uy, (2.0f * (((float) M_PI) * yi)), xi) * sqrtf(fmaf(t_0, (-ux * (ux * maxCos)), 1.0f))));
} else {
tmp = fmaf(uy, fmaf(2.0f, (((float) M_PI) * yi), (uy * fmaf(-1.3333333333333333f, ((uy * yi) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))), (-2.0f * (xi * (((float) M_PI) * ((float) M_PI))))))), xi);
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(Float32(1.0) - ux) * maxCos) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(7.000000186963007e-5)) tmp = fma(Float32(zi * t_0), ux, Float32(fma(uy, Float32(Float32(2.0) * Float32(Float32(pi) * yi)), xi) * sqrt(fma(t_0, Float32(Float32(-ux) * Float32(ux * maxCos)), Float32(1.0))))); else tmp = fma(uy, fma(Float32(2.0), Float32(Float32(pi) * yi), Float32(uy * fma(Float32(-1.3333333333333333), Float32(Float32(uy * yi) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))), Float32(Float32(-2.0) * Float32(xi * Float32(Float32(pi) * Float32(pi))))))), xi); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(1 - ux\right) \cdot maxCos\\
\mathbf{if}\;uy \cdot 2 \leq 7.000000186963007 \cdot 10^{-5}:\\
\;\;\;\;\mathsf{fma}\left(zi \cdot t\_0, ux, \mathsf{fma}\left(uy, 2 \cdot \left(\pi \cdot yi\right), xi\right) \cdot \sqrt{\mathsf{fma}\left(t\_0, \left(-ux\right) \cdot \left(ux \cdot maxCos\right), 1\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(uy, \mathsf{fma}\left(2, \pi \cdot yi, uy \cdot \mathsf{fma}\left(-1.3333333333333333, \left(uy \cdot yi\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), -2 \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right)\right)\right), xi\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 7.00000019e-5Initial program 99.4%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3299.4
Applied rewrites99.4%
Taylor expanded in uy around 0
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites99.3%
Applied rewrites99.5%
Taylor expanded in ux around 0
Applied rewrites99.1%
if 7.00000019e-5 < (*.f32 uy #s(literal 2 binary32)) Initial program 98.0%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3297.7
Applied rewrites97.7%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3293.8
Applied rewrites93.8%
Taylor expanded in uy around 0
Applied rewrites55.6%
Taylor expanded in uy around 0
Applied rewrites72.6%
Final simplification86.9%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(if (<= (* uy 2.0) 7.000000186963007e-5)
(+
(* (* ux (* (- 1.0 ux) maxCos)) zi)
(*
(fma 2.0 (* uy (* PI yi)) xi)
(sqrt (fma (* maxCos maxCos) (* ux (- ux)) 1.0))))
(fma
uy
(fma
2.0
(* PI yi)
(*
uy
(fma
-1.3333333333333333
(* (* uy yi) (* PI (* PI PI)))
(* -2.0 (* xi (* PI PI))))))
xi)))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float tmp;
if ((uy * 2.0f) <= 7.000000186963007e-5f) {
tmp = ((ux * ((1.0f - ux) * maxCos)) * zi) + (fmaf(2.0f, (uy * (((float) M_PI) * yi)), xi) * sqrtf(fmaf((maxCos * maxCos), (ux * -ux), 1.0f)));
} else {
tmp = fmaf(uy, fmaf(2.0f, (((float) M_PI) * yi), (uy * fmaf(-1.3333333333333333f, ((uy * yi) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))), (-2.0f * (xi * (((float) M_PI) * ((float) M_PI))))))), xi);
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(7.000000186963007e-5)) tmp = Float32(Float32(Float32(ux * Float32(Float32(Float32(1.0) - ux) * maxCos)) * zi) + Float32(fma(Float32(2.0), Float32(uy * Float32(Float32(pi) * yi)), xi) * sqrt(fma(Float32(maxCos * maxCos), Float32(ux * Float32(-ux)), Float32(1.0))))); else tmp = fma(uy, fma(Float32(2.0), Float32(Float32(pi) * yi), Float32(uy * fma(Float32(-1.3333333333333333), Float32(Float32(uy * yi) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))), Float32(Float32(-2.0) * Float32(xi * Float32(Float32(pi) * Float32(pi))))))), xi); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;uy \cdot 2 \leq 7.000000186963007 \cdot 10^{-5}:\\
\;\;\;\;\left(ux \cdot \left(\left(1 - ux\right) \cdot maxCos\right)\right) \cdot zi + \mathsf{fma}\left(2, uy \cdot \left(\pi \cdot yi\right), xi\right) \cdot \sqrt{\mathsf{fma}\left(maxCos \cdot maxCos, ux \cdot \left(-ux\right), 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(uy, \mathsf{fma}\left(2, \pi \cdot yi, uy \cdot \mathsf{fma}\left(-1.3333333333333333, \left(uy \cdot yi\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), -2 \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right)\right)\right), xi\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 7.00000019e-5Initial program 99.4%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3299.4
Applied rewrites99.4%
Taylor expanded in uy around 0
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites99.3%
Taylor expanded in ux around 0
Applied rewrites98.8%
if 7.00000019e-5 < (*.f32 uy #s(literal 2 binary32)) Initial program 98.0%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3297.7
Applied rewrites97.7%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3293.8
Applied rewrites93.8%
Taylor expanded in uy around 0
Applied rewrites55.6%
Taylor expanded in uy around 0
Applied rewrites72.6%
Final simplification86.7%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(if (<= (* uy 2.0) 7.000000186963007e-5)
(+ (* (* ux (* (- 1.0 ux) maxCos)) zi) (fma 2.0 (* uy (* PI yi)) xi))
(fma
uy
(fma
2.0
(* PI yi)
(*
uy
(fma
-1.3333333333333333
(* (* uy yi) (* PI (* PI PI)))
(* -2.0 (* xi (* PI PI))))))
xi)))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float tmp;
if ((uy * 2.0f) <= 7.000000186963007e-5f) {
tmp = ((ux * ((1.0f - ux) * maxCos)) * zi) + fmaf(2.0f, (uy * (((float) M_PI) * yi)), xi);
} else {
tmp = fmaf(uy, fmaf(2.0f, (((float) M_PI) * yi), (uy * fmaf(-1.3333333333333333f, ((uy * yi) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))), (-2.0f * (xi * (((float) M_PI) * ((float) M_PI))))))), xi);
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) tmp = Float32(0.0) if (Float32(uy * Float32(2.0)) <= Float32(7.000000186963007e-5)) tmp = Float32(Float32(Float32(ux * Float32(Float32(Float32(1.0) - ux) * maxCos)) * zi) + fma(Float32(2.0), Float32(uy * Float32(Float32(pi) * yi)), xi)); else tmp = fma(uy, fma(Float32(2.0), Float32(Float32(pi) * yi), Float32(uy * fma(Float32(-1.3333333333333333), Float32(Float32(uy * yi) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))), Float32(Float32(-2.0) * Float32(xi * Float32(Float32(pi) * Float32(pi))))))), xi); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;uy \cdot 2 \leq 7.000000186963007 \cdot 10^{-5}:\\
\;\;\;\;\left(ux \cdot \left(\left(1 - ux\right) \cdot maxCos\right)\right) \cdot zi + \mathsf{fma}\left(2, uy \cdot \left(\pi \cdot yi\right), xi\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(uy, \mathsf{fma}\left(2, \pi \cdot yi, uy \cdot \mathsf{fma}\left(-1.3333333333333333, \left(uy \cdot yi\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right), -2 \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right)\right)\right), xi\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 7.00000019e-5Initial program 99.4%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3299.4
Applied rewrites99.4%
Taylor expanded in uy around 0
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites99.3%
Taylor expanded in maxCos around 0
Applied rewrites98.7%
if 7.00000019e-5 < (*.f32 uy #s(literal 2 binary32)) Initial program 98.0%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3297.7
Applied rewrites97.7%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3293.8
Applied rewrites93.8%
Taylor expanded in uy around 0
Applied rewrites55.6%
Taylor expanded in uy around 0
Applied rewrites72.6%
Final simplification86.7%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (+ (* (* ux (* (- 1.0 ux) maxCos)) zi) (fma 2.0 (* uy (* PI yi)) xi)))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return ((ux * ((1.0f - ux) * maxCos)) * zi) + fmaf(2.0f, (uy * (((float) M_PI) * yi)), xi);
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(Float32(ux * Float32(Float32(Float32(1.0) - ux) * maxCos)) * zi) + fma(Float32(2.0), Float32(uy * Float32(Float32(pi) * yi)), xi)) end
\begin{array}{l}
\\
\left(ux \cdot \left(\left(1 - ux\right) \cdot maxCos\right)\right) \cdot zi + \mathsf{fma}\left(2, uy \cdot \left(\pi \cdot yi\right), xi\right)
\end{array}
Initial program 98.7%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3298.6
Applied rewrites98.6%
Taylor expanded in uy around 0
associate-*r*N/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites80.3%
Taylor expanded in maxCos around 0
Applied rewrites80.0%
Final simplification80.0%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (fma uy (fma 2.0 (* PI yi) (* (* xi (* PI PI)) (* uy -2.0))) xi))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf(uy, fmaf(2.0f, (((float) M_PI) * yi), ((xi * (((float) M_PI) * ((float) M_PI))) * (uy * -2.0f))), xi);
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(uy, fma(Float32(2.0), Float32(Float32(pi) * yi), Float32(Float32(xi * Float32(Float32(pi) * Float32(pi))) * Float32(uy * Float32(-2.0)))), xi) end
\begin{array}{l}
\\
\mathsf{fma}\left(uy, \mathsf{fma}\left(2, \pi \cdot yi, \left(xi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(uy \cdot -2\right)\right), xi\right)
\end{array}
Initial program 98.7%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3298.6
Applied rewrites98.6%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3287.6
Applied rewrites87.6%
Taylor expanded in uy around 0
Applied rewrites74.2%
Final simplification74.2%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (fma (* uy (* 2.0 PI)) yi xi))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf((uy * (2.0f * ((float) M_PI))), yi, xi);
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(Float32(uy * Float32(Float32(2.0) * Float32(pi))), yi, xi) end
\begin{array}{l}
\\
\mathsf{fma}\left(uy \cdot \left(2 \cdot \pi\right), yi, xi\right)
\end{array}
Initial program 98.7%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3298.6
Applied rewrites98.6%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3287.6
Applied rewrites87.6%
Taylor expanded in uy around 0
Applied rewrites69.9%
Applied rewrites70.0%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (fma 2.0 (* uy (* PI yi)) xi))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf(2.0f, (uy * (((float) M_PI) * yi)), xi);
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(Float32(2.0), Float32(uy * Float32(Float32(pi) * yi)), xi) end
\begin{array}{l}
\\
\mathsf{fma}\left(2, uy \cdot \left(\pi \cdot yi\right), xi\right)
\end{array}
Initial program 98.7%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3298.6
Applied rewrites98.6%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3287.6
Applied rewrites87.6%
Taylor expanded in uy around 0
Applied rewrites69.9%
Final simplification69.9%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (* 2.0 (* uy (* PI yi))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return 2.0f * (uy * (((float) M_PI) * yi));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * yi))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = single(2.0) * (uy * (single(pi) * yi)); end
\begin{array}{l}
\\
2 \cdot \left(uy \cdot \left(\pi \cdot yi\right)\right)
\end{array}
Initial program 98.7%
lift-*.f32N/A
lift-PI.f32N/A
add-sqr-sqrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lift-PI.f32N/A
lower-sqrt.f32N/A
lift-PI.f32N/A
lower-sqrt.f3298.6
Applied rewrites98.6%
Taylor expanded in ux around 0
lower-fma.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3287.6
Applied rewrites87.6%
Taylor expanded in uy around 0
Applied rewrites69.9%
Taylor expanded in uy around inf
Applied rewrites30.6%
Final simplification30.6%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (* maxCos (* ux zi)))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return maxCos * (ux * zi);
}
real(4) function code(xi, yi, zi, ux, uy, maxcos)
real(4), intent (in) :: xi
real(4), intent (in) :: yi
real(4), intent (in) :: zi
real(4), intent (in) :: ux
real(4), intent (in) :: uy
real(4), intent (in) :: maxcos
code = maxcos * (ux * zi)
end function
function code(xi, yi, zi, ux, uy, maxCos) return Float32(maxCos * Float32(ux * zi)) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = maxCos * (ux * zi); end
\begin{array}{l}
\\
maxCos \cdot \left(ux \cdot zi\right)
\end{array}
Initial program 98.7%
Taylor expanded in zi around inf
lower-*.f32N/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f3214.7
Applied rewrites14.7%
Taylor expanded in ux around 0
Applied rewrites12.6%
herbie shell --seed 2024219
(FPCore (xi yi zi ux uy maxCos)
:name "UniformSampleCone 2"
:precision binary32
:pre (and (and (and (and (and (and (<= -10000.0 xi) (<= xi 10000.0)) (and (<= -10000.0 yi) (<= yi 10000.0))) (and (<= -10000.0 zi) (<= zi 10000.0))) (and (<= 2.328306437e-10 ux) (<= ux 1.0))) (and (<= 2.328306437e-10 uy) (<= uy 1.0))) (and (<= 0.0 maxCos) (<= maxCos 1.0)))
(+ (+ (* (* (cos (* (* uy 2.0) PI)) (sqrt (- 1.0 (* (* (* (- 1.0 ux) maxCos) ux) (* (* (- 1.0 ux) maxCos) ux))))) xi) (* (* (sin (* (* uy 2.0) PI)) (sqrt (- 1.0 (* (* (* (- 1.0 ux) maxCos) ux) (* (* (- 1.0 ux) maxCos) ux))))) yi)) (* (* (* (- 1.0 ux) maxCos) ux) zi)))