
(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 29 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))) (t_1 (* 2.0 (* PI uy))))
(fma
(* zi t_0)
maxCos
(*
(pow (+ 1.0 (* t_0 (* maxCos (* ux (* maxCos (+ ux -1.0)))))) 0.5)
(+ (* xi (cos t_1)) (* yi (sin t_1)))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ux * (1.0f - ux);
float t_1 = 2.0f * (((float) M_PI) * uy);
return fmaf((zi * t_0), maxCos, (powf((1.0f + (t_0 * (maxCos * (ux * (maxCos * (ux + -1.0f)))))), 0.5f) * ((xi * cosf(t_1)) + (yi * sinf(t_1)))));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(ux * Float32(Float32(1.0) - ux)) t_1 = Float32(Float32(2.0) * Float32(Float32(pi) * uy)) return fma(Float32(zi * t_0), maxCos, Float32((Float32(Float32(1.0) + Float32(t_0 * Float32(maxCos * Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0))))))) ^ Float32(0.5)) * Float32(Float32(xi * cos(t_1)) + Float32(yi * sin(t_1))))) end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := ux \cdot \left(1 - ux\right)\\
t_1 := 2 \cdot \left(\pi \cdot uy\right)\\
\mathsf{fma}\left(zi \cdot t\_0, maxCos, {\left(1 + t\_0 \cdot \left(maxCos \cdot \left(ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\right)\right)\right)}^{0.5} \cdot \left(xi \cdot \cos t\_1 + yi \cdot \sin t\_1\right)\right)
\end{array}
\end{array}
Initial program 99.1%
Simplified99.1%
Applied egg-rr99.1%
Final simplification99.1%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* uy (* 2.0 PI))))
(+
(*
(sqrt
(+ 1.0 (* (- 1.0 ux) (* ux (* maxCos (* ux (* maxCos (+ ux -1.0))))))))
(+ (* yi (sin t_0)) (* xi (cos t_0))))
(* zi (* ux (* ux (- (/ maxCos ux) maxCos)))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = uy * (2.0f * ((float) M_PI));
return (sqrtf((1.0f + ((1.0f - ux) * (ux * (maxCos * (ux * (maxCos * (ux + -1.0f)))))))) * ((yi * sinf(t_0)) + (xi * cosf(t_0)))) + (zi * (ux * (ux * ((maxCos / ux) - maxCos))));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(uy * Float32(Float32(2.0) * Float32(pi))) return Float32(Float32(sqrt(Float32(Float32(1.0) + Float32(Float32(Float32(1.0) - ux) * Float32(ux * Float32(maxCos * Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0))))))))) * Float32(Float32(yi * sin(t_0)) + Float32(xi * cos(t_0)))) + Float32(zi * Float32(ux * Float32(ux * Float32(Float32(maxCos / ux) - maxCos))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = uy * (single(2.0) * single(pi)); tmp = (sqrt((single(1.0) + ((single(1.0) - ux) * (ux * (maxCos * (ux * (maxCos * (ux + single(-1.0))))))))) * ((yi * sin(t_0)) + (xi * cos(t_0)))) + (zi * (ux * (ux * ((maxCos / ux) - maxCos)))); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := uy \cdot \left(2 \cdot \pi\right)\\
\sqrt{1 + \left(1 - ux\right) \cdot \left(ux \cdot \left(maxCos \cdot \left(ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\right)\right)\right)} \cdot \left(yi \cdot \sin t\_0 + xi \cdot \cos t\_0\right) + zi \cdot \left(ux \cdot \left(ux \cdot \left(\frac{maxCos}{ux} - maxCos\right)\right)\right)
\end{array}
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in ux around inf
unpow2N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
/-lowering-/.f3299.1%
Simplified99.1%
Final simplification99.1%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* uy (* 2.0 PI))))
(+
(*
(sqrt
(+ 1.0 (* (- 1.0 ux) (* ux (* maxCos (* ux (* maxCos (+ ux -1.0))))))))
(+ (* yi (sin t_0)) (* xi (cos t_0))))
(* zi (* ux (* (- 1.0 ux) maxCos))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = uy * (2.0f * ((float) M_PI));
return (sqrtf((1.0f + ((1.0f - ux) * (ux * (maxCos * (ux * (maxCos * (ux + -1.0f)))))))) * ((yi * sinf(t_0)) + (xi * cosf(t_0)))) + (zi * (ux * ((1.0f - ux) * maxCos)));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(uy * Float32(Float32(2.0) * Float32(pi))) return Float32(Float32(sqrt(Float32(Float32(1.0) + Float32(Float32(Float32(1.0) - ux) * Float32(ux * Float32(maxCos * Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0))))))))) * Float32(Float32(yi * sin(t_0)) + Float32(xi * cos(t_0)))) + Float32(zi * Float32(ux * Float32(Float32(Float32(1.0) - ux) * maxCos)))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = uy * (single(2.0) * single(pi)); tmp = (sqrt((single(1.0) + ((single(1.0) - ux) * (ux * (maxCos * (ux * (maxCos * (ux + single(-1.0))))))))) * ((yi * sin(t_0)) + (xi * cos(t_0)))) + (zi * (ux * ((single(1.0) - ux) * maxCos))); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := uy \cdot \left(2 \cdot \pi\right)\\
\sqrt{1 + \left(1 - ux\right) \cdot \left(ux \cdot \left(maxCos \cdot \left(ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\right)\right)\right)} \cdot \left(yi \cdot \sin t\_0 + xi \cdot \cos t\_0\right) + zi \cdot \left(ux \cdot \left(\left(1 - ux\right) \cdot maxCos\right)\right)
\end{array}
\end{array}
Initial program 99.1%
Simplified99.1%
Final simplification99.1%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* 2.0 (* PI uy))))
(fma
(* zi (* ux (- 1.0 ux)))
maxCos
(+ (* 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 * (((float) M_PI) * uy);
return fmaf((zi * (ux * (1.0f - ux))), maxCos, ((xi * cosf(t_0)) + (yi * sinf(t_0))));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(2.0) * Float32(Float32(pi) * uy)) return fma(Float32(zi * Float32(ux * Float32(Float32(1.0) - ux))), maxCos, Float32(Float32(xi * cos(t_0)) + Float32(yi * sin(t_0)))) end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 2 \cdot \left(\pi \cdot uy\right)\\
\mathsf{fma}\left(zi \cdot \left(ux \cdot \left(1 - ux\right)\right), maxCos, xi \cdot \cos t\_0 + yi \cdot \sin t\_0\right)
\end{array}
\end{array}
Initial program 99.1%
Simplified99.1%
Applied egg-rr99.1%
Taylor expanded in maxCos around 0
Simplified99.0%
Final simplification99.0%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* 2.0 (* PI uy))))
(+
(* xi (cos t_0))
(+ (* yi (sin t_0)) (* zi (* (* ux (- 1.0 ux)) maxCos))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = 2.0f * (((float) M_PI) * uy);
return (xi * cosf(t_0)) + ((yi * sinf(t_0)) + (zi * ((ux * (1.0f - ux)) * maxCos)));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(2.0) * Float32(Float32(pi) * uy)) return Float32(Float32(xi * cos(t_0)) + Float32(Float32(yi * sin(t_0)) + Float32(zi * Float32(Float32(ux * Float32(Float32(1.0) - ux)) * maxCos)))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = single(2.0) * (single(pi) * uy); tmp = (xi * cos(t_0)) + ((yi * sin(t_0)) + (zi * ((ux * (single(1.0) - ux)) * maxCos))); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 2 \cdot \left(\pi \cdot uy\right)\\
xi \cdot \cos t\_0 + \left(yi \cdot \sin t\_0 + zi \cdot \left(\left(ux \cdot \left(1 - ux\right)\right) \cdot maxCos\right)\right)
\end{array}
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.9%
Final simplification98.9%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(+
(+ (* yi (sin (* 2.0 (* PI uy)))) (* zi (* (* ux (- 1.0 ux)) maxCos)))
(+
xi
(*
uy
(*
uy
(+
(* xi (* (* PI PI) -2.0))
(* xi (* (* PI (* PI (* PI PI))) (* uy (* uy 0.6666666666666666))))))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return ((yi * sinf((2.0f * (((float) M_PI) * uy)))) + (zi * ((ux * (1.0f - ux)) * maxCos))) + (xi + (uy * (uy * ((xi * ((((float) M_PI) * ((float) M_PI)) * -2.0f)) + (xi * ((((float) M_PI) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))) * (uy * (uy * 0.6666666666666666f))))))));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(Float32(yi * sin(Float32(Float32(2.0) * Float32(Float32(pi) * uy)))) + Float32(zi * Float32(Float32(ux * Float32(Float32(1.0) - ux)) * maxCos))) + Float32(xi + Float32(uy * Float32(uy * Float32(Float32(xi * Float32(Float32(Float32(pi) * Float32(pi)) * Float32(-2.0))) + Float32(xi * Float32(Float32(Float32(pi) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))) * Float32(uy * Float32(uy * Float32(0.6666666666666666)))))))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = ((yi * sin((single(2.0) * (single(pi) * uy)))) + (zi * ((ux * (single(1.0) - ux)) * maxCos))) + (xi + (uy * (uy * ((xi * ((single(pi) * single(pi)) * single(-2.0))) + (xi * ((single(pi) * (single(pi) * (single(pi) * single(pi)))) * (uy * (uy * single(0.6666666666666666))))))))); end
\begin{array}{l}
\\
\left(yi \cdot \sin \left(2 \cdot \left(\pi \cdot uy\right)\right) + zi \cdot \left(\left(ux \cdot \left(1 - ux\right)\right) \cdot maxCos\right)\right) + \left(xi + uy \cdot \left(uy \cdot \left(xi \cdot \left(\left(\pi \cdot \pi\right) \cdot -2\right) + xi \cdot \left(\left(\pi \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right) \cdot \left(uy \cdot \left(uy \cdot 0.6666666666666666\right)\right)\right)\right)\right)\right)
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.9%
Taylor expanded in uy around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
pow-lowering-pow.f32N/A
Simplified96.5%
+-commutativeN/A
+-lowering-+.f32N/A
Applied egg-rr96.5%
Final simplification96.5%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(+
(+
xi
(+
(* yi (sin (* 2.0 (* PI uy))))
(*
uy
(*
uy
(+
(* xi (* (* PI PI) -2.0))
(* xi (* (* PI (* PI (* PI PI))) (* uy (* uy 0.6666666666666666)))))))))
(* (* (- 1.0 ux) maxCos) (* zi ux))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return (xi + ((yi * sinf((2.0f * (((float) M_PI) * uy)))) + (uy * (uy * ((xi * ((((float) M_PI) * ((float) M_PI)) * -2.0f)) + (xi * ((((float) M_PI) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))) * (uy * (uy * 0.6666666666666666f))))))))) + (((1.0f - ux) * maxCos) * (zi * ux));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(xi + Float32(Float32(yi * sin(Float32(Float32(2.0) * Float32(Float32(pi) * uy)))) + Float32(uy * Float32(uy * Float32(Float32(xi * Float32(Float32(Float32(pi) * Float32(pi)) * Float32(-2.0))) + Float32(xi * Float32(Float32(Float32(pi) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))) * Float32(uy * Float32(uy * Float32(0.6666666666666666)))))))))) + Float32(Float32(Float32(Float32(1.0) - ux) * maxCos) * Float32(zi * ux))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = (xi + ((yi * sin((single(2.0) * (single(pi) * uy)))) + (uy * (uy * ((xi * ((single(pi) * single(pi)) * single(-2.0))) + (xi * ((single(pi) * (single(pi) * (single(pi) * single(pi)))) * (uy * (uy * single(0.6666666666666666)))))))))) + (((single(1.0) - ux) * maxCos) * (zi * ux)); end
\begin{array}{l}
\\
\left(xi + \left(yi \cdot \sin \left(2 \cdot \left(\pi \cdot uy\right)\right) + uy \cdot \left(uy \cdot \left(xi \cdot \left(\left(\pi \cdot \pi\right) \cdot -2\right) + xi \cdot \left(\left(\pi \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right) \cdot \left(uy \cdot \left(uy \cdot 0.6666666666666666\right)\right)\right)\right)\right)\right)\right) + \left(\left(1 - ux\right) \cdot maxCos\right) \cdot \left(zi \cdot ux\right)
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.9%
Taylor expanded in uy around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
pow-lowering-pow.f32N/A
Simplified96.5%
Applied egg-rr96.5%
Final simplification96.5%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(if (<= uy 0.011420000344514847)
(+
(+ xi (* (* ux maxCos) (* zi (- 1.0 ux))))
(*
uy
(+
(* 2.0 (* PI yi))
(*
uy
(+
(* -2.0 (* xi (* PI PI)))
(* -1.3333333333333333 (* (* PI (* PI PI)) (* uy yi))))))))
(+
(+
(* yi (sin (* 2.0 (* PI uy))))
(* xi (+ 1.0 (* uy (* -2.0 (* uy (* PI PI)))))))
(* maxCos (* zi ux)))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float tmp;
if (uy <= 0.011420000344514847f) {
tmp = (xi + ((ux * maxCos) * (zi * (1.0f - ux)))) + (uy * ((2.0f * (((float) M_PI) * yi)) + (uy * ((-2.0f * (xi * (((float) M_PI) * ((float) M_PI)))) + (-1.3333333333333333f * ((((float) M_PI) * (((float) M_PI) * ((float) M_PI))) * (uy * yi)))))));
} else {
tmp = ((yi * sinf((2.0f * (((float) M_PI) * uy)))) + (xi * (1.0f + (uy * (-2.0f * (uy * (((float) M_PI) * ((float) M_PI)))))))) + (maxCos * (zi * ux));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) tmp = Float32(0.0) if (uy <= Float32(0.011420000344514847)) tmp = Float32(Float32(xi + Float32(Float32(ux * maxCos) * Float32(zi * Float32(Float32(1.0) - ux)))) + Float32(uy * Float32(Float32(Float32(2.0) * Float32(Float32(pi) * yi)) + Float32(uy * Float32(Float32(Float32(-2.0) * Float32(xi * Float32(Float32(pi) * Float32(pi)))) + Float32(Float32(-1.3333333333333333) * Float32(Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))) * Float32(uy * yi)))))))); else tmp = Float32(Float32(Float32(yi * sin(Float32(Float32(2.0) * Float32(Float32(pi) * uy)))) + Float32(xi * Float32(Float32(1.0) + Float32(uy * Float32(Float32(-2.0) * Float32(uy * Float32(Float32(pi) * Float32(pi)))))))) + Float32(maxCos * Float32(zi * ux))); end return tmp end
function tmp_2 = code(xi, yi, zi, ux, uy, maxCos) tmp = single(0.0); if (uy <= single(0.011420000344514847)) tmp = (xi + ((ux * maxCos) * (zi * (single(1.0) - ux)))) + (uy * ((single(2.0) * (single(pi) * yi)) + (uy * ((single(-2.0) * (xi * (single(pi) * single(pi)))) + (single(-1.3333333333333333) * ((single(pi) * (single(pi) * single(pi))) * (uy * yi))))))); else tmp = ((yi * sin((single(2.0) * (single(pi) * uy)))) + (xi * (single(1.0) + (uy * (single(-2.0) * (uy * (single(pi) * single(pi)))))))) + (maxCos * (zi * ux)); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;uy \leq 0.011420000344514847:\\
\;\;\;\;\left(xi + \left(ux \cdot maxCos\right) \cdot \left(zi \cdot \left(1 - ux\right)\right)\right) + uy \cdot \left(2 \cdot \left(\pi \cdot yi\right) + uy \cdot \left(-2 \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right) + -1.3333333333333333 \cdot \left(\left(\pi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(uy \cdot yi\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(yi \cdot \sin \left(2 \cdot \left(\pi \cdot uy\right)\right) + xi \cdot \left(1 + uy \cdot \left(-2 \cdot \left(uy \cdot \left(\pi \cdot \pi\right)\right)\right)\right)\right) + maxCos \cdot \left(zi \cdot ux\right)\\
\end{array}
\end{array}
if uy < 0.0114200003Initial program 99.3%
Simplified99.3%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified99.3%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified99.1%
if 0.0114200003 < uy Initial program 98.2%
Simplified98.2%
Taylor expanded in ux around 0
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
+-lowering-+.f32N/A
Simplified95.6%
Taylor expanded in uy around 0
*-lft-identityN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
distribute-rgt-inN/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
Simplified76.9%
Final simplification94.3%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (+ (+ (* yi (sin (* 2.0 (* PI uy)))) (* zi (* (* ux (- 1.0 ux)) maxCos))) (+ xi (* (* -2.0 (* uy uy)) (* xi (* PI PI))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return ((yi * sinf((2.0f * (((float) M_PI) * uy)))) + (zi * ((ux * (1.0f - ux)) * maxCos))) + (xi + ((-2.0f * (uy * uy)) * (xi * (((float) M_PI) * ((float) M_PI)))));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(Float32(yi * sin(Float32(Float32(2.0) * Float32(Float32(pi) * uy)))) + Float32(zi * Float32(Float32(ux * Float32(Float32(1.0) - ux)) * maxCos))) + Float32(xi + Float32(Float32(Float32(-2.0) * Float32(uy * uy)) * Float32(xi * Float32(Float32(pi) * Float32(pi)))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = ((yi * sin((single(2.0) * (single(pi) * uy)))) + (zi * ((ux * (single(1.0) - ux)) * maxCos))) + (xi + ((single(-2.0) * (uy * uy)) * (xi * (single(pi) * single(pi))))); end
\begin{array}{l}
\\
\left(yi \cdot \sin \left(2 \cdot \left(\pi \cdot uy\right)\right) + zi \cdot \left(\left(ux \cdot \left(1 - ux\right)\right) \cdot maxCos\right)\right) + \left(xi + \left(-2 \cdot \left(uy \cdot uy\right)\right) \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right)\right)
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.9%
Taylor expanded in uy around 0
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3294.8%
Simplified94.8%
Final simplification94.8%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(if (<= uy 0.05000000074505806)
(+
(+ xi (* (* ux maxCos) (* zi (- 1.0 ux))))
(*
uy
(+
(* 2.0 (* PI yi))
(*
uy
(+
(* -2.0 (* xi (* PI PI)))
(* -1.3333333333333333 (* (* PI (* PI PI)) (* uy yi))))))))
(+
xi
(+ (* yi (sin (* 2.0 (* PI uy)))) (* zi (* (* ux (- 1.0 ux)) maxCos))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float tmp;
if (uy <= 0.05000000074505806f) {
tmp = (xi + ((ux * maxCos) * (zi * (1.0f - ux)))) + (uy * ((2.0f * (((float) M_PI) * yi)) + (uy * ((-2.0f * (xi * (((float) M_PI) * ((float) M_PI)))) + (-1.3333333333333333f * ((((float) M_PI) * (((float) M_PI) * ((float) M_PI))) * (uy * yi)))))));
} else {
tmp = xi + ((yi * sinf((2.0f * (((float) M_PI) * uy)))) + (zi * ((ux * (1.0f - ux)) * maxCos)));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) tmp = Float32(0.0) if (uy <= Float32(0.05000000074505806)) tmp = Float32(Float32(xi + Float32(Float32(ux * maxCos) * Float32(zi * Float32(Float32(1.0) - ux)))) + Float32(uy * Float32(Float32(Float32(2.0) * Float32(Float32(pi) * yi)) + Float32(uy * Float32(Float32(Float32(-2.0) * Float32(xi * Float32(Float32(pi) * Float32(pi)))) + Float32(Float32(-1.3333333333333333) * Float32(Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))) * Float32(uy * yi)))))))); else tmp = Float32(xi + Float32(Float32(yi * sin(Float32(Float32(2.0) * Float32(Float32(pi) * uy)))) + Float32(zi * Float32(Float32(ux * Float32(Float32(1.0) - ux)) * maxCos)))); end return tmp end
function tmp_2 = code(xi, yi, zi, ux, uy, maxCos) tmp = single(0.0); if (uy <= single(0.05000000074505806)) tmp = (xi + ((ux * maxCos) * (zi * (single(1.0) - ux)))) + (uy * ((single(2.0) * (single(pi) * yi)) + (uy * ((single(-2.0) * (xi * (single(pi) * single(pi)))) + (single(-1.3333333333333333) * ((single(pi) * (single(pi) * single(pi))) * (uy * yi))))))); else tmp = xi + ((yi * sin((single(2.0) * (single(pi) * uy)))) + (zi * ((ux * (single(1.0) - ux)) * maxCos))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;uy \leq 0.05000000074505806:\\
\;\;\;\;\left(xi + \left(ux \cdot maxCos\right) \cdot \left(zi \cdot \left(1 - ux\right)\right)\right) + uy \cdot \left(2 \cdot \left(\pi \cdot yi\right) + uy \cdot \left(-2 \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right) + -1.3333333333333333 \cdot \left(\left(\pi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(uy \cdot yi\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;xi + \left(yi \cdot \sin \left(2 \cdot \left(\pi \cdot uy\right)\right) + zi \cdot \left(\left(ux \cdot \left(1 - ux\right)\right) \cdot maxCos\right)\right)\\
\end{array}
\end{array}
if uy < 0.0500000007Initial program 99.3%
Simplified99.3%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified99.1%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified97.9%
if 0.0500000007 < uy Initial program 98.0%
Simplified98.0%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified97.9%
Taylor expanded in uy around 0
Simplified71.5%
Final simplification93.9%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(if (<= uy 0.05000000074505806)
(+
(+ xi (* (* ux maxCos) (* zi (- 1.0 ux))))
(*
uy
(+
(* 2.0 (* PI yi))
(*
uy
(+
(* -2.0 (* xi (* PI PI)))
(* -1.3333333333333333 (* (* PI (* PI PI)) (* uy yi))))))))
(+ (* maxCos (* zi ux)) (+ xi (* yi (sin (* 2.0 (* PI uy))))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float tmp;
if (uy <= 0.05000000074505806f) {
tmp = (xi + ((ux * maxCos) * (zi * (1.0f - ux)))) + (uy * ((2.0f * (((float) M_PI) * yi)) + (uy * ((-2.0f * (xi * (((float) M_PI) * ((float) M_PI)))) + (-1.3333333333333333f * ((((float) M_PI) * (((float) M_PI) * ((float) M_PI))) * (uy * yi)))))));
} else {
tmp = (maxCos * (zi * ux)) + (xi + (yi * sinf((2.0f * (((float) M_PI) * uy)))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) tmp = Float32(0.0) if (uy <= Float32(0.05000000074505806)) tmp = Float32(Float32(xi + Float32(Float32(ux * maxCos) * Float32(zi * Float32(Float32(1.0) - ux)))) + Float32(uy * Float32(Float32(Float32(2.0) * Float32(Float32(pi) * yi)) + Float32(uy * Float32(Float32(Float32(-2.0) * Float32(xi * Float32(Float32(pi) * Float32(pi)))) + Float32(Float32(-1.3333333333333333) * Float32(Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))) * Float32(uy * yi)))))))); else tmp = Float32(Float32(maxCos * Float32(zi * ux)) + Float32(xi + Float32(yi * sin(Float32(Float32(2.0) * Float32(Float32(pi) * uy)))))); end return tmp end
function tmp_2 = code(xi, yi, zi, ux, uy, maxCos) tmp = single(0.0); if (uy <= single(0.05000000074505806)) tmp = (xi + ((ux * maxCos) * (zi * (single(1.0) - ux)))) + (uy * ((single(2.0) * (single(pi) * yi)) + (uy * ((single(-2.0) * (xi * (single(pi) * single(pi)))) + (single(-1.3333333333333333) * ((single(pi) * (single(pi) * single(pi))) * (uy * yi))))))); else tmp = (maxCos * (zi * ux)) + (xi + (yi * sin((single(2.0) * (single(pi) * uy))))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;uy \leq 0.05000000074505806:\\
\;\;\;\;\left(xi + \left(ux \cdot maxCos\right) \cdot \left(zi \cdot \left(1 - ux\right)\right)\right) + uy \cdot \left(2 \cdot \left(\pi \cdot yi\right) + uy \cdot \left(-2 \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right) + -1.3333333333333333 \cdot \left(\left(\pi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(uy \cdot yi\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;maxCos \cdot \left(zi \cdot ux\right) + \left(xi + yi \cdot \sin \left(2 \cdot \left(\pi \cdot uy\right)\right)\right)\\
\end{array}
\end{array}
if uy < 0.0500000007Initial program 99.3%
Simplified99.3%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified99.1%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified97.9%
if 0.0500000007 < uy Initial program 98.0%
Simplified98.0%
Taylor expanded in ux around 0
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
+-lowering-+.f32N/A
Simplified95.2%
Taylor expanded in uy around 0
Simplified69.9%
Final simplification93.6%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(if (<= uy 0.08749999850988388)
(+
(+ xi (* (* ux maxCos) (* zi (- 1.0 ux))))
(*
uy
(+
(* 2.0 (* PI yi))
(*
uy
(+
(* -2.0 (* xi (* PI PI)))
(* -1.3333333333333333 (* (* PI (* PI PI)) (* uy yi))))))))
(+ (* maxCos (* zi ux)) (* yi (sin (* PI (* 2.0 uy)))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float tmp;
if (uy <= 0.08749999850988388f) {
tmp = (xi + ((ux * maxCos) * (zi * (1.0f - ux)))) + (uy * ((2.0f * (((float) M_PI) * yi)) + (uy * ((-2.0f * (xi * (((float) M_PI) * ((float) M_PI)))) + (-1.3333333333333333f * ((((float) M_PI) * (((float) M_PI) * ((float) M_PI))) * (uy * yi)))))));
} else {
tmp = (maxCos * (zi * ux)) + (yi * sinf((((float) M_PI) * (2.0f * uy))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) tmp = Float32(0.0) if (uy <= Float32(0.08749999850988388)) tmp = Float32(Float32(xi + Float32(Float32(ux * maxCos) * Float32(zi * Float32(Float32(1.0) - ux)))) + Float32(uy * Float32(Float32(Float32(2.0) * Float32(Float32(pi) * yi)) + Float32(uy * Float32(Float32(Float32(-2.0) * Float32(xi * Float32(Float32(pi) * Float32(pi)))) + Float32(Float32(-1.3333333333333333) * Float32(Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))) * Float32(uy * yi)))))))); else tmp = Float32(Float32(maxCos * Float32(zi * ux)) + Float32(yi * sin(Float32(Float32(pi) * Float32(Float32(2.0) * uy))))); end return tmp end
function tmp_2 = code(xi, yi, zi, ux, uy, maxCos) tmp = single(0.0); if (uy <= single(0.08749999850988388)) tmp = (xi + ((ux * maxCos) * (zi * (single(1.0) - ux)))) + (uy * ((single(2.0) * (single(pi) * yi)) + (uy * ((single(-2.0) * (xi * (single(pi) * single(pi)))) + (single(-1.3333333333333333) * ((single(pi) * (single(pi) * single(pi))) * (uy * yi))))))); else tmp = (maxCos * (zi * ux)) + (yi * sin((single(pi) * (single(2.0) * uy)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;uy \leq 0.08749999850988388:\\
\;\;\;\;\left(xi + \left(ux \cdot maxCos\right) \cdot \left(zi \cdot \left(1 - ux\right)\right)\right) + uy \cdot \left(2 \cdot \left(\pi \cdot yi\right) + uy \cdot \left(-2 \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right) + -1.3333333333333333 \cdot \left(\left(\pi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(uy \cdot yi\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;maxCos \cdot \left(zi \cdot ux\right) + yi \cdot \sin \left(\pi \cdot \left(2 \cdot uy\right)\right)\\
\end{array}
\end{array}
if uy < 0.0874999985Initial program 99.2%
Simplified99.3%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified99.1%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified96.9%
if 0.0874999985 < uy Initial program 97.9%
Simplified97.9%
Taylor expanded in ux around 0
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
+-lowering-+.f32N/A
Simplified94.6%
Taylor expanded in xi around 0
+-commutativeN/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
sin-lowering-sin.f32N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f3265.2%
Simplified65.2%
Final simplification92.8%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(if (<= uy 0.0006600000197067857)
(+
(+ xi (* (* ux maxCos) (* zi (- 1.0 ux))))
(*
yi
(+ (* PI (* 2.0 uy)) (/ (* (* PI PI) (* -2.0 (* xi (* uy uy)))) yi))))
(+
(* maxCos (* zi ux))
(+
xi
(*
uy
(+
(* 2.0 (* PI yi))
(*
uy
(+
(* (* PI PI) (* xi -2.0))
(* (* PI (* yi (* PI PI))) (* uy -1.3333333333333333))))))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float tmp;
if (uy <= 0.0006600000197067857f) {
tmp = (xi + ((ux * maxCos) * (zi * (1.0f - ux)))) + (yi * ((((float) M_PI) * (2.0f * uy)) + (((((float) M_PI) * ((float) M_PI)) * (-2.0f * (xi * (uy * uy)))) / yi)));
} else {
tmp = (maxCos * (zi * ux)) + (xi + (uy * ((2.0f * (((float) M_PI) * yi)) + (uy * (((((float) M_PI) * ((float) M_PI)) * (xi * -2.0f)) + ((((float) M_PI) * (yi * (((float) M_PI) * ((float) M_PI)))) * (uy * -1.3333333333333333f)))))));
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) tmp = Float32(0.0) if (uy <= Float32(0.0006600000197067857)) tmp = Float32(Float32(xi + Float32(Float32(ux * maxCos) * Float32(zi * Float32(Float32(1.0) - ux)))) + Float32(yi * Float32(Float32(Float32(pi) * Float32(Float32(2.0) * uy)) + Float32(Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(-2.0) * Float32(xi * Float32(uy * uy)))) / yi)))); else tmp = Float32(Float32(maxCos * Float32(zi * ux)) + Float32(xi + Float32(uy * Float32(Float32(Float32(2.0) * Float32(Float32(pi) * yi)) + Float32(uy * Float32(Float32(Float32(Float32(pi) * Float32(pi)) * Float32(xi * Float32(-2.0))) + Float32(Float32(Float32(pi) * Float32(yi * Float32(Float32(pi) * Float32(pi)))) * Float32(uy * Float32(-1.3333333333333333))))))))); end return tmp end
function tmp_2 = code(xi, yi, zi, ux, uy, maxCos) tmp = single(0.0); if (uy <= single(0.0006600000197067857)) tmp = (xi + ((ux * maxCos) * (zi * (single(1.0) - ux)))) + (yi * ((single(pi) * (single(2.0) * uy)) + (((single(pi) * single(pi)) * (single(-2.0) * (xi * (uy * uy)))) / yi))); else tmp = (maxCos * (zi * ux)) + (xi + (uy * ((single(2.0) * (single(pi) * yi)) + (uy * (((single(pi) * single(pi)) * (xi * single(-2.0))) + ((single(pi) * (yi * (single(pi) * single(pi)))) * (uy * single(-1.3333333333333333)))))))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;uy \leq 0.0006600000197067857:\\
\;\;\;\;\left(xi + \left(ux \cdot maxCos\right) \cdot \left(zi \cdot \left(1 - ux\right)\right)\right) + yi \cdot \left(\pi \cdot \left(2 \cdot uy\right) + \frac{\left(\pi \cdot \pi\right) \cdot \left(-2 \cdot \left(xi \cdot \left(uy \cdot uy\right)\right)\right)}{yi}\right)\\
\mathbf{else}:\\
\;\;\;\;maxCos \cdot \left(zi \cdot ux\right) + \left(xi + uy \cdot \left(2 \cdot \left(\pi \cdot yi\right) + uy \cdot \left(\left(\pi \cdot \pi\right) \cdot \left(xi \cdot -2\right) + \left(\pi \cdot \left(yi \cdot \left(\pi \cdot \pi\right)\right)\right) \cdot \left(uy \cdot -1.3333333333333333\right)\right)\right)\right)\\
\end{array}
\end{array}
if uy < 6.6000002e-4Initial program 99.4%
Simplified99.5%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified99.4%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified98.9%
Taylor expanded in yi around inf
*-lowering-*.f32N/A
+-commutativeN/A
+-lowering-+.f32N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
associate-*l/N/A
associate-*r*N/A
*-commutativeN/A
/-lowering-/.f32N/A
Simplified99.0%
if 6.6000002e-4 < uy Initial program 98.2%
Simplified98.2%
Taylor expanded in ux around 0
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
+-lowering-+.f32N/A
Simplified96.2%
Taylor expanded in uy around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified66.4%
Final simplification89.2%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(+
(+ xi (* (* ux maxCos) (* zi (- 1.0 ux))))
(*
uy
(+
(* 2.0 (* PI yi))
(*
uy
(+
(* -2.0 (* xi (* PI PI)))
(* -1.3333333333333333 (* (* PI (* PI PI)) (* uy yi)))))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return (xi + ((ux * maxCos) * (zi * (1.0f - ux)))) + (uy * ((2.0f * (((float) M_PI) * yi)) + (uy * ((-2.0f * (xi * (((float) M_PI) * ((float) M_PI)))) + (-1.3333333333333333f * ((((float) M_PI) * (((float) M_PI) * ((float) M_PI))) * (uy * yi)))))));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(xi + Float32(Float32(ux * maxCos) * Float32(zi * Float32(Float32(1.0) - ux)))) + Float32(uy * Float32(Float32(Float32(2.0) * Float32(Float32(pi) * yi)) + Float32(uy * Float32(Float32(Float32(-2.0) * Float32(xi * Float32(Float32(pi) * Float32(pi)))) + Float32(Float32(-1.3333333333333333) * Float32(Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))) * Float32(uy * yi)))))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = (xi + ((ux * maxCos) * (zi * (single(1.0) - ux)))) + (uy * ((single(2.0) * (single(pi) * yi)) + (uy * ((single(-2.0) * (xi * (single(pi) * single(pi)))) + (single(-1.3333333333333333) * ((single(pi) * (single(pi) * single(pi))) * (uy * yi))))))); end
\begin{array}{l}
\\
\left(xi + \left(ux \cdot maxCos\right) \cdot \left(zi \cdot \left(1 - ux\right)\right)\right) + uy \cdot \left(2 \cdot \left(\pi \cdot yi\right) + uy \cdot \left(-2 \cdot \left(xi \cdot \left(\pi \cdot \pi\right)\right) + -1.3333333333333333 \cdot \left(\left(\pi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(uy \cdot yi\right)\right)\right)\right)
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.9%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified89.6%
Final simplification89.6%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (+ (+ xi (* (* ux maxCos) (* zi (- 1.0 ux)))) (* yi (+ (* PI (* 2.0 uy)) (/ (* (* PI PI) (* -2.0 (* xi (* uy uy)))) yi)))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return (xi + ((ux * maxCos) * (zi * (1.0f - ux)))) + (yi * ((((float) M_PI) * (2.0f * uy)) + (((((float) M_PI) * ((float) M_PI)) * (-2.0f * (xi * (uy * uy)))) / yi)));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(xi + Float32(Float32(ux * maxCos) * Float32(zi * Float32(Float32(1.0) - ux)))) + Float32(yi * Float32(Float32(Float32(pi) * Float32(Float32(2.0) * uy)) + Float32(Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(-2.0) * Float32(xi * Float32(uy * uy)))) / yi)))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = (xi + ((ux * maxCos) * (zi * (single(1.0) - ux)))) + (yi * ((single(pi) * (single(2.0) * uy)) + (((single(pi) * single(pi)) * (single(-2.0) * (xi * (uy * uy)))) / yi))); end
\begin{array}{l}
\\
\left(xi + \left(ux \cdot maxCos\right) \cdot \left(zi \cdot \left(1 - ux\right)\right)\right) + yi \cdot \left(\pi \cdot \left(2 \cdot uy\right) + \frac{\left(\pi \cdot \pi\right) \cdot \left(-2 \cdot \left(xi \cdot \left(uy \cdot uy\right)\right)\right)}{yi}\right)
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.9%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified85.9%
Taylor expanded in yi around inf
*-lowering-*.f32N/A
+-commutativeN/A
+-lowering-+.f32N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
associate-*l/N/A
associate-*r*N/A
*-commutativeN/A
/-lowering-/.f32N/A
Simplified86.0%
Final simplification86.0%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (+ (+ xi (* ux (- (* zi maxCos) (* maxCos (* zi ux))))) (* uy (+ (* 2.0 (* PI yi)) (* (* xi (* PI PI)) (* uy -2.0))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return (xi + (ux * ((zi * maxCos) - (maxCos * (zi * ux))))) + (uy * ((2.0f * (((float) M_PI) * yi)) + ((xi * (((float) M_PI) * ((float) M_PI))) * (uy * -2.0f))));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(xi + Float32(ux * Float32(Float32(zi * maxCos) - Float32(maxCos * Float32(zi * ux))))) + Float32(uy * Float32(Float32(Float32(2.0) * Float32(Float32(pi) * yi)) + Float32(Float32(xi * Float32(Float32(pi) * Float32(pi))) * Float32(uy * Float32(-2.0)))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = (xi + (ux * ((zi * maxCos) - (maxCos * (zi * ux))))) + (uy * ((single(2.0) * (single(pi) * yi)) + ((xi * (single(pi) * single(pi))) * (uy * single(-2.0))))); end
\begin{array}{l}
\\
\left(xi + ux \cdot \left(zi \cdot maxCos - maxCos \cdot \left(zi \cdot ux\right)\right)\right) + uy \cdot \left(2 \cdot \left(\pi \cdot yi\right) + \left(xi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(uy \cdot -2\right)\right)
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.9%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified85.9%
Taylor expanded in ux around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f3285.9%
Simplified85.9%
Final simplification85.9%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (+ (+ xi (* (* ux maxCos) (* zi (- 1.0 ux)))) (* uy (+ (* 2.0 (* PI yi)) (* (* xi (* PI PI)) (* uy -2.0))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return (xi + ((ux * maxCos) * (zi * (1.0f - ux)))) + (uy * ((2.0f * (((float) M_PI) * yi)) + ((xi * (((float) M_PI) * ((float) M_PI))) * (uy * -2.0f))));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(xi + Float32(Float32(ux * maxCos) * Float32(zi * Float32(Float32(1.0) - ux)))) + Float32(uy * Float32(Float32(Float32(2.0) * Float32(Float32(pi) * yi)) + Float32(Float32(xi * Float32(Float32(pi) * Float32(pi))) * Float32(uy * Float32(-2.0)))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = (xi + ((ux * maxCos) * (zi * (single(1.0) - ux)))) + (uy * ((single(2.0) * (single(pi) * yi)) + ((xi * (single(pi) * single(pi))) * (uy * single(-2.0))))); end
\begin{array}{l}
\\
\left(xi + \left(ux \cdot maxCos\right) \cdot \left(zi \cdot \left(1 - ux\right)\right)\right) + uy \cdot \left(2 \cdot \left(\pi \cdot yi\right) + \left(xi \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(uy \cdot -2\right)\right)
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.9%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified85.9%
Final simplification85.9%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* (* ux ux) (/ (* 2.0 (* uy (* PI yi))) (* ux ux)))))
(if (<= yi -3.99999992980668e-13)
t_0
(if (<= yi 1.000000013351432e-10)
(+ xi (* (- 1.0 ux) (* maxCos (* zi ux))))
t_0))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = (ux * ux) * ((2.0f * (uy * (((float) M_PI) * yi))) / (ux * ux));
float tmp;
if (yi <= -3.99999992980668e-13f) {
tmp = t_0;
} else if (yi <= 1.000000013351432e-10f) {
tmp = xi + ((1.0f - ux) * (maxCos * (zi * ux)));
} else {
tmp = t_0;
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(ux * ux) * Float32(Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * yi))) / Float32(ux * ux))) tmp = Float32(0.0) if (yi <= Float32(-3.99999992980668e-13)) tmp = t_0; elseif (yi <= Float32(1.000000013351432e-10)) tmp = Float32(xi + Float32(Float32(Float32(1.0) - ux) * Float32(maxCos * Float32(zi * ux)))); else tmp = t_0; end return tmp end
function tmp_2 = code(xi, yi, zi, ux, uy, maxCos) t_0 = (ux * ux) * ((single(2.0) * (uy * (single(pi) * yi))) / (ux * ux)); tmp = single(0.0); if (yi <= single(-3.99999992980668e-13)) tmp = t_0; elseif (yi <= single(1.000000013351432e-10)) tmp = xi + ((single(1.0) - ux) * (maxCos * (zi * ux))); else tmp = t_0; end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(ux \cdot ux\right) \cdot \frac{2 \cdot \left(uy \cdot \left(\pi \cdot yi\right)\right)}{ux \cdot ux}\\
\mathbf{if}\;yi \leq -3.99999992980668 \cdot 10^{-13}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;yi \leq 1.000000013351432 \cdot 10^{-10}:\\
\;\;\;\;xi + \left(1 - ux\right) \cdot \left(maxCos \cdot \left(zi \cdot ux\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if yi < -3.99999993e-13 or 1.00000001e-10 < yi Initial program 98.8%
Simplified98.8%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.5%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified82.5%
Taylor expanded in ux around -inf
Simplified82.1%
Taylor expanded in yi around inf
associate-*r/N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
/-lowering-/.f32N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
unpow2N/A
*-lowering-*.f3256.4%
Simplified56.4%
if -3.99999993e-13 < yi < 1.00000001e-10Initial program 99.2%
Simplified99.2%
Taylor expanded in ux around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3299.0%
Simplified99.0%
Taylor expanded in uy around 0
+-commutativeN/A
+-lowering-+.f32N/A
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
--lowering--.f3268.2%
Simplified68.2%
Final simplification64.5%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (+ (* maxCos (* zi ux)) (+ xi (* uy (+ (* 2.0 (* PI yi)) (* xi (* -2.0 (* uy (* PI PI)))))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return (maxCos * (zi * ux)) + (xi + (uy * ((2.0f * (((float) M_PI) * yi)) + (xi * (-2.0f * (uy * (((float) M_PI) * ((float) M_PI))))))));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(maxCos * Float32(zi * ux)) + Float32(xi + Float32(uy * Float32(Float32(Float32(2.0) * Float32(Float32(pi) * yi)) + Float32(xi * Float32(Float32(-2.0) * Float32(uy * Float32(Float32(pi) * Float32(pi))))))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = (maxCos * (zi * ux)) + (xi + (uy * ((single(2.0) * (single(pi) * yi)) + (xi * (single(-2.0) * (uy * (single(pi) * single(pi)))))))); end
\begin{array}{l}
\\
maxCos \cdot \left(zi \cdot ux\right) + \left(xi + uy \cdot \left(2 \cdot \left(\pi \cdot yi\right) + xi \cdot \left(-2 \cdot \left(uy \cdot \left(\pi \cdot \pi\right)\right)\right)\right)\right)
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in ux around 0
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
+-lowering-+.f32N/A
Simplified95.2%
Taylor expanded in uy around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
Simplified82.6%
Final simplification82.6%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* 2.0 (* uy (* PI yi)))))
(if (<= yi -3.99999992980668e-13)
t_0
(if (<= yi 1.000000013351432e-10)
(+ xi (* (- 1.0 ux) (* maxCos (* zi ux))))
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) * yi));
float tmp;
if (yi <= -3.99999992980668e-13f) {
tmp = t_0;
} else if (yi <= 1.000000013351432e-10f) {
tmp = xi + ((1.0f - ux) * (maxCos * (zi * ux)));
} else {
tmp = t_0;
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * yi))) tmp = Float32(0.0) if (yi <= Float32(-3.99999992980668e-13)) tmp = t_0; elseif (yi <= Float32(1.000000013351432e-10)) tmp = Float32(xi + Float32(Float32(Float32(1.0) - ux) * Float32(maxCos * Float32(zi * ux)))); else tmp = t_0; end return tmp end
function tmp_2 = code(xi, yi, zi, ux, uy, maxCos) t_0 = single(2.0) * (uy * (single(pi) * yi)); tmp = single(0.0); if (yi <= single(-3.99999992980668e-13)) tmp = t_0; elseif (yi <= single(1.000000013351432e-10)) tmp = xi + ((single(1.0) - ux) * (maxCos * (zi * ux))); else tmp = t_0; end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 2 \cdot \left(uy \cdot \left(\pi \cdot yi\right)\right)\\
\mathbf{if}\;yi \leq -3.99999992980668 \cdot 10^{-13}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;yi \leq 1.000000013351432 \cdot 10^{-10}:\\
\;\;\;\;xi + \left(1 - ux\right) \cdot \left(maxCos \cdot \left(zi \cdot ux\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if yi < -3.99999993e-13 or 1.00000001e-10 < yi Initial program 98.8%
Simplified98.8%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.5%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified82.5%
Taylor expanded in yi around inf
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
PI-lowering-PI.f3256.4%
Simplified56.4%
if -3.99999993e-13 < yi < 1.00000001e-10Initial program 99.2%
Simplified99.2%
Taylor expanded in ux around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3299.0%
Simplified99.0%
Taylor expanded in uy around 0
+-commutativeN/A
+-lowering-+.f32N/A
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
--lowering--.f3268.2%
Simplified68.2%
Final simplification64.5%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* 2.0 (* uy (* PI yi)))))
(if (<= yi -3.99999992980668e-13)
t_0
(if (<= yi 1.000000013351432e-10)
(+ xi (* zi (* ux (* (- 1.0 ux) maxCos))))
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) * yi));
float tmp;
if (yi <= -3.99999992980668e-13f) {
tmp = t_0;
} else if (yi <= 1.000000013351432e-10f) {
tmp = xi + (zi * (ux * ((1.0f - ux) * maxCos)));
} else {
tmp = t_0;
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * yi))) tmp = Float32(0.0) if (yi <= Float32(-3.99999992980668e-13)) tmp = t_0; elseif (yi <= Float32(1.000000013351432e-10)) tmp = Float32(xi + Float32(zi * Float32(ux * Float32(Float32(Float32(1.0) - ux) * maxCos)))); else tmp = t_0; end return tmp end
function tmp_2 = code(xi, yi, zi, ux, uy, maxCos) t_0 = single(2.0) * (uy * (single(pi) * yi)); tmp = single(0.0); if (yi <= single(-3.99999992980668e-13)) tmp = t_0; elseif (yi <= single(1.000000013351432e-10)) tmp = xi + (zi * (ux * ((single(1.0) - ux) * maxCos))); else tmp = t_0; end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 2 \cdot \left(uy \cdot \left(\pi \cdot yi\right)\right)\\
\mathbf{if}\;yi \leq -3.99999992980668 \cdot 10^{-13}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;yi \leq 1.000000013351432 \cdot 10^{-10}:\\
\;\;\;\;xi + zi \cdot \left(ux \cdot \left(\left(1 - ux\right) \cdot maxCos\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if yi < -3.99999993e-13 or 1.00000001e-10 < yi Initial program 98.8%
Simplified98.8%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.5%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified82.5%
Taylor expanded in yi around inf
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
PI-lowering-PI.f3256.4%
Simplified56.4%
if -3.99999993e-13 < yi < 1.00000001e-10Initial program 99.2%
Simplified99.2%
Taylor expanded in ux around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3299.0%
Simplified99.0%
Taylor expanded in uy around 0
Simplified68.2%
Final simplification64.4%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* 2.0 (* uy (* PI yi)))))
(if (<= yi -3.99999992980668e-13)
t_0
(if (<= yi 1.000000013351432e-10)
(+ xi (* (* ux maxCos) (* zi (- 1.0 ux))))
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) * yi));
float tmp;
if (yi <= -3.99999992980668e-13f) {
tmp = t_0;
} else if (yi <= 1.000000013351432e-10f) {
tmp = xi + ((ux * maxCos) * (zi * (1.0f - ux)));
} else {
tmp = t_0;
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * yi))) tmp = Float32(0.0) if (yi <= Float32(-3.99999992980668e-13)) tmp = t_0; elseif (yi <= Float32(1.000000013351432e-10)) tmp = Float32(xi + Float32(Float32(ux * maxCos) * Float32(zi * Float32(Float32(1.0) - ux)))); else tmp = t_0; end return tmp end
function tmp_2 = code(xi, yi, zi, ux, uy, maxCos) t_0 = single(2.0) * (uy * (single(pi) * yi)); tmp = single(0.0); if (yi <= single(-3.99999992980668e-13)) tmp = t_0; elseif (yi <= single(1.000000013351432e-10)) tmp = xi + ((ux * maxCos) * (zi * (single(1.0) - ux))); else tmp = t_0; end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 2 \cdot \left(uy \cdot \left(\pi \cdot yi\right)\right)\\
\mathbf{if}\;yi \leq -3.99999992980668 \cdot 10^{-13}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;yi \leq 1.000000013351432 \cdot 10^{-10}:\\
\;\;\;\;xi + \left(ux \cdot maxCos\right) \cdot \left(zi \cdot \left(1 - ux\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if yi < -3.99999993e-13 or 1.00000001e-10 < yi Initial program 98.8%
Simplified98.8%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.5%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified82.5%
Taylor expanded in yi around inf
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
PI-lowering-PI.f3256.4%
Simplified56.4%
if -3.99999993e-13 < yi < 1.00000001e-10Initial program 99.2%
Simplified99.2%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified99.1%
Taylor expanded in uy around 0
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f3268.1%
Simplified68.1%
Final simplification64.4%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (+ (* (- 1.0 ux) (* maxCos (* zi ux))) (+ xi (* 2.0 (* uy (* PI yi))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return ((1.0f - ux) * (maxCos * (zi * ux))) + (xi + (2.0f * (uy * (((float) M_PI) * yi))));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(Float32(Float32(1.0) - ux) * Float32(maxCos * Float32(zi * ux))) + Float32(xi + Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * yi))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = ((single(1.0) - ux) * (maxCos * (zi * ux))) + (xi + (single(2.0) * (uy * (single(pi) * yi)))); end
\begin{array}{l}
\\
\left(1 - ux\right) \cdot \left(maxCos \cdot \left(zi \cdot ux\right)\right) + \left(xi + 2 \cdot \left(uy \cdot \left(\pi \cdot yi\right)\right)\right)
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in ux around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.9%
Simplified98.9%
Taylor expanded in uy around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
+-lowering-+.f32N/A
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
--lowering--.f32N/A
Simplified82.4%
Final simplification82.4%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (+ xi (+ (* (* ux maxCos) (* zi (- 1.0 ux))) (* 2.0 (* uy (* PI yi))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return xi + (((ux * maxCos) * (zi * (1.0f - ux))) + (2.0f * (uy * (((float) M_PI) * yi))));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(xi + Float32(Float32(Float32(ux * maxCos) * Float32(zi * Float32(Float32(1.0) - ux))) + Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * yi))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = xi + (((ux * maxCos) * (zi * (single(1.0) - ux))) + (single(2.0) * (uy * (single(pi) * yi)))); end
\begin{array}{l}
\\
xi + \left(\left(ux \cdot maxCos\right) \cdot \left(zi \cdot \left(1 - ux\right)\right) + 2 \cdot \left(uy \cdot \left(\pi \cdot yi\right)\right)\right)
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.9%
Taylor expanded in uy around 0
+-lowering-+.f32N/A
+-commutativeN/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
PI-lowering-PI.f3282.4%
Simplified82.4%
Final simplification82.4%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* 2.0 (* uy (* PI yi)))))
(if (<= yi -3.99999992980668e-13)
t_0
(if (<= yi 1.000000013351432e-10) (+ xi (* maxCos (* zi ux))) 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) * yi));
float tmp;
if (yi <= -3.99999992980668e-13f) {
tmp = t_0;
} else if (yi <= 1.000000013351432e-10f) {
tmp = xi + (maxCos * (zi * ux));
} else {
tmp = t_0;
}
return tmp;
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * yi))) tmp = Float32(0.0) if (yi <= Float32(-3.99999992980668e-13)) tmp = t_0; elseif (yi <= Float32(1.000000013351432e-10)) tmp = Float32(xi + Float32(maxCos * Float32(zi * ux))); else tmp = t_0; end return tmp end
function tmp_2 = code(xi, yi, zi, ux, uy, maxCos) t_0 = single(2.0) * (uy * (single(pi) * yi)); tmp = single(0.0); if (yi <= single(-3.99999992980668e-13)) tmp = t_0; elseif (yi <= single(1.000000013351432e-10)) tmp = xi + (maxCos * (zi * ux)); else tmp = t_0; end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 2 \cdot \left(uy \cdot \left(\pi \cdot yi\right)\right)\\
\mathbf{if}\;yi \leq -3.99999992980668 \cdot 10^{-13}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;yi \leq 1.000000013351432 \cdot 10^{-10}:\\
\;\;\;\;xi + maxCos \cdot \left(zi \cdot ux\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if yi < -3.99999993e-13 or 1.00000001e-10 < yi Initial program 98.8%
Simplified98.8%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.5%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified82.5%
Taylor expanded in yi around inf
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
PI-lowering-PI.f3256.4%
Simplified56.4%
if -3.99999993e-13 < yi < 1.00000001e-10Initial program 99.2%
Simplified99.2%
Taylor expanded in ux around 0
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
+-lowering-+.f32N/A
Simplified94.3%
Taylor expanded in uy around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f3265.4%
Simplified65.4%
Final simplification62.5%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (+ xi (+ (* maxCos (* zi ux)) (* 2.0 (* uy (* PI yi))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return xi + ((maxCos * (zi * ux)) + (2.0f * (uy * (((float) M_PI) * yi))));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(xi + Float32(Float32(maxCos * Float32(zi * ux)) + Float32(Float32(2.0) * Float32(uy * Float32(Float32(pi) * yi))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = xi + ((maxCos * (zi * ux)) + (single(2.0) * (uy * (single(pi) * yi)))); end
\begin{array}{l}
\\
xi + \left(maxCos \cdot \left(zi \cdot ux\right) + 2 \cdot \left(uy \cdot \left(\pi \cdot yi\right)\right)\right)
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in ux around 0
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
+-lowering-+.f32N/A
Simplified95.2%
Taylor expanded in uy around 0
+-lowering-+.f32N/A
+-commutativeN/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
PI-lowering-PI.f3279.1%
Simplified79.1%
Final simplification79.1%
(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 99.1%
Simplified99.1%
Taylor expanded in maxCos around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cos-lowering-cos.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
Simplified98.9%
Taylor expanded in uy around 0
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified85.9%
Taylor expanded in yi around inf
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
PI-lowering-PI.f3229.9%
Simplified29.9%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (* ux (* zi maxCos)))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return ux * (zi * maxCos);
}
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 = ux * (zi * maxcos)
end function
function code(xi, yi, zi, ux, uy, maxCos) return Float32(ux * Float32(zi * maxCos)) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = ux * (zi * maxCos); end
\begin{array}{l}
\\
ux \cdot \left(zi \cdot maxCos\right)
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in zi around inf
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
--lowering--.f3214.5%
Simplified14.5%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
*-lowering-*.f3212.7%
Simplified12.7%
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f3212.7%
Applied egg-rr12.7%
Final simplification12.7%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (* 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)
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(maxCos * Float32(zi * ux)) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = maxCos * (zi * ux); end
\begin{array}{l}
\\
maxCos \cdot \left(zi \cdot ux\right)
\end{array}
Initial program 99.1%
Simplified99.1%
Taylor expanded in zi around inf
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
--lowering--.f3214.5%
Simplified14.5%
Taylor expanded in ux around 0
*-lowering-*.f32N/A
*-lowering-*.f3212.7%
Simplified12.7%
Final simplification12.7%
herbie shell --seed 2024141
(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)))