
(FPCore (cosTheta_i u1 u2) :precision binary32 (* (sqrt (- (log (- 1.0 u1)))) (cos (* (* 2.0 PI) u2))))
float code(float cosTheta_i, float u1, float u2) {
return sqrtf(-logf((1.0f - u1))) * cosf(((2.0f * ((float) M_PI)) * u2));
}
function code(cosTheta_i, u1, u2) return Float32(sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) end
function tmp = code(cosTheta_i, u1, u2) tmp = sqrt(-log((single(1.0) - u1))) * cos(((single(2.0) * single(pi)) * u2)); end
\begin{array}{l}
\\
\sqrt{-\log \left(1 - u1\right)} \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right)
\end{array}
Herbie found 22 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (cosTheta_i u1 u2) :precision binary32 (* (sqrt (- (log (- 1.0 u1)))) (cos (* (* 2.0 PI) u2))))
float code(float cosTheta_i, float u1, float u2) {
return sqrtf(-logf((1.0f - u1))) * cosf(((2.0f * ((float) M_PI)) * u2));
}
function code(cosTheta_i, u1, u2) return Float32(sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) end
function tmp = code(cosTheta_i, u1, u2) tmp = sqrt(-log((single(1.0) - u1))) * cos(((single(2.0) * single(pi)) * u2)); end
\begin{array}{l}
\\
\sqrt{-\log \left(1 - u1\right)} \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right)
\end{array}
(FPCore (cosTheta_i u1 u2)
:precision binary32
(if (<= u1 0.01899999938905239)
(*
(fma
(fma
(/ (fma 0.5 (* (- 0.25 (/ 0.0625 u1)) u1) 0.16666666666666666) (sqrt u1))
u1
(/ 0.25 (sqrt u1)))
(* u1 u1)
(sqrt u1))
(sin (fma (+ PI PI) (- u2) (/ PI 2.0))))
(/ (cos (* (+ PI PI) u2)) (pow (- (log (- 1.0 u1))) -0.5))))
float code(float cosTheta_i, float u1, float u2) {
float tmp;
if (u1 <= 0.01899999938905239f) {
tmp = fmaf(fmaf((fmaf(0.5f, ((0.25f - (0.0625f / u1)) * u1), 0.16666666666666666f) / sqrtf(u1)), u1, (0.25f / sqrtf(u1))), (u1 * u1), sqrtf(u1)) * sinf(fmaf((((float) M_PI) + ((float) M_PI)), -u2, (((float) M_PI) / 2.0f)));
} else {
tmp = cosf(((((float) M_PI) + ((float) M_PI)) * u2)) / powf(-logf((1.0f - u1)), -0.5f);
}
return tmp;
}
function code(cosTheta_i, u1, u2) tmp = Float32(0.0) if (u1 <= Float32(0.01899999938905239)) tmp = Float32(fma(fma(Float32(fma(Float32(0.5), Float32(Float32(Float32(0.25) - Float32(Float32(0.0625) / u1)) * u1), Float32(0.16666666666666666)) / sqrt(u1)), u1, Float32(Float32(0.25) / sqrt(u1))), Float32(u1 * u1), sqrt(u1)) * sin(fma(Float32(Float32(pi) + Float32(pi)), Float32(-u2), Float32(Float32(pi) / Float32(2.0))))); else tmp = Float32(cos(Float32(Float32(Float32(pi) + Float32(pi)) * u2)) / (Float32(-log(Float32(Float32(1.0) - u1))) ^ Float32(-0.5))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;u1 \leq 0.01899999938905239:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\frac{\mathsf{fma}\left(0.5, \left(0.25 - \frac{0.0625}{u1}\right) \cdot u1, 0.16666666666666666\right)}{\sqrt{u1}}, u1, \frac{0.25}{\sqrt{u1}}\right), u1 \cdot u1, \sqrt{u1}\right) \cdot \sin \left(\mathsf{fma}\left(\pi + \pi, -u2, \frac{\pi}{2}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\cos \left(\left(\pi + \pi\right) \cdot u2\right)}{{\left(-\log \left(1 - u1\right)\right)}^{-0.5}}\\
\end{array}
\end{array}
if u1 < 0.0189999994Initial program 57.5%
Taylor expanded in u1 around 0
Applied rewrites93.3%
lift-cos.f32N/A
lift-*.f32N/A
lift-PI.f32N/A
lift-*.f32N/A
cos-neg-revN/A
associate-*l*N/A
*-commutativeN/A
sin-+PI/2-revN/A
lower-sin.f32N/A
lift-/.f32N/A
lift-PI.f32N/A
lower-+.f32N/A
Applied rewrites93.4%
lift-+.f32N/A
lift-neg.f32N/A
lift-*.f32N/A
distribute-rgt-neg-inN/A
lift-PI.f32N/A
lift-PI.f32N/A
lift-+.f32N/A
count-2-revN/A
lower-fma.f32N/A
count-2-revN/A
lift-+.f32N/A
lift-PI.f32N/A
lift-PI.f32N/A
lower-neg.f3293.5
Applied rewrites93.5%
if 0.0189999994 < u1 Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.1
Applied rewrites88.1%
lift-sqrt.f32N/A
pow1/2N/A
metadata-evalN/A
pow-negN/A
lower-/.f32N/A
lower-pow.f3287.8
Applied rewrites87.8%
Taylor expanded in u2 around inf
pow-flipN/A
metadata-evalN/A
pow1/2N/A
lower-/.f32N/A
lower-cos.f32N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f32N/A
count-2-revN/A
lift-+.f32N/A
lift-PI.f32N/A
lift-PI.f32N/A
lower-pow.f32N/A
Applied rewrites57.5%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(if (<= u1 0.012000000104308128)
(*
(fma
(fma (/ 0.13541666666666666 (sqrt u1)) u1 (/ 0.25 (sqrt u1)))
(* u1 u1)
(sqrt u1))
(sin (fma (+ PI PI) (- u2) (/ PI 2.0))))
(/ (cos (* (+ PI PI) u2)) (pow (- (log (- 1.0 u1))) -0.5))))
float code(float cosTheta_i, float u1, float u2) {
float tmp;
if (u1 <= 0.012000000104308128f) {
tmp = fmaf(fmaf((0.13541666666666666f / sqrtf(u1)), u1, (0.25f / sqrtf(u1))), (u1 * u1), sqrtf(u1)) * sinf(fmaf((((float) M_PI) + ((float) M_PI)), -u2, (((float) M_PI) / 2.0f)));
} else {
tmp = cosf(((((float) M_PI) + ((float) M_PI)) * u2)) / powf(-logf((1.0f - u1)), -0.5f);
}
return tmp;
}
function code(cosTheta_i, u1, u2) tmp = Float32(0.0) if (u1 <= Float32(0.012000000104308128)) tmp = Float32(fma(fma(Float32(Float32(0.13541666666666666) / sqrt(u1)), u1, Float32(Float32(0.25) / sqrt(u1))), Float32(u1 * u1), sqrt(u1)) * sin(fma(Float32(Float32(pi) + Float32(pi)), Float32(-u2), Float32(Float32(pi) / Float32(2.0))))); else tmp = Float32(cos(Float32(Float32(Float32(pi) + Float32(pi)) * u2)) / (Float32(-log(Float32(Float32(1.0) - u1))) ^ Float32(-0.5))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;u1 \leq 0.012000000104308128:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\frac{0.13541666666666666}{\sqrt{u1}}, u1, \frac{0.25}{\sqrt{u1}}\right), u1 \cdot u1, \sqrt{u1}\right) \cdot \sin \left(\mathsf{fma}\left(\pi + \pi, -u2, \frac{\pi}{2}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\cos \left(\left(\pi + \pi\right) \cdot u2\right)}{{\left(-\log \left(1 - u1\right)\right)}^{-0.5}}\\
\end{array}
\end{array}
if u1 < 0.0120000001Initial program 57.5%
Taylor expanded in u1 around 0
Applied rewrites93.3%
lift-cos.f32N/A
lift-*.f32N/A
lift-PI.f32N/A
lift-*.f32N/A
cos-neg-revN/A
associate-*l*N/A
*-commutativeN/A
sin-+PI/2-revN/A
lower-sin.f32N/A
lift-/.f32N/A
lift-PI.f32N/A
lower-+.f32N/A
Applied rewrites93.4%
lift-+.f32N/A
lift-neg.f32N/A
lift-*.f32N/A
distribute-rgt-neg-inN/A
lift-PI.f32N/A
lift-PI.f32N/A
lift-+.f32N/A
count-2-revN/A
lower-fma.f32N/A
count-2-revN/A
lift-+.f32N/A
lift-PI.f32N/A
lift-PI.f32N/A
lower-neg.f3293.5
Applied rewrites93.5%
Taylor expanded in u1 around 0
Applied rewrites92.1%
if 0.0120000001 < u1 Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.1
Applied rewrites88.1%
lift-sqrt.f32N/A
pow1/2N/A
metadata-evalN/A
pow-negN/A
lower-/.f32N/A
lower-pow.f3287.8
Applied rewrites87.8%
Taylor expanded in u2 around inf
pow-flipN/A
metadata-evalN/A
pow1/2N/A
lower-/.f32N/A
lower-cos.f32N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f32N/A
count-2-revN/A
lift-+.f32N/A
lift-PI.f32N/A
lift-PI.f32N/A
lower-pow.f32N/A
Applied rewrites57.5%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(if (<= u1 0.029999999329447746)
(*
(sqrt (* (fma (fma (fma 0.25 u1 0.3333333333333333) u1 0.5) u1 1.0) u1))
(cos (* (* 2.0 PI) u2)))
(/ (cos (* (+ PI PI) u2)) (pow (- (log (- 1.0 u1))) -0.5))))
float code(float cosTheta_i, float u1, float u2) {
float tmp;
if (u1 <= 0.029999999329447746f) {
tmp = sqrtf((fmaf(fmaf(fmaf(0.25f, u1, 0.3333333333333333f), u1, 0.5f), u1, 1.0f) * u1)) * cosf(((2.0f * ((float) M_PI)) * u2));
} else {
tmp = cosf(((((float) M_PI) + ((float) M_PI)) * u2)) / powf(-logf((1.0f - u1)), -0.5f);
}
return tmp;
}
function code(cosTheta_i, u1, u2) tmp = Float32(0.0) if (u1 <= Float32(0.029999999329447746)) tmp = Float32(sqrt(Float32(fma(fma(fma(Float32(0.25), u1, Float32(0.3333333333333333)), u1, Float32(0.5)), u1, Float32(1.0)) * u1)) * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))); else tmp = Float32(cos(Float32(Float32(Float32(pi) + Float32(pi)) * u2)) / (Float32(-log(Float32(Float32(1.0) - u1))) ^ Float32(-0.5))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;u1 \leq 0.029999999329447746:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.25, u1, 0.3333333333333333\right), u1, 0.5\right), u1, 1\right) \cdot u1} \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\cos \left(\left(\pi + \pi\right) \cdot u2\right)}{{\left(-\log \left(1 - u1\right)\right)}^{-0.5}}\\
\end{array}
\end{array}
if u1 < 0.0299999993Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
+-commutativeN/A
lower-fma.f3293.5
Applied rewrites93.5%
if 0.0299999993 < u1 Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.1
Applied rewrites88.1%
lift-sqrt.f32N/A
pow1/2N/A
metadata-evalN/A
pow-negN/A
lower-/.f32N/A
lower-pow.f3287.8
Applied rewrites87.8%
Taylor expanded in u2 around inf
pow-flipN/A
metadata-evalN/A
pow1/2N/A
lower-/.f32N/A
lower-cos.f32N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f32N/A
count-2-revN/A
lift-+.f32N/A
lift-PI.f32N/A
lift-PI.f32N/A
lower-pow.f32N/A
Applied rewrites57.5%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(if (<= u1 0.029999999329447746)
(*
(sqrt (* (fma (fma (fma 0.25 u1 0.3333333333333333) u1 0.5) u1 1.0) u1))
(cos (* (* 2.0 PI) u2)))
(* (sqrt (- (log (- 1.0 u1)))) (cos (* (+ PI PI) u2)))))
float code(float cosTheta_i, float u1, float u2) {
float tmp;
if (u1 <= 0.029999999329447746f) {
tmp = sqrtf((fmaf(fmaf(fmaf(0.25f, u1, 0.3333333333333333f), u1, 0.5f), u1, 1.0f) * u1)) * cosf(((2.0f * ((float) M_PI)) * u2));
} else {
tmp = sqrtf(-logf((1.0f - u1))) * cosf(((((float) M_PI) + ((float) M_PI)) * u2));
}
return tmp;
}
function code(cosTheta_i, u1, u2) tmp = Float32(0.0) if (u1 <= Float32(0.029999999329447746)) tmp = Float32(sqrt(Float32(fma(fma(fma(Float32(0.25), u1, Float32(0.3333333333333333)), u1, Float32(0.5)), u1, Float32(1.0)) * u1)) * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))); else tmp = Float32(sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) * cos(Float32(Float32(Float32(pi) + Float32(pi)) * u2))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;u1 \leq 0.029999999329447746:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.25, u1, 0.3333333333333333\right), u1, 0.5\right), u1, 1\right) \cdot u1} \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{-\log \left(1 - u1\right)} \cdot \cos \left(\left(\pi + \pi\right) \cdot u2\right)\\
\end{array}
\end{array}
if u1 < 0.0299999993Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
+-commutativeN/A
lower-fma.f3293.5
Applied rewrites93.5%
if 0.0299999993 < u1 Initial program 57.5%
lift-PI.f32N/A
lift-*.f32N/A
count-2-revN/A
lower-+.f32N/A
lift-PI.f32N/A
lift-PI.f3257.5
Applied rewrites57.5%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (cos (* (+ PI PI) u2))))
(if (<= u1 0.012000000104308128)
(* (sqrt (* (fma (fma 0.3333333333333333 u1 0.5) u1 1.0) u1)) t_0)
(* (sqrt (- (log (- 1.0 u1)))) t_0))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = cosf(((((float) M_PI) + ((float) M_PI)) * u2));
float tmp;
if (u1 <= 0.012000000104308128f) {
tmp = sqrtf((fmaf(fmaf(0.3333333333333333f, u1, 0.5f), u1, 1.0f) * u1)) * t_0;
} else {
tmp = sqrtf(-logf((1.0f - u1))) * t_0;
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = cos(Float32(Float32(Float32(pi) + Float32(pi)) * u2)) tmp = Float32(0.0) if (u1 <= Float32(0.012000000104308128)) tmp = Float32(sqrt(Float32(fma(fma(Float32(0.3333333333333333), u1, Float32(0.5)), u1, Float32(1.0)) * u1)) * t_0); else tmp = Float32(sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) * t_0); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\left(\pi + \pi\right) \cdot u2\right)\\
\mathbf{if}\;u1 \leq 0.012000000104308128:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(\mathsf{fma}\left(0.3333333333333333, u1, 0.5\right), u1, 1\right) \cdot u1} \cdot t\_0\\
\mathbf{else}:\\
\;\;\;\;\sqrt{-\log \left(1 - u1\right)} \cdot t\_0\\
\end{array}
\end{array}
if u1 < 0.0120000001Initial program 57.5%
lift-cos.f32N/A
sin-+PI/2-revN/A
lower-sin.f32N/A
lift-*.f32N/A
*-commutativeN/A
lower-fma.f32N/A
lift-PI.f32N/A
lift-*.f32N/A
count-2-revN/A
lower-+.f32N/A
lift-PI.f32N/A
lift-PI.f32N/A
lower-/.f32N/A
lift-PI.f3257.5
Applied rewrites57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
+-commutativeN/A
lower-fma.f3291.6
Applied rewrites91.6%
Applied rewrites91.8%
if 0.0120000001 < u1 Initial program 57.5%
lift-PI.f32N/A
lift-*.f32N/A
count-2-revN/A
lower-+.f32N/A
lift-PI.f32N/A
lift-PI.f3257.5
Applied rewrites57.5%
(FPCore (cosTheta_i u1 u2) :precision binary32 (if (<= u1 0.0026000000070780516) (* (sqrt (fma (* u1 u1) 0.5 u1)) (cos (* (* 2.0 PI) u2))) (* (sqrt (- (log (- 1.0 u1)))) (cos (* (+ PI PI) u2)))))
float code(float cosTheta_i, float u1, float u2) {
float tmp;
if (u1 <= 0.0026000000070780516f) {
tmp = sqrtf(fmaf((u1 * u1), 0.5f, u1)) * cosf(((2.0f * ((float) M_PI)) * u2));
} else {
tmp = sqrtf(-logf((1.0f - u1))) * cosf(((((float) M_PI) + ((float) M_PI)) * u2));
}
return tmp;
}
function code(cosTheta_i, u1, u2) tmp = Float32(0.0) if (u1 <= Float32(0.0026000000070780516)) tmp = Float32(sqrt(fma(Float32(u1 * u1), Float32(0.5), u1)) * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))); else tmp = Float32(sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) * cos(Float32(Float32(Float32(pi) + Float32(pi)) * u2))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;u1 \leq 0.0026000000070780516:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(u1 \cdot u1, 0.5, u1\right)} \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{-\log \left(1 - u1\right)} \cdot \cos \left(\left(\pi + \pi\right) \cdot u2\right)\\
\end{array}
\end{array}
if u1 < 0.00260000001Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.1
Applied rewrites88.1%
Taylor expanded in u1 around inf
distribute-lft-inN/A
inv-powN/A
pow-prod-upN/A
metadata-evalN/A
unpow1N/A
lower-fma.f32N/A
pow2N/A
lift-*.f3288.2
Applied rewrites88.2%
if 0.00260000001 < u1 Initial program 57.5%
lift-PI.f32N/A
lift-*.f32N/A
count-2-revN/A
lower-+.f32N/A
lift-PI.f32N/A
lift-PI.f3257.5
Applied rewrites57.5%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (sqrt (- (log (- 1.0 u1))))) (t_1 (cos (* (* 2.0 PI) u2))))
(if (<= (* t_0 t_1) 0.09399999678134918)
(* (sqrt (fma (* u1 u1) 0.5 u1)) t_1)
(*
t_0
(fma
(fma
(* 0.6666666666666666 (* u2 u2))
(* (* PI PI) (* PI PI))
(* (* PI PI) -2.0))
(* u2 u2)
1.0)))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = sqrtf(-logf((1.0f - u1)));
float t_1 = cosf(((2.0f * ((float) M_PI)) * u2));
float tmp;
if ((t_0 * t_1) <= 0.09399999678134918f) {
tmp = sqrtf(fmaf((u1 * u1), 0.5f, u1)) * t_1;
} else {
tmp = t_0 * fmaf(fmaf((0.6666666666666666f * (u2 * u2)), ((((float) M_PI) * ((float) M_PI)) * (((float) M_PI) * ((float) M_PI))), ((((float) M_PI) * ((float) M_PI)) * -2.0f)), (u2 * u2), 1.0f);
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) t_1 = cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2)) tmp = Float32(0.0) if (Float32(t_0 * t_1) <= Float32(0.09399999678134918)) tmp = Float32(sqrt(fma(Float32(u1 * u1), Float32(0.5), u1)) * t_1); else tmp = Float32(t_0 * fma(fma(Float32(Float32(0.6666666666666666) * Float32(u2 * u2)), Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(pi) * Float32(pi))), Float32(Float32(Float32(pi) * Float32(pi)) * Float32(-2.0))), Float32(u2 * u2), Float32(1.0))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{-\log \left(1 - u1\right)}\\
t_1 := \cos \left(\left(2 \cdot \pi\right) \cdot u2\right)\\
\mathbf{if}\;t\_0 \cdot t\_1 \leq 0.09399999678134918:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(u1 \cdot u1, 0.5, u1\right)} \cdot t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_0 \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.6666666666666666 \cdot \left(u2 \cdot u2\right), \left(\pi \cdot \pi\right) \cdot \left(\pi \cdot \pi\right), \left(\pi \cdot \pi\right) \cdot -2\right), u2 \cdot u2, 1\right)\\
\end{array}
\end{array}
if (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.0939999968Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.1
Applied rewrites88.1%
Taylor expanded in u1 around inf
distribute-lft-inN/A
inv-powN/A
pow-prod-upN/A
metadata-evalN/A
unpow1N/A
lower-fma.f32N/A
pow2N/A
lift-*.f3288.2
Applied rewrites88.2%
if 0.0939999968 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) Initial program 57.5%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites54.3%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (sqrt (- (log (- 1.0 u1))))))
(if (<= (* t_0 (cos (* (* 2.0 PI) u2))) 0.09399999678134918)
(* (sqrt (* (fma 0.5 u1 1.0) u1)) (cos (* (+ PI PI) u2)))
(*
t_0
(fma
(fma
(* 0.6666666666666666 (* u2 u2))
(* (* PI PI) (* PI PI))
(* (* PI PI) -2.0))
(* u2 u2)
1.0)))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = sqrtf(-logf((1.0f - u1)));
float tmp;
if ((t_0 * cosf(((2.0f * ((float) M_PI)) * u2))) <= 0.09399999678134918f) {
tmp = sqrtf((fmaf(0.5f, u1, 1.0f) * u1)) * cosf(((((float) M_PI) + ((float) M_PI)) * u2));
} else {
tmp = t_0 * fmaf(fmaf((0.6666666666666666f * (u2 * u2)), ((((float) M_PI) * ((float) M_PI)) * (((float) M_PI) * ((float) M_PI))), ((((float) M_PI) * ((float) M_PI)) * -2.0f)), (u2 * u2), 1.0f);
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) tmp = Float32(0.0) if (Float32(t_0 * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) <= Float32(0.09399999678134918)) tmp = Float32(sqrt(Float32(fma(Float32(0.5), u1, Float32(1.0)) * u1)) * cos(Float32(Float32(Float32(pi) + Float32(pi)) * u2))); else tmp = Float32(t_0 * fma(fma(Float32(Float32(0.6666666666666666) * Float32(u2 * u2)), Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(pi) * Float32(pi))), Float32(Float32(Float32(pi) * Float32(pi)) * Float32(-2.0))), Float32(u2 * u2), Float32(1.0))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{-\log \left(1 - u1\right)}\\
\mathbf{if}\;t\_0 \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right) \leq 0.09399999678134918:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(0.5, u1, 1\right) \cdot u1} \cdot \cos \left(\left(\pi + \pi\right) \cdot u2\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0 \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.6666666666666666 \cdot \left(u2 \cdot u2\right), \left(\pi \cdot \pi\right) \cdot \left(\pi \cdot \pi\right), \left(\pi \cdot \pi\right) \cdot -2\right), u2 \cdot u2, 1\right)\\
\end{array}
\end{array}
if (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.0939999968Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.1
Applied rewrites88.1%
lift-PI.f32N/A
lift-*.f32N/A
count-2-revN/A
lift-+.f32N/A
lift-PI.f32N/A
lift-PI.f3288.1
Applied rewrites88.1%
if 0.0939999968 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) Initial program 57.5%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites54.3%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0
(fma
(fma
(* 0.6666666666666666 (* u2 u2))
(* (* PI PI) (* PI PI))
(* (* PI PI) -2.0))
(* u2 u2)
1.0))
(t_1 (sqrt (- (log (- 1.0 u1)))))
(t_2 (* t_1 (cos (* (* 2.0 PI) u2)))))
(if (<= t_2 0.0019000000320374966)
(* (sqrt u1) (sin (fma (+ PI PI) (- u2) (/ PI 2.0))))
(if (<= t_2 0.052000001072883606)
(* (sqrt (* (fma 0.5 u1 1.0) u1)) t_0)
(* t_1 t_0)))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = fmaf(fmaf((0.6666666666666666f * (u2 * u2)), ((((float) M_PI) * ((float) M_PI)) * (((float) M_PI) * ((float) M_PI))), ((((float) M_PI) * ((float) M_PI)) * -2.0f)), (u2 * u2), 1.0f);
float t_1 = sqrtf(-logf((1.0f - u1)));
float t_2 = t_1 * cosf(((2.0f * ((float) M_PI)) * u2));
float tmp;
if (t_2 <= 0.0019000000320374966f) {
tmp = sqrtf(u1) * sinf(fmaf((((float) M_PI) + ((float) M_PI)), -u2, (((float) M_PI) / 2.0f)));
} else if (t_2 <= 0.052000001072883606f) {
tmp = sqrtf((fmaf(0.5f, u1, 1.0f) * u1)) * t_0;
} else {
tmp = t_1 * t_0;
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = fma(fma(Float32(Float32(0.6666666666666666) * Float32(u2 * u2)), Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(pi) * Float32(pi))), Float32(Float32(Float32(pi) * Float32(pi)) * Float32(-2.0))), Float32(u2 * u2), Float32(1.0)) t_1 = sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) t_2 = Float32(t_1 * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) tmp = Float32(0.0) if (t_2 <= Float32(0.0019000000320374966)) tmp = Float32(sqrt(u1) * sin(fma(Float32(Float32(pi) + Float32(pi)), Float32(-u2), Float32(Float32(pi) / Float32(2.0))))); elseif (t_2 <= Float32(0.052000001072883606)) tmp = Float32(sqrt(Float32(fma(Float32(0.5), u1, Float32(1.0)) * u1)) * t_0); else tmp = Float32(t_1 * t_0); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(\mathsf{fma}\left(0.6666666666666666 \cdot \left(u2 \cdot u2\right), \left(\pi \cdot \pi\right) \cdot \left(\pi \cdot \pi\right), \left(\pi \cdot \pi\right) \cdot -2\right), u2 \cdot u2, 1\right)\\
t_1 := \sqrt{-\log \left(1 - u1\right)}\\
t_2 := t\_1 \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right)\\
\mathbf{if}\;t\_2 \leq 0.0019000000320374966:\\
\;\;\;\;\sqrt{u1} \cdot \sin \left(\mathsf{fma}\left(\pi + \pi, -u2, \frac{\pi}{2}\right)\right)\\
\mathbf{elif}\;t\_2 \leq 0.052000001072883606:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(0.5, u1, 1\right) \cdot u1} \cdot t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot t\_0\\
\end{array}
\end{array}
if (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.00190000003Initial program 57.5%
Taylor expanded in u1 around 0
Applied rewrites93.3%
lift-cos.f32N/A
lift-*.f32N/A
lift-PI.f32N/A
lift-*.f32N/A
cos-neg-revN/A
associate-*l*N/A
*-commutativeN/A
sin-+PI/2-revN/A
lower-sin.f32N/A
lift-/.f32N/A
lift-PI.f32N/A
lower-+.f32N/A
Applied rewrites93.4%
lift-+.f32N/A
lift-neg.f32N/A
lift-*.f32N/A
distribute-rgt-neg-inN/A
lift-PI.f32N/A
lift-PI.f32N/A
lift-+.f32N/A
count-2-revN/A
lower-fma.f32N/A
count-2-revN/A
lift-+.f32N/A
lift-PI.f32N/A
lift-PI.f32N/A
lower-neg.f3293.5
Applied rewrites93.5%
Taylor expanded in u1 around 0
flip3--N/A
metadata-evalN/A
pow3N/A
metadata-evalN/A
diff-logN/A
lift-sqrt.f3276.6
Applied rewrites76.6%
if 0.00190000003 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.0520000011Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.1
Applied rewrites88.1%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites81.9%
if 0.0520000011 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) Initial program 57.5%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites54.3%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0
(fma
(fma
(* 0.6666666666666666 (* u2 u2))
(* (* PI PI) (* PI PI))
(* (* PI PI) -2.0))
(* u2 u2)
1.0))
(t_1 (sqrt (- (log (- 1.0 u1)))))
(t_2 (* t_1 (cos (* (* 2.0 PI) u2)))))
(if (<= t_2 0.0019000000320374966)
(* (sqrt u1) (sin (fma -2.0 (* PI u2) (* PI 0.5))))
(if (<= t_2 0.052000001072883606)
(* (sqrt (* (fma 0.5 u1 1.0) u1)) t_0)
(* t_1 t_0)))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = fmaf(fmaf((0.6666666666666666f * (u2 * u2)), ((((float) M_PI) * ((float) M_PI)) * (((float) M_PI) * ((float) M_PI))), ((((float) M_PI) * ((float) M_PI)) * -2.0f)), (u2 * u2), 1.0f);
float t_1 = sqrtf(-logf((1.0f - u1)));
float t_2 = t_1 * cosf(((2.0f * ((float) M_PI)) * u2));
float tmp;
if (t_2 <= 0.0019000000320374966f) {
tmp = sqrtf(u1) * sinf(fmaf(-2.0f, (((float) M_PI) * u2), (((float) M_PI) * 0.5f)));
} else if (t_2 <= 0.052000001072883606f) {
tmp = sqrtf((fmaf(0.5f, u1, 1.0f) * u1)) * t_0;
} else {
tmp = t_1 * t_0;
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = fma(fma(Float32(Float32(0.6666666666666666) * Float32(u2 * u2)), Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(pi) * Float32(pi))), Float32(Float32(Float32(pi) * Float32(pi)) * Float32(-2.0))), Float32(u2 * u2), Float32(1.0)) t_1 = sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) t_2 = Float32(t_1 * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) tmp = Float32(0.0) if (t_2 <= Float32(0.0019000000320374966)) tmp = Float32(sqrt(u1) * sin(fma(Float32(-2.0), Float32(Float32(pi) * u2), Float32(Float32(pi) * Float32(0.5))))); elseif (t_2 <= Float32(0.052000001072883606)) tmp = Float32(sqrt(Float32(fma(Float32(0.5), u1, Float32(1.0)) * u1)) * t_0); else tmp = Float32(t_1 * t_0); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(\mathsf{fma}\left(0.6666666666666666 \cdot \left(u2 \cdot u2\right), \left(\pi \cdot \pi\right) \cdot \left(\pi \cdot \pi\right), \left(\pi \cdot \pi\right) \cdot -2\right), u2 \cdot u2, 1\right)\\
t_1 := \sqrt{-\log \left(1 - u1\right)}\\
t_2 := t\_1 \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right)\\
\mathbf{if}\;t\_2 \leq 0.0019000000320374966:\\
\;\;\;\;\sqrt{u1} \cdot \sin \left(\mathsf{fma}\left(-2, \pi \cdot u2, \pi \cdot 0.5\right)\right)\\
\mathbf{elif}\;t\_2 \leq 0.052000001072883606:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(0.5, u1, 1\right) \cdot u1} \cdot t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot t\_0\\
\end{array}
\end{array}
if (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.00190000003Initial program 57.5%
Taylor expanded in u1 around 0
Applied rewrites93.3%
lift-cos.f32N/A
lift-*.f32N/A
lift-PI.f32N/A
lift-*.f32N/A
cos-neg-revN/A
associate-*l*N/A
*-commutativeN/A
sin-+PI/2-revN/A
lower-sin.f32N/A
lift-/.f32N/A
lift-PI.f32N/A
lower-+.f32N/A
Applied rewrites93.4%
Taylor expanded in u1 around 0
Applied rewrites76.5%
if 0.00190000003 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.0520000011Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.1
Applied rewrites88.1%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites81.9%
if 0.0520000011 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) Initial program 57.5%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites54.3%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0
(fma
(fma
(* 0.6666666666666666 (* u2 u2))
(* (* PI PI) (* PI PI))
(* (* PI PI) -2.0))
(* u2 u2)
1.0))
(t_1 (sqrt (- (log (- 1.0 u1)))))
(t_2 (* t_1 (cos (* (* 2.0 PI) u2)))))
(if (<= t_2 0.0019000000320374966)
(* (sqrt u1) (cos (* u2 (+ PI PI))))
(if (<= t_2 0.052000001072883606)
(* (sqrt (* (fma 0.5 u1 1.0) u1)) t_0)
(* t_1 t_0)))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = fmaf(fmaf((0.6666666666666666f * (u2 * u2)), ((((float) M_PI) * ((float) M_PI)) * (((float) M_PI) * ((float) M_PI))), ((((float) M_PI) * ((float) M_PI)) * -2.0f)), (u2 * u2), 1.0f);
float t_1 = sqrtf(-logf((1.0f - u1)));
float t_2 = t_1 * cosf(((2.0f * ((float) M_PI)) * u2));
float tmp;
if (t_2 <= 0.0019000000320374966f) {
tmp = sqrtf(u1) * cosf((u2 * (((float) M_PI) + ((float) M_PI))));
} else if (t_2 <= 0.052000001072883606f) {
tmp = sqrtf((fmaf(0.5f, u1, 1.0f) * u1)) * t_0;
} else {
tmp = t_1 * t_0;
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = fma(fma(Float32(Float32(0.6666666666666666) * Float32(u2 * u2)), Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(pi) * Float32(pi))), Float32(Float32(Float32(pi) * Float32(pi)) * Float32(-2.0))), Float32(u2 * u2), Float32(1.0)) t_1 = sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) t_2 = Float32(t_1 * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) tmp = Float32(0.0) if (t_2 <= Float32(0.0019000000320374966)) tmp = Float32(sqrt(u1) * cos(Float32(u2 * Float32(Float32(pi) + Float32(pi))))); elseif (t_2 <= Float32(0.052000001072883606)) tmp = Float32(sqrt(Float32(fma(Float32(0.5), u1, Float32(1.0)) * u1)) * t_0); else tmp = Float32(t_1 * t_0); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(\mathsf{fma}\left(0.6666666666666666 \cdot \left(u2 \cdot u2\right), \left(\pi \cdot \pi\right) \cdot \left(\pi \cdot \pi\right), \left(\pi \cdot \pi\right) \cdot -2\right), u2 \cdot u2, 1\right)\\
t_1 := \sqrt{-\log \left(1 - u1\right)}\\
t_2 := t\_1 \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right)\\
\mathbf{if}\;t\_2 \leq 0.0019000000320374966:\\
\;\;\;\;\sqrt{u1} \cdot \cos \left(u2 \cdot \left(\pi + \pi\right)\right)\\
\mathbf{elif}\;t\_2 \leq 0.052000001072883606:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(0.5, u1, 1\right) \cdot u1} \cdot t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot t\_0\\
\end{array}
\end{array}
if (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.00190000003Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
lower-sqrt.f32N/A
*-commutativeN/A
associate-*l*N/A
lift-*.f32N/A
lift-PI.f32N/A
lift-*.f32N/A
lift-cos.f3276.5
lift-*.f32N/A
*-commutativeN/A
lower-*.f3276.5
lift-PI.f32N/A
lift-*.f32N/A
count-2-revN/A
lower-+.f32N/A
lift-PI.f32N/A
lift-PI.f3276.5
Applied rewrites76.5%
if 0.00190000003 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.0520000011Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.1
Applied rewrites88.1%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites81.9%
if 0.0520000011 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) Initial program 57.5%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites54.3%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (- (log (- 1.0 u1))))
(t_1 (* (sqrt t_0) (cos (* (* 2.0 PI) u2)))))
(if (<= t_1 0.0019000000320374966)
(* (sqrt u1) (cos (* u2 (+ PI PI))))
(if (<= t_1 0.09399999678134918)
(*
(sqrt (* (fma 0.5 u1 1.0) u1))
(fma
(fma
(* 0.6666666666666666 (* u2 u2))
(* (* PI PI) (* PI PI))
(* (* PI PI) -2.0))
(* u2 u2)
1.0))
(* (exp (* (log t_0) 0.5)) (fma (* -2.0 (* u2 u2)) (* PI PI) 1.0))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = -logf((1.0f - u1));
float t_1 = sqrtf(t_0) * cosf(((2.0f * ((float) M_PI)) * u2));
float tmp;
if (t_1 <= 0.0019000000320374966f) {
tmp = sqrtf(u1) * cosf((u2 * (((float) M_PI) + ((float) M_PI))));
} else if (t_1 <= 0.09399999678134918f) {
tmp = sqrtf((fmaf(0.5f, u1, 1.0f) * u1)) * fmaf(fmaf((0.6666666666666666f * (u2 * u2)), ((((float) M_PI) * ((float) M_PI)) * (((float) M_PI) * ((float) M_PI))), ((((float) M_PI) * ((float) M_PI)) * -2.0f)), (u2 * u2), 1.0f);
} else {
tmp = expf((logf(t_0) * 0.5f)) * fmaf((-2.0f * (u2 * u2)), (((float) M_PI) * ((float) M_PI)), 1.0f);
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(-log(Float32(Float32(1.0) - u1))) t_1 = Float32(sqrt(t_0) * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) tmp = Float32(0.0) if (t_1 <= Float32(0.0019000000320374966)) tmp = Float32(sqrt(u1) * cos(Float32(u2 * Float32(Float32(pi) + Float32(pi))))); elseif (t_1 <= Float32(0.09399999678134918)) tmp = Float32(sqrt(Float32(fma(Float32(0.5), u1, Float32(1.0)) * u1)) * fma(fma(Float32(Float32(0.6666666666666666) * Float32(u2 * u2)), Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(pi) * Float32(pi))), Float32(Float32(Float32(pi) * Float32(pi)) * Float32(-2.0))), Float32(u2 * u2), Float32(1.0))); else tmp = Float32(exp(Float32(log(t_0) * Float32(0.5))) * fma(Float32(Float32(-2.0) * Float32(u2 * u2)), Float32(Float32(pi) * Float32(pi)), Float32(1.0))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -\log \left(1 - u1\right)\\
t_1 := \sqrt{t\_0} \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right)\\
\mathbf{if}\;t\_1 \leq 0.0019000000320374966:\\
\;\;\;\;\sqrt{u1} \cdot \cos \left(u2 \cdot \left(\pi + \pi\right)\right)\\
\mathbf{elif}\;t\_1 \leq 0.09399999678134918:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(0.5, u1, 1\right) \cdot u1} \cdot \mathsf{fma}\left(\mathsf{fma}\left(0.6666666666666666 \cdot \left(u2 \cdot u2\right), \left(\pi \cdot \pi\right) \cdot \left(\pi \cdot \pi\right), \left(\pi \cdot \pi\right) \cdot -2\right), u2 \cdot u2, 1\right)\\
\mathbf{else}:\\
\;\;\;\;e^{\log t\_0 \cdot 0.5} \cdot \mathsf{fma}\left(-2 \cdot \left(u2 \cdot u2\right), \pi \cdot \pi, 1\right)\\
\end{array}
\end{array}
if (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.00190000003Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
lower-sqrt.f32N/A
*-commutativeN/A
associate-*l*N/A
lift-*.f32N/A
lift-PI.f32N/A
lift-*.f32N/A
lift-cos.f3276.5
lift-*.f32N/A
*-commutativeN/A
lower-*.f3276.5
lift-PI.f32N/A
lift-*.f32N/A
count-2-revN/A
lower-+.f32N/A
lift-PI.f32N/A
lift-PI.f3276.5
Applied rewrites76.5%
if 0.00190000003 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.0939999968Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.1
Applied rewrites88.1%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites81.9%
if 0.0939999968 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) Initial program 57.5%
lift-sqrt.f32N/A
lift-neg.f32N/A
lift--.f32N/A
lift-log.f32N/A
pow1/2N/A
pow-to-expN/A
lower-exp.f32N/A
lower-*.f32N/A
lower-log.f32N/A
lift-log.f32N/A
lift--.f32N/A
lift-neg.f3257.5
Applied rewrites57.5%
Taylor expanded in u2 around 0
+-commutativeN/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lift-PI.f32N/A
lift-PI.f3253.1
Applied rewrites53.1%
(FPCore (cosTheta_i u1 u2) :precision binary32 (if (<= (cos (* (* 2.0 PI) u2)) 0.9999849796295166) (* (sqrt u1) (cos (* u2 (+ PI PI)))) (/ 1.0 (pow (- (log1p (- u1))) -0.5))))
float code(float cosTheta_i, float u1, float u2) {
float tmp;
if (cosf(((2.0f * ((float) M_PI)) * u2)) <= 0.9999849796295166f) {
tmp = sqrtf(u1) * cosf((u2 * (((float) M_PI) + ((float) M_PI))));
} else {
tmp = 1.0f / powf(-log1pf(-u1), -0.5f);
}
return tmp;
}
function code(cosTheta_i, u1, u2) tmp = Float32(0.0) if (cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2)) <= Float32(0.9999849796295166)) tmp = Float32(sqrt(u1) * cos(Float32(u2 * Float32(Float32(pi) + Float32(pi))))); else tmp = Float32(Float32(1.0) / (Float32(-log1p(Float32(-u1))) ^ Float32(-0.5))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\left(2 \cdot \pi\right) \cdot u2\right) \leq 0.9999849796295166:\\
\;\;\;\;\sqrt{u1} \cdot \cos \left(u2 \cdot \left(\pi + \pi\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{{\left(-\mathsf{log1p}\left(-u1\right)\right)}^{-0.5}}\\
\end{array}
\end{array}
if (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2)) < 0.99998498Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
lower-sqrt.f32N/A
*-commutativeN/A
associate-*l*N/A
lift-*.f32N/A
lift-PI.f32N/A
lift-*.f32N/A
lift-cos.f3276.5
lift-*.f32N/A
*-commutativeN/A
lower-*.f3276.5
lift-PI.f32N/A
lift-*.f32N/A
count-2-revN/A
lower-+.f32N/A
lift-PI.f32N/A
lift-PI.f3276.5
Applied rewrites76.5%
if 0.99998498 < (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2)) Initial program 57.5%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.8
Applied rewrites49.8%
lift-sqrt.f32N/A
lift-neg.f32N/A
lift--.f32N/A
lift-log.f32N/A
pow1/2N/A
metadata-evalN/A
pow-negN/A
lower-/.f32N/A
lower-pow.f32N/A
lift-log.f32N/A
lift--.f32N/A
lift-neg.f3249.8
Applied rewrites49.8%
lift--.f32N/A
lift-log.f32N/A
*-lft-identityN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
mul-1-negN/A
lower-log1p.f32N/A
lower-neg.f3280.0
Applied rewrites80.0%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (sqrt (- (log (- 1.0 u1)))))
(t_1 (fma (* -2.0 (* u2 u2)) (* PI PI) 1.0)))
(if (<= (* t_0 (cos (* (* 2.0 PI) u2))) 0.052000001072883606)
(* (sqrt (* (* (+ (/ 1.0 u1) 0.5) u1) u1)) t_1)
(* t_0 t_1))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = sqrtf(-logf((1.0f - u1)));
float t_1 = fmaf((-2.0f * (u2 * u2)), (((float) M_PI) * ((float) M_PI)), 1.0f);
float tmp;
if ((t_0 * cosf(((2.0f * ((float) M_PI)) * u2))) <= 0.052000001072883606f) {
tmp = sqrtf(((((1.0f / u1) + 0.5f) * u1) * u1)) * t_1;
} else {
tmp = t_0 * t_1;
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) t_1 = fma(Float32(Float32(-2.0) * Float32(u2 * u2)), Float32(Float32(pi) * Float32(pi)), Float32(1.0)) tmp = Float32(0.0) if (Float32(t_0 * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) <= Float32(0.052000001072883606)) tmp = Float32(sqrt(Float32(Float32(Float32(Float32(Float32(1.0) / u1) + Float32(0.5)) * u1) * u1)) * t_1); else tmp = Float32(t_0 * t_1); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{-\log \left(1 - u1\right)}\\
t_1 := \mathsf{fma}\left(-2 \cdot \left(u2 \cdot u2\right), \pi \cdot \pi, 1\right)\\
\mathbf{if}\;t\_0 \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right) \leq 0.052000001072883606:\\
\;\;\;\;\sqrt{\left(\left(\frac{1}{u1} + 0.5\right) \cdot u1\right) \cdot u1} \cdot t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_0 \cdot t\_1\\
\end{array}
\end{array}
if (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.0520000011Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.1
Applied rewrites88.1%
lift-fma.f32N/A
metadata-evalN/A
metadata-evalN/A
pow-plusN/A
inv-powN/A
distribute-rgt-inN/A
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-+.f32N/A
lower-/.f3288.0
Applied rewrites88.0%
Taylor expanded in u2 around 0
+-commutativeN/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lift-PI.f32N/A
lift-PI.f3279.2
Applied rewrites79.2%
if 0.0520000011 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) Initial program 57.5%
Taylor expanded in u2 around 0
+-commutativeN/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lift-PI.f32N/A
lift-PI.f3253.1
Applied rewrites53.1%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (sqrt (- (log (- 1.0 u1)))))
(t_1 (fma (* -2.0 (* u2 u2)) (* PI PI) 1.0)))
(if (<= (* t_0 (cos (* (* 2.0 PI) u2))) 0.052000001072883606)
(* (sqrt (* (fma 0.5 u1 1.0) u1)) t_1)
(* t_0 t_1))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = sqrtf(-logf((1.0f - u1)));
float t_1 = fmaf((-2.0f * (u2 * u2)), (((float) M_PI) * ((float) M_PI)), 1.0f);
float tmp;
if ((t_0 * cosf(((2.0f * ((float) M_PI)) * u2))) <= 0.052000001072883606f) {
tmp = sqrtf((fmaf(0.5f, u1, 1.0f) * u1)) * t_1;
} else {
tmp = t_0 * t_1;
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) t_1 = fma(Float32(Float32(-2.0) * Float32(u2 * u2)), Float32(Float32(pi) * Float32(pi)), Float32(1.0)) tmp = Float32(0.0) if (Float32(t_0 * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) <= Float32(0.052000001072883606)) tmp = Float32(sqrt(Float32(fma(Float32(0.5), u1, Float32(1.0)) * u1)) * t_1); else tmp = Float32(t_0 * t_1); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{-\log \left(1 - u1\right)}\\
t_1 := \mathsf{fma}\left(-2 \cdot \left(u2 \cdot u2\right), \pi \cdot \pi, 1\right)\\
\mathbf{if}\;t\_0 \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right) \leq 0.052000001072883606:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(0.5, u1, 1\right) \cdot u1} \cdot t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_0 \cdot t\_1\\
\end{array}
\end{array}
if (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.0520000011Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.1
Applied rewrites88.1%
Taylor expanded in u2 around 0
+-commutativeN/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lift-PI.f32N/A
lift-PI.f3279.3
Applied rewrites79.3%
if 0.0520000011 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) Initial program 57.5%
Taylor expanded in u2 around 0
+-commutativeN/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lift-PI.f32N/A
lift-PI.f3253.1
Applied rewrites53.1%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (- (log (- 1.0 u1)))))
(if (<= (* (sqrt t_0) (cos (* (* 2.0 PI) u2))) 0.07599999755620956)
(* (sqrt (* (fma 0.5 u1 1.0) u1)) (fma (* -2.0 (* u2 u2)) (* PI PI) 1.0))
(sqrt (exp (log t_0))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = -logf((1.0f - u1));
float tmp;
if ((sqrtf(t_0) * cosf(((2.0f * ((float) M_PI)) * u2))) <= 0.07599999755620956f) {
tmp = sqrtf((fmaf(0.5f, u1, 1.0f) * u1)) * fmaf((-2.0f * (u2 * u2)), (((float) M_PI) * ((float) M_PI)), 1.0f);
} else {
tmp = sqrtf(expf(logf(t_0)));
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(-log(Float32(Float32(1.0) - u1))) tmp = Float32(0.0) if (Float32(sqrt(t_0) * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) <= Float32(0.07599999755620956)) tmp = Float32(sqrt(Float32(fma(Float32(0.5), u1, Float32(1.0)) * u1)) * fma(Float32(Float32(-2.0) * Float32(u2 * u2)), Float32(Float32(pi) * Float32(pi)), Float32(1.0))); else tmp = sqrt(exp(log(t_0))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -\log \left(1 - u1\right)\\
\mathbf{if}\;\sqrt{t\_0} \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right) \leq 0.07599999755620956:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(0.5, u1, 1\right) \cdot u1} \cdot \mathsf{fma}\left(-2 \cdot \left(u2 \cdot u2\right), \pi \cdot \pi, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{e^{\log t\_0}}\\
\end{array}
\end{array}
if (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.0759999976Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.1
Applied rewrites88.1%
Taylor expanded in u2 around 0
+-commutativeN/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lift-PI.f32N/A
lift-PI.f3279.3
Applied rewrites79.3%
if 0.0759999976 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) Initial program 57.5%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.8
Applied rewrites49.8%
lift-neg.f32N/A
lift--.f32N/A
lift-log.f32N/A
rem-exp-logN/A
lower-exp.f32N/A
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-log.f3249.8
Applied rewrites49.8%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (sqrt (- (log (- 1.0 u1))))))
(if (<= (* t_0 (cos (* (* 2.0 PI) u2))) 0.07599999755620956)
(* (sqrt (* (fma 0.5 u1 1.0) u1)) (fma (* -2.0 (* u2 u2)) (* PI PI) 1.0))
t_0)))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = sqrtf(-logf((1.0f - u1)));
float tmp;
if ((t_0 * cosf(((2.0f * ((float) M_PI)) * u2))) <= 0.07599999755620956f) {
tmp = sqrtf((fmaf(0.5f, u1, 1.0f) * u1)) * fmaf((-2.0f * (u2 * u2)), (((float) M_PI) * ((float) M_PI)), 1.0f);
} else {
tmp = t_0;
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) tmp = Float32(0.0) if (Float32(t_0 * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) <= Float32(0.07599999755620956)) tmp = Float32(sqrt(Float32(fma(Float32(0.5), u1, Float32(1.0)) * u1)) * fma(Float32(Float32(-2.0) * Float32(u2 * u2)), Float32(Float32(pi) * Float32(pi)), Float32(1.0))); else tmp = t_0; end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{-\log \left(1 - u1\right)}\\
\mathbf{if}\;t\_0 \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right) \leq 0.07599999755620956:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(0.5, u1, 1\right) \cdot u1} \cdot \mathsf{fma}\left(-2 \cdot \left(u2 \cdot u2\right), \pi \cdot \pi, 1\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.0759999976Initial program 57.5%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.1
Applied rewrites88.1%
Taylor expanded in u2 around 0
+-commutativeN/A
associate-*r*N/A
lower-fma.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
lift-PI.f32N/A
lift-PI.f3279.3
Applied rewrites79.3%
if 0.0759999976 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) Initial program 57.5%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.8
Applied rewrites49.8%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (sqrt (- (log (- 1.0 u1))))))
(if (<= (* t_0 (cos (* (* 2.0 PI) u2))) 0.20000000298023224)
(sqrt (* (fma (fma (fma 0.25 u1 0.3333333333333333) u1 0.5) u1 1.0) u1))
t_0)))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = sqrtf(-logf((1.0f - u1)));
float tmp;
if ((t_0 * cosf(((2.0f * ((float) M_PI)) * u2))) <= 0.20000000298023224f) {
tmp = sqrtf((fmaf(fmaf(fmaf(0.25f, u1, 0.3333333333333333f), u1, 0.5f), u1, 1.0f) * u1));
} else {
tmp = t_0;
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) tmp = Float32(0.0) if (Float32(t_0 * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) <= Float32(0.20000000298023224)) tmp = sqrt(Float32(fma(fma(fma(Float32(0.25), u1, Float32(0.3333333333333333)), u1, Float32(0.5)), u1, Float32(1.0)) * u1)); else tmp = t_0; end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{-\log \left(1 - u1\right)}\\
\mathbf{if}\;t\_0 \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right) \leq 0.20000000298023224:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.25, u1, 0.3333333333333333\right), u1, 0.5\right), u1, 1\right) \cdot u1}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.200000003Initial program 57.5%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.8
Applied rewrites49.8%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
+-commutativeN/A
lower-fma.f3276.5
Applied rewrites76.5%
if 0.200000003 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) Initial program 57.5%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.8
Applied rewrites49.8%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (sqrt (- (log (- 1.0 u1))))))
(if (<= (* t_0 (cos (* (* 2.0 PI) u2))) 0.10199999809265137)
(sqrt (* (fma (fma 0.3333333333333333 u1 0.5) u1 1.0) u1))
t_0)))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = sqrtf(-logf((1.0f - u1)));
float tmp;
if ((t_0 * cosf(((2.0f * ((float) M_PI)) * u2))) <= 0.10199999809265137f) {
tmp = sqrtf((fmaf(fmaf(0.3333333333333333f, u1, 0.5f), u1, 1.0f) * u1));
} else {
tmp = t_0;
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) tmp = Float32(0.0) if (Float32(t_0 * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) <= Float32(0.10199999809265137)) tmp = sqrt(Float32(fma(fma(Float32(0.3333333333333333), u1, Float32(0.5)), u1, Float32(1.0)) * u1)); else tmp = t_0; end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{-\log \left(1 - u1\right)}\\
\mathbf{if}\;t\_0 \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right) \leq 0.10199999809265137:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(\mathsf{fma}\left(0.3333333333333333, u1, 0.5\right), u1, 1\right) \cdot u1}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.102Initial program 57.5%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.8
Applied rewrites49.8%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
+-commutativeN/A
lower-fma.f3275.3
Applied rewrites75.3%
if 0.102 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) Initial program 57.5%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.8
Applied rewrites49.8%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (sqrt (- (log (- 1.0 u1))))))
(if (<= (* t_0 (cos (* (* 2.0 PI) u2))) 0.052000001072883606)
(sqrt (* (fma 0.5 u1 1.0) u1))
t_0)))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = sqrtf(-logf((1.0f - u1)));
float tmp;
if ((t_0 * cosf(((2.0f * ((float) M_PI)) * u2))) <= 0.052000001072883606f) {
tmp = sqrtf((fmaf(0.5f, u1, 1.0f) * u1));
} else {
tmp = t_0;
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) tmp = Float32(0.0) if (Float32(t_0 * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) <= Float32(0.052000001072883606)) tmp = sqrt(Float32(fma(Float32(0.5), u1, Float32(1.0)) * u1)); else tmp = t_0; end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{-\log \left(1 - u1\right)}\\
\mathbf{if}\;t\_0 \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right) \leq 0.052000001072883606:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(0.5, u1, 1\right) \cdot u1}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) < 0.0520000011Initial program 57.5%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.8
Applied rewrites49.8%
Taylor expanded in u1 around 0
Applied rewrites6.6%
Taylor expanded in u1 around 0
flip3--N/A
metadata-evalN/A
pow3N/A
metadata-evalN/A
diff-logN/A
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3272.7
Applied rewrites72.7%
if 0.0520000011 < (*.f32 (sqrt.f32 (neg.f32 (log.f32 (-.f32 #s(literal 1 binary32) u1)))) (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2))) Initial program 57.5%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.8
Applied rewrites49.8%
(FPCore (cosTheta_i u1 u2) :precision binary32 (sqrt (* (fma 0.5 u1 1.0) u1)))
float code(float cosTheta_i, float u1, float u2) {
return sqrtf((fmaf(0.5f, u1, 1.0f) * u1));
}
function code(cosTheta_i, u1, u2) return sqrt(Float32(fma(Float32(0.5), u1, Float32(1.0)) * u1)) end
\begin{array}{l}
\\
\sqrt{\mathsf{fma}\left(0.5, u1, 1\right) \cdot u1}
\end{array}
Initial program 57.5%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.8
Applied rewrites49.8%
Taylor expanded in u1 around 0
Applied rewrites6.6%
Taylor expanded in u1 around 0
flip3--N/A
metadata-evalN/A
pow3N/A
metadata-evalN/A
diff-logN/A
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3272.7
Applied rewrites72.7%
(FPCore (cosTheta_i u1 u2) :precision binary32 (sqrt (- (- u1))))
float code(float cosTheta_i, float u1, float u2) {
return sqrtf(-(-u1));
}
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(4) function code(costheta_i, u1, u2)
use fmin_fmax_functions
real(4), intent (in) :: costheta_i
real(4), intent (in) :: u1
real(4), intent (in) :: u2
code = sqrt(-(-u1))
end function
function code(cosTheta_i, u1, u2) return sqrt(Float32(-Float32(-u1))) end
function tmp = code(cosTheta_i, u1, u2) tmp = sqrt(-(-u1)); end
\begin{array}{l}
\\
\sqrt{-\left(-u1\right)}
\end{array}
Initial program 57.5%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.8
Applied rewrites49.8%
Taylor expanded in u1 around 0
mul-1-negN/A
lower-neg.f3264.8
Applied rewrites64.8%
herbie shell --seed 2025138
(FPCore (cosTheta_i u1 u2)
:name "Beckmann Sample, near normal, slope_x"
:precision binary32
:pre (and (and (and (> cosTheta_i 0.9999) (<= cosTheta_i 1.0)) (and (<= 2.328306437e-10 u1) (<= u1 1.0))) (and (<= 2.328306437e-10 u2) (<= u2 1.0)))
(* (sqrt (- (log (- 1.0 u1)))) (cos (* (* 2.0 PI) u2))))