
(FPCore (ux uy maxCos) :precision binary32 (let* ((t_0 (+ (- 1.0 ux) (* ux maxCos)))) (* (sin (* (* uy 2.0) (PI))) (sqrt (- 1.0 (* t_0 t_0))))))
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(1 - ux\right) + ux \cdot maxCos\\
\sin \left(\left(uy \cdot 2\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \sqrt{1 - t\_0 \cdot t\_0}
\end{array}
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 17 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (ux uy maxCos) :precision binary32 (let* ((t_0 (+ (- 1.0 ux) (* ux maxCos)))) (* (sin (* (* uy 2.0) (PI))) (sqrt (- 1.0 (* t_0 t_0))))))
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(1 - ux\right) + ux \cdot maxCos\\
\sin \left(\left(uy \cdot 2\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \sqrt{1 - t\_0 \cdot t\_0}
\end{array}
\end{array}
(FPCore (ux uy maxCos)
:precision binary32
(*
(sqrt
(*
(* ux ux)
(-
(/ (- (* (/ ux maxCos) (- 2.0 maxCos)) ux) (* (/ ux maxCos) ux))
(pow (- 1.0 maxCos) 2.0))))
(sin (* (PI) (* 2.0 uy)))))\begin{array}{l}
\\
\sqrt{\left(ux \cdot ux\right) \cdot \left(\frac{\frac{ux}{maxCos} \cdot \left(2 - maxCos\right) - ux}{\frac{ux}{maxCos} \cdot ux} - {\left(1 - maxCos\right)}^{2}\right)} \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(2 \cdot uy\right)\right)
\end{array}
Initial program 55.5%
lift-*.f32N/A
lift-+.f32N/A
lift--.f32N/A
sub-negN/A
associate-+l+N/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-+.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lower-*.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
neg-mul-1N/A
lift-*.f32N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites18.6%
Taylor expanded in ux around -inf
*-commutativeN/A
lower-*.f32N/A
Applied rewrites98.2%
Applied rewrites98.2%
Final simplification98.2%
(FPCore (ux uy maxCos)
:precision binary32
(*
(sqrt
(*
(-
(/ (- (* ux (- 2.0 maxCos)) (* ux maxCos)) (* ux ux))
(pow (- 1.0 maxCos) 2.0))
(* ux ux)))
(sin (* (PI) (* 2.0 uy)))))\begin{array}{l}
\\
\sqrt{\left(\frac{ux \cdot \left(2 - maxCos\right) - ux \cdot maxCos}{ux \cdot ux} - {\left(1 - maxCos\right)}^{2}\right) \cdot \left(ux \cdot ux\right)} \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(2 \cdot uy\right)\right)
\end{array}
Initial program 55.5%
lift-*.f32N/A
lift-+.f32N/A
lift--.f32N/A
sub-negN/A
associate-+l+N/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-+.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lower-*.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
neg-mul-1N/A
lift-*.f32N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites17.2%
Taylor expanded in ux around -inf
*-commutativeN/A
lower-*.f32N/A
Applied rewrites98.2%
Applied rewrites98.2%
Final simplification98.2%
(FPCore (ux uy maxCos) :precision binary32 (* (sqrt (* (* (- (/ (- (- 2.0 maxCos) maxCos) ux) (pow (- 1.0 maxCos) 2.0)) ux) ux)) (sin (* (PI) (* 2.0 uy)))))
\begin{array}{l}
\\
\sqrt{\left(\left(\frac{\left(2 - maxCos\right) - maxCos}{ux} - {\left(1 - maxCos\right)}^{2}\right) \cdot ux\right) \cdot ux} \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(2 \cdot uy\right)\right)
\end{array}
Initial program 55.5%
lift-*.f32N/A
lift-+.f32N/A
lift--.f32N/A
sub-negN/A
associate-+l+N/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-+.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lower-*.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
neg-mul-1N/A
lift-*.f32N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites18.6%
Taylor expanded in ux around -inf
*-commutativeN/A
lower-*.f32N/A
Applied rewrites98.2%
Applied rewrites98.2%
Final simplification98.2%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= (* 2.0 uy) 0.000699999975040555)
(*
(sqrt
(-
(- ux (* ux maxCos))
(* (* (fma maxCos ux (- 1.0 ux)) ux) (- maxCos 1.0))))
(* (* (PI) 2.0) uy))
(* (sqrt (* (* (- (/ 2.0 ux) 1.0) ux) ux)) (sin (* (PI) (* 2.0 uy))))))\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;2 \cdot uy \leq 0.000699999975040555:\\
\;\;\;\;\sqrt{\left(ux - ux \cdot maxCos\right) - \left(\mathsf{fma}\left(maxCos, ux, 1 - ux\right) \cdot ux\right) \cdot \left(maxCos - 1\right)} \cdot \left(\left(\mathsf{PI}\left(\right) \cdot 2\right) \cdot uy\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(\left(\frac{2}{ux} - 1\right) \cdot ux\right) \cdot ux} \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(2 \cdot uy\right)\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 6.99999975e-4Initial program 51.4%
lift-*.f32N/A
lift-+.f32N/A
lift--.f32N/A
sub-negN/A
associate-+l+N/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-+.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lower-*.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
neg-mul-1N/A
lift-*.f32N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites21.9%
lift--.f32N/A
lift-+.f32N/A
associate--r+N/A
lower--.f32N/A
lift-fma.f32N/A
*-commutativeN/A
lift-*.f32N/A
+-commutativeN/A
lift--.f32N/A
associate-+l-N/A
associate--r-N/A
metadata-evalN/A
lower-+.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
lower--.f3271.2
lift-*.f32N/A
lift-*.f32N/A
associate-*r*N/A
lift-+.f32N/A
Applied rewrites70.0%
Taylor expanded in uy around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-PI.f3296.1
Applied rewrites95.7%
if 6.99999975e-4 < (*.f32 uy #s(literal 2 binary32)) Initial program 62.2%
lift-*.f32N/A
lift-+.f32N/A
lift--.f32N/A
sub-negN/A
associate-+l+N/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-+.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lower-*.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
neg-mul-1N/A
lift-*.f32N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites13.5%
Taylor expanded in ux around -inf
*-commutativeN/A
lower-*.f32N/A
Applied rewrites97.8%
Applied rewrites97.8%
Taylor expanded in maxCos around 0
Applied rewrites91.3%
Final simplification68.0%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= (* 2.0 uy) 0.000699999975040555)
(*
(sqrt
(-
(- ux (* ux maxCos))
(* (* (fma maxCos ux (- 1.0 ux)) ux) (- maxCos 1.0))))
(* (* (PI) 2.0) uy))
(* (sqrt (* (- (/ 2.0 ux) 1.0) (* ux ux))) (sin (* (PI) (* 2.0 uy))))))\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;2 \cdot uy \leq 0.000699999975040555:\\
\;\;\;\;\sqrt{\left(ux - ux \cdot maxCos\right) - \left(\mathsf{fma}\left(maxCos, ux, 1 - ux\right) \cdot ux\right) \cdot \left(maxCos - 1\right)} \cdot \left(\left(\mathsf{PI}\left(\right) \cdot 2\right) \cdot uy\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(\frac{2}{ux} - 1\right) \cdot \left(ux \cdot ux\right)} \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(2 \cdot uy\right)\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 6.99999975e-4Initial program 51.4%
lift-*.f32N/A
lift-+.f32N/A
lift--.f32N/A
sub-negN/A
associate-+l+N/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-+.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lower-*.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
neg-mul-1N/A
lift-*.f32N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites23.3%
lift--.f32N/A
lift-+.f32N/A
associate--r+N/A
lower--.f32N/A
lift-fma.f32N/A
*-commutativeN/A
lift-*.f32N/A
+-commutativeN/A
lift--.f32N/A
associate-+l-N/A
associate--r-N/A
metadata-evalN/A
lower-+.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
lower--.f3272.9
lift-*.f32N/A
lift-*.f32N/A
associate-*r*N/A
lift-+.f32N/A
Applied rewrites70.6%
Taylor expanded in uy around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-PI.f3293.2
Applied rewrites95.7%
if 6.99999975e-4 < (*.f32 uy #s(literal 2 binary32)) Initial program 62.2%
lift-*.f32N/A
lift-+.f32N/A
lift--.f32N/A
sub-negN/A
associate-+l+N/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-+.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lower-*.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
neg-mul-1N/A
lift-*.f32N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites13.5%
Taylor expanded in ux around -inf
*-commutativeN/A
lower-*.f32N/A
Applied rewrites97.8%
Taylor expanded in maxCos around 0
Applied rewrites91.2%
Final simplification67.7%
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0 (- (* ux maxCos) ux)))
(*
(sqrt (- (* (* (- -1.0 t_0) ux) (- maxCos 1.0)) t_0))
(sin (* (PI) (* 2.0 uy))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := ux \cdot maxCos - ux\\
\sqrt{\left(\left(-1 - t\_0\right) \cdot ux\right) \cdot \left(maxCos - 1\right) - t\_0} \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(2 \cdot uy\right)\right)
\end{array}
\end{array}
Initial program 55.5%
lift-*.f32N/A
lift-+.f32N/A
lift--.f32N/A
sub-negN/A
associate-+l+N/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-+.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lower-*.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
neg-mul-1N/A
lift-*.f32N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites17.3%
lift--.f32N/A
lift-+.f32N/A
associate--r+N/A
lower--.f32N/A
lift-fma.f32N/A
*-commutativeN/A
lift-*.f32N/A
+-commutativeN/A
lift--.f32N/A
associate-+l-N/A
associate--r-N/A
metadata-evalN/A
lower-+.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
lower--.f3248.9
lift-*.f32N/A
lift-*.f32N/A
associate-*r*N/A
lift-+.f32N/A
Applied rewrites50.4%
lift-fma.f32N/A
lift-*.f32N/A
+-commutativeN/A
lift--.f32N/A
associate-+l-N/A
lift--.f32N/A
+-lft-identityN/A
lift-+.f32N/A
lower--.f3298.2
lift-+.f32N/A
+-lft-identity98.2
lift-*.f32N/A
*-commutativeN/A
lift-*.f3298.2
Applied rewrites98.2%
Final simplification98.2%
(FPCore (ux uy maxCos)
:precision binary32
(*
(sqrt
(-
(* (- 1.0 maxCos) ux)
(* (* (fma maxCos ux (- 1.0 ux)) ux) (- maxCos 1.0))))
(sin (* (PI) (* 2.0 uy)))))\begin{array}{l}
\\
\sqrt{\left(1 - maxCos\right) \cdot ux - \left(\mathsf{fma}\left(maxCos, ux, 1 - ux\right) \cdot ux\right) \cdot \left(maxCos - 1\right)} \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(2 \cdot uy\right)\right)
\end{array}
Initial program 55.5%
lift-*.f32N/A
lift-+.f32N/A
lift--.f32N/A
sub-negN/A
associate-+l+N/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-+.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lower-*.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
neg-mul-1N/A
lift-*.f32N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites17.3%
lift--.f32N/A
lift-+.f32N/A
associate--r+N/A
lower--.f32N/A
lift-fma.f32N/A
*-commutativeN/A
lift-*.f32N/A
+-commutativeN/A
lift--.f32N/A
associate-+l-N/A
associate--r-N/A
metadata-evalN/A
lower-+.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
lower--.f3249.8
lift-*.f32N/A
lift-*.f32N/A
associate-*r*N/A
lift-+.f32N/A
Applied rewrites48.6%
lift-+.f32N/A
+-lft-identity48.6
lift--.f32N/A
lift-*.f32N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lift--.f32N/A
lower-*.f32N/A
lift--.f32N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lift--.f3249.0
Applied rewrites49.1%
Final simplification49.4%
(FPCore (ux uy maxCos) :precision binary32 (* (sqrt (- (- ux (* ux maxCos)) (* (* (- 1.0 ux) ux) (- maxCos 1.0)))) (sin (* (PI) (* 2.0 uy)))))
\begin{array}{l}
\\
\sqrt{\left(ux - ux \cdot maxCos\right) - \left(\left(1 - ux\right) \cdot ux\right) \cdot \left(maxCos - 1\right)} \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(2 \cdot uy\right)\right)
\end{array}
Initial program 55.5%
lift-*.f32N/A
lift-+.f32N/A
lift--.f32N/A
sub-negN/A
associate-+l+N/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-+.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lower-*.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
neg-mul-1N/A
lift-*.f32N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites17.4%
lift--.f32N/A
lift-+.f32N/A
associate--r+N/A
lower--.f32N/A
lift-fma.f32N/A
*-commutativeN/A
lift-*.f32N/A
+-commutativeN/A
lift--.f32N/A
associate-+l-N/A
associate--r-N/A
metadata-evalN/A
lower-+.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
lower--.f3250.7
lift-*.f32N/A
lift-*.f32N/A
associate-*r*N/A
lift-+.f32N/A
Applied rewrites48.4%
Taylor expanded in maxCos around 0
lower--.f3296.1
Applied rewrites96.1%
Final simplification96.1%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= (* 2.0 uy) 0.000699999975040555)
(*
(sqrt
(-
(- ux (* ux maxCos))
(* (* (fma maxCos ux (- 1.0 ux)) ux) (- maxCos 1.0))))
(* (* (PI) 2.0) uy))
(* (sqrt (+ (* (- 1.0 ux) ux) ux)) (sin (* (PI) (* 2.0 uy))))))\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;2 \cdot uy \leq 0.000699999975040555:\\
\;\;\;\;\sqrt{\left(ux - ux \cdot maxCos\right) - \left(\mathsf{fma}\left(maxCos, ux, 1 - ux\right) \cdot ux\right) \cdot \left(maxCos - 1\right)} \cdot \left(\left(\mathsf{PI}\left(\right) \cdot 2\right) \cdot uy\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(1 - ux\right) \cdot ux + ux} \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(2 \cdot uy\right)\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 6.99999975e-4Initial program 51.4%
lift-*.f32N/A
lift-+.f32N/A
lift--.f32N/A
sub-negN/A
associate-+l+N/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-+.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lower-*.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
neg-mul-1N/A
lift-*.f32N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites20.9%
lift--.f32N/A
lift-+.f32N/A
associate--r+N/A
lower--.f32N/A
lift-fma.f32N/A
*-commutativeN/A
lift-*.f32N/A
+-commutativeN/A
lift--.f32N/A
associate-+l-N/A
associate--r-N/A
metadata-evalN/A
lower-+.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
lower--.f3271.4
lift-*.f32N/A
lift-*.f32N/A
associate-*r*N/A
lift-+.f32N/A
Applied rewrites71.5%
Taylor expanded in uy around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-PI.f3296.1
Applied rewrites95.7%
if 6.99999975e-4 < (*.f32 uy #s(literal 2 binary32)) Initial program 62.2%
lift-*.f32N/A
lift-+.f32N/A
lift--.f32N/A
sub-negN/A
associate-+l+N/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-+.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lower-*.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
neg-mul-1N/A
lift-*.f32N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites13.5%
Taylor expanded in maxCos around 0
associate-*r*N/A
mul-1-negN/A
cancel-sign-subN/A
lower-+.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower--.f3291.0
Applied rewrites91.0%
Final simplification67.5%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= (* 2.0 uy) 0.007499999832361937)
(*
(sqrt
(-
(- ux (* ux maxCos))
(* (* (fma maxCos ux (- 1.0 ux)) ux) (- maxCos 1.0))))
(* (* (PI) 2.0) uy))
(* (sqrt (* (fma maxCos -2.0 2.0) ux)) (sin (* (PI) (* 2.0 uy))))))\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;2 \cdot uy \leq 0.007499999832361937:\\
\;\;\;\;\sqrt{\left(ux - ux \cdot maxCos\right) - \left(\mathsf{fma}\left(maxCos, ux, 1 - ux\right) \cdot ux\right) \cdot \left(maxCos - 1\right)} \cdot \left(\left(\mathsf{PI}\left(\right) \cdot 2\right) \cdot uy\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(maxCos, -2, 2\right) \cdot ux} \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \left(2 \cdot uy\right)\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.00749999983Initial program 53.5%
lift-*.f32N/A
lift-+.f32N/A
lift--.f32N/A
sub-negN/A
associate-+l+N/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-+.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lower-*.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
neg-mul-1N/A
lift-*.f32N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites22.0%
lift--.f32N/A
lift-+.f32N/A
associate--r+N/A
lower--.f32N/A
lift-fma.f32N/A
*-commutativeN/A
lift-*.f32N/A
+-commutativeN/A
lift--.f32N/A
associate-+l-N/A
associate--r-N/A
metadata-evalN/A
lower-+.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
lower--.f3263.8
lift-*.f32N/A
lift-*.f32N/A
associate-*r*N/A
lift-+.f32N/A
Applied rewrites62.4%
Taylor expanded in uy around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-PI.f3291.2
Applied rewrites92.4%
if 0.00749999983 < (*.f32 uy #s(literal 2 binary32)) Initial program 61.6%
Taylor expanded in ux around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f327.3
Applied rewrites7.3%
Final simplification54.0%
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0 (+ (- 1.0 ux) (* ux maxCos))))
(if (<= (- 1.0 (* t_0 t_0)) 0.0003000000142492354)
(* (* (+ uy uy) (PI)) (sqrt (* (fma maxCos -2.0 2.0) ux)))
(* (sqrt (- 1.0 (* (- 1.0 ux) t_0))) (* (* (PI) 2.0) uy)))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(1 - ux\right) + ux \cdot maxCos\\
\mathbf{if}\;1 - t\_0 \cdot t\_0 \leq 0.0003000000142492354:\\
\;\;\;\;\left(\left(uy + uy\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \sqrt{\mathsf{fma}\left(maxCos, -2, 2\right) \cdot ux}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{1 - \left(1 - ux\right) \cdot t\_0} \cdot \left(\left(\mathsf{PI}\left(\right) \cdot 2\right) \cdot uy\right)\\
\end{array}
\end{array}
if (-.f32 #s(literal 1 binary32) (*.f32 (+.f32 (-.f32 #s(literal 1 binary32) ux) (*.f32 ux maxCos)) (+.f32 (-.f32 #s(literal 1 binary32) ux) (*.f32 ux maxCos)))) < 3.00000014e-4Initial program 34.6%
Taylor expanded in uy around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-PI.f3231.0
Applied rewrites31.0%
Taylor expanded in ux around 0
Applied rewrites7.3%
Applied rewrites7.3%
Taylor expanded in ux around 0
*-commutativeN/A
lower-*.f32N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f3273.5
Applied rewrites72.8%
if 3.00000014e-4 < (-.f32 #s(literal 1 binary32) (*.f32 (+.f32 (-.f32 #s(literal 1 binary32) ux) (*.f32 ux maxCos)) (+.f32 (-.f32 #s(literal 1 binary32) ux) (*.f32 ux maxCos)))) Initial program 89.8%
Taylor expanded in uy around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-PI.f3274.0
Applied rewrites74.0%
Taylor expanded in maxCos around 0
lower--.f3269.2
Applied rewrites69.2%
Final simplification71.7%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= ux 5.999999848427251e-5)
(* (* (+ uy uy) (PI)) (sqrt (* (fma maxCos -2.0 2.0) ux)))
(*
(sqrt
(- 1.0 (* (+ (- 1.0 ux) (* ux maxCos)) (+ (- (* ux maxCos) ux) 1.0))))
(* (* (PI) 2.0) uy))))\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;ux \leq 5.999999848427251 \cdot 10^{-5}:\\
\;\;\;\;\left(\left(uy + uy\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \sqrt{\mathsf{fma}\left(maxCos, -2, 2\right) \cdot ux}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{1 - \left(\left(1 - ux\right) + ux \cdot maxCos\right) \cdot \left(\left(ux \cdot maxCos - ux\right) + 1\right)} \cdot \left(\left(\mathsf{PI}\left(\right) \cdot 2\right) \cdot uy\right)\\
\end{array}
\end{array}
if ux < 5.99999985e-5Initial program 33.0%
Taylor expanded in uy around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-PI.f3230.6
Applied rewrites30.6%
Taylor expanded in ux around 0
Applied rewrites7.3%
Applied rewrites7.3%
Taylor expanded in ux around 0
*-commutativeN/A
lower-*.f32N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f3275.2
Applied rewrites74.5%
if 5.99999985e-5 < ux Initial program 88.4%
Taylor expanded in uy around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-PI.f3271.6
Applied rewrites71.6%
lift-+.f32N/A
lift--.f32N/A
associate-+l-N/A
lower--.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
lower--.f3271.7
Applied rewrites71.7%
Final simplification72.1%
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0 (+ (- 1.0 ux) (* ux maxCos))))
(if (<= ux 5.999999848427251e-5)
(* (* (+ uy uy) (PI)) (sqrt (* (fma maxCos -2.0 2.0) ux)))
(* (sqrt (- 1.0 (* t_0 t_0))) (* (* (PI) 2.0) uy)))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(1 - ux\right) + ux \cdot maxCos\\
\mathbf{if}\;ux \leq 5.999999848427251 \cdot 10^{-5}:\\
\;\;\;\;\left(\left(uy + uy\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \sqrt{\mathsf{fma}\left(maxCos, -2, 2\right) \cdot ux}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{1 - t\_0 \cdot t\_0} \cdot \left(\left(\mathsf{PI}\left(\right) \cdot 2\right) \cdot uy\right)\\
\end{array}
\end{array}
if ux < 5.99999985e-5Initial program 33.0%
Taylor expanded in uy around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-PI.f3230.6
Applied rewrites30.6%
Taylor expanded in ux around 0
Applied rewrites7.3%
Applied rewrites7.3%
Taylor expanded in ux around 0
*-commutativeN/A
lower-*.f32N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f3271.5
Applied rewrites74.5%
if 5.99999985e-5 < ux Initial program 88.4%
Taylor expanded in uy around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-PI.f3271.6
Applied rewrites71.6%
Final simplification73.8%
(FPCore (ux uy maxCos)
:precision binary32
(*
(sqrt
(-
(- ux (* ux maxCos))
(* (* (fma maxCos ux (- 1.0 ux)) ux) (- maxCos 1.0))))
(* (* (PI) 2.0) uy)))\begin{array}{l}
\\
\sqrt{\left(ux - ux \cdot maxCos\right) - \left(\mathsf{fma}\left(maxCos, ux, 1 - ux\right) \cdot ux\right) \cdot \left(maxCos - 1\right)} \cdot \left(\left(\mathsf{PI}\left(\right) \cdot 2\right) \cdot uy\right)
\end{array}
Initial program 55.5%
lift-*.f32N/A
lift-+.f32N/A
lift--.f32N/A
sub-negN/A
associate-+l+N/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-+.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lower-*.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
neg-mul-1N/A
lift-*.f32N/A
*-commutativeN/A
distribute-rgt-outN/A
lower-*.f32N/A
Applied rewrites17.3%
lift--.f32N/A
lift-+.f32N/A
associate--r+N/A
lower--.f32N/A
lift-fma.f32N/A
*-commutativeN/A
lift-*.f32N/A
+-commutativeN/A
lift--.f32N/A
associate-+l-N/A
associate--r-N/A
metadata-evalN/A
lower-+.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
lower--.f3249.3
lift-*.f32N/A
lift-*.f32N/A
associate-*r*N/A
lift-+.f32N/A
Applied rewrites48.5%
Taylor expanded in uy around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-PI.f3279.5
Applied rewrites79.3%
Final simplification79.5%
(FPCore (ux uy maxCos) :precision binary32 (* (* (+ uy uy) (PI)) (sqrt (* (fma maxCos -2.0 2.0) ux))))
\begin{array}{l}
\\
\left(\left(uy + uy\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \sqrt{\mathsf{fma}\left(maxCos, -2, 2\right) \cdot ux}
\end{array}
Initial program 55.5%
Taylor expanded in uy around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-PI.f3247.3
Applied rewrites47.3%
Taylor expanded in ux around 0
Applied rewrites7.1%
Applied rewrites7.1%
Taylor expanded in ux around 0
*-commutativeN/A
lower-*.f32N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f3263.4
Applied rewrites63.3%
Final simplification63.4%
(FPCore (ux uy maxCos) :precision binary32 (* (sqrt (fma -1.0 1.0 1.0)) (* (+ uy uy) (PI))))
\begin{array}{l}
\\
\sqrt{\mathsf{fma}\left(-1, 1, 1\right)} \cdot \left(\left(uy + uy\right) \cdot \mathsf{PI}\left(\right)\right)
\end{array}
Initial program 55.5%
Taylor expanded in uy around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-PI.f3247.3
Applied rewrites47.3%
Taylor expanded in ux around 0
Applied rewrites7.1%
Applied rewrites7.1%
lift--.f32N/A
sub-negN/A
+-commutativeN/A
neg-mul-1N/A
lower-fma.f3220.3
Applied rewrites20.3%
Final simplification20.3%
(FPCore (ux uy maxCos) :precision binary32 (* (sqrt (- 1.0 1.0)) (+ uy uy)))
float code(float ux, float uy, float maxCos) {
return sqrtf((1.0f - 1.0f)) * (uy + uy);
}
real(4) function code(ux, uy, maxcos)
real(4), intent (in) :: ux
real(4), intent (in) :: uy
real(4), intent (in) :: maxcos
code = sqrt((1.0e0 - 1.0e0)) * (uy + uy)
end function
function code(ux, uy, maxCos) return Float32(sqrt(Float32(Float32(1.0) - Float32(1.0))) * Float32(uy + uy)) end
function tmp = code(ux, uy, maxCos) tmp = sqrt((single(1.0) - single(1.0))) * (uy + uy); end
\begin{array}{l}
\\
\sqrt{1 - 1} \cdot \left(uy + uy\right)
\end{array}
Initial program 55.5%
Taylor expanded in uy around 0
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-PI.f3247.3
Applied rewrites47.3%
Taylor expanded in ux around 0
Applied rewrites7.1%
Applied rewrites7.1%
Applied rewrites7.1%
Final simplification7.1%
herbie shell --seed 2024271
(FPCore (ux uy maxCos)
:name "UniformSampleCone, y"
:precision binary32
:pre (and (and (and (<= 2.328306437e-10 ux) (<= ux 1.0)) (and (<= 2.328306437e-10 uy) (<= uy 1.0))) (and (<= 0.0 maxCos) (<= maxCos 1.0)))
(* (sin (* (* uy 2.0) (PI))) (sqrt (- 1.0 (* (+ (- 1.0 ux) (* ux maxCos)) (+ (- 1.0 ux) (* ux maxCos)))))))