
(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 10 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
(*
(sin (* (* uy 2.0) (PI)))
(sqrt
(*
(* ux ux)
(+
(/ (+ 2.0 (* maxCos -2.0)) ux)
(* (- 1.0 maxCos) (* maxCos (- (/ -1.0 maxCos) -1.0))))))))\begin{array}{l}
\\
\sin \left(\left(uy \cdot 2\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \sqrt{\left(ux \cdot ux\right) \cdot \left(\frac{2 + maxCos \cdot -2}{ux} + \left(1 - maxCos\right) \cdot \left(maxCos \cdot \left(\frac{-1}{maxCos} - -1\right)\right)\right)}
\end{array}
Initial program 57.7%
Taylor expanded in maxCos around inf
lower-*.f32N/A
associate--l+N/A
div-subN/A
lower-+.f32N/A
lower-/.f32N/A
lower--.f3257.9
Applied rewrites57.9%
Taylor expanded in ux around -inf
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower--.f32N/A
lower-/.f32N/A
lower-*.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-+.f32N/A
mul-1-negN/A
lower-neg.f32N/A
sub-negN/A
Applied rewrites98.0%
lift-/.f32N/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
lift-/.f32N/A
div-invN/A
lift-/.f32N/A
associate-*l*N/A
lift-/.f32N/A
inv-powN/A
pow-plusN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
lower-+.f32N/A
lower-*.f3298.2
Applied rewrites98.2%
lift-neg.f32N/A
lift-+.f32N/A
lift-/.f32N/A
lift-+.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
associate-*l*N/A
lower-*.f32N/A
lift-+.f32N/A
lift-neg.f32N/A
unsub-negN/A
lower--.f32N/A
lower-*.f3298.3
Applied rewrites98.3%
Final simplification98.3%
(FPCore (ux uy maxCos) :precision binary32 (* (sin (* (* uy 2.0) (PI))) (sqrt (* (* ux ux) (+ (/ (+ 2.0 (* maxCos -2.0)) ux) -1.0)))))
\begin{array}{l}
\\
\sin \left(\left(uy \cdot 2\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \sqrt{\left(ux \cdot ux\right) \cdot \left(\frac{2 + maxCos \cdot -2}{ux} + -1\right)}
\end{array}
Initial program 57.7%
Taylor expanded in maxCos around inf
lower-*.f32N/A
associate--l+N/A
div-subN/A
lower-+.f32N/A
lower-/.f32N/A
lower--.f3257.9
Applied rewrites57.9%
Taylor expanded in ux around -inf
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower--.f32N/A
lower-/.f32N/A
lower-*.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-+.f32N/A
mul-1-negN/A
lower-neg.f32N/A
sub-negN/A
Applied rewrites98.0%
lift-/.f32N/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
lift-/.f32N/A
div-invN/A
lift-/.f32N/A
associate-*l*N/A
lift-/.f32N/A
inv-powN/A
pow-plusN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
lower-+.f32N/A
lower-*.f3298.2
Applied rewrites98.2%
Taylor expanded in maxCos around 0
Applied rewrites96.6%
Final simplification96.6%
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0 (* 2.0 (* uy (PI)))))
(if (<= (* uy 2.0) 0.0008299999753944576)
(*
t_0
(sqrt
(*
(* ux ux)
(+
(/ (+ 2.0 (* maxCos -2.0)) ux)
(* maxCos (* (- 1.0 maxCos) (- (/ -1.0 maxCos) -1.0)))))))
(* (* ux (sin t_0)) (sqrt (+ -1.0 (/ 2.0 ux)))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := 2 \cdot \left(uy \cdot \mathsf{PI}\left(\right)\right)\\
\mathbf{if}\;uy \cdot 2 \leq 0.0008299999753944576:\\
\;\;\;\;t\_0 \cdot \sqrt{\left(ux \cdot ux\right) \cdot \left(\frac{2 + maxCos \cdot -2}{ux} + maxCos \cdot \left(\left(1 - maxCos\right) \cdot \left(\frac{-1}{maxCos} - -1\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(ux \cdot \sin t\_0\right) \cdot \sqrt{-1 + \frac{2}{ux}}\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 8.29999975e-4Initial program 58.6%
Taylor expanded in maxCos around inf
lower-*.f32N/A
associate--l+N/A
div-subN/A
lower-+.f32N/A
lower-/.f32N/A
lower--.f3258.7
Applied rewrites58.7%
Taylor expanded in ux around -inf
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower--.f32N/A
lower-/.f32N/A
lower-*.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-+.f32N/A
mul-1-negN/A
lower-neg.f32N/A
sub-negN/A
Applied rewrites97.9%
lift-/.f32N/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
lift-/.f32N/A
div-invN/A
lift-/.f32N/A
associate-*l*N/A
lift-/.f32N/A
inv-powN/A
pow-plusN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
lower-+.f32N/A
lower-*.f3298.4
Applied rewrites98.4%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3297.8
Applied rewrites97.8%
if 8.29999975e-4 < (*.f32 uy #s(literal 2 binary32)) Initial program 55.8%
Taylor expanded in maxCos around inf
lower-*.f32N/A
associate--l+N/A
div-subN/A
lower-+.f32N/A
lower-/.f32N/A
lower--.f3256.2
Applied rewrites56.2%
Taylor expanded in ux around -inf
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower--.f32N/A
lower-/.f32N/A
lower-*.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-+.f32N/A
mul-1-negN/A
lower-neg.f32N/A
sub-negN/A
Applied rewrites98.2%
Taylor expanded in maxCos around 0
Applied rewrites97.7%
Taylor expanded in maxCos around 0
lower-*.f32N/A
lower-*.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-sqrt.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f3293.9
Applied rewrites93.9%
Final simplification96.5%
(FPCore (ux uy maxCos)
:precision binary32
(if (<= (* uy 2.0) 0.0008299999753944576)
(*
(* 2.0 (* uy (PI)))
(sqrt
(*
(* ux ux)
(+
(/ (+ 2.0 (* maxCos -2.0)) ux)
(* maxCos (* (- 1.0 maxCos) (- (/ -1.0 maxCos) -1.0)))))))
(* (sin (* (* uy 2.0) (PI))) (* ux (sqrt (+ -1.0 (/ 2.0 ux)))))))\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;uy \cdot 2 \leq 0.0008299999753944576:\\
\;\;\;\;\left(2 \cdot \left(uy \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sqrt{\left(ux \cdot ux\right) \cdot \left(\frac{2 + maxCos \cdot -2}{ux} + maxCos \cdot \left(\left(1 - maxCos\right) \cdot \left(\frac{-1}{maxCos} - -1\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\sin \left(\left(uy \cdot 2\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(ux \cdot \sqrt{-1 + \frac{2}{ux}}\right)\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 8.29999975e-4Initial program 58.6%
Taylor expanded in maxCos around inf
lower-*.f32N/A
associate--l+N/A
div-subN/A
lower-+.f32N/A
lower-/.f32N/A
lower--.f3258.7
Applied rewrites58.7%
Taylor expanded in ux around -inf
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower--.f32N/A
lower-/.f32N/A
lower-*.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-+.f32N/A
mul-1-negN/A
lower-neg.f32N/A
sub-negN/A
Applied rewrites97.9%
lift-/.f32N/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
lift-/.f32N/A
div-invN/A
lift-/.f32N/A
associate-*l*N/A
lift-/.f32N/A
inv-powN/A
pow-plusN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
lower-+.f32N/A
lower-*.f3298.4
Applied rewrites98.4%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3297.8
Applied rewrites97.8%
if 8.29999975e-4 < (*.f32 uy #s(literal 2 binary32)) Initial program 55.8%
Taylor expanded in maxCos around inf
lower-*.f32N/A
associate--l+N/A
div-subN/A
lower-+.f32N/A
lower-/.f32N/A
lower--.f3256.2
Applied rewrites56.2%
Taylor expanded in ux around -inf
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower--.f32N/A
lower-/.f32N/A
lower-*.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-+.f32N/A
mul-1-negN/A
lower-neg.f32N/A
sub-negN/A
Applied rewrites98.2%
Taylor expanded in maxCos around 0
lower-*.f32N/A
lower-sqrt.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f3293.8
Applied rewrites93.8%
Final simplification96.5%
(FPCore (ux uy maxCos)
:precision binary32
(let* ((t_0 (* 2.0 (* uy (PI)))))
(if (<= (* uy 2.0) 0.006000000052154064)
(*
t_0
(sqrt
(*
(* ux ux)
(+
(/ (+ 2.0 (* maxCos -2.0)) ux)
(* maxCos (* (- 1.0 maxCos) (- (/ -1.0 maxCos) -1.0)))))))
(* (sin t_0) (sqrt (* ux (- 2.0 maxCos)))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := 2 \cdot \left(uy \cdot \mathsf{PI}\left(\right)\right)\\
\mathbf{if}\;uy \cdot 2 \leq 0.006000000052154064:\\
\;\;\;\;t\_0 \cdot \sqrt{\left(ux \cdot ux\right) \cdot \left(\frac{2 + maxCos \cdot -2}{ux} + maxCos \cdot \left(\left(1 - maxCos\right) \cdot \left(\frac{-1}{maxCos} - -1\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\sin t\_0 \cdot \sqrt{ux \cdot \left(2 - maxCos\right)}\\
\end{array}
\end{array}
if (*.f32 uy #s(literal 2 binary32)) < 0.00600000005Initial program 59.5%
Taylor expanded in maxCos around inf
lower-*.f32N/A
associate--l+N/A
div-subN/A
lower-+.f32N/A
lower-/.f32N/A
lower--.f3259.8
Applied rewrites59.8%
Taylor expanded in ux around -inf
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower--.f32N/A
lower-/.f32N/A
lower-*.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-+.f32N/A
mul-1-negN/A
lower-neg.f32N/A
sub-negN/A
Applied rewrites98.0%
lift-/.f32N/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
lift-/.f32N/A
div-invN/A
lift-/.f32N/A
associate-*l*N/A
lift-/.f32N/A
inv-powN/A
pow-plusN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
lower-+.f32N/A
lower-*.f3298.3
Applied rewrites98.3%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3295.9
Applied rewrites95.9%
if 0.00600000005 < (*.f32 uy #s(literal 2 binary32)) Initial program 52.4%
Taylor expanded in maxCos around 0
lower--.f3251.4
Applied rewrites51.4%
Taylor expanded in ux around 0
lower-*.f32N/A
lower-sqrt.f32N/A
sub-negN/A
mul-1-negN/A
lower-*.f32N/A
mul-1-negN/A
sub-negN/A
lower--.f32N/A
lower-sin.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3278.9
Applied rewrites78.9%
Final simplification91.5%
(FPCore (ux uy maxCos)
:precision binary32
(*
(* 2.0 (* uy (PI)))
(sqrt
(*
(* ux ux)
(+
(/ (+ 2.0 (* maxCos -2.0)) ux)
(* maxCos (* (- 1.0 maxCos) (- (/ -1.0 maxCos) -1.0))))))))\begin{array}{l}
\\
\left(2 \cdot \left(uy \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sqrt{\left(ux \cdot ux\right) \cdot \left(\frac{2 + maxCos \cdot -2}{ux} + maxCos \cdot \left(\left(1 - maxCos\right) \cdot \left(\frac{-1}{maxCos} - -1\right)\right)\right)}
\end{array}
Initial program 57.7%
Taylor expanded in maxCos around inf
lower-*.f32N/A
associate--l+N/A
div-subN/A
lower-+.f32N/A
lower-/.f32N/A
lower--.f3257.9
Applied rewrites57.9%
Taylor expanded in ux around -inf
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower--.f32N/A
lower-/.f32N/A
lower-*.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-+.f32N/A
mul-1-negN/A
lower-neg.f32N/A
sub-negN/A
Applied rewrites98.0%
lift-/.f32N/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
lift-/.f32N/A
div-invN/A
lift-/.f32N/A
associate-*l*N/A
lift-/.f32N/A
inv-powN/A
pow-plusN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
lower-+.f32N/A
lower-*.f3298.2
Applied rewrites98.2%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3283.2
Applied rewrites83.2%
Final simplification83.2%
(FPCore (ux uy maxCos) :precision binary32 (* (* 2.0 (* uy (PI))) (sqrt (* (* ux ux) (+ -1.0 (/ (* maxCos (+ -2.0 (/ 2.0 maxCos))) ux))))))
\begin{array}{l}
\\
\left(2 \cdot \left(uy \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sqrt{\left(ux \cdot ux\right) \cdot \left(-1 + \frac{maxCos \cdot \left(-2 + \frac{2}{maxCos}\right)}{ux}\right)}
\end{array}
Initial program 57.7%
Taylor expanded in maxCos around inf
lower-*.f32N/A
associate--l+N/A
div-subN/A
lower-+.f32N/A
lower-/.f32N/A
lower--.f3257.9
Applied rewrites57.9%
Taylor expanded in ux around -inf
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lower--.f32N/A
lower-/.f32N/A
lower-*.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-+.f32N/A
mul-1-negN/A
lower-neg.f32N/A
sub-negN/A
Applied rewrites98.0%
Taylor expanded in maxCos around 0
Applied rewrites96.3%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f3281.5
Applied rewrites81.5%
Final simplification81.5%
(FPCore (ux uy maxCos) :precision binary32 (* (* 2.0 (* uy (PI))) (sqrt (* (* ux ux) (+ -1.0 (/ 2.0 ux))))))
\begin{array}{l}
\\
\left(2 \cdot \left(uy \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sqrt{\left(ux \cdot ux\right) \cdot \left(-1 + \frac{2}{ux}\right)}
\end{array}
Initial program 57.7%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-sqrt.f32N/A
sub-negN/A
+-commutativeN/A
unpow2N/A
distribute-rgt-neg-inN/A
lower-fma.f32N/A
Applied rewrites31.2%
lift--.f32N/A
lift-*.f32N/A
+-commutativeN/A
lift-+.f3220.5
lift-*.f32N/A
*-commutativeN/A
lower-*.f3220.5
Applied rewrites20.5%
Taylor expanded in maxCos around 0
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-sqrt.f32N/A
+-commutativeN/A
lower-fma.f32N/A
lower--.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f3220.4
Applied rewrites20.4%
Taylor expanded in ux around inf
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f3278.6
Applied rewrites78.6%
Final simplification78.6%
(FPCore (ux uy maxCos) :precision binary32 (* (* 2.0 (* uy (PI))) (sqrt (* ux (- 2.0 ux)))))
\begin{array}{l}
\\
\left(2 \cdot \left(uy \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sqrt{ux \cdot \left(2 - ux\right)}
\end{array}
Initial program 57.7%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-sqrt.f32N/A
sub-negN/A
+-commutativeN/A
unpow2N/A
distribute-rgt-neg-inN/A
lower-fma.f32N/A
Applied rewrites31.2%
lift--.f32N/A
lift-*.f32N/A
+-commutativeN/A
lift-+.f3220.5
lift-*.f32N/A
*-commutativeN/A
lower-*.f3220.5
Applied rewrites20.5%
Taylor expanded in maxCos around 0
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-sqrt.f32N/A
+-commutativeN/A
lower-fma.f32N/A
lower--.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f3220.4
Applied rewrites20.4%
Taylor expanded in ux around 0
lower-*.f32N/A
mul-1-negN/A
unsub-negN/A
lower--.f3278.6
Applied rewrites78.6%
(FPCore (ux uy maxCos) :precision binary32 (* (* 2.0 (* uy (PI))) (sqrt (* 2.0 ux))))
\begin{array}{l}
\\
\left(2 \cdot \left(uy \cdot \mathsf{PI}\left(\right)\right)\right) \cdot \sqrt{2 \cdot ux}
\end{array}
Initial program 57.7%
Taylor expanded in uy around 0
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-sqrt.f32N/A
sub-negN/A
+-commutativeN/A
unpow2N/A
distribute-rgt-neg-inN/A
lower-fma.f32N/A
Applied rewrites31.2%
lift--.f32N/A
lift-*.f32N/A
+-commutativeN/A
lift-+.f3220.5
lift-*.f32N/A
*-commutativeN/A
lower-*.f3220.5
Applied rewrites20.5%
Taylor expanded in maxCos around 0
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-sqrt.f32N/A
+-commutativeN/A
lower-fma.f32N/A
lower--.f32N/A
sub-negN/A
metadata-evalN/A
lower-+.f3220.4
Applied rewrites20.4%
Taylor expanded in ux around 0
lower-*.f3264.1
Applied rewrites64.1%
herbie shell --seed 2024219
(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)))))))