
(FPCore (xi yi zi ux uy maxCos)
: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)))
(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}
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}
Herbie found 26 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (xi yi zi ux uy maxCos)
: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)))
(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}
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}
(FPCore (xi yi zi ux uy maxCos)
: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)))
(let* ((t_0 (* (* maxCos (- 1.0 ux)) ux))
(t_1 (* PI (+ uy uy)))
(t_2 (sqrt (fma (* (- ux 1.0) maxCos) (* t_0 ux) 1.0))))
(fma (* yi t_2) (sin t_1) (fma (* xi t_2) (cos t_1) (* zi t_0)))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = (maxCos * (1.0f - ux)) * ux;
float t_1 = ((float) M_PI) * (uy + uy);
float t_2 = sqrtf(fmaf(((ux - 1.0f) * maxCos), (t_0 * ux), 1.0f));
return fmaf((yi * t_2), sinf(t_1), fmaf((xi * t_2), cosf(t_1), (zi * t_0)));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(maxCos * Float32(Float32(1.0) - ux)) * ux) t_1 = Float32(Float32(pi) * Float32(uy + uy)) t_2 = sqrt(fma(Float32(Float32(ux - Float32(1.0)) * maxCos), Float32(t_0 * ux), Float32(1.0))) return fma(Float32(yi * t_2), sin(t_1), fma(Float32(xi * t_2), cos(t_1), Float32(zi * t_0))) end
\begin{array}{l}
t_0 := \left(maxCos \cdot \left(1 - ux\right)\right) \cdot ux\\
t_1 := \pi \cdot \left(uy + uy\right)\\
t_2 := \sqrt{\mathsf{fma}\left(\left(ux - 1\right) \cdot maxCos, t\_0 \cdot ux, 1\right)}\\
\mathsf{fma}\left(yi \cdot t\_2, \sin t\_1, \mathsf{fma}\left(xi \cdot t\_2, \cos t\_1, zi \cdot t\_0\right)\right)
\end{array}
Initial program 98.9%
Applied rewrites99.0%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(fma
(- 1.0 ux)
(* (* zi ux) maxCos)
(*
(sqrt
(fma
(* (- ux 1.0) maxCos)
(* (* (* maxCos (- 1.0 ux)) ux) ux)
1.0))
(fma
(sin (* PI (+ uy uy)))
yi
(*
(cos
(+ (- (fma (+ uy uy) PI 1.5707963705062866)) 1.5707963705062866))
xi)))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf((1.0f - ux), ((zi * ux) * maxCos), (sqrtf(fmaf(((ux - 1.0f) * maxCos), (((maxCos * (1.0f - ux)) * ux) * ux), 1.0f)) * fmaf(sinf((((float) M_PI) * (uy + uy))), yi, (cosf((-fmaf((uy + uy), ((float) M_PI), 1.5707963705062866f) + 1.5707963705062866f)) * xi))));
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(Float32(Float32(1.0) - ux), Float32(Float32(zi * ux) * maxCos), Float32(sqrt(fma(Float32(Float32(ux - Float32(1.0)) * maxCos), Float32(Float32(Float32(maxCos * Float32(Float32(1.0) - ux)) * ux) * ux), Float32(1.0))) * fma(sin(Float32(Float32(pi) * Float32(uy + uy))), yi, Float32(cos(Float32(Float32(-fma(Float32(uy + uy), Float32(pi), Float32(1.5707963705062866))) + Float32(1.5707963705062866))) * xi)))) end
\mathsf{fma}\left(1 - ux, \left(zi \cdot ux\right) \cdot maxCos, \sqrt{\mathsf{fma}\left(\left(ux - 1\right) \cdot maxCos, \left(\left(maxCos \cdot \left(1 - ux\right)\right) \cdot ux\right) \cdot ux, 1\right)} \cdot \mathsf{fma}\left(\sin \left(\pi \cdot \left(uy + uy\right)\right), yi, \cos \left(\left(-\mathsf{fma}\left(uy + uy, \pi, 1.5707963705062866\right)\right) + 1.5707963705062866\right) \cdot xi\right)\right)
Initial program 98.9%
Applied rewrites99.0%
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f3299.0%
Applied rewrites99.0%
lift-cos.f32N/A
lift-*.f32N/A
lift-+.f32N/A
distribute-rgt-inN/A
lift-*.f32N/A
lift-*.f32N/A
cos-sum-revN/A
sub-negate-revN/A
sub-negateN/A
cos-sum-revN/A
lift-*.f32N/A
lift-*.f32N/A
distribute-rgt-inN/A
lift-+.f32N/A
lift-*.f32N/A
sin-+PI/2-revN/A
Applied rewrites98.9%
Evaluated real constant98.9%
Evaluated real constant98.9%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(let* ((t_0 (* PI (+ uy uy))))
(fma
(- 1.0 ux)
(* maxCos (* zi ux))
(*
(sqrt
(fma
(* (- ux 1.0) maxCos)
(* (* (* maxCos (- 1.0 ux)) ux) ux)
1.0))
(fma (sin t_0) yi (* (cos t_0) xi))))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ((float) M_PI) * (uy + uy);
return fmaf((1.0f - ux), (maxCos * (zi * ux)), (sqrtf(fmaf(((ux - 1.0f) * maxCos), (((maxCos * (1.0f - ux)) * ux) * ux), 1.0f)) * fmaf(sinf(t_0), yi, (cosf(t_0) * xi))));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(pi) * Float32(uy + uy)) return fma(Float32(Float32(1.0) - ux), Float32(maxCos * Float32(zi * ux)), Float32(sqrt(fma(Float32(Float32(ux - Float32(1.0)) * maxCos), Float32(Float32(Float32(maxCos * Float32(Float32(1.0) - ux)) * ux) * ux), Float32(1.0))) * fma(sin(t_0), yi, Float32(cos(t_0) * xi)))) end
\begin{array}{l}
t_0 := \pi \cdot \left(uy + uy\right)\\
\mathsf{fma}\left(1 - ux, maxCos \cdot \left(zi \cdot ux\right), \sqrt{\mathsf{fma}\left(\left(ux - 1\right) \cdot maxCos, \left(\left(maxCos \cdot \left(1 - ux\right)\right) \cdot ux\right) \cdot ux, 1\right)} \cdot \mathsf{fma}\left(\sin t\_0, yi, \cos t\_0 \cdot xi\right)\right)
\end{array}
Initial program 98.9%
Applied rewrites99.0%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(let* ((t_0 (* PI (+ uy uy))))
(fma
(- 1.0 ux)
(* (* zi ux) maxCos)
(*
(sqrt (fma (* (- ux 1.0) maxCos) (* (* maxCos ux) ux) 1.0))
(fma (sin t_0) yi (* (cos t_0) xi))))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ((float) M_PI) * (uy + uy);
return fmaf((1.0f - ux), ((zi * ux) * maxCos), (sqrtf(fmaf(((ux - 1.0f) * maxCos), ((maxCos * ux) * ux), 1.0f)) * fmaf(sinf(t_0), yi, (cosf(t_0) * xi))));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(pi) * Float32(uy + uy)) return fma(Float32(Float32(1.0) - ux), Float32(Float32(zi * ux) * maxCos), Float32(sqrt(fma(Float32(Float32(ux - Float32(1.0)) * maxCos), Float32(Float32(maxCos * ux) * ux), Float32(1.0))) * fma(sin(t_0), yi, Float32(cos(t_0) * xi)))) end
\begin{array}{l}
t_0 := \pi \cdot \left(uy + uy\right)\\
\mathsf{fma}\left(1 - ux, \left(zi \cdot ux\right) \cdot maxCos, \sqrt{\mathsf{fma}\left(\left(ux - 1\right) \cdot maxCos, \left(maxCos \cdot ux\right) \cdot ux, 1\right)} \cdot \mathsf{fma}\left(\sin t\_0, yi, \cos t\_0 \cdot xi\right)\right)
\end{array}
Initial program 98.9%
Applied rewrites99.0%
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f3299.0%
Applied rewrites99.0%
Taylor expanded in ux around 0
Applied rewrites98.8%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(let* ((t_0 (* (+ PI PI) uy)))
(fma
(- 1.0 ux)
(* (* zi ux) maxCos)
(fma (sin t_0) yi (* (cos t_0) xi)))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = (((float) M_PI) + ((float) M_PI)) * uy;
return fmaf((1.0f - ux), ((zi * ux) * maxCos), fmaf(sinf(t_0), yi, (cosf(t_0) * xi)));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(Float32(pi) + Float32(pi)) * uy) return fma(Float32(Float32(1.0) - ux), Float32(Float32(zi * ux) * maxCos), fma(sin(t_0), yi, Float32(cos(t_0) * xi))) end
\begin{array}{l}
t_0 := \left(\pi + \pi\right) \cdot uy\\
\mathsf{fma}\left(1 - ux, \left(zi \cdot ux\right) \cdot maxCos, \mathsf{fma}\left(\sin t\_0, yi, \cos t\_0 \cdot xi\right)\right)
\end{array}
Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in maxCos around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-fma.f32N/A
Applied rewrites98.7%
lift-fma.f32N/A
*-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
associate-*l*N/A
lift-*.f32N/A
associate-*r*N/A
lift-*.f32N/A
lower-fma.f3298.7%
lift-fma.f32N/A
lift-*.f32N/A
Applied rewrites98.7%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(let* ((t_0 (* 2.0 (* uy PI))))
(if (<= uy 0.008500000461935997)
(+
xi
(fma
maxCos
(* ux (* zi (- 1.0 ux)))
(*
uy
(fma
2.0
(* yi PI)
(*
uy
(fma
-2.0
(* xi (pow PI 2.0))
(* -1.3333333333333333 (* uy (* yi (pow PI 3.0))))))))))
(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 tmp;
if (uy <= 0.008500000461935997f) {
tmp = xi + fmaf(maxCos, (ux * (zi * (1.0f - ux))), (uy * fmaf(2.0f, (yi * ((float) M_PI)), (uy * fmaf(-2.0f, (xi * powf(((float) M_PI), 2.0f)), (-1.3333333333333333f * (uy * (yi * powf(((float) M_PI), 3.0f)))))))));
} 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))) tmp = Float32(0.0) if (uy <= Float32(0.008500000461935997)) tmp = Float32(xi + fma(maxCos, Float32(ux * Float32(zi * Float32(Float32(1.0) - ux))), Float32(uy * fma(Float32(2.0), Float32(yi * Float32(pi)), Float32(uy * fma(Float32(-2.0), Float32(xi * (Float32(pi) ^ Float32(2.0))), Float32(Float32(-1.3333333333333333) * Float32(uy * Float32(yi * (Float32(pi) ^ Float32(3.0))))))))))); else tmp = fma(maxCos, Float32(ux * zi), fma(xi, cos(t_0), Float32(yi * sin(t_0)))); end return tmp end
\begin{array}{l}
t_0 := 2 \cdot \left(uy \cdot \pi\right)\\
\mathbf{if}\;uy \leq 0.008500000461935997:\\
\;\;\;\;xi + \mathsf{fma}\left(maxCos, ux \cdot \left(zi \cdot \left(1 - ux\right)\right), uy \cdot \mathsf{fma}\left(2, yi \cdot \pi, uy \cdot \mathsf{fma}\left(-2, xi \cdot {\pi}^{2}, -1.3333333333333333 \cdot \left(uy \cdot \left(yi \cdot {\pi}^{3}\right)\right)\right)\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}
if uy < 0.00850000046Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in maxCos around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-fma.f32N/A
Applied rewrites98.7%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
Applied rewrites89.6%
if 0.00850000046 < uy Initial program 98.9%
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
Applied rewrites95.6%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(if (<= uy 0.10849999636411667)
(+
xi
(fma
maxCos
(* ux (* zi (- 1.0 ux)))
(*
uy
(fma
2.0
(* yi PI)
(*
uy
(fma
-2.0
(* xi (pow PI 2.0))
(* -1.3333333333333333 (* uy (* yi (pow PI 3.0))))))))))
(fma
(- 1.0 ux)
(* (* zi ux) maxCos)
(*
(sqrt
(fma
(* (- ux 1.0) maxCos)
(* (* (* maxCos (- 1.0 ux)) ux) ux)
1.0))
(* yi (+ (sin (* 2.0 (* uy PI))) (/ xi yi)))))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float tmp;
if (uy <= 0.10849999636411667f) {
tmp = xi + fmaf(maxCos, (ux * (zi * (1.0f - ux))), (uy * fmaf(2.0f, (yi * ((float) M_PI)), (uy * fmaf(-2.0f, (xi * powf(((float) M_PI), 2.0f)), (-1.3333333333333333f * (uy * (yi * powf(((float) M_PI), 3.0f)))))))));
} else {
tmp = fmaf((1.0f - ux), ((zi * ux) * maxCos), (sqrtf(fmaf(((ux - 1.0f) * maxCos), (((maxCos * (1.0f - ux)) * ux) * ux), 1.0f)) * (yi * (sinf((2.0f * (uy * ((float) M_PI)))) + (xi / yi)))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) tmp = Float32(0.0) if (uy <= Float32(0.10849999636411667)) tmp = Float32(xi + fma(maxCos, Float32(ux * Float32(zi * Float32(Float32(1.0) - ux))), Float32(uy * fma(Float32(2.0), Float32(yi * Float32(pi)), Float32(uy * fma(Float32(-2.0), Float32(xi * (Float32(pi) ^ Float32(2.0))), Float32(Float32(-1.3333333333333333) * Float32(uy * Float32(yi * (Float32(pi) ^ Float32(3.0))))))))))); else tmp = fma(Float32(Float32(1.0) - ux), Float32(Float32(zi * ux) * maxCos), Float32(sqrt(fma(Float32(Float32(ux - Float32(1.0)) * maxCos), Float32(Float32(Float32(maxCos * Float32(Float32(1.0) - ux)) * ux) * ux), Float32(1.0))) * Float32(yi * Float32(sin(Float32(Float32(2.0) * Float32(uy * Float32(pi)))) + Float32(xi / yi))))); end return tmp end
\begin{array}{l}
\mathbf{if}\;uy \leq 0.10849999636411667:\\
\;\;\;\;xi + \mathsf{fma}\left(maxCos, ux \cdot \left(zi \cdot \left(1 - ux\right)\right), uy \cdot \mathsf{fma}\left(2, yi \cdot \pi, uy \cdot \mathsf{fma}\left(-2, xi \cdot {\pi}^{2}, -1.3333333333333333 \cdot \left(uy \cdot \left(yi \cdot {\pi}^{3}\right)\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 - ux, \left(zi \cdot ux\right) \cdot maxCos, \sqrt{\mathsf{fma}\left(\left(ux - 1\right) \cdot maxCos, \left(\left(maxCos \cdot \left(1 - ux\right)\right) \cdot ux\right) \cdot ux, 1\right)} \cdot \left(yi \cdot \left(\sin \left(2 \cdot \left(uy \cdot \pi\right)\right) + \frac{xi}{yi}\right)\right)\right)\\
\end{array}
if uy < 0.108499996Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in maxCos around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-fma.f32N/A
Applied rewrites98.7%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
Applied rewrites89.6%
if 0.108499996 < uy Initial program 98.9%
Applied rewrites99.0%
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f3299.0%
Applied rewrites99.0%
Taylor expanded in yi around inf
lower-*.f32N/A
lower-+.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-/.f32N/A
Applied rewrites98.3%
Taylor expanded in uy around 0
lower-/.f3288.3%
Applied rewrites88.3%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(if (<= uy 0.10849999636411667)
(+
xi
(fma
maxCos
(* ux (* zi (- 1.0 ux)))
(*
uy
(fma
2.0
(* yi PI)
(*
uy
(fma
-2.0
(* xi (pow PI 2.0))
(* -1.3333333333333333 (* uy (* yi (pow PI 3.0))))))))))
(fma
(* (* maxCos (- 1.0 ux)) zi)
ux
(*
(* (sin (* (+ uy uy) PI)) yi)
(sqrt
(fma
(* (* (* (- ux 1.0) maxCos) ux) maxCos)
(* (- 1.0 ux) ux)
1.0))))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float tmp;
if (uy <= 0.10849999636411667f) {
tmp = xi + fmaf(maxCos, (ux * (zi * (1.0f - ux))), (uy * fmaf(2.0f, (yi * ((float) M_PI)), (uy * fmaf(-2.0f, (xi * powf(((float) M_PI), 2.0f)), (-1.3333333333333333f * (uy * (yi * powf(((float) M_PI), 3.0f)))))))));
} else {
tmp = fmaf(((maxCos * (1.0f - ux)) * zi), ux, ((sinf(((uy + uy) * ((float) M_PI))) * yi) * sqrtf(fmaf(((((ux - 1.0f) * maxCos) * ux) * maxCos), ((1.0f - ux) * ux), 1.0f))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) tmp = Float32(0.0) if (uy <= Float32(0.10849999636411667)) tmp = Float32(xi + fma(maxCos, Float32(ux * Float32(zi * Float32(Float32(1.0) - ux))), Float32(uy * fma(Float32(2.0), Float32(yi * Float32(pi)), Float32(uy * fma(Float32(-2.0), Float32(xi * (Float32(pi) ^ Float32(2.0))), Float32(Float32(-1.3333333333333333) * Float32(uy * Float32(yi * (Float32(pi) ^ Float32(3.0))))))))))); else tmp = fma(Float32(Float32(maxCos * Float32(Float32(1.0) - ux)) * zi), ux, Float32(Float32(sin(Float32(Float32(uy + uy) * Float32(pi))) * yi) * sqrt(fma(Float32(Float32(Float32(Float32(ux - Float32(1.0)) * maxCos) * ux) * maxCos), Float32(Float32(Float32(1.0) - ux) * ux), Float32(1.0))))); end return tmp end
\begin{array}{l}
\mathbf{if}\;uy \leq 0.10849999636411667:\\
\;\;\;\;xi + \mathsf{fma}\left(maxCos, ux \cdot \left(zi \cdot \left(1 - ux\right)\right), uy \cdot \mathsf{fma}\left(2, yi \cdot \pi, uy \cdot \mathsf{fma}\left(-2, xi \cdot {\pi}^{2}, -1.3333333333333333 \cdot \left(uy \cdot \left(yi \cdot {\pi}^{3}\right)\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(maxCos \cdot \left(1 - ux\right)\right) \cdot zi, ux, \left(\sin \left(\left(uy + uy\right) \cdot \pi\right) \cdot yi\right) \cdot \sqrt{\mathsf{fma}\left(\left(\left(\left(ux - 1\right) \cdot maxCos\right) \cdot ux\right) \cdot maxCos, \left(1 - ux\right) \cdot ux, 1\right)}\right)\\
\end{array}
if uy < 0.108499996Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in maxCos around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-fma.f32N/A
Applied rewrites98.7%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
Applied rewrites89.6%
if 0.108499996 < uy Initial program 98.9%
Applied rewrites99.0%
Taylor expanded in xi around 0
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3244.8%
Applied rewrites44.8%
Applied rewrites44.8%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(if (<= uy 0.10849999636411667)
(+
xi
(fma
maxCos
(* ux (* zi (- 1.0 ux)))
(*
uy
(fma
2.0
(* yi PI)
(*
uy
(fma
-2.0
(* xi (pow PI 2.0))
(* -1.3333333333333333 (* uy (* yi (pow PI 3.0))))))))))
(fma
(- 1.0 ux)
(* (* maxCos ux) zi)
(*
(sqrt (fma (* (- ux 1.0) maxCos) (* (* maxCos ux) ux) 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 <= 0.10849999636411667f) {
tmp = xi + fmaf(maxCos, (ux * (zi * (1.0f - ux))), (uy * fmaf(2.0f, (yi * ((float) M_PI)), (uy * fmaf(-2.0f, (xi * powf(((float) M_PI), 2.0f)), (-1.3333333333333333f * (uy * (yi * powf(((float) M_PI), 3.0f)))))))));
} else {
tmp = fmaf((1.0f - ux), ((maxCos * ux) * zi), (sqrtf(fmaf(((ux - 1.0f) * maxCos), ((maxCos * ux) * ux), 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 (uy <= Float32(0.10849999636411667)) tmp = Float32(xi + fma(maxCos, Float32(ux * Float32(zi * Float32(Float32(1.0) - ux))), Float32(uy * fma(Float32(2.0), Float32(yi * Float32(pi)), Float32(uy * fma(Float32(-2.0), Float32(xi * (Float32(pi) ^ Float32(2.0))), Float32(Float32(-1.3333333333333333) * Float32(uy * Float32(yi * (Float32(pi) ^ Float32(3.0))))))))))); else tmp = fma(Float32(Float32(1.0) - ux), Float32(Float32(maxCos * ux) * zi), Float32(sqrt(fma(Float32(Float32(ux - Float32(1.0)) * maxCos), Float32(Float32(maxCos * ux) * ux), Float32(1.0))) * Float32(yi * sin(Float32(Float32(2.0) * Float32(uy * Float32(pi))))))); end return tmp end
\begin{array}{l}
\mathbf{if}\;uy \leq 0.10849999636411667:\\
\;\;\;\;xi + \mathsf{fma}\left(maxCos, ux \cdot \left(zi \cdot \left(1 - ux\right)\right), uy \cdot \mathsf{fma}\left(2, yi \cdot \pi, uy \cdot \mathsf{fma}\left(-2, xi \cdot {\pi}^{2}, -1.3333333333333333 \cdot \left(uy \cdot \left(yi \cdot {\pi}^{3}\right)\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 - ux, \left(maxCos \cdot ux\right) \cdot zi, \sqrt{\mathsf{fma}\left(\left(ux - 1\right) \cdot maxCos, \left(maxCos \cdot ux\right) \cdot ux, 1\right)} \cdot \left(yi \cdot \sin \left(2 \cdot \left(uy \cdot \pi\right)\right)\right)\right)\\
\end{array}
if uy < 0.108499996Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in maxCos around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-fma.f32N/A
Applied rewrites98.7%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
Applied rewrites89.6%
if 0.108499996 < uy Initial program 98.9%
Applied rewrites99.0%
Taylor expanded in xi around 0
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3244.8%
Applied rewrites44.8%
Taylor expanded in ux around 0
Applied rewrites44.7%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(let* ((t_0 (* (- ux 1.0) maxCos)))
(if (<= uy 0.009999999776482582)
(fma
(- 1.0 ux)
(* (* zi ux) maxCos)
(*
(sqrt (fma t_0 (* (* (* maxCos (- 1.0 ux)) ux) ux) 1.0))
(+
xi
(*
uy
(fma -2.0 (* uy (* xi (pow PI 2.0))) (* 2.0 (* yi PI)))))))
(fma
(- 1.0 ux)
(* (* maxCos ux) zi)
(*
(sqrt (fma t_0 (* (* maxCos ux) ux) 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 = (ux - 1.0f) * maxCos;
float tmp;
if (uy <= 0.009999999776482582f) {
tmp = fmaf((1.0f - ux), ((zi * ux) * maxCos), (sqrtf(fmaf(t_0, (((maxCos * (1.0f - ux)) * ux) * ux), 1.0f)) * (xi + (uy * fmaf(-2.0f, (uy * (xi * powf(((float) M_PI), 2.0f))), (2.0f * (yi * ((float) M_PI))))))));
} else {
tmp = fmaf((1.0f - ux), ((maxCos * ux) * zi), (sqrtf(fmaf(t_0, ((maxCos * ux) * ux), 1.0f)) * (yi * sinf((2.0f * (uy * ((float) M_PI)))))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(ux - Float32(1.0)) * maxCos) tmp = Float32(0.0) if (uy <= Float32(0.009999999776482582)) tmp = fma(Float32(Float32(1.0) - ux), Float32(Float32(zi * ux) * maxCos), Float32(sqrt(fma(t_0, Float32(Float32(Float32(maxCos * Float32(Float32(1.0) - ux)) * ux) * ux), Float32(1.0))) * Float32(xi + Float32(uy * fma(Float32(-2.0), Float32(uy * Float32(xi * (Float32(pi) ^ Float32(2.0)))), Float32(Float32(2.0) * Float32(yi * Float32(pi)))))))); else tmp = fma(Float32(Float32(1.0) - ux), Float32(Float32(maxCos * ux) * zi), Float32(sqrt(fma(t_0, Float32(Float32(maxCos * ux) * ux), Float32(1.0))) * Float32(yi * sin(Float32(Float32(2.0) * Float32(uy * Float32(pi))))))); end return tmp end
\begin{array}{l}
t_0 := \left(ux - 1\right) \cdot maxCos\\
\mathbf{if}\;uy \leq 0.009999999776482582:\\
\;\;\;\;\mathsf{fma}\left(1 - ux, \left(zi \cdot ux\right) \cdot maxCos, \sqrt{\mathsf{fma}\left(t\_0, \left(\left(maxCos \cdot \left(1 - ux\right)\right) \cdot ux\right) \cdot ux, 1\right)} \cdot \left(xi + uy \cdot \mathsf{fma}\left(-2, uy \cdot \left(xi \cdot {\pi}^{2}\right), 2 \cdot \left(yi \cdot \pi\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 - ux, \left(maxCos \cdot ux\right) \cdot zi, \sqrt{\mathsf{fma}\left(t\_0, \left(maxCos \cdot ux\right) \cdot ux, 1\right)} \cdot \left(yi \cdot \sin \left(2 \cdot \left(uy \cdot \pi\right)\right)\right)\right)\\
\end{array}
if uy < 0.00999999978Initial program 98.9%
Applied rewrites99.0%
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f3299.0%
Applied rewrites99.0%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-*.f32N/A
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-pow.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3286.1%
Applied rewrites86.1%
if 0.00999999978 < uy Initial program 98.9%
Applied rewrites99.0%
Taylor expanded in xi around 0
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3244.8%
Applied rewrites44.8%
Taylor expanded in ux around 0
Applied rewrites44.7%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(let* ((t_0 (* (- ux 1.0) maxCos)) (t_1 (* (* maxCos ux) zi)))
(if (<= uy 0.009999999776482582)
(fma
(- 1.0 ux)
t_1
(*
(sqrt (fma t_0 (* (* (* maxCos (- 1.0 ux)) ux) ux) 1.0))
(+
xi
(*
uy
(fma -2.0 (* uy (* xi (pow PI 2.0))) (* 2.0 (* yi PI)))))))
(fma
(- 1.0 ux)
t_1
(*
(sqrt (fma t_0 (* (* maxCos ux) ux) 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 = (ux - 1.0f) * maxCos;
float t_1 = (maxCos * ux) * zi;
float tmp;
if (uy <= 0.009999999776482582f) {
tmp = fmaf((1.0f - ux), t_1, (sqrtf(fmaf(t_0, (((maxCos * (1.0f - ux)) * ux) * ux), 1.0f)) * (xi + (uy * fmaf(-2.0f, (uy * (xi * powf(((float) M_PI), 2.0f))), (2.0f * (yi * ((float) M_PI))))))));
} else {
tmp = fmaf((1.0f - ux), t_1, (sqrtf(fmaf(t_0, ((maxCos * ux) * ux), 1.0f)) * (yi * sinf((2.0f * (uy * ((float) M_PI)))))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(ux - Float32(1.0)) * maxCos) t_1 = Float32(Float32(maxCos * ux) * zi) tmp = Float32(0.0) if (uy <= Float32(0.009999999776482582)) tmp = fma(Float32(Float32(1.0) - ux), t_1, Float32(sqrt(fma(t_0, Float32(Float32(Float32(maxCos * Float32(Float32(1.0) - ux)) * ux) * ux), Float32(1.0))) * Float32(xi + Float32(uy * fma(Float32(-2.0), Float32(uy * Float32(xi * (Float32(pi) ^ Float32(2.0)))), Float32(Float32(2.0) * Float32(yi * Float32(pi)))))))); else tmp = fma(Float32(Float32(1.0) - ux), t_1, Float32(sqrt(fma(t_0, Float32(Float32(maxCos * ux) * ux), Float32(1.0))) * Float32(yi * sin(Float32(Float32(2.0) * Float32(uy * Float32(pi))))))); end return tmp end
\begin{array}{l}
t_0 := \left(ux - 1\right) \cdot maxCos\\
t_1 := \left(maxCos \cdot ux\right) \cdot zi\\
\mathbf{if}\;uy \leq 0.009999999776482582:\\
\;\;\;\;\mathsf{fma}\left(1 - ux, t\_1, \sqrt{\mathsf{fma}\left(t\_0, \left(\left(maxCos \cdot \left(1 - ux\right)\right) \cdot ux\right) \cdot ux, 1\right)} \cdot \left(xi + uy \cdot \mathsf{fma}\left(-2, uy \cdot \left(xi \cdot {\pi}^{2}\right), 2 \cdot \left(yi \cdot \pi\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 - ux, t\_1, \sqrt{\mathsf{fma}\left(t\_0, \left(maxCos \cdot ux\right) \cdot ux, 1\right)} \cdot \left(yi \cdot \sin \left(2 \cdot \left(uy \cdot \pi\right)\right)\right)\right)\\
\end{array}
if uy < 0.00999999978Initial program 98.9%
Applied rewrites99.0%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-*.f32N/A
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-pow.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3286.1%
Applied rewrites86.1%
if 0.00999999978 < uy Initial program 98.9%
Applied rewrites99.0%
Taylor expanded in xi around 0
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3244.8%
Applied rewrites44.8%
Taylor expanded in ux around 0
Applied rewrites44.7%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(if (<= uy 0.009999999776482582)
(fma
maxCos
(* ux (* zi (- 1.0 ux)))
(+
xi
(* uy (fma -2.0 (* uy (* xi (pow PI 2.0))) (* 2.0 (* yi PI))))))
(fma
(- 1.0 ux)
(* (* maxCos ux) zi)
(*
(sqrt (fma (* (- ux 1.0) maxCos) (* (* maxCos ux) ux) 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 <= 0.009999999776482582f) {
tmp = fmaf(maxCos, (ux * (zi * (1.0f - ux))), (xi + (uy * fmaf(-2.0f, (uy * (xi * powf(((float) M_PI), 2.0f))), (2.0f * (yi * ((float) M_PI)))))));
} else {
tmp = fmaf((1.0f - ux), ((maxCos * ux) * zi), (sqrtf(fmaf(((ux - 1.0f) * maxCos), ((maxCos * ux) * ux), 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 (uy <= Float32(0.009999999776482582)) tmp = fma(maxCos, Float32(ux * Float32(zi * Float32(Float32(1.0) - ux))), Float32(xi + Float32(uy * fma(Float32(-2.0), Float32(uy * Float32(xi * (Float32(pi) ^ Float32(2.0)))), Float32(Float32(2.0) * Float32(yi * Float32(pi))))))); else tmp = fma(Float32(Float32(1.0) - ux), Float32(Float32(maxCos * ux) * zi), Float32(sqrt(fma(Float32(Float32(ux - Float32(1.0)) * maxCos), Float32(Float32(maxCos * ux) * ux), Float32(1.0))) * Float32(yi * sin(Float32(Float32(2.0) * Float32(uy * Float32(pi))))))); end return tmp end
\begin{array}{l}
\mathbf{if}\;uy \leq 0.009999999776482582:\\
\;\;\;\;\mathsf{fma}\left(maxCos, ux \cdot \left(zi \cdot \left(1 - ux\right)\right), xi + uy \cdot \mathsf{fma}\left(-2, uy \cdot \left(xi \cdot {\pi}^{2}\right), 2 \cdot \left(yi \cdot \pi\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 - ux, \left(maxCos \cdot ux\right) \cdot zi, \sqrt{\mathsf{fma}\left(\left(ux - 1\right) \cdot maxCos, \left(maxCos \cdot ux\right) \cdot ux, 1\right)} \cdot \left(yi \cdot \sin \left(2 \cdot \left(uy \cdot \pi\right)\right)\right)\right)\\
\end{array}
if uy < 0.00999999978Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in maxCos around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-fma.f32N/A
Applied rewrites98.7%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-*.f32N/A
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-pow.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3285.9%
Applied rewrites85.9%
if 0.00999999978 < uy Initial program 98.9%
Applied rewrites99.0%
Taylor expanded in xi around 0
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3244.8%
Applied rewrites44.8%
Taylor expanded in ux around 0
Applied rewrites44.7%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(let* ((t_0 (* ux (* zi (- 1.0 ux)))))
(if (<= uy 0.009999999776482582)
(fma
maxCos
t_0
(+
xi
(* uy (fma -2.0 (* uy (* xi (pow PI 2.0))) (* 2.0 (* yi PI))))))
(fma maxCos t_0 (* yi (sin (* 2.0 (* uy PI))))))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ux * (zi * (1.0f - ux));
float tmp;
if (uy <= 0.009999999776482582f) {
tmp = fmaf(maxCos, t_0, (xi + (uy * fmaf(-2.0f, (uy * (xi * powf(((float) M_PI), 2.0f))), (2.0f * (yi * ((float) M_PI)))))));
} else {
tmp = fmaf(maxCos, t_0, (yi * sinf((2.0f * (uy * ((float) M_PI))))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(ux * Float32(zi * Float32(Float32(1.0) - ux))) tmp = Float32(0.0) if (uy <= Float32(0.009999999776482582)) tmp = fma(maxCos, t_0, Float32(xi + Float32(uy * fma(Float32(-2.0), Float32(uy * Float32(xi * (Float32(pi) ^ Float32(2.0)))), Float32(Float32(2.0) * Float32(yi * Float32(pi))))))); else tmp = fma(maxCos, t_0, Float32(yi * sin(Float32(Float32(2.0) * Float32(uy * Float32(pi)))))); end return tmp end
\begin{array}{l}
t_0 := ux \cdot \left(zi \cdot \left(1 - ux\right)\right)\\
\mathbf{if}\;uy \leq 0.009999999776482582:\\
\;\;\;\;\mathsf{fma}\left(maxCos, t\_0, xi + uy \cdot \mathsf{fma}\left(-2, uy \cdot \left(xi \cdot {\pi}^{2}\right), 2 \cdot \left(yi \cdot \pi\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(maxCos, t\_0, yi \cdot \sin \left(2 \cdot \left(uy \cdot \pi\right)\right)\right)\\
\end{array}
if uy < 0.00999999978Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in maxCos around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-fma.f32N/A
Applied rewrites98.7%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-*.f32N/A
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-pow.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3285.9%
Applied rewrites85.9%
if 0.00999999978 < uy Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in maxCos around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-fma.f32N/A
Applied rewrites98.7%
Taylor expanded in xi around 0
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3244.7%
Applied rewrites44.7%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(let* ((t_0 (* ux (* zi (- 1.0 ux)))))
(if (<= uy 0.30000001192092896)
(fma
maxCos
t_0
(+
xi
(* uy (fma -2.0 (* uy (* xi (pow PI 2.0))) (* 2.0 (* yi PI))))))
(fma maxCos t_0 (* xi (cos (* 2.0 (* uy PI))))))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ux * (zi * (1.0f - ux));
float tmp;
if (uy <= 0.30000001192092896f) {
tmp = fmaf(maxCos, t_0, (xi + (uy * fmaf(-2.0f, (uy * (xi * powf(((float) M_PI), 2.0f))), (2.0f * (yi * ((float) M_PI)))))));
} else {
tmp = fmaf(maxCos, t_0, (xi * cosf((2.0f * (uy * ((float) M_PI))))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(ux * Float32(zi * Float32(Float32(1.0) - ux))) tmp = Float32(0.0) if (uy <= Float32(0.30000001192092896)) tmp = fma(maxCos, t_0, Float32(xi + Float32(uy * fma(Float32(-2.0), Float32(uy * Float32(xi * (Float32(pi) ^ Float32(2.0)))), Float32(Float32(2.0) * Float32(yi * Float32(pi))))))); else tmp = fma(maxCos, t_0, Float32(xi * cos(Float32(Float32(2.0) * Float32(uy * Float32(pi)))))); end return tmp end
\begin{array}{l}
t_0 := ux \cdot \left(zi \cdot \left(1 - ux\right)\right)\\
\mathbf{if}\;uy \leq 0.30000001192092896:\\
\;\;\;\;\mathsf{fma}\left(maxCos, t\_0, xi + uy \cdot \mathsf{fma}\left(-2, uy \cdot \left(xi \cdot {\pi}^{2}\right), 2 \cdot \left(yi \cdot \pi\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(maxCos, t\_0, xi \cdot \cos \left(2 \cdot \left(uy \cdot \pi\right)\right)\right)\\
\end{array}
if uy < 0.300000012Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in maxCos around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-fma.f32N/A
Applied rewrites98.7%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-*.f32N/A
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-pow.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3285.9%
Applied rewrites85.9%
if 0.300000012 < uy Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in maxCos around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-fma.f32N/A
Applied rewrites98.7%
Taylor expanded in yi around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-cos.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3259.1%
Applied rewrites59.1%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(fma
maxCos
(* ux (* zi (- 1.0 ux)))
(+
xi
(* uy (fma -2.0 (* uy (* xi (pow PI 2.0))) (* 2.0 (* yi PI)))))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf(maxCos, (ux * (zi * (1.0f - ux))), (xi + (uy * fmaf(-2.0f, (uy * (xi * powf(((float) M_PI), 2.0f))), (2.0f * (yi * ((float) M_PI)))))));
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(maxCos, Float32(ux * Float32(zi * Float32(Float32(1.0) - ux))), Float32(xi + Float32(uy * fma(Float32(-2.0), Float32(uy * Float32(xi * (Float32(pi) ^ Float32(2.0)))), Float32(Float32(2.0) * Float32(yi * Float32(pi))))))) end
\mathsf{fma}\left(maxCos, ux \cdot \left(zi \cdot \left(1 - ux\right)\right), xi + uy \cdot \mathsf{fma}\left(-2, uy \cdot \left(xi \cdot {\pi}^{2}\right), 2 \cdot \left(yi \cdot \pi\right)\right)\right)
Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in maxCos around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-fma.f32N/A
Applied rewrites98.7%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-*.f32N/A
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-pow.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3285.9%
Applied rewrites85.9%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(fma
(- 1.0 ux)
(* (* zi ux) maxCos)
(*
(sqrt
(fma
(* (- ux 1.0) maxCos)
(* (* (* maxCos (- 1.0 ux)) ux) ux)
1.0))
(fma (+ uy uy) (* yi PI) xi))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf((1.0f - ux), ((zi * ux) * maxCos), (sqrtf(fmaf(((ux - 1.0f) * maxCos), (((maxCos * (1.0f - ux)) * ux) * ux), 1.0f)) * fmaf((uy + uy), (yi * ((float) M_PI)), xi)));
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(Float32(Float32(1.0) - ux), Float32(Float32(zi * ux) * maxCos), Float32(sqrt(fma(Float32(Float32(ux - Float32(1.0)) * maxCos), Float32(Float32(Float32(maxCos * Float32(Float32(1.0) - ux)) * ux) * ux), Float32(1.0))) * fma(Float32(uy + uy), Float32(yi * Float32(pi)), xi))) end
\mathsf{fma}\left(1 - ux, \left(zi \cdot ux\right) \cdot maxCos, \sqrt{\mathsf{fma}\left(\left(ux - 1\right) \cdot maxCos, \left(\left(maxCos \cdot \left(1 - ux\right)\right) \cdot ux\right) \cdot ux, 1\right)} \cdot \mathsf{fma}\left(uy + uy, yi \cdot \pi, xi\right)\right)
Initial program 98.9%
Applied rewrites99.0%
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f3299.0%
Applied rewrites99.0%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3281.9%
Applied rewrites81.9%
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
lift-*.f32N/A
associate-*r*N/A
count-2N/A
lift-+.f32N/A
lower-fma.f3282.0%
Applied rewrites82.0%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(fma
(- 1.0 ux)
(* (* zi ux) maxCos)
(*
(sqrt
(fma
(* (- ux 1.0) maxCos)
(* (* (* maxCos (- 1.0 ux)) ux) ux)
1.0))
(fma (* (+ PI PI) uy) yi xi))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf((1.0f - ux), ((zi * ux) * maxCos), (sqrtf(fmaf(((ux - 1.0f) * maxCos), (((maxCos * (1.0f - ux)) * ux) * ux), 1.0f)) * fmaf(((((float) M_PI) + ((float) M_PI)) * uy), yi, xi)));
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(Float32(Float32(1.0) - ux), Float32(Float32(zi * ux) * maxCos), Float32(sqrt(fma(Float32(Float32(ux - Float32(1.0)) * maxCos), Float32(Float32(Float32(maxCos * Float32(Float32(1.0) - ux)) * ux) * ux), Float32(1.0))) * fma(Float32(Float32(Float32(pi) + Float32(pi)) * uy), yi, xi))) end
\mathsf{fma}\left(1 - ux, \left(zi \cdot ux\right) \cdot maxCos, \sqrt{\mathsf{fma}\left(\left(ux - 1\right) \cdot maxCos, \left(\left(maxCos \cdot \left(1 - ux\right)\right) \cdot ux\right) \cdot ux, 1\right)} \cdot \mathsf{fma}\left(\left(\pi + \pi\right) \cdot uy, yi, xi\right)\right)
Initial program 98.9%
Applied rewrites99.0%
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f3299.0%
Applied rewrites99.0%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3281.9%
Applied rewrites81.9%
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
lift-*.f32N/A
associate-*r*N/A
count-2N/A
lift-+.f32N/A
lift-*.f32N/A
*-commutativeN/A
associate-*r*N/A
lift-+.f32N/A
count-2N/A
associate-*r*N/A
lift-*.f32N/A
lift-*.f32N/A
lower-fma.f3282.0%
lift-*.f32N/A
lift-*.f32N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
count-2-revN/A
lower-+.f3282.0%
Applied rewrites82.0%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(fma
(- 1.0 ux)
(* (* zi ux) maxCos)
(*
(sqrt
(fma
(* (- ux 1.0) maxCos)
(* (* (* maxCos (- 1.0 ux)) ux) ux)
1.0))
(fma (* (+ uy uy) yi) PI xi))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf((1.0f - ux), ((zi * ux) * maxCos), (sqrtf(fmaf(((ux - 1.0f) * maxCos), (((maxCos * (1.0f - ux)) * ux) * ux), 1.0f)) * fmaf(((uy + uy) * yi), ((float) M_PI), xi)));
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(Float32(Float32(1.0) - ux), Float32(Float32(zi * ux) * maxCos), Float32(sqrt(fma(Float32(Float32(ux - Float32(1.0)) * maxCos), Float32(Float32(Float32(maxCos * Float32(Float32(1.0) - ux)) * ux) * ux), Float32(1.0))) * fma(Float32(Float32(uy + uy) * yi), Float32(pi), xi))) end
\mathsf{fma}\left(1 - ux, \left(zi \cdot ux\right) \cdot maxCos, \sqrt{\mathsf{fma}\left(\left(ux - 1\right) \cdot maxCos, \left(\left(maxCos \cdot \left(1 - ux\right)\right) \cdot ux\right) \cdot ux, 1\right)} \cdot \mathsf{fma}\left(\left(uy + uy\right) \cdot yi, \pi, xi\right)\right)
Initial program 98.9%
Applied rewrites99.0%
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f3299.0%
Applied rewrites99.0%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3281.9%
Applied rewrites81.9%
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
lift-*.f32N/A
associate-*r*N/A
count-2N/A
lift-+.f32N/A
lift-*.f32N/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f3281.9%
Applied rewrites81.9%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(fma
(- 1.0 ux)
(* (* zi ux) maxCos)
(*
(sqrt (fma (* (- ux 1.0) maxCos) (* (* maxCos ux) ux) 1.0))
(+ xi (* 2.0 (* uy (* yi PI)))))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf((1.0f - ux), ((zi * ux) * maxCos), (sqrtf(fmaf(((ux - 1.0f) * maxCos), ((maxCos * ux) * ux), 1.0f)) * (xi + (2.0f * (uy * (yi * ((float) M_PI)))))));
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(Float32(Float32(1.0) - ux), Float32(Float32(zi * ux) * maxCos), Float32(sqrt(fma(Float32(Float32(ux - Float32(1.0)) * maxCos), Float32(Float32(maxCos * ux) * ux), Float32(1.0))) * Float32(xi + Float32(Float32(2.0) * Float32(uy * Float32(yi * Float32(pi))))))) end
\mathsf{fma}\left(1 - ux, \left(zi \cdot ux\right) \cdot maxCos, \sqrt{\mathsf{fma}\left(\left(ux - 1\right) \cdot maxCos, \left(maxCos \cdot ux\right) \cdot ux, 1\right)} \cdot \left(xi + 2 \cdot \left(uy \cdot \left(yi \cdot \pi\right)\right)\right)\right)
Initial program 98.9%
Applied rewrites99.0%
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f3299.0%
Applied rewrites99.0%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3281.9%
Applied rewrites81.9%
Taylor expanded in ux around 0
Applied rewrites81.8%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(+ xi (fma 2.0 (* uy (* yi PI)) (* maxCos (* ux (* zi (- 1.0 ux)))))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return xi + fmaf(2.0f, (uy * (yi * ((float) M_PI))), (maxCos * (ux * (zi * (1.0f - ux)))));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(xi + fma(Float32(2.0), Float32(uy * Float32(yi * Float32(pi))), Float32(maxCos * Float32(ux * Float32(zi * Float32(Float32(1.0) - ux)))))) end
xi + \mathsf{fma}\left(2, uy \cdot \left(yi \cdot \pi\right), maxCos \cdot \left(ux \cdot \left(zi \cdot \left(1 - ux\right)\right)\right)\right)
Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in maxCos around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-fma.f32N/A
Applied rewrites98.7%
Taylor expanded in uy around 0
lower-+.f32N/A
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f3281.8%
Applied rewrites81.8%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(+ xi (* maxCos (* ux (* zi (- 1.0 ux))))))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return xi + (maxCos * (ux * (zi * (1.0f - ux))));
}
real(4) function code(xi, yi, zi, ux, uy, maxcos)
use fmin_fmax_functions
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 = xi + (maxcos * (ux * (zi * (1.0e0 - ux))))
end function
function code(xi, yi, zi, ux, uy, maxCos) return Float32(xi + Float32(maxCos * Float32(ux * Float32(zi * Float32(Float32(1.0) - ux))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = xi + (maxCos * (ux * (zi * (single(1.0) - ux)))); end
xi + maxCos \cdot \left(ux \cdot \left(zi \cdot \left(1 - ux\right)\right)\right)
Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in ux around 0
lower-+.f32N/A
lower-*.f32N/A
lower-*.f3249.4%
Applied rewrites49.4%
Taylor expanded in maxCos around 0
lower-+.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f3251.6%
Applied rewrites51.6%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(+ xi (* (* zi maxCos) ux)))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return xi + ((zi * maxCos) * ux);
}
real(4) function code(xi, yi, zi, ux, uy, maxcos)
use fmin_fmax_functions
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 = xi + ((zi * maxcos) * ux)
end function
function code(xi, yi, zi, ux, uy, maxCos) return Float32(xi + Float32(Float32(zi * maxCos) * ux)) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = xi + ((zi * maxCos) * ux); end
xi + \left(zi \cdot maxCos\right) \cdot ux
Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in ux around 0
lower-+.f32N/A
lower-*.f32N/A
lower-*.f3249.4%
Applied rewrites49.4%
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f3249.4%
Applied rewrites49.4%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(fma (* zi maxCos) ux xi))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf((zi * maxCos), ux, xi);
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(Float32(zi * maxCos), ux, xi) end
\mathsf{fma}\left(zi \cdot maxCos, ux, xi\right)
Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in ux around 0
lower-+.f32N/A
lower-*.f32N/A
lower-*.f3249.4%
Applied rewrites49.4%
lift-+.f32N/A
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
lift-*.f32N/A
lift-*.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f3249.4%
Applied rewrites49.4%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(fma (* maxCos ux) zi xi))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf((maxCos * ux), zi, xi);
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(Float32(maxCos * ux), zi, xi) end
\mathsf{fma}\left(maxCos \cdot ux, zi, xi\right)
Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in ux around 0
lower-+.f32N/A
lower-*.f32N/A
lower-*.f3249.4%
Applied rewrites49.4%
lift-+.f32N/A
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
lift-*.f32N/A
lift-*.f32N/A
+-commutativeN/A
lift-*.f32N/A
lower-fma.f3249.4%
Applied rewrites49.4%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(* (* maxCos zi) ux))float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return (maxCos * zi) * ux;
}
real(4) function code(xi, yi, zi, ux, uy, maxcos)
use fmin_fmax_functions
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 * zi) * ux
end function
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(maxCos * zi) * ux) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = (maxCos * zi) * ux; end
\left(maxCos \cdot zi\right) \cdot ux
Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in ux around 0
lower-+.f32N/A
lower-*.f32N/A
lower-*.f3249.4%
Applied rewrites49.4%
Taylor expanded in xi around 0
lower-*.f32N/A
lower-*.f3212.1%
Applied rewrites12.1%
lift-*.f32N/A
lift-*.f32N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f3212.1%
Applied rewrites12.1%
(FPCore (xi yi zi ux uy maxCos)
: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)))
(* 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)
use fmin_fmax_functions
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
maxCos \cdot \left(ux \cdot zi\right)
Initial program 98.9%
Taylor expanded in uy around 0
lower-fma.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower--.f32N/A
Applied rewrites51.6%
Taylor expanded in ux around 0
lower-+.f32N/A
lower-*.f32N/A
lower-*.f3249.4%
Applied rewrites49.4%
Taylor expanded in xi around 0
lower-*.f32N/A
lower-*.f3212.1%
Applied rewrites12.1%
herbie shell --seed 2025356
(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)))