
(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 25 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.026000000536441803)
(*
(sqrt (fma (* (fma (fma 0.25 u1 0.3333333333333333) u1 0.5) u1) u1 u1))
(cos (* u2 (+ PI PI))))
(* (/ 1.0 (pow (- (log (- 1.0 u1))) -0.5)) (cos (* (* 2.0 PI) u2)))))
float code(float cosTheta_i, float u1, float u2) {
float tmp;
if (u1 <= 0.026000000536441803f) {
tmp = sqrtf(fmaf((fmaf(fmaf(0.25f, u1, 0.3333333333333333f), u1, 0.5f) * u1), u1, u1)) * cosf((u2 * (((float) M_PI) + ((float) M_PI))));
} else {
tmp = (1.0f / powf(-logf((1.0f - u1)), -0.5f)) * cosf(((2.0f * ((float) M_PI)) * u2));
}
return tmp;
}
function code(cosTheta_i, u1, u2) tmp = Float32(0.0) if (u1 <= Float32(0.026000000536441803)) tmp = Float32(sqrt(fma(Float32(fma(fma(Float32(0.25), u1, Float32(0.3333333333333333)), u1, Float32(0.5)) * u1), u1, u1)) * cos(Float32(u2 * Float32(Float32(pi) + Float32(pi))))); else tmp = Float32(Float32(Float32(1.0) / (Float32(-log(Float32(Float32(1.0) - u1))) ^ Float32(-0.5))) * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;u1 \leq 0.026000000536441803:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.25, u1, 0.3333333333333333\right), u1, 0.5\right) \cdot u1, u1, u1\right)} \cdot \cos \left(u2 \cdot \left(\pi + \pi\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{{\left(-\log \left(1 - u1\right)\right)}^{-0.5}} \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right)\\
\end{array}
\end{array}
if u1 < 0.0260000005Initial program 57.7%
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%
lift-*.f32N/A
lift-fma.f32N/A
lift-fma.f32N/A
lift-fma.f32N/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
*-commutativeN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-inN/A
lower-fma.f32N/A
lower-*.f32N/A
*-commutativeN/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Applied rewrites93.7%
Applied rewrites93.7%
if 0.0260000005 < u1 Initial program 57.7%
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.f3257.7
Applied rewrites57.7%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(if (<= u1 0.019999999552965164)
(*
(sqrt (fma (* (fma (fma 0.25 u1 0.3333333333333333) u1 0.5) u1) u1 u1))
(cos (* u2 (+ PI PI))))
(* (exp (* (log (- (log (- 1.0 u1)))) 0.5)) (cos (* (* 2.0 PI) u2)))))
float code(float cosTheta_i, float u1, float u2) {
float tmp;
if (u1 <= 0.019999999552965164f) {
tmp = sqrtf(fmaf((fmaf(fmaf(0.25f, u1, 0.3333333333333333f), u1, 0.5f) * u1), u1, u1)) * cosf((u2 * (((float) M_PI) + ((float) M_PI))));
} else {
tmp = expf((logf(-logf((1.0f - u1))) * 0.5f)) * cosf(((2.0f * ((float) M_PI)) * u2));
}
return tmp;
}
function code(cosTheta_i, u1, u2) tmp = Float32(0.0) if (u1 <= Float32(0.019999999552965164)) tmp = Float32(sqrt(fma(Float32(fma(fma(Float32(0.25), u1, Float32(0.3333333333333333)), u1, Float32(0.5)) * u1), u1, u1)) * cos(Float32(u2 * Float32(Float32(pi) + Float32(pi))))); else tmp = Float32(exp(Float32(log(Float32(-log(Float32(Float32(1.0) - u1)))) * Float32(0.5))) * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;u1 \leq 0.019999999552965164:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.25, u1, 0.3333333333333333\right), u1, 0.5\right) \cdot u1, u1, u1\right)} \cdot \cos \left(u2 \cdot \left(\pi + \pi\right)\right)\\
\mathbf{else}:\\
\;\;\;\;e^{\log \left(-\log \left(1 - u1\right)\right) \cdot 0.5} \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right)\\
\end{array}
\end{array}
if u1 < 0.0199999996Initial program 57.7%
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%
lift-*.f32N/A
lift-fma.f32N/A
lift-fma.f32N/A
lift-fma.f32N/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
*-commutativeN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-inN/A
lower-fma.f32N/A
lower-*.f32N/A
*-commutativeN/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Applied rewrites93.7%
Applied rewrites93.7%
if 0.0199999996 < u1 Initial program 57.7%
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.7
Applied rewrites57.7%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(if (<= u1 0.026000000536441803)
(*
(sqrt (fma (* (fma (fma 0.25 u1 0.3333333333333333) u1 0.5) u1) u1 u1))
(cos (* u2 (+ PI PI))))
(*
(/ 1.0 (/ 1.0 (sqrt (- (log (- 1.0 u1))))))
(sin (fma (+ PI PI) u2 (* PI 0.5))))))
float code(float cosTheta_i, float u1, float u2) {
float tmp;
if (u1 <= 0.026000000536441803f) {
tmp = sqrtf(fmaf((fmaf(fmaf(0.25f, u1, 0.3333333333333333f), u1, 0.5f) * u1), u1, u1)) * cosf((u2 * (((float) M_PI) + ((float) M_PI))));
} else {
tmp = (1.0f / (1.0f / sqrtf(-logf((1.0f - u1))))) * sinf(fmaf((((float) M_PI) + ((float) M_PI)), u2, (((float) M_PI) * 0.5f)));
}
return tmp;
}
function code(cosTheta_i, u1, u2) tmp = Float32(0.0) if (u1 <= Float32(0.026000000536441803)) tmp = Float32(sqrt(fma(Float32(fma(fma(Float32(0.25), u1, Float32(0.3333333333333333)), u1, Float32(0.5)) * u1), u1, u1)) * cos(Float32(u2 * Float32(Float32(pi) + Float32(pi))))); else tmp = Float32(Float32(Float32(1.0) / Float32(Float32(1.0) / sqrt(Float32(-log(Float32(Float32(1.0) - u1)))))) * sin(fma(Float32(Float32(pi) + Float32(pi)), u2, Float32(Float32(pi) * Float32(0.5))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;u1 \leq 0.026000000536441803:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.25, u1, 0.3333333333333333\right), u1, 0.5\right) \cdot u1, u1, u1\right)} \cdot \cos \left(u2 \cdot \left(\pi + \pi\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{1}{\sqrt{-\log \left(1 - u1\right)}}} \cdot \sin \left(\mathsf{fma}\left(\pi + \pi, u2, \pi \cdot 0.5\right)\right)\\
\end{array}
\end{array}
if u1 < 0.0260000005Initial program 57.7%
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%
lift-*.f32N/A
lift-fma.f32N/A
lift-fma.f32N/A
lift-fma.f32N/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
*-commutativeN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-inN/A
lower-fma.f32N/A
lower-*.f32N/A
*-commutativeN/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Applied rewrites93.7%
Applied rewrites93.7%
if 0.0260000005 < u1 Initial program 57.7%
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.7
Applied rewrites57.7%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
associate-*l*N/A
lower-fma.f32N/A
count-2-revN/A
lift-+.f32N/A
lift-PI.f32N/A
lift-PI.f32N/A
*-commutativeN/A
lower-*.f32N/A
lift-PI.f3257.7
Applied rewrites57.7%
lift-sqrt.f32N/A
lift-neg.f32N/A
lift-log.f32N/A
lift--.f32N/A
pow1/2N/A
metadata-evalN/A
pow-flipN/A
lower-/.f32N/A
metadata-evalN/A
pow-negN/A
pow1/2N/A
lower-/.f32N/A
lift--.f32N/A
lift-log.f32N/A
lift-neg.f32N/A
lift-sqrt.f3257.6
Applied rewrites57.6%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (* u2 (+ PI PI))))
(if (<= u1 0.026000000536441803)
(*
(sqrt (fma (* (fma (fma 0.25 u1 0.3333333333333333) u1 0.5) u1) u1 u1))
(cos t_0))
(* (sqrt (- (log (- 1.0 u1)))) (sin (+ (- t_0) (/ PI 2.0)))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = u2 * (((float) M_PI) + ((float) M_PI));
float tmp;
if (u1 <= 0.026000000536441803f) {
tmp = sqrtf(fmaf((fmaf(fmaf(0.25f, u1, 0.3333333333333333f), u1, 0.5f) * u1), u1, u1)) * cosf(t_0);
} else {
tmp = sqrtf(-logf((1.0f - u1))) * sinf((-t_0 + (((float) M_PI) / 2.0f)));
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(u2 * Float32(Float32(pi) + Float32(pi))) tmp = Float32(0.0) if (u1 <= Float32(0.026000000536441803)) tmp = Float32(sqrt(fma(Float32(fma(fma(Float32(0.25), u1, Float32(0.3333333333333333)), u1, Float32(0.5)) * u1), u1, u1)) * cos(t_0)); else tmp = Float32(sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) * sin(Float32(Float32(-t_0) + Float32(Float32(pi) / Float32(2.0))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := u2 \cdot \left(\pi + \pi\right)\\
\mathbf{if}\;u1 \leq 0.026000000536441803:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.25, u1, 0.3333333333333333\right), u1, 0.5\right) \cdot u1, u1, u1\right)} \cdot \cos t\_0\\
\mathbf{else}:\\
\;\;\;\;\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(-t\_0\right) + \frac{\pi}{2}\right)\\
\end{array}
\end{array}
if u1 < 0.0260000005Initial program 57.7%
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%
lift-*.f32N/A
lift-fma.f32N/A
lift-fma.f32N/A
lift-fma.f32N/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
*-commutativeN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-inN/A
lower-fma.f32N/A
lower-*.f32N/A
*-commutativeN/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Applied rewrites93.7%
Applied rewrites93.7%
if 0.0260000005 < u1 Initial program 57.7%
lift-cos.f32N/A
cos-neg-revN/A
sin-+PI/2-revN/A
lower-sin.f32N/A
lower-+.f32N/A
lower-neg.f32N/A
lift-*.f32N/A
*-commutativeN/A
lower-*.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.7
Applied rewrites57.7%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (* u2 (+ PI PI))) (t_1 (log (- 1.0 u1))))
(if (<= t_1 -0.014000000432133675)
(* (sqrt (- t_1)) (sin (+ (- t_0) (/ PI 2.0))))
(* (sqrt (fma (* (fma 0.3333333333333333 u1 0.5) u1) u1 u1)) (cos t_0)))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = u2 * (((float) M_PI) + ((float) M_PI));
float t_1 = logf((1.0f - u1));
float tmp;
if (t_1 <= -0.014000000432133675f) {
tmp = sqrtf(-t_1) * sinf((-t_0 + (((float) M_PI) / 2.0f)));
} else {
tmp = sqrtf(fmaf((fmaf(0.3333333333333333f, u1, 0.5f) * u1), u1, u1)) * cosf(t_0);
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(u2 * Float32(Float32(pi) + Float32(pi))) t_1 = log(Float32(Float32(1.0) - u1)) tmp = Float32(0.0) if (t_1 <= Float32(-0.014000000432133675)) tmp = Float32(sqrt(Float32(-t_1)) * sin(Float32(Float32(-t_0) + Float32(Float32(pi) / Float32(2.0))))); else tmp = Float32(sqrt(fma(Float32(fma(Float32(0.3333333333333333), u1, Float32(0.5)) * u1), u1, u1)) * cos(t_0)); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := u2 \cdot \left(\pi + \pi\right)\\
t_1 := \log \left(1 - u1\right)\\
\mathbf{if}\;t\_1 \leq -0.014000000432133675:\\
\;\;\;\;\sqrt{-t\_1} \cdot \sin \left(\left(-t\_0\right) + \frac{\pi}{2}\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(\mathsf{fma}\left(0.3333333333333333, u1, 0.5\right) \cdot u1, u1, u1\right)} \cdot \cos t\_0\\
\end{array}
\end{array}
if (log.f32 (-.f32 #s(literal 1 binary32) u1)) < -0.0140000004Initial program 57.7%
lift-cos.f32N/A
cos-neg-revN/A
sin-+PI/2-revN/A
lower-sin.f32N/A
lower-+.f32N/A
lower-neg.f32N/A
lift-*.f32N/A
*-commutativeN/A
lower-*.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.7
Applied rewrites57.7%
if -0.0140000004 < (log.f32 (-.f32 #s(literal 1 binary32) u1)) Initial program 57.7%
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%
lift-*.f32N/A
lift-fma.f32N/A
lift-fma.f32N/A
lift-fma.f32N/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
*-commutativeN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-inN/A
lower-fma.f32N/A
lower-*.f32N/A
*-commutativeN/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Applied rewrites93.7%
Applied rewrites93.7%
Taylor expanded in u1 around 0
Applied rewrites91.8%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (log (- 1.0 u1))))
(if (<= t_0 -0.014000000432133675)
(* (sqrt (- t_0)) (cos (* (+ PI PI) u2)))
(*
(sqrt (fma (* (fma 0.3333333333333333 u1 0.5) u1) u1 u1))
(cos (* u2 (+ PI PI)))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = logf((1.0f - u1));
float tmp;
if (t_0 <= -0.014000000432133675f) {
tmp = sqrtf(-t_0) * cosf(((((float) M_PI) + ((float) M_PI)) * u2));
} else {
tmp = sqrtf(fmaf((fmaf(0.3333333333333333f, u1, 0.5f) * u1), u1, u1)) * cosf((u2 * (((float) M_PI) + ((float) M_PI))));
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = log(Float32(Float32(1.0) - u1)) tmp = Float32(0.0) if (t_0 <= Float32(-0.014000000432133675)) tmp = Float32(sqrt(Float32(-t_0)) * cos(Float32(Float32(Float32(pi) + Float32(pi)) * u2))); else tmp = Float32(sqrt(fma(Float32(fma(Float32(0.3333333333333333), u1, Float32(0.5)) * u1), u1, u1)) * cos(Float32(u2 * Float32(Float32(pi) + Float32(pi))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \log \left(1 - u1\right)\\
\mathbf{if}\;t\_0 \leq -0.014000000432133675:\\
\;\;\;\;\sqrt{-t\_0} \cdot \cos \left(\left(\pi + \pi\right) \cdot u2\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(\mathsf{fma}\left(0.3333333333333333, u1, 0.5\right) \cdot u1, u1, u1\right)} \cdot \cos \left(u2 \cdot \left(\pi + \pi\right)\right)\\
\end{array}
\end{array}
if (log.f32 (-.f32 #s(literal 1 binary32) u1)) < -0.0140000004Initial program 57.7%
lift-PI.f32N/A
lift-*.f32N/A
count-2-revN/A
lower-+.f32N/A
lift-PI.f32N/A
lift-PI.f3257.7
Applied rewrites57.7%
if -0.0140000004 < (log.f32 (-.f32 #s(literal 1 binary32) u1)) Initial program 57.7%
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%
lift-*.f32N/A
lift-fma.f32N/A
lift-fma.f32N/A
lift-fma.f32N/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
*-commutativeN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-inN/A
lower-fma.f32N/A
lower-*.f32N/A
*-commutativeN/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Applied rewrites93.7%
Applied rewrites93.7%
Taylor expanded in u1 around 0
Applied rewrites91.8%
(FPCore (cosTheta_i u1 u2) :precision binary32 (if (<= u1 0.006000000052154064) (* (sqrt (* (/ -1.0 (- (* 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.006000000052154064f) {
tmp = sqrtf(((-1.0f / ((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.006000000052154064)) tmp = Float32(sqrt(Float32(Float32(Float32(-1.0) / Float32(Float32(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
function tmp_2 = code(cosTheta_i, u1, u2) tmp = single(0.0); if (u1 <= single(0.006000000052154064)) tmp = sqrt(((single(-1.0) / ((single(0.5) * u1) - single(1.0))) * u1)) * cos(((single(2.0) * single(pi)) * u2)); else tmp = sqrt(-log((single(1.0) - u1))) * cos(((single(pi) + single(pi)) * u2)); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;u1 \leq 0.006000000052154064:\\
\;\;\;\;\sqrt{\frac{-1}{0.5 \cdot u1 - 1} \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.00600000005Initial program 57.7%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.0
Applied rewrites88.0%
lift-fma.f32N/A
flip-+N/A
lower-/.f32N/A
metadata-evalN/A
lower--.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f3288.0
Applied rewrites88.0%
Taylor expanded in u1 around 0
Applied rewrites89.7%
if 0.00600000005 < u1 Initial program 57.7%
lift-PI.f32N/A
lift-*.f32N/A
count-2-revN/A
lower-+.f32N/A
lift-PI.f32N/A
lift-PI.f3257.7
Applied rewrites57.7%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (log (- 1.0 u1))))
(if (<= t_0 -0.0035000001080334187)
(* (sqrt (- t_0)) (cos (* (+ PI PI) u2)))
(* (fma (* (sqrt u1) u1) 0.25 (sqrt u1)) (cos (* (* 2.0 PI) u2))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = logf((1.0f - u1));
float tmp;
if (t_0 <= -0.0035000001080334187f) {
tmp = sqrtf(-t_0) * cosf(((((float) M_PI) + ((float) M_PI)) * u2));
} else {
tmp = fmaf((sqrtf(u1) * u1), 0.25f, sqrtf(u1)) * cosf(((2.0f * ((float) M_PI)) * u2));
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = log(Float32(Float32(1.0) - u1)) tmp = Float32(0.0) if (t_0 <= Float32(-0.0035000001080334187)) tmp = Float32(sqrt(Float32(-t_0)) * cos(Float32(Float32(Float32(pi) + Float32(pi)) * u2))); else tmp = Float32(fma(Float32(sqrt(u1) * u1), Float32(0.25), sqrt(u1)) * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \log \left(1 - u1\right)\\
\mathbf{if}\;t\_0 \leq -0.0035000001080334187:\\
\;\;\;\;\sqrt{-t\_0} \cdot \cos \left(\left(\pi + \pi\right) \cdot u2\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\sqrt{u1} \cdot u1, 0.25, \sqrt{u1}\right) \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right)\\
\end{array}
\end{array}
if (log.f32 (-.f32 #s(literal 1 binary32) u1)) < -0.00350000011Initial program 57.7%
lift-PI.f32N/A
lift-*.f32N/A
count-2-revN/A
lower-+.f32N/A
lift-PI.f32N/A
lift-PI.f3257.7
Applied rewrites57.7%
if -0.00350000011 < (log.f32 (-.f32 #s(literal 1 binary32) u1)) Initial program 57.7%
Taylor expanded in u1 around 0
+-commutativeN/A
*-commutativeN/A
pow1/2N/A
pow-divN/A
metadata-evalN/A
metadata-evalN/A
sqrt-pow2N/A
lower-fma.f32N/A
cube-multN/A
rem-square-sqrtN/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower-sqrt.f3288.2
Applied rewrites88.2%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (log (- 1.0 u1))))
(if (<= t_0 -0.0035000001080334187)
(* (sqrt (- t_0)) (cos (* (+ PI PI) u2)))
(* (sqrt (fma u1 (* 0.5 u1) u1)) (cos (* u2 (+ PI PI)))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = logf((1.0f - u1));
float tmp;
if (t_0 <= -0.0035000001080334187f) {
tmp = sqrtf(-t_0) * cosf(((((float) M_PI) + ((float) M_PI)) * u2));
} else {
tmp = sqrtf(fmaf(u1, (0.5f * u1), u1)) * cosf((u2 * (((float) M_PI) + ((float) M_PI))));
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = log(Float32(Float32(1.0) - u1)) tmp = Float32(0.0) if (t_0 <= Float32(-0.0035000001080334187)) tmp = Float32(sqrt(Float32(-t_0)) * cos(Float32(Float32(Float32(pi) + Float32(pi)) * u2))); else tmp = Float32(sqrt(fma(u1, Float32(Float32(0.5) * u1), u1)) * cos(Float32(u2 * Float32(Float32(pi) + Float32(pi))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \log \left(1 - u1\right)\\
\mathbf{if}\;t\_0 \leq -0.0035000001080334187:\\
\;\;\;\;\sqrt{-t\_0} \cdot \cos \left(\left(\pi + \pi\right) \cdot u2\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(u1, 0.5 \cdot u1, u1\right)} \cdot \cos \left(u2 \cdot \left(\pi + \pi\right)\right)\\
\end{array}
\end{array}
if (log.f32 (-.f32 #s(literal 1 binary32) u1)) < -0.00350000011Initial program 57.7%
lift-PI.f32N/A
lift-*.f32N/A
count-2-revN/A
lower-+.f32N/A
lift-PI.f32N/A
lift-PI.f3257.7
Applied rewrites57.7%
if -0.00350000011 < (log.f32 (-.f32 #s(literal 1 binary32) u1)) Initial program 57.7%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.0
Applied rewrites88.0%
lift-fma.f32N/A
flip-+N/A
lower-/.f32N/A
metadata-evalN/A
lower--.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f3288.0
Applied rewrites88.0%
Applied rewrites88.1%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (sqrt (- (log (- 1.0 u1))))))
(if (<= (* t_0 (cos (* (* 2.0 PI) u2))) 0.06499999761581421)
(* (sqrt (fma u1 (* 0.5 u1) u1)) (cos (* u2 (+ PI PI))))
(*
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.06499999761581421f) {
tmp = sqrtf(fmaf(u1, (0.5f * u1), u1)) * cosf((u2 * (((float) M_PI) + ((float) M_PI))));
} 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.06499999761581421)) tmp = Float32(sqrt(fma(u1, Float32(Float32(0.5) * u1), u1)) * cos(Float32(u2 * Float32(Float32(pi) + Float32(pi))))); 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.06499999761581421:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(u1, 0.5 \cdot u1, u1\right)} \cdot \cos \left(u2 \cdot \left(\pi + \pi\right)\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.0649999976Initial program 57.7%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.0
Applied rewrites88.0%
lift-fma.f32N/A
flip-+N/A
lower-/.f32N/A
metadata-evalN/A
lower--.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f3288.0
Applied rewrites88.0%
Applied rewrites88.1%
if 0.0649999976 < (*.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.7%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites54.2%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (sqrt (- (log (- 1.0 u1))))))
(if (<= (* t_0 (cos (* (* 2.0 PI) u2))) 0.06499999761581421)
(* (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.06499999761581421f) {
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.06499999761581421)) 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.06499999761581421:\\
\;\;\;\;\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.0649999976Initial program 57.7%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.0
Applied rewrites88.0%
lift-PI.f32N/A
lift-*.f32N/A
count-2-revN/A
lift-+.f32N/A
lift-PI.f32N/A
lift-PI.f3288.0
Applied rewrites88.0%
if 0.0649999976 < (*.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.7%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites54.2%
(FPCore (cosTheta_i u1 u2) :precision binary32 (if (<= (cos (* (* 2.0 PI) u2)) 0.999970018863678) (* (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.999970018863678f) {
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.999970018863678)) 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.999970018863678:\\
\;\;\;\;\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.999970019Initial program 57.7%
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.999970019 < (cos.f32 (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2)) Initial program 57.7%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.3
Applied rewrites49.3%
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.3
Applied rewrites49.3%
lift--.f32N/A
lift-log.f32N/A
*-rgt-identityN/A
*-commutativeN/A
metadata-evalN/A
fp-cancel-sign-sub-invN/A
mul-1-negN/A
lift-neg.f32N/A
lower-log1p.f3279.7
Applied rewrites79.7%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (sqrt (- (log (- 1.0 u1))))) (t_1 (* -2.0 (* u2 u2))))
(if (<= (* t_0 (cos (* (* 2.0 PI) u2))) 0.16200000047683716)
(*
(sqrt
(fma
u1
1.0
(* u1 (* (fma (fma 0.25 u1 0.3333333333333333) u1 0.5) u1))))
(fma t_1 (* PI PI) 1.0))
(fma t_1 (* (* PI PI) t_0) t_0))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = sqrtf(-logf((1.0f - u1)));
float t_1 = -2.0f * (u2 * u2);
float tmp;
if ((t_0 * cosf(((2.0f * ((float) M_PI)) * u2))) <= 0.16200000047683716f) {
tmp = sqrtf(fmaf(u1, 1.0f, (u1 * (fmaf(fmaf(0.25f, u1, 0.3333333333333333f), u1, 0.5f) * u1)))) * fmaf(t_1, (((float) M_PI) * ((float) M_PI)), 1.0f);
} else {
tmp = fmaf(t_1, ((((float) M_PI) * ((float) M_PI)) * t_0), t_0);
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) t_1 = Float32(Float32(-2.0) * Float32(u2 * u2)) tmp = Float32(0.0) if (Float32(t_0 * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) <= Float32(0.16200000047683716)) tmp = Float32(sqrt(fma(u1, Float32(1.0), Float32(u1 * Float32(fma(fma(Float32(0.25), u1, Float32(0.3333333333333333)), u1, Float32(0.5)) * u1)))) * fma(t_1, Float32(Float32(pi) * Float32(pi)), Float32(1.0))); else tmp = fma(t_1, Float32(Float32(Float32(pi) * Float32(pi)) * t_0), t_0); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{-\log \left(1 - u1\right)}\\
t_1 := -2 \cdot \left(u2 \cdot u2\right)\\
\mathbf{if}\;t\_0 \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right) \leq 0.16200000047683716:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(u1, 1, u1 \cdot \left(\mathsf{fma}\left(\mathsf{fma}\left(0.25, u1, 0.3333333333333333\right), u1, 0.5\right) \cdot u1\right)\right)} \cdot \mathsf{fma}\left(t\_1, \pi \cdot \pi, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t\_1, \left(\pi \cdot \pi\right) \cdot t\_0, t\_0\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.162Initial program 57.7%
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%
lift-*.f32N/A
lift-fma.f32N/A
lift-fma.f32N/A
lift-fma.f32N/A
+-commutativeN/A
+-commutativeN/A
*-commutativeN/A
*-commutativeN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-inN/A
lower-fma.f32N/A
lower-*.f32N/A
*-commutativeN/A
*-commutativeN/A
+-commutativeN/A
+-commutativeN/A
Applied rewrites93.7%
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.f3283.8
Applied rewrites83.8%
if 0.162 < (*.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.7%
Taylor expanded in u2 around 0
+-commutativeN/A
associate-*r*N/A
lower-fma.f32N/A
Applied rewrites52.8%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (sqrt (- (log (- 1.0 u1))))) (t_1 (* -2.0 (* u2 u2))))
(if (<= (* t_0 (cos (* (* 2.0 PI) u2))) 0.16200000047683716)
(*
(sqrt (* (fma (fma (fma 0.25 u1 0.3333333333333333) u1 0.5) u1 1.0) u1))
(fma t_1 (* PI PI) 1.0))
(fma t_1 (* (* PI PI) t_0) t_0))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = sqrtf(-logf((1.0f - u1)));
float t_1 = -2.0f * (u2 * u2);
float tmp;
if ((t_0 * cosf(((2.0f * ((float) M_PI)) * u2))) <= 0.16200000047683716f) {
tmp = sqrtf((fmaf(fmaf(fmaf(0.25f, u1, 0.3333333333333333f), u1, 0.5f), u1, 1.0f) * u1)) * fmaf(t_1, (((float) M_PI) * ((float) M_PI)), 1.0f);
} else {
tmp = fmaf(t_1, ((((float) M_PI) * ((float) M_PI)) * t_0), t_0);
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) t_1 = Float32(Float32(-2.0) * Float32(u2 * u2)) tmp = Float32(0.0) if (Float32(t_0 * cos(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) <= Float32(0.16200000047683716)) tmp = Float32(sqrt(Float32(fma(fma(fma(Float32(0.25), u1, Float32(0.3333333333333333)), u1, Float32(0.5)), u1, Float32(1.0)) * u1)) * fma(t_1, Float32(Float32(pi) * Float32(pi)), Float32(1.0))); else tmp = fma(t_1, Float32(Float32(Float32(pi) * Float32(pi)) * t_0), t_0); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{-\log \left(1 - u1\right)}\\
t_1 := -2 \cdot \left(u2 \cdot u2\right)\\
\mathbf{if}\;t\_0 \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right) \leq 0.16200000047683716:\\
\;\;\;\;\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 \mathsf{fma}\left(t\_1, \pi \cdot \pi, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t\_1, \left(\pi \cdot \pi\right) \cdot t\_0, t\_0\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.162Initial program 57.7%
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%
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.f3283.7
Applied rewrites83.7%
if 0.162 < (*.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.7%
Taylor expanded in u2 around 0
+-commutativeN/A
associate-*r*N/A
lower-fma.f32N/A
Applied rewrites52.8%
(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.16200000047683716)
(*
(sqrt (* (fma (fma (fma 0.25 u1 0.3333333333333333) u1 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.16200000047683716f) {
tmp = sqrtf((fmaf(fmaf(fmaf(0.25f, u1, 0.3333333333333333f), u1, 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.16200000047683716)) tmp = Float32(sqrt(Float32(fma(fma(fma(Float32(0.25), u1, Float32(0.3333333333333333)), u1, 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.16200000047683716:\\
\;\;\;\;\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 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.162Initial program 57.7%
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%
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.f3283.7
Applied rewrites83.7%
if 0.162 < (*.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.7%
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.f3252.8
Applied rewrites52.8%
(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.07999999821186066)
(* (sqrt (* (/ -1.0 (- (* 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.07999999821186066f) {
tmp = sqrtf(((-1.0f / ((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.07999999821186066)) tmp = Float32(sqrt(Float32(Float32(Float32(-1.0) / Float32(Float32(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.07999999821186066:\\
\;\;\;\;\sqrt{\frac{-1}{0.5 \cdot u1 - 1} \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.0799999982Initial program 57.7%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.0
Applied rewrites88.0%
lift-fma.f32N/A
flip-+N/A
lower-/.f32N/A
metadata-evalN/A
lower--.f32N/A
lower-*.f32N/A
lower-*.f32N/A
lower-*.f32N/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.1
Applied rewrites79.1%
Taylor expanded in u1 around 0
Applied rewrites80.6%
if 0.0799999982 < (*.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.7%
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.f3252.8
Applied rewrites52.8%
(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.06499999761581421)
(* (sqrt (fma (* u1 u1) 0.5 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.06499999761581421f) {
tmp = sqrtf(fmaf((u1 * u1), 0.5f, 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.06499999761581421)) tmp = Float32(sqrt(fma(Float32(u1 * u1), Float32(0.5), 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.06499999761581421:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(u1 \cdot u1, 0.5, u1\right)} \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.0649999976Initial program 57.7%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.0
Applied rewrites88.0%
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.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.2
Applied rewrites79.2%
if 0.0649999976 < (*.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.7%
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.f3252.8
Applied rewrites52.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.07699999958276749)
(* (sqrt (fma (* u1 u1) 0.5 u1)) (fma (* -2.0 (* u2 u2)) (* PI PI) 1.0))
(/ 1.0 (/ 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.07699999958276749f) {
tmp = sqrtf(fmaf((u1 * u1), 0.5f, u1)) * fmaf((-2.0f * (u2 * u2)), (((float) M_PI) * ((float) M_PI)), 1.0f);
} else {
tmp = 1.0f / (1.0f / 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.07699999958276749)) tmp = Float32(sqrt(fma(Float32(u1 * u1), Float32(0.5), u1)) * fma(Float32(Float32(-2.0) * Float32(u2 * u2)), Float32(Float32(pi) * Float32(pi)), Float32(1.0))); else tmp = Float32(Float32(1.0) / Float32(Float32(1.0) / 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.07699999958276749:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(u1 \cdot u1, 0.5, u1\right)} \cdot \mathsf{fma}\left(-2 \cdot \left(u2 \cdot u2\right), \pi \cdot \pi, 1\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{1}{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.0769999996Initial program 57.7%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
+-commutativeN/A
lower-fma.f3288.0
Applied rewrites88.0%
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.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.2
Applied rewrites79.2%
if 0.0769999996 < (*.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.7%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.3
Applied rewrites49.3%
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.3
Applied rewrites49.3%
lift-pow.f32N/A
lift-neg.f32N/A
lift--.f32N/A
lift-log.f32N/A
metadata-evalN/A
pow-negN/A
pow1/2N/A
lower-/.f32N/A
lower-sqrt.f32N/A
lift-log.f32N/A
lift--.f32N/A
lift-neg.f3249.3
Applied rewrites49.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.15000000596046448)
(sqrt
(-
(*
(- (* (- (* (- (* -0.25 u1) 0.3333333333333333) u1) 0.5) u1) 1.0)
u1)))
(/ 1.0 (/ 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.15000000596046448f) {
tmp = sqrtf(-(((((((-0.25f * u1) - 0.3333333333333333f) * u1) - 0.5f) * u1) - 1.0f) * u1));
} else {
tmp = 1.0f / (1.0f / 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.15000000596046448)) tmp = sqrt(Float32(-Float32(Float32(Float32(Float32(Float32(Float32(Float32(Float32(-0.25) * u1) - Float32(0.3333333333333333)) * u1) - Float32(0.5)) * u1) - Float32(1.0)) * u1))); else tmp = Float32(Float32(1.0) / Float32(Float32(1.0) / t_0)); end return tmp end
function tmp_2 = code(cosTheta_i, u1, u2) t_0 = sqrt(-log((single(1.0) - u1))); tmp = single(0.0); if ((t_0 * cos(((single(2.0) * single(pi)) * u2))) <= single(0.15000000596046448)) tmp = sqrt(-(((((((single(-0.25) * u1) - single(0.3333333333333333)) * u1) - single(0.5)) * u1) - single(1.0)) * u1)); else tmp = single(1.0) / (single(1.0) / t_0); end tmp_2 = 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.15000000596046448:\\
\;\;\;\;\sqrt{-\left(\left(\left(-0.25 \cdot u1 - 0.3333333333333333\right) \cdot u1 - 0.5\right) \cdot u1 - 1\right) \cdot u1}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{1}{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.150000006Initial program 57.7%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.3
Applied rewrites49.3%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
lower--.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower--.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f3276.3
Applied rewrites76.3%
if 0.150000006 < (*.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.7%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.3
Applied rewrites49.3%
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.3
Applied rewrites49.3%
lift-pow.f32N/A
lift-neg.f32N/A
lift--.f32N/A
lift-log.f32N/A
metadata-evalN/A
pow-negN/A
pow1/2N/A
lower-/.f32N/A
lower-sqrt.f32N/A
lift-log.f32N/A
lift--.f32N/A
lift-neg.f3249.3
Applied rewrites49.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.11500000208616257)
(sqrt (- (* (- (* (- (* -0.3333333333333333 u1) 0.5) u1) 1.0) u1)))
(/ 1.0 (/ 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.11500000208616257f) {
tmp = sqrtf(-(((((-0.3333333333333333f * u1) - 0.5f) * u1) - 1.0f) * u1));
} else {
tmp = 1.0f / (1.0f / 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.11500000208616257)) tmp = sqrt(Float32(-Float32(Float32(Float32(Float32(Float32(Float32(-0.3333333333333333) * u1) - Float32(0.5)) * u1) - Float32(1.0)) * u1))); else tmp = Float32(Float32(1.0) / Float32(Float32(1.0) / t_0)); end return tmp end
function tmp_2 = code(cosTheta_i, u1, u2) t_0 = sqrt(-log((single(1.0) - u1))); tmp = single(0.0); if ((t_0 * cos(((single(2.0) * single(pi)) * u2))) <= single(0.11500000208616257)) tmp = sqrt(-(((((single(-0.3333333333333333) * u1) - single(0.5)) * u1) - single(1.0)) * u1)); else tmp = single(1.0) / (single(1.0) / t_0); end tmp_2 = 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.11500000208616257:\\
\;\;\;\;\sqrt{-\left(\left(-0.3333333333333333 \cdot u1 - 0.5\right) \cdot u1 - 1\right) \cdot u1}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{1}{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.115000002Initial program 57.7%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.3
Applied rewrites49.3%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
lower--.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f3275.1
Applied rewrites75.1%
if 0.115000002 < (*.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.7%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.3
Applied rewrites49.3%
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.3
Applied rewrites49.3%
lift-pow.f32N/A
lift-neg.f32N/A
lift--.f32N/A
lift-log.f32N/A
metadata-evalN/A
pow-negN/A
pow1/2N/A
lower-/.f32N/A
lower-sqrt.f32N/A
lift-log.f32N/A
lift--.f32N/A
lift-neg.f3249.3
Applied rewrites49.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.11500000208616257)
(sqrt (- (* (- (* (- (* -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.11500000208616257f) {
tmp = sqrtf(-(((((-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.11500000208616257)) tmp = sqrt(Float32(-Float32(Float32(Float32(Float32(Float32(Float32(-0.3333333333333333) * u1) - Float32(0.5)) * u1) - Float32(1.0)) * u1))); else tmp = t_0; end return tmp end
function tmp_2 = code(cosTheta_i, u1, u2) t_0 = sqrt(-log((single(1.0) - u1))); tmp = single(0.0); if ((t_0 * cos(((single(2.0) * single(pi)) * u2))) <= single(0.11500000208616257)) tmp = sqrt(-(((((single(-0.3333333333333333) * u1) - single(0.5)) * u1) - single(1.0)) * u1)); else tmp = t_0; end tmp_2 = 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.11500000208616257:\\
\;\;\;\;\sqrt{-\left(\left(-0.3333333333333333 \cdot u1 - 0.5\right) \cdot 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.115000002Initial program 57.7%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.3
Applied rewrites49.3%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
lower--.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f3275.1
Applied rewrites75.1%
if 0.115000002 < (*.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.7%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.3
Applied rewrites49.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.061000000685453415)
(fma (/ u1 (/ 1.0 (sqrt u1))) 0.25 (sqrt 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.061000000685453415f) {
tmp = fmaf((u1 / (1.0f / sqrtf(u1))), 0.25f, sqrtf(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.061000000685453415)) tmp = fma(Float32(u1 / Float32(Float32(1.0) / sqrt(u1))), Float32(0.25), sqrt(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.061000000685453415:\\
\;\;\;\;\mathsf{fma}\left(\frac{u1}{\frac{1}{\sqrt{u1}}}, 0.25, \sqrt{u1}\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.0610000007Initial program 57.7%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.3
Applied rewrites49.3%
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.3
Applied rewrites49.3%
Taylor expanded in u1 around 0
*-commutativeN/A
pow-flipN/A
metadata-evalN/A
pow1/2N/A
lower-fma.f32N/A
lower-/.f32N/A
metadata-evalN/A
pow-negN/A
pow1/2N/A
lower-/.f32N/A
lower-sqrt.f32N/A
lower-sqrt.f3272.8
Applied rewrites72.8%
if 0.0610000007 < (*.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.7%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.3
Applied rewrites49.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.05999999865889549)
(sqrt (- (* (- (* -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.05999999865889549f) {
tmp = sqrtf(-(((-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.05999999865889549)) tmp = sqrt(Float32(-Float32(Float32(Float32(Float32(-0.5) * u1) - Float32(1.0)) * u1))); else tmp = t_0; end return tmp end
function tmp_2 = code(cosTheta_i, u1, u2) t_0 = sqrt(-log((single(1.0) - u1))); tmp = single(0.0); if ((t_0 * cos(((single(2.0) * single(pi)) * u2))) <= single(0.05999999865889549)) tmp = sqrt(-(((single(-0.5) * u1) - single(1.0)) * u1)); else tmp = t_0; end tmp_2 = 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.05999999865889549:\\
\;\;\;\;\sqrt{-\left(-0.5 \cdot 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.0599999987Initial program 57.7%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.3
Applied rewrites49.3%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
lower--.f32N/A
lower-*.f3272.7
Applied rewrites72.7%
if 0.0599999987 < (*.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.7%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.3
Applied rewrites49.3%
(FPCore (cosTheta_i u1 u2) :precision binary32 (sqrt (- (* (- (* -0.5 u1) 1.0) u1))))
float code(float cosTheta_i, float u1, float u2) {
return sqrtf(-(((-0.5f * u1) - 1.0f) * 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(-((((-0.5e0) * u1) - 1.0e0) * u1))
end function
function code(cosTheta_i, u1, u2) return sqrt(Float32(-Float32(Float32(Float32(Float32(-0.5) * u1) - Float32(1.0)) * u1))) end
function tmp = code(cosTheta_i, u1, u2) tmp = sqrt(-(((single(-0.5) * u1) - single(1.0)) * u1)); end
\begin{array}{l}
\\
\sqrt{-\left(-0.5 \cdot u1 - 1\right) \cdot u1}
\end{array}
Initial program 57.7%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.3
Applied rewrites49.3%
Taylor expanded in u1 around 0
*-commutativeN/A
lower-*.f32N/A
lower--.f32N/A
lower-*.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.7%
Taylor expanded in u2 around 0
lift-log.f32N/A
lift--.f32N/A
lift-neg.f32N/A
lift-sqrt.f3249.3
Applied rewrites49.3%
Taylor expanded in u1 around 0
mul-1-negN/A
lower-neg.f3264.9
Applied rewrites64.9%
herbie shell --seed 2025139
(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))))