
(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 22 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 (* uy (* 2.0 PI)))
(t_1
(sqrt
(+
1.0
(* (- 1.0 ux) (* maxCos (* (* ux maxCos) (* ux (+ ux -1.0)))))))))
(fma
(cos t_0)
(* t_1 xi)
(fma (* t_1 (sin t_0)) yi (* ux (* (* (- 1.0 ux) maxCos) zi))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = uy * (2.0f * ((float) M_PI));
float t_1 = sqrtf((1.0f + ((1.0f - ux) * (maxCos * ((ux * maxCos) * (ux * (ux + -1.0f)))))));
return fmaf(cosf(t_0), (t_1 * xi), fmaf((t_1 * sinf(t_0)), yi, (ux * (((1.0f - ux) * maxCos) * zi))));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(uy * Float32(Float32(2.0) * Float32(pi))) t_1 = sqrt(Float32(Float32(1.0) + Float32(Float32(Float32(1.0) - ux) * Float32(maxCos * Float32(Float32(ux * maxCos) * Float32(ux * Float32(ux + Float32(-1.0)))))))) return fma(cos(t_0), Float32(t_1 * xi), fma(Float32(t_1 * sin(t_0)), yi, Float32(ux * Float32(Float32(Float32(Float32(1.0) - ux) * maxCos) * zi)))) end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := uy \cdot \left(2 \cdot \pi\right)\\
t_1 := \sqrt{1 + \left(1 - ux\right) \cdot \left(maxCos \cdot \left(\left(ux \cdot maxCos\right) \cdot \left(ux \cdot \left(ux + -1\right)\right)\right)\right)}\\
\mathsf{fma}\left(\cos t_0, t_1 \cdot xi, \mathsf{fma}\left(t_1 \cdot \sin t_0, yi, ux \cdot \left(\left(\left(1 - ux\right) \cdot maxCos\right) \cdot zi\right)\right)\right)
\end{array}
\end{array}
Initial program 98.9%
Simplified98.9%
Final simplification98.9%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* ux (* maxCos (+ ux -1.0))))
(t_1 (sqrt (- 1.0 (* t_0 t_0))))
(t_2 (* PI (* uy 2.0))))
(+
(+ (* xi (* (cos t_2) t_1)) (* yi (* t_1 (sin t_2))))
(* zi (* ux (* (- 1.0 ux) maxCos))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ux * (maxCos * (ux + -1.0f));
float t_1 = sqrtf((1.0f - (t_0 * t_0)));
float t_2 = ((float) M_PI) * (uy * 2.0f);
return ((xi * (cosf(t_2) * t_1)) + (yi * (t_1 * sinf(t_2)))) + (zi * (ux * ((1.0f - ux) * maxCos)));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0)))) t_1 = sqrt(Float32(Float32(1.0) - Float32(t_0 * t_0))) t_2 = Float32(Float32(pi) * Float32(uy * Float32(2.0))) return Float32(Float32(Float32(xi * Float32(cos(t_2) * t_1)) + Float32(yi * Float32(t_1 * sin(t_2)))) + Float32(zi * Float32(ux * Float32(Float32(Float32(1.0) - ux) * maxCos)))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = ux * (maxCos * (ux + single(-1.0))); t_1 = sqrt((single(1.0) - (t_0 * t_0))); t_2 = single(pi) * (uy * single(2.0)); tmp = ((xi * (cos(t_2) * t_1)) + (yi * (t_1 * sin(t_2)))) + (zi * (ux * ((single(1.0) - ux) * maxCos))); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\\
t_1 := \sqrt{1 - t_0 \cdot t_0}\\
t_2 := \pi \cdot \left(uy \cdot 2\right)\\
\left(xi \cdot \left(\cos t_2 \cdot t_1\right) + yi \cdot \left(t_1 \cdot \sin t_2\right)\right) + zi \cdot \left(ux \cdot \left(\left(1 - ux\right) \cdot maxCos\right)\right)
\end{array}
\end{array}
Initial program 98.9%
Final simplification98.9%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* ux (* maxCos (+ ux -1.0))))
(t_1 (sqrt (- 1.0 (* t_0 t_0))))
(t_2 (* PI (* uy 2.0))))
(+
(+ (* xi (* (cos t_2) t_1)) (* yi (* t_1 (sin t_2))))
(* zi (* (- 1.0 ux) (* ux maxCos))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ux * (maxCos * (ux + -1.0f));
float t_1 = sqrtf((1.0f - (t_0 * t_0)));
float t_2 = ((float) M_PI) * (uy * 2.0f);
return ((xi * (cosf(t_2) * t_1)) + (yi * (t_1 * sinf(t_2)))) + (zi * ((1.0f - ux) * (ux * maxCos)));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0)))) t_1 = sqrt(Float32(Float32(1.0) - Float32(t_0 * t_0))) t_2 = Float32(Float32(pi) * Float32(uy * Float32(2.0))) return Float32(Float32(Float32(xi * Float32(cos(t_2) * t_1)) + Float32(yi * Float32(t_1 * sin(t_2)))) + Float32(zi * Float32(Float32(Float32(1.0) - ux) * Float32(ux * maxCos)))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = ux * (maxCos * (ux + single(-1.0))); t_1 = sqrt((single(1.0) - (t_0 * t_0))); t_2 = single(pi) * (uy * single(2.0)); tmp = ((xi * (cos(t_2) * t_1)) + (yi * (t_1 * sin(t_2)))) + (zi * ((single(1.0) - ux) * (ux * maxCos))); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\\
t_1 := \sqrt{1 - t_0 \cdot t_0}\\
t_2 := \pi \cdot \left(uy \cdot 2\right)\\
\left(xi \cdot \left(\cos t_2 \cdot t_1\right) + yi \cdot \left(t_1 \cdot \sin t_2\right)\right) + zi \cdot \left(\left(1 - ux\right) \cdot \left(ux \cdot maxCos\right)\right)
\end{array}
\end{array}
Initial program 98.9%
Taylor expanded in ux around 0 98.8%
associate-*r*98.8%
neg-mul-198.8%
unpow298.8%
associate-*r*98.8%
neg-mul-198.8%
*-commutative98.8%
associate-*r*98.8%
distribute-rgt-in98.9%
*-rgt-identity98.9%
distribute-lft-in98.9%
*-commutative98.9%
+-commutative98.9%
mul-1-neg98.9%
sub-neg98.9%
*-commutative98.9%
associate-*r*98.9%
*-commutative98.9%
Simplified98.9%
Final simplification98.9%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(fma
(- maxCos (* ux maxCos))
(* ux zi)
(*
(sqrt
(+ 1.0 (* ux (* ux (* maxCos (* (- 1.0 ux) (- (* ux maxCos) maxCos)))))))
(+ (* xi (cos (* PI (* uy -2.0)))) (* (sin (* uy (* 2.0 PI))) yi)))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf((maxCos - (ux * maxCos)), (ux * zi), (sqrtf((1.0f + (ux * (ux * (maxCos * ((1.0f - ux) * ((ux * maxCos) - maxCos))))))) * ((xi * cosf((((float) M_PI) * (uy * -2.0f)))) + (sinf((uy * (2.0f * ((float) M_PI)))) * yi))));
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(Float32(maxCos - Float32(ux * maxCos)), Float32(ux * zi), Float32(sqrt(Float32(Float32(1.0) + Float32(ux * Float32(ux * Float32(maxCos * Float32(Float32(Float32(1.0) - ux) * Float32(Float32(ux * maxCos) - maxCos))))))) * Float32(Float32(xi * cos(Float32(Float32(pi) * Float32(uy * Float32(-2.0))))) + Float32(sin(Float32(uy * Float32(Float32(2.0) * Float32(pi)))) * yi)))) end
\begin{array}{l}
\\
\mathsf{fma}\left(maxCos - ux \cdot maxCos, ux \cdot zi, \sqrt{1 + ux \cdot \left(ux \cdot \left(maxCos \cdot \left(\left(1 - ux\right) \cdot \left(ux \cdot maxCos - maxCos\right)\right)\right)\right)} \cdot \left(xi \cdot \cos \left(\pi \cdot \left(uy \cdot -2\right)\right) + \sin \left(uy \cdot \left(2 \cdot \pi\right)\right) \cdot yi\right)\right)
\end{array}
Initial program 98.9%
Simplified98.8%
Final simplification98.8%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* PI (* uy 2.0))) (t_1 (* ux (* maxCos (+ ux -1.0)))))
(+
(* zi (* ux (* (- 1.0 ux) maxCos)))
(+ (* xi (* (cos t_0) (sqrt (- 1.0 (* t_1 t_1))))) (* yi (sin t_0))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ((float) M_PI) * (uy * 2.0f);
float t_1 = ux * (maxCos * (ux + -1.0f));
return (zi * (ux * ((1.0f - ux) * maxCos))) + ((xi * (cosf(t_0) * sqrtf((1.0f - (t_1 * t_1))))) + (yi * sinf(t_0)));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(pi) * Float32(uy * Float32(2.0))) t_1 = Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0)))) return Float32(Float32(zi * Float32(ux * Float32(Float32(Float32(1.0) - ux) * maxCos))) + Float32(Float32(xi * Float32(cos(t_0) * sqrt(Float32(Float32(1.0) - Float32(t_1 * t_1))))) + Float32(yi * sin(t_0)))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = single(pi) * (uy * single(2.0)); t_1 = ux * (maxCos * (ux + single(-1.0))); tmp = (zi * (ux * ((single(1.0) - ux) * maxCos))) + ((xi * (cos(t_0) * sqrt((single(1.0) - (t_1 * t_1))))) + (yi * sin(t_0))); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \pi \cdot \left(uy \cdot 2\right)\\
t_1 := ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\\
zi \cdot \left(ux \cdot \left(\left(1 - ux\right) \cdot maxCos\right)\right) + \left(xi \cdot \left(\cos t_0 \cdot \sqrt{1 - t_1 \cdot t_1}\right) + yi \cdot \sin t_0\right)
\end{array}
\end{array}
Initial program 98.9%
Taylor expanded in ux around 0 98.7%
associate-*r*98.7%
Simplified98.7%
Final simplification98.7%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* PI (* uy 2.0))) (t_1 (* ux (* maxCos (+ ux -1.0)))))
(+
(* zi (* (- 1.0 ux) (* ux maxCos)))
(+ (* xi (* (cos t_0) (sqrt (- 1.0 (* t_1 t_1))))) (* yi (sin t_0))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ((float) M_PI) * (uy * 2.0f);
float t_1 = ux * (maxCos * (ux + -1.0f));
return (zi * ((1.0f - ux) * (ux * maxCos))) + ((xi * (cosf(t_0) * sqrtf((1.0f - (t_1 * t_1))))) + (yi * sinf(t_0)));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(Float32(pi) * Float32(uy * Float32(2.0))) t_1 = Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0)))) return Float32(Float32(zi * Float32(Float32(Float32(1.0) - ux) * Float32(ux * maxCos))) + Float32(Float32(xi * Float32(cos(t_0) * sqrt(Float32(Float32(1.0) - Float32(t_1 * t_1))))) + Float32(yi * sin(t_0)))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = single(pi) * (uy * single(2.0)); t_1 = ux * (maxCos * (ux + single(-1.0))); tmp = (zi * ((single(1.0) - ux) * (ux * maxCos))) + ((xi * (cos(t_0) * sqrt((single(1.0) - (t_1 * t_1))))) + (yi * sin(t_0))); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \pi \cdot \left(uy \cdot 2\right)\\
t_1 := ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\\
zi \cdot \left(\left(1 - ux\right) \cdot \left(ux \cdot maxCos\right)\right) + \left(xi \cdot \left(\cos t_0 \cdot \sqrt{1 - t_1 \cdot t_1}\right) + yi \cdot \sin t_0\right)
\end{array}
\end{array}
Initial program 98.9%
Taylor expanded in ux around 0 98.8%
associate-*r*98.8%
neg-mul-198.8%
unpow298.8%
associate-*r*98.8%
neg-mul-198.8%
*-commutative98.8%
associate-*r*98.8%
distribute-rgt-in98.9%
*-rgt-identity98.9%
distribute-lft-in98.9%
*-commutative98.9%
+-commutative98.9%
mul-1-neg98.9%
sub-neg98.9%
*-commutative98.9%
associate-*r*98.9%
*-commutative98.9%
Simplified98.9%
Taylor expanded in ux around 0 98.8%
associate-*r*98.7%
Simplified98.8%
Final simplification98.8%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* ux (* maxCos (+ ux -1.0)))) (t_1 (* PI (* uy 2.0))))
(+
(+ (* xi (* (cos t_1) (sqrt (- 1.0 (* t_0 t_0))))) (* yi (sin t_1)))
(* (* ux maxCos) zi))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ux * (maxCos * (ux + -1.0f));
float t_1 = ((float) M_PI) * (uy * 2.0f);
return ((xi * (cosf(t_1) * sqrtf((1.0f - (t_0 * t_0))))) + (yi * sinf(t_1))) + ((ux * maxCos) * zi);
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0)))) t_1 = Float32(Float32(pi) * Float32(uy * Float32(2.0))) return Float32(Float32(Float32(xi * Float32(cos(t_1) * sqrt(Float32(Float32(1.0) - Float32(t_0 * t_0))))) + Float32(yi * sin(t_1))) + Float32(Float32(ux * maxCos) * zi)) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = ux * (maxCos * (ux + single(-1.0))); t_1 = single(pi) * (uy * single(2.0)); tmp = ((xi * (cos(t_1) * sqrt((single(1.0) - (t_0 * t_0))))) + (yi * sin(t_1))) + ((ux * maxCos) * zi); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\\
t_1 := \pi \cdot \left(uy \cdot 2\right)\\
\left(xi \cdot \left(\cos t_1 \cdot \sqrt{1 - t_0 \cdot t_0}\right) + yi \cdot \sin t_1\right) + \left(ux \cdot maxCos\right) \cdot zi
\end{array}
\end{array}
Initial program 98.9%
Taylor expanded in ux around 0 98.8%
associate-*r*98.8%
neg-mul-198.8%
unpow298.8%
associate-*r*98.8%
neg-mul-198.8%
*-commutative98.8%
associate-*r*98.8%
distribute-rgt-in98.9%
*-rgt-identity98.9%
distribute-lft-in98.9%
*-commutative98.9%
+-commutative98.9%
mul-1-neg98.9%
sub-neg98.9%
*-commutative98.9%
associate-*r*98.9%
*-commutative98.9%
Simplified98.9%
Taylor expanded in ux around 0 98.8%
associate-*r*98.7%
Simplified98.8%
Taylor expanded in ux around 0 94.5%
Final simplification94.5%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* ux (* maxCos (+ ux -1.0)))))
(+
(* zi (* ux (* (- 1.0 ux) maxCos)))
(+
(* xi (* (cos (* PI (* uy 2.0))) (sqrt (- 1.0 (* t_0 t_0)))))
(* 2.0 (* PI (* uy yi)))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ux * (maxCos * (ux + -1.0f));
return (zi * (ux * ((1.0f - ux) * maxCos))) + ((xi * (cosf((((float) M_PI) * (uy * 2.0f))) * sqrtf((1.0f - (t_0 * t_0))))) + (2.0f * (((float) M_PI) * (uy * yi))));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0)))) return Float32(Float32(zi * Float32(ux * Float32(Float32(Float32(1.0) - ux) * maxCos))) + Float32(Float32(xi * Float32(cos(Float32(Float32(pi) * Float32(uy * Float32(2.0)))) * sqrt(Float32(Float32(1.0) - Float32(t_0 * t_0))))) + Float32(Float32(2.0) * Float32(Float32(pi) * Float32(uy * yi))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = ux * (maxCos * (ux + single(-1.0))); tmp = (zi * (ux * ((single(1.0) - ux) * maxCos))) + ((xi * (cos((single(pi) * (uy * single(2.0)))) * sqrt((single(1.0) - (t_0 * t_0))))) + (single(2.0) * (single(pi) * (uy * yi)))); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\\
zi \cdot \left(ux \cdot \left(\left(1 - ux\right) \cdot maxCos\right)\right) + \left(xi \cdot \left(\cos \left(\pi \cdot \left(uy \cdot 2\right)\right) \cdot \sqrt{1 - t_0 \cdot t_0}\right) + 2 \cdot \left(\pi \cdot \left(uy \cdot yi\right)\right)\right)
\end{array}
\end{array}
Initial program 98.9%
Taylor expanded in ux around 0 98.7%
associate-*r*98.7%
Simplified98.7%
*-commutative98.7%
*-commutative98.7%
associate-*r*98.7%
add-cbrt-cube98.7%
add-cbrt-cube89.0%
cbrt-unprod83.2%
pow383.3%
associate-*r*83.3%
*-commutative83.3%
*-commutative83.3%
*-commutative83.3%
pow383.3%
Applied egg-rr83.3%
Taylor expanded in uy around 0 88.5%
associate-*r*88.5%
Simplified88.5%
Final simplification88.5%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* ux (* maxCos (+ ux -1.0)))))
(+
(* zi (* ux (* (- 1.0 ux) maxCos)))
(+
(* xi (* (cos (* PI (* uy 2.0))) (sqrt (- 1.0 (* t_0 t_0)))))
(* (* uy (* 2.0 PI)) yi)))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ux * (maxCos * (ux + -1.0f));
return (zi * (ux * ((1.0f - ux) * maxCos))) + ((xi * (cosf((((float) M_PI) * (uy * 2.0f))) * sqrtf((1.0f - (t_0 * t_0))))) + ((uy * (2.0f * ((float) M_PI))) * yi));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0)))) return Float32(Float32(zi * Float32(ux * Float32(Float32(Float32(1.0) - ux) * maxCos))) + Float32(Float32(xi * Float32(cos(Float32(Float32(pi) * Float32(uy * Float32(2.0)))) * sqrt(Float32(Float32(1.0) - Float32(t_0 * t_0))))) + Float32(Float32(uy * Float32(Float32(2.0) * Float32(pi))) * yi))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = ux * (maxCos * (ux + single(-1.0))); tmp = (zi * (ux * ((single(1.0) - ux) * maxCos))) + ((xi * (cos((single(pi) * (uy * single(2.0)))) * sqrt((single(1.0) - (t_0 * t_0))))) + ((uy * (single(2.0) * single(pi))) * yi)); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\\
zi \cdot \left(ux \cdot \left(\left(1 - ux\right) \cdot maxCos\right)\right) + \left(xi \cdot \left(\cos \left(\pi \cdot \left(uy \cdot 2\right)\right) \cdot \sqrt{1 - t_0 \cdot t_0}\right) + \left(uy \cdot \left(2 \cdot \pi\right)\right) \cdot yi\right)
\end{array}
\end{array}
Initial program 98.9%
Taylor expanded in ux around 0 98.7%
associate-*r*98.7%
Simplified98.7%
Taylor expanded in uy around 0 88.5%
associate-*r*88.5%
*-commutative88.5%
*-commutative88.5%
associate-*r*88.5%
associate-*l*88.5%
Simplified88.5%
Final simplification88.5%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(let* ((t_0 (* ux (* maxCos (+ ux -1.0)))))
(+
(* zi (* (- 1.0 ux) (* ux maxCos)))
(+
(* xi (* (cos (* PI (* uy 2.0))) (sqrt (- 1.0 (* t_0 t_0)))))
(* (* uy (* 2.0 PI)) yi)))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
float t_0 = ux * (maxCos * (ux + -1.0f));
return (zi * ((1.0f - ux) * (ux * maxCos))) + ((xi * (cosf((((float) M_PI) * (uy * 2.0f))) * sqrtf((1.0f - (t_0 * t_0))))) + ((uy * (2.0f * ((float) M_PI))) * yi));
}
function code(xi, yi, zi, ux, uy, maxCos) t_0 = Float32(ux * Float32(maxCos * Float32(ux + Float32(-1.0)))) return Float32(Float32(zi * Float32(Float32(Float32(1.0) - ux) * Float32(ux * maxCos))) + Float32(Float32(xi * Float32(cos(Float32(Float32(pi) * Float32(uy * Float32(2.0)))) * sqrt(Float32(Float32(1.0) - Float32(t_0 * t_0))))) + Float32(Float32(uy * Float32(Float32(2.0) * Float32(pi))) * yi))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) t_0 = ux * (maxCos * (ux + single(-1.0))); tmp = (zi * ((single(1.0) - ux) * (ux * maxCos))) + ((xi * (cos((single(pi) * (uy * single(2.0)))) * sqrt((single(1.0) - (t_0 * t_0))))) + ((uy * (single(2.0) * single(pi))) * yi)); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := ux \cdot \left(maxCos \cdot \left(ux + -1\right)\right)\\
zi \cdot \left(\left(1 - ux\right) \cdot \left(ux \cdot maxCos\right)\right) + \left(xi \cdot \left(\cos \left(\pi \cdot \left(uy \cdot 2\right)\right) \cdot \sqrt{1 - t_0 \cdot t_0}\right) + \left(uy \cdot \left(2 \cdot \pi\right)\right) \cdot yi\right)
\end{array}
\end{array}
Initial program 98.9%
Taylor expanded in ux around 0 98.8%
associate-*r*98.8%
neg-mul-198.8%
unpow298.8%
associate-*r*98.8%
neg-mul-198.8%
*-commutative98.8%
associate-*r*98.8%
distribute-rgt-in98.9%
*-rgt-identity98.9%
distribute-lft-in98.9%
*-commutative98.9%
+-commutative98.9%
mul-1-neg98.9%
sub-neg98.9%
*-commutative98.9%
associate-*r*98.9%
*-commutative98.9%
Simplified98.9%
Taylor expanded in ux around 0 98.8%
associate-*r*98.7%
Simplified98.8%
Taylor expanded in uy around 0 88.5%
associate-*r*88.5%
*-commutative88.5%
*-commutative88.5%
associate-*r*88.5%
associate-*l*88.5%
Simplified88.5%
Final simplification88.5%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (fma (cos (* uy (* 2.0 PI))) (* xi (sqrt (+ 1.0 (* (* maxCos (* ux (* ux maxCos))) (+ ux -1.0))))) (- (* maxCos (* ux zi)) (* maxCos (* ux (* ux zi))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf(cosf((uy * (2.0f * ((float) M_PI)))), (xi * sqrtf((1.0f + ((maxCos * (ux * (ux * maxCos))) * (ux + -1.0f))))), ((maxCos * (ux * zi)) - (maxCos * (ux * (ux * zi)))));
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(cos(Float32(uy * Float32(Float32(2.0) * Float32(pi)))), Float32(xi * sqrt(Float32(Float32(1.0) + Float32(Float32(maxCos * Float32(ux * Float32(ux * maxCos))) * Float32(ux + Float32(-1.0)))))), Float32(Float32(maxCos * Float32(ux * zi)) - Float32(maxCos * Float32(ux * Float32(ux * zi))))) end
\begin{array}{l}
\\
\mathsf{fma}\left(\cos \left(uy \cdot \left(2 \cdot \pi\right)\right), xi \cdot \sqrt{1 + \left(maxCos \cdot \left(ux \cdot \left(ux \cdot maxCos\right)\right)\right) \cdot \left(ux + -1\right)}, maxCos \cdot \left(ux \cdot zi\right) - maxCos \cdot \left(ux \cdot \left(ux \cdot zi\right)\right)\right)
\end{array}
Initial program 98.9%
Simplified98.9%
Taylor expanded in uy around 0 54.4%
*-commutative54.4%
associate-*r*54.3%
*-commutative54.3%
associate-*r*54.3%
associate-*r*54.3%
*-commutative54.3%
associate-*l*54.3%
sub-neg54.3%
mul-1-neg54.3%
distribute-lft-in54.3%
*-rgt-identity54.3%
mul-1-neg54.3%
distribute-rgt-neg-out54.3%
unsub-neg54.3%
*-commutative54.3%
Simplified54.3%
Taylor expanded in ux around 0 54.0%
Taylor expanded in ux around 0 54.0%
+-commutative54.0%
mul-1-neg54.0%
unsub-neg54.0%
unpow254.0%
associate-*l*54.1%
Simplified54.1%
Final simplification54.1%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(fma
(cos (* uy (* 2.0 PI)))
(*
(sqrt
(+ 1.0 (* (- 1.0 ux) (* maxCos (* (* ux maxCos) (* ux (+ ux -1.0)))))))
xi)
(* ux (* zi (- maxCos (* ux maxCos))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf(cosf((uy * (2.0f * ((float) M_PI)))), (sqrtf((1.0f + ((1.0f - ux) * (maxCos * ((ux * maxCos) * (ux * (ux + -1.0f))))))) * xi), (ux * (zi * (maxCos - (ux * maxCos)))));
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(cos(Float32(uy * Float32(Float32(2.0) * Float32(pi)))), Float32(sqrt(Float32(Float32(1.0) + Float32(Float32(Float32(1.0) - ux) * Float32(maxCos * Float32(Float32(ux * maxCos) * Float32(ux * Float32(ux + Float32(-1.0)))))))) * xi), Float32(ux * Float32(zi * Float32(maxCos - Float32(ux * maxCos))))) end
\begin{array}{l}
\\
\mathsf{fma}\left(\cos \left(uy \cdot \left(2 \cdot \pi\right)\right), \sqrt{1 + \left(1 - ux\right) \cdot \left(maxCos \cdot \left(\left(ux \cdot maxCos\right) \cdot \left(ux \cdot \left(ux + -1\right)\right)\right)\right)} \cdot xi, ux \cdot \left(zi \cdot \left(maxCos - ux \cdot maxCos\right)\right)\right)
\end{array}
Initial program 98.9%
Simplified98.9%
Taylor expanded in uy around 0 54.4%
*-commutative54.4%
associate-*r*54.3%
*-commutative54.3%
associate-*r*54.3%
associate-*r*54.3%
*-commutative54.3%
associate-*l*54.3%
sub-neg54.3%
mul-1-neg54.3%
distribute-lft-in54.3%
*-rgt-identity54.3%
mul-1-neg54.3%
distribute-rgt-neg-out54.3%
unsub-neg54.3%
*-commutative54.3%
Simplified54.3%
Final simplification54.3%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(fma
(cos (* uy (* 2.0 PI)))
(*
(sqrt
(+ 1.0 (* (- 1.0 ux) (* maxCos (* (* ux maxCos) (* ux (+ ux -1.0)))))))
xi)
(* (* (- 1.0 ux) maxCos) (* ux zi))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf(cosf((uy * (2.0f * ((float) M_PI)))), (sqrtf((1.0f + ((1.0f - ux) * (maxCos * ((ux * maxCos) * (ux * (ux + -1.0f))))))) * xi), (((1.0f - ux) * maxCos) * (ux * zi)));
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(cos(Float32(uy * Float32(Float32(2.0) * Float32(pi)))), Float32(sqrt(Float32(Float32(1.0) + Float32(Float32(Float32(1.0) - ux) * Float32(maxCos * Float32(Float32(ux * maxCos) * Float32(ux * Float32(ux + Float32(-1.0)))))))) * xi), Float32(Float32(Float32(Float32(1.0) - ux) * maxCos) * Float32(ux * zi))) end
\begin{array}{l}
\\
\mathsf{fma}\left(\cos \left(uy \cdot \left(2 \cdot \pi\right)\right), \sqrt{1 + \left(1 - ux\right) \cdot \left(maxCos \cdot \left(\left(ux \cdot maxCos\right) \cdot \left(ux \cdot \left(ux + -1\right)\right)\right)\right)} \cdot xi, \left(\left(1 - ux\right) \cdot maxCos\right) \cdot \left(ux \cdot zi\right)\right)
\end{array}
Initial program 98.9%
Simplified98.9%
Taylor expanded in uy around 0 54.4%
associate-*r*54.3%
Simplified54.3%
Final simplification54.3%
(FPCore (xi yi zi ux uy maxCos)
:precision binary32
(fma
(cos (* uy (* 2.0 PI)))
(*
(sqrt
(+ 1.0 (* (- 1.0 ux) (* maxCos (* (* ux maxCos) (* ux (+ ux -1.0)))))))
xi)
(* (- 1.0 ux) (* maxCos (* ux zi)))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf(cosf((uy * (2.0f * ((float) M_PI)))), (sqrtf((1.0f + ((1.0f - ux) * (maxCos * ((ux * maxCos) * (ux * (ux + -1.0f))))))) * xi), ((1.0f - ux) * (maxCos * (ux * zi))));
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(cos(Float32(uy * Float32(Float32(2.0) * Float32(pi)))), Float32(sqrt(Float32(Float32(1.0) + Float32(Float32(Float32(1.0) - ux) * Float32(maxCos * Float32(Float32(ux * maxCos) * Float32(ux * Float32(ux + Float32(-1.0)))))))) * xi), Float32(Float32(Float32(1.0) - ux) * Float32(maxCos * Float32(ux * zi)))) end
\begin{array}{l}
\\
\mathsf{fma}\left(\cos \left(uy \cdot \left(2 \cdot \pi\right)\right), \sqrt{1 + \left(1 - ux\right) \cdot \left(maxCos \cdot \left(\left(ux \cdot maxCos\right) \cdot \left(ux \cdot \left(ux + -1\right)\right)\right)\right)} \cdot xi, \left(1 - ux\right) \cdot \left(maxCos \cdot \left(ux \cdot zi\right)\right)\right)
\end{array}
Initial program 98.9%
Simplified98.9%
Taylor expanded in uy around 0 54.4%
Final simplification54.4%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (fma (cos (* uy (* 2.0 PI))) (* xi (sqrt (+ 1.0 (* (* maxCos (* ux (* ux maxCos))) (+ ux -1.0))))) (* ux (- (* maxCos zi) (* maxCos (* ux zi))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf(cosf((uy * (2.0f * ((float) M_PI)))), (xi * sqrtf((1.0f + ((maxCos * (ux * (ux * maxCos))) * (ux + -1.0f))))), (ux * ((maxCos * zi) - (maxCos * (ux * zi)))));
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(cos(Float32(uy * Float32(Float32(2.0) * Float32(pi)))), Float32(xi * sqrt(Float32(Float32(1.0) + Float32(Float32(maxCos * Float32(ux * Float32(ux * maxCos))) * Float32(ux + Float32(-1.0)))))), Float32(ux * Float32(Float32(maxCos * zi) - Float32(maxCos * Float32(ux * zi))))) end
\begin{array}{l}
\\
\mathsf{fma}\left(\cos \left(uy \cdot \left(2 \cdot \pi\right)\right), xi \cdot \sqrt{1 + \left(maxCos \cdot \left(ux \cdot \left(ux \cdot maxCos\right)\right)\right) \cdot \left(ux + -1\right)}, ux \cdot \left(maxCos \cdot zi - maxCos \cdot \left(ux \cdot zi\right)\right)\right)
\end{array}
Initial program 98.9%
Simplified98.9%
Taylor expanded in uy around 0 54.4%
*-commutative54.4%
associate-*r*54.3%
*-commutative54.3%
associate-*r*54.3%
associate-*r*54.3%
*-commutative54.3%
associate-*l*54.3%
sub-neg54.3%
mul-1-neg54.3%
distribute-lft-in54.3%
*-rgt-identity54.3%
mul-1-neg54.3%
distribute-rgt-neg-out54.3%
unsub-neg54.3%
*-commutative54.3%
Simplified54.3%
Taylor expanded in ux around 0 54.0%
Taylor expanded in ux around 0 54.0%
Final simplification54.0%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (fma (cos (* uy (* 2.0 PI))) (* xi (sqrt (+ 1.0 (* (* maxCos (* ux (* ux maxCos))) (+ ux -1.0))))) (* maxCos (* zi (- ux (* ux ux))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf(cosf((uy * (2.0f * ((float) M_PI)))), (xi * sqrtf((1.0f + ((maxCos * (ux * (ux * maxCos))) * (ux + -1.0f))))), (maxCos * (zi * (ux - (ux * ux)))));
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(cos(Float32(uy * Float32(Float32(2.0) * Float32(pi)))), Float32(xi * sqrt(Float32(Float32(1.0) + Float32(Float32(maxCos * Float32(ux * Float32(ux * maxCos))) * Float32(ux + Float32(-1.0)))))), Float32(maxCos * Float32(zi * Float32(ux - Float32(ux * ux))))) end
\begin{array}{l}
\\
\mathsf{fma}\left(\cos \left(uy \cdot \left(2 \cdot \pi\right)\right), xi \cdot \sqrt{1 + \left(maxCos \cdot \left(ux \cdot \left(ux \cdot maxCos\right)\right)\right) \cdot \left(ux + -1\right)}, maxCos \cdot \left(zi \cdot \left(ux - ux \cdot ux\right)\right)\right)
\end{array}
Initial program 98.9%
Simplified98.9%
Taylor expanded in uy around 0 54.4%
*-commutative54.4%
associate-*r*54.3%
*-commutative54.3%
associate-*r*54.3%
associate-*r*54.3%
*-commutative54.3%
associate-*l*54.3%
sub-neg54.3%
mul-1-neg54.3%
distribute-lft-in54.3%
*-rgt-identity54.3%
mul-1-neg54.3%
distribute-rgt-neg-out54.3%
unsub-neg54.3%
*-commutative54.3%
Simplified54.3%
Taylor expanded in ux around 0 54.0%
Taylor expanded in ux around 0 54.0%
+-commutative54.0%
mul-1-neg54.0%
unsub-neg54.0%
unpow254.0%
associate-*l*54.1%
Simplified54.1%
Taylor expanded in maxCos around 0 54.0%
*-commutative54.0%
distribute-rgt-out--54.0%
unpow254.0%
Simplified54.0%
Final simplification54.0%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (fma (cos (* uy (* 2.0 PI))) (* xi (sqrt (+ 1.0 (* (* maxCos (* ux (* ux maxCos))) (+ ux -1.0))))) (* ux (* zi (- maxCos (* ux maxCos))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return fmaf(cosf((uy * (2.0f * ((float) M_PI)))), (xi * sqrtf((1.0f + ((maxCos * (ux * (ux * maxCos))) * (ux + -1.0f))))), (ux * (zi * (maxCos - (ux * maxCos)))));
}
function code(xi, yi, zi, ux, uy, maxCos) return fma(cos(Float32(uy * Float32(Float32(2.0) * Float32(pi)))), Float32(xi * sqrt(Float32(Float32(1.0) + Float32(Float32(maxCos * Float32(ux * Float32(ux * maxCos))) * Float32(ux + Float32(-1.0)))))), Float32(ux * Float32(zi * Float32(maxCos - Float32(ux * maxCos))))) end
\begin{array}{l}
\\
\mathsf{fma}\left(\cos \left(uy \cdot \left(2 \cdot \pi\right)\right), xi \cdot \sqrt{1 + \left(maxCos \cdot \left(ux \cdot \left(ux \cdot maxCos\right)\right)\right) \cdot \left(ux + -1\right)}, ux \cdot \left(zi \cdot \left(maxCos - ux \cdot maxCos\right)\right)\right)
\end{array}
Initial program 98.9%
Simplified98.9%
Taylor expanded in uy around 0 54.4%
*-commutative54.4%
associate-*r*54.3%
*-commutative54.3%
associate-*r*54.3%
associate-*r*54.3%
*-commutative54.3%
associate-*l*54.3%
sub-neg54.3%
mul-1-neg54.3%
distribute-lft-in54.3%
*-rgt-identity54.3%
mul-1-neg54.3%
distribute-rgt-neg-out54.3%
unsub-neg54.3%
*-commutative54.3%
Simplified54.3%
Taylor expanded in ux around 0 54.0%
Final simplification54.0%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (+ (* maxCos (* ux zi)) (* (cos (* 2.0 (* uy PI))) (* xi (sqrt (+ 1.0 (* (pow (* ux maxCos) 2.0) (+ ux -1.0))))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return (maxCos * (ux * zi)) + (cosf((2.0f * (uy * ((float) M_PI)))) * (xi * sqrtf((1.0f + (powf((ux * maxCos), 2.0f) * (ux + -1.0f))))));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(maxCos * Float32(ux * zi)) + Float32(cos(Float32(Float32(2.0) * Float32(uy * Float32(pi)))) * Float32(xi * sqrt(Float32(Float32(1.0) + Float32((Float32(ux * maxCos) ^ Float32(2.0)) * Float32(ux + Float32(-1.0)))))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = (maxCos * (ux * zi)) + (cos((single(2.0) * (uy * single(pi)))) * (xi * sqrt((single(1.0) + (((ux * maxCos) ^ single(2.0)) * (ux + single(-1.0))))))); end
\begin{array}{l}
\\
maxCos \cdot \left(ux \cdot zi\right) + \cos \left(2 \cdot \left(uy \cdot \pi\right)\right) \cdot \left(xi \cdot \sqrt{1 + {\left(ux \cdot maxCos\right)}^{2} \cdot \left(ux + -1\right)}\right)
\end{array}
Initial program 98.9%
Simplified98.9%
Taylor expanded in uy around 0 54.4%
*-commutative54.4%
associate-*r*54.3%
*-commutative54.3%
associate-*r*54.3%
associate-*r*54.3%
*-commutative54.3%
associate-*l*54.3%
sub-neg54.3%
mul-1-neg54.3%
distribute-lft-in54.3%
*-rgt-identity54.3%
mul-1-neg54.3%
distribute-rgt-neg-out54.3%
unsub-neg54.3%
*-commutative54.3%
Simplified54.3%
Taylor expanded in ux around 0 54.0%
Taylor expanded in ux around 0 50.9%
fma-udef50.9%
associate-*r*50.9%
*-commutative50.9%
associate-*l*50.9%
associate-*r*50.9%
*-commutative50.9%
pow250.9%
Applied egg-rr50.9%
Final simplification50.9%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (+ (* maxCos (* ux zi)) (* xi (cos (* 2.0 (* uy PI))))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return (maxCos * (ux * zi)) + (xi * cosf((2.0f * (uy * ((float) M_PI)))));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(Float32(maxCos * Float32(ux * zi)) + Float32(xi * cos(Float32(Float32(2.0) * Float32(uy * Float32(pi)))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = (maxCos * (ux * zi)) + (xi * cos((single(2.0) * (uy * single(pi))))); end
\begin{array}{l}
\\
maxCos \cdot \left(ux \cdot zi\right) + xi \cdot \cos \left(2 \cdot \left(uy \cdot \pi\right)\right)
\end{array}
Initial program 98.9%
Simplified98.9%
Taylor expanded in uy around 0 54.4%
*-commutative54.4%
associate-*r*54.3%
*-commutative54.3%
associate-*r*54.3%
associate-*r*54.3%
*-commutative54.3%
associate-*l*54.3%
sub-neg54.3%
mul-1-neg54.3%
distribute-lft-in54.3%
*-rgt-identity54.3%
mul-1-neg54.3%
distribute-rgt-neg-out54.3%
unsub-neg54.3%
*-commutative54.3%
Simplified54.3%
Taylor expanded in ux around 0 54.0%
Taylor expanded in ux around 0 50.9%
Taylor expanded in ux around 0 50.9%
Final simplification50.9%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (* xi (cos (* 2.0 (* uy PI)))))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return xi * cosf((2.0f * (uy * ((float) M_PI))));
}
function code(xi, yi, zi, ux, uy, maxCos) return Float32(xi * cos(Float32(Float32(2.0) * Float32(uy * Float32(pi))))) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = xi * cos((single(2.0) * (uy * single(pi)))); end
\begin{array}{l}
\\
xi \cdot \cos \left(2 \cdot \left(uy \cdot \pi\right)\right)
\end{array}
Initial program 98.9%
Simplified98.9%
Taylor expanded in uy around 0 54.4%
*-commutative54.4%
associate-*r*54.3%
*-commutative54.3%
associate-*r*54.3%
associate-*r*54.3%
*-commutative54.3%
associate-*l*54.3%
sub-neg54.3%
mul-1-neg54.3%
distribute-lft-in54.3%
*-rgt-identity54.3%
mul-1-neg54.3%
distribute-rgt-neg-out54.3%
unsub-neg54.3%
*-commutative54.3%
Simplified54.3%
Taylor expanded in ux around 0 54.0%
Taylor expanded in ux around 0 50.9%
Taylor expanded in ux around 0 44.5%
Final simplification44.5%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (* maxCos (* ux zi)))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return maxCos * (ux * zi);
}
real(4) function code(xi, yi, zi, ux, uy, maxcos)
real(4), intent (in) :: xi
real(4), intent (in) :: yi
real(4), intent (in) :: zi
real(4), intent (in) :: ux
real(4), intent (in) :: uy
real(4), intent (in) :: maxcos
code = maxcos * (ux * zi)
end function
function code(xi, yi, zi, ux, uy, maxCos) return Float32(maxCos * Float32(ux * zi)) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = maxCos * (ux * zi); end
\begin{array}{l}
\\
maxCos \cdot \left(ux \cdot zi\right)
\end{array}
Initial program 98.9%
Simplified98.9%
Taylor expanded in uy around 0 54.4%
*-commutative54.4%
associate-*r*54.3%
*-commutative54.3%
associate-*r*54.3%
associate-*r*54.3%
*-commutative54.3%
associate-*l*54.3%
sub-neg54.3%
mul-1-neg54.3%
distribute-lft-in54.3%
*-rgt-identity54.3%
mul-1-neg54.3%
distribute-rgt-neg-out54.3%
unsub-neg54.3%
*-commutative54.3%
Simplified54.3%
Taylor expanded in ux around 0 54.0%
Taylor expanded in ux around 0 50.9%
Taylor expanded in xi around 0 13.3%
Final simplification13.3%
(FPCore (xi yi zi ux uy maxCos) :precision binary32 (* ux (* maxCos zi)))
float code(float xi, float yi, float zi, float ux, float uy, float maxCos) {
return ux * (maxCos * zi);
}
real(4) function code(xi, yi, zi, ux, uy, maxcos)
real(4), intent (in) :: xi
real(4), intent (in) :: yi
real(4), intent (in) :: zi
real(4), intent (in) :: ux
real(4), intent (in) :: uy
real(4), intent (in) :: maxcos
code = ux * (maxcos * zi)
end function
function code(xi, yi, zi, ux, uy, maxCos) return Float32(ux * Float32(maxCos * zi)) end
function tmp = code(xi, yi, zi, ux, uy, maxCos) tmp = ux * (maxCos * zi); end
\begin{array}{l}
\\
ux \cdot \left(maxCos \cdot zi\right)
\end{array}
Initial program 98.9%
Simplified98.9%
Taylor expanded in uy around 0 54.4%
*-commutative54.4%
associate-*r*54.3%
*-commutative54.3%
associate-*r*54.3%
associate-*r*54.3%
*-commutative54.3%
associate-*l*54.3%
sub-neg54.3%
mul-1-neg54.3%
distribute-lft-in54.3%
*-rgt-identity54.3%
mul-1-neg54.3%
distribute-rgt-neg-out54.3%
unsub-neg54.3%
*-commutative54.3%
Simplified54.3%
Taylor expanded in ux around 0 54.0%
Taylor expanded in ux around 0 50.9%
Taylor expanded in xi around 0 13.3%
*-commutative13.3%
associate-*l*13.3%
Simplified13.3%
Final simplification13.3%
herbie shell --seed 2023274
(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)))