
(FPCore (cosTheta_i u1 u2) :precision binary32 (* (sqrt (- (log (- 1.0 u1)))) (sin (* (* 2.0 PI) u2))))
float code(float cosTheta_i, float u1, float u2) {
return sqrtf(-logf((1.0f - u1))) * sinf(((2.0f * ((float) M_PI)) * u2));
}
function code(cosTheta_i, u1, u2) return Float32(sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) * sin(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) end
function tmp = code(cosTheta_i, u1, u2) tmp = sqrt(-log((single(1.0) - u1))) * sin(((single(2.0) * single(pi)) * u2)); end
\begin{array}{l}
\\
\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right)
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 21 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (cosTheta_i u1 u2) :precision binary32 (* (sqrt (- (log (- 1.0 u1)))) (sin (* (* 2.0 PI) u2))))
float code(float cosTheta_i, float u1, float u2) {
return sqrtf(-logf((1.0f - u1))) * sinf(((2.0f * ((float) M_PI)) * u2));
}
function code(cosTheta_i, u1, u2) return Float32(sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) * sin(Float32(Float32(Float32(2.0) * Float32(pi)) * u2))) end
function tmp = code(cosTheta_i, u1, u2) tmp = sqrt(-log((single(1.0) - u1))) * sin(((single(2.0) * single(pi)) * u2)); end
\begin{array}{l}
\\
\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right)
\end{array}
(FPCore (cosTheta_i u1 u2) :precision binary32 (* (sqrt (- (log1p (- u1)))) (sin (* 2.0 (* PI u2)))))
float code(float cosTheta_i, float u1, float u2) {
return sqrtf(-log1pf(-u1)) * sinf((2.0f * (((float) M_PI) * u2)));
}
function code(cosTheta_i, u1, u2) return Float32(sqrt(Float32(-log1p(Float32(-u1)))) * sin(Float32(Float32(2.0) * Float32(Float32(pi) * u2)))) end
\begin{array}{l}
\\
\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(2 \cdot \left(\pi \cdot u2\right)\right)
\end{array}
Initial program 59.1%
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
neg-lowering-neg.f32N/A
sub-negN/A
log1p-defineN/A
log1p-lowering-log1p.f32N/A
neg-lowering-neg.f32N/A
sin-lowering-sin.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.3%
Simplified98.3%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (* u2 (* 2.0 PI)))
(t_1 (+ 1.0 (* u1 (+ 0.5 (* u1 (+ 0.3333333333333333 (* u1 0.25))))))))
(if (<= t_0 0.09000000357627869)
(*
(sqrt (- (log1p (- u1))))
(* u2 (* PI (+ 2.0 (* (* -1.3333333333333333 (* u2 u2)) (* PI PI))))))
(* (sin t_0) (pow (* t_1 (* (* u1 u1) t_1)) 0.25)))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = u2 * (2.0f * ((float) M_PI));
float t_1 = 1.0f + (u1 * (0.5f + (u1 * (0.3333333333333333f + (u1 * 0.25f)))));
float tmp;
if (t_0 <= 0.09000000357627869f) {
tmp = sqrtf(-log1pf(-u1)) * (u2 * (((float) M_PI) * (2.0f + ((-1.3333333333333333f * (u2 * u2)) * (((float) M_PI) * ((float) M_PI))))));
} else {
tmp = sinf(t_0) * powf((t_1 * ((u1 * u1) * t_1)), 0.25f);
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi))) t_1 = Float32(Float32(1.0) + Float32(u1 * Float32(Float32(0.5) + Float32(u1 * Float32(Float32(0.3333333333333333) + Float32(u1 * Float32(0.25))))))) tmp = Float32(0.0) if (t_0 <= Float32(0.09000000357627869)) tmp = Float32(sqrt(Float32(-log1p(Float32(-u1)))) * Float32(u2 * Float32(Float32(pi) * Float32(Float32(2.0) + Float32(Float32(Float32(-1.3333333333333333) * Float32(u2 * u2)) * Float32(Float32(pi) * Float32(pi))))))); else tmp = Float32(sin(t_0) * (Float32(t_1 * Float32(Float32(u1 * u1) * t_1)) ^ Float32(0.25))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := u2 \cdot \left(2 \cdot \pi\right)\\
t_1 := 1 + u1 \cdot \left(0.5 + u1 \cdot \left(0.3333333333333333 + u1 \cdot 0.25\right)\right)\\
\mathbf{if}\;t\_0 \leq 0.09000000357627869:\\
\;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(u2 \cdot \left(\pi \cdot \left(2 + \left(-1.3333333333333333 \cdot \left(u2 \cdot u2\right)\right) \cdot \left(\pi \cdot \pi\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin t\_0 \cdot {\left(t\_1 \cdot \left(\left(u1 \cdot u1\right) \cdot t\_1\right)\right)}^{0.25}\\
\end{array}
\end{array}
if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 0.0900000036Initial program 59.6%
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
neg-lowering-neg.f32N/A
sub-negN/A
log1p-defineN/A
log1p-lowering-log1p.f32N/A
neg-lowering-neg.f32N/A
sin-lowering-sin.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.4%
Simplified98.4%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f32N/A
associate-*r*N/A
*-commutativeN/A
unpow3N/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
Simplified98.3%
if 0.0900000036 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) Initial program 57.3%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3295.5%
Simplified95.5%
pow1/2N/A
metadata-evalN/A
pow-prod-upN/A
pow-prod-downN/A
pow-lowering-pow.f32N/A
Applied egg-rr95.6%
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f32N/A
Applied egg-rr95.7%
Final simplification97.7%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (* u2 (* 2.0 PI))))
(if (<= t_0 0.09000000357627869)
(*
(sqrt (- (log1p (- u1))))
(* u2 (* PI (+ 2.0 (* (* -1.3333333333333333 (* u2 u2)) (* PI PI))))))
(*
(sin t_0)
(sqrt
(/
(*
u1
(-
1.0
(*
(* u1 u1)
(+ 0.25 (* u1 (+ 0.3333333333333333 (* u1 0.3611111111111111)))))))
(- 1.0 (* u1 (+ 0.5 (* u1 (+ 0.3333333333333333 (* u1 0.25))))))))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = u2 * (2.0f * ((float) M_PI));
float tmp;
if (t_0 <= 0.09000000357627869f) {
tmp = sqrtf(-log1pf(-u1)) * (u2 * (((float) M_PI) * (2.0f + ((-1.3333333333333333f * (u2 * u2)) * (((float) M_PI) * ((float) M_PI))))));
} else {
tmp = sinf(t_0) * sqrtf(((u1 * (1.0f - ((u1 * u1) * (0.25f + (u1 * (0.3333333333333333f + (u1 * 0.3611111111111111f))))))) / (1.0f - (u1 * (0.5f + (u1 * (0.3333333333333333f + (u1 * 0.25f))))))));
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi))) tmp = Float32(0.0) if (t_0 <= Float32(0.09000000357627869)) tmp = Float32(sqrt(Float32(-log1p(Float32(-u1)))) * Float32(u2 * Float32(Float32(pi) * Float32(Float32(2.0) + Float32(Float32(Float32(-1.3333333333333333) * Float32(u2 * u2)) * Float32(Float32(pi) * Float32(pi))))))); else tmp = Float32(sin(t_0) * sqrt(Float32(Float32(u1 * Float32(Float32(1.0) - Float32(Float32(u1 * u1) * Float32(Float32(0.25) + Float32(u1 * Float32(Float32(0.3333333333333333) + Float32(u1 * Float32(0.3611111111111111)))))))) / Float32(Float32(1.0) - Float32(u1 * Float32(Float32(0.5) + Float32(u1 * Float32(Float32(0.3333333333333333) + Float32(u1 * Float32(0.25)))))))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := u2 \cdot \left(2 \cdot \pi\right)\\
\mathbf{if}\;t\_0 \leq 0.09000000357627869:\\
\;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(u2 \cdot \left(\pi \cdot \left(2 + \left(-1.3333333333333333 \cdot \left(u2 \cdot u2\right)\right) \cdot \left(\pi \cdot \pi\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin t\_0 \cdot \sqrt{\frac{u1 \cdot \left(1 - \left(u1 \cdot u1\right) \cdot \left(0.25 + u1 \cdot \left(0.3333333333333333 + u1 \cdot 0.3611111111111111\right)\right)\right)}{1 - u1 \cdot \left(0.5 + u1 \cdot \left(0.3333333333333333 + u1 \cdot 0.25\right)\right)}}\\
\end{array}
\end{array}
if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 0.0900000036Initial program 59.6%
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
neg-lowering-neg.f32N/A
sub-negN/A
log1p-defineN/A
log1p-lowering-log1p.f32N/A
neg-lowering-neg.f32N/A
sin-lowering-sin.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.4%
Simplified98.4%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f32N/A
associate-*r*N/A
*-commutativeN/A
unpow3N/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
Simplified98.3%
if 0.0900000036 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) Initial program 57.3%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3295.5%
Simplified95.5%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
/-lowering-/.f32N/A
Applied egg-rr95.7%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3295.7%
Simplified95.7%
Final simplification97.7%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (+ 0.5 (* u1 (+ 0.3333333333333333 (* u1 0.25))))))
(*
(sqrt
(/
(* u1 (- 1.0 (* t_0 (* (+ 0.5 (* u1 0.3333333333333333)) (* u1 u1)))))
(- 1.0 (* u1 t_0))))
(sin (* u2 (* 2.0 PI))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = 0.5f + (u1 * (0.3333333333333333f + (u1 * 0.25f)));
return sqrtf(((u1 * (1.0f - (t_0 * ((0.5f + (u1 * 0.3333333333333333f)) * (u1 * u1))))) / (1.0f - (u1 * t_0)))) * sinf((u2 * (2.0f * ((float) M_PI))));
}
function code(cosTheta_i, u1, u2) t_0 = Float32(Float32(0.5) + Float32(u1 * Float32(Float32(0.3333333333333333) + Float32(u1 * Float32(0.25))))) return Float32(sqrt(Float32(Float32(u1 * Float32(Float32(1.0) - Float32(t_0 * Float32(Float32(Float32(0.5) + Float32(u1 * Float32(0.3333333333333333))) * Float32(u1 * u1))))) / Float32(Float32(1.0) - Float32(u1 * t_0)))) * sin(Float32(u2 * Float32(Float32(2.0) * Float32(pi))))) end
function tmp = code(cosTheta_i, u1, u2) t_0 = single(0.5) + (u1 * (single(0.3333333333333333) + (u1 * single(0.25)))); tmp = sqrt(((u1 * (single(1.0) - (t_0 * ((single(0.5) + (u1 * single(0.3333333333333333))) * (u1 * u1))))) / (single(1.0) - (u1 * t_0)))) * sin((u2 * (single(2.0) * single(pi)))); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.5 + u1 \cdot \left(0.3333333333333333 + u1 \cdot 0.25\right)\\
\sqrt{\frac{u1 \cdot \left(1 - t\_0 \cdot \left(\left(0.5 + u1 \cdot 0.3333333333333333\right) \cdot \left(u1 \cdot u1\right)\right)\right)}{1 - u1 \cdot t\_0}} \cdot \sin \left(u2 \cdot \left(2 \cdot \pi\right)\right)
\end{array}
\end{array}
Initial program 59.1%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3294.0%
Simplified94.0%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
/-lowering-/.f32N/A
Applied egg-rr94.1%
Taylor expanded in u1 around 0
Simplified94.7%
Final simplification94.7%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (* u2 (* 2.0 PI))))
(if (<= t_0 0.09000000357627869)
(*
(sqrt (- (log1p (- u1))))
(* u2 (* PI (+ 2.0 (* (* -1.3333333333333333 (* u2 u2)) (* PI PI))))))
(*
(sin t_0)
(sqrt
(+
u1
(* (+ 0.5 (* u1 (+ 0.3333333333333333 (* u1 0.25)))) (* u1 u1))))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = u2 * (2.0f * ((float) M_PI));
float tmp;
if (t_0 <= 0.09000000357627869f) {
tmp = sqrtf(-log1pf(-u1)) * (u2 * (((float) M_PI) * (2.0f + ((-1.3333333333333333f * (u2 * u2)) * (((float) M_PI) * ((float) M_PI))))));
} else {
tmp = sinf(t_0) * sqrtf((u1 + ((0.5f + (u1 * (0.3333333333333333f + (u1 * 0.25f)))) * (u1 * u1))));
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi))) tmp = Float32(0.0) if (t_0 <= Float32(0.09000000357627869)) tmp = Float32(sqrt(Float32(-log1p(Float32(-u1)))) * Float32(u2 * Float32(Float32(pi) * Float32(Float32(2.0) + Float32(Float32(Float32(-1.3333333333333333) * Float32(u2 * u2)) * Float32(Float32(pi) * Float32(pi))))))); else tmp = Float32(sin(t_0) * sqrt(Float32(u1 + Float32(Float32(Float32(0.5) + Float32(u1 * Float32(Float32(0.3333333333333333) + Float32(u1 * Float32(0.25))))) * Float32(u1 * u1))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := u2 \cdot \left(2 \cdot \pi\right)\\
\mathbf{if}\;t\_0 \leq 0.09000000357627869:\\
\;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(u2 \cdot \left(\pi \cdot \left(2 + \left(-1.3333333333333333 \cdot \left(u2 \cdot u2\right)\right) \cdot \left(\pi \cdot \pi\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin t\_0 \cdot \sqrt{u1 + \left(0.5 + u1 \cdot \left(0.3333333333333333 + u1 \cdot 0.25\right)\right) \cdot \left(u1 \cdot u1\right)}\\
\end{array}
\end{array}
if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 0.0900000036Initial program 59.6%
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
neg-lowering-neg.f32N/A
sub-negN/A
log1p-defineN/A
log1p-lowering-log1p.f32N/A
neg-lowering-neg.f32N/A
sin-lowering-sin.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.4%
Simplified98.4%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f32N/A
associate-*r*N/A
*-commutativeN/A
unpow3N/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
Simplified98.3%
if 0.0900000036 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) Initial program 57.3%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3295.5%
Simplified95.5%
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
+-lowering-+.f32N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f3295.7%
Applied egg-rr95.7%
Final simplification97.7%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (* u2 (* 2.0 PI))))
(if (<= t_0 0.09000000357627869)
(*
(sqrt (- (log1p (- u1))))
(* u2 (* PI (+ 2.0 (* (* -1.3333333333333333 (* u2 u2)) (* PI PI))))))
(*
(sin t_0)
(sqrt
(*
u1
(+ 1.0 (* u1 (+ 0.5 (* u1 (+ 0.3333333333333333 (* u1 0.25))))))))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = u2 * (2.0f * ((float) M_PI));
float tmp;
if (t_0 <= 0.09000000357627869f) {
tmp = sqrtf(-log1pf(-u1)) * (u2 * (((float) M_PI) * (2.0f + ((-1.3333333333333333f * (u2 * u2)) * (((float) M_PI) * ((float) M_PI))))));
} else {
tmp = sinf(t_0) * sqrtf((u1 * (1.0f + (u1 * (0.5f + (u1 * (0.3333333333333333f + (u1 * 0.25f))))))));
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi))) tmp = Float32(0.0) if (t_0 <= Float32(0.09000000357627869)) tmp = Float32(sqrt(Float32(-log1p(Float32(-u1)))) * Float32(u2 * Float32(Float32(pi) * Float32(Float32(2.0) + Float32(Float32(Float32(-1.3333333333333333) * Float32(u2 * u2)) * Float32(Float32(pi) * Float32(pi))))))); else tmp = Float32(sin(t_0) * sqrt(Float32(u1 * Float32(Float32(1.0) + Float32(u1 * Float32(Float32(0.5) + Float32(u1 * Float32(Float32(0.3333333333333333) + Float32(u1 * Float32(0.25)))))))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := u2 \cdot \left(2 \cdot \pi\right)\\
\mathbf{if}\;t\_0 \leq 0.09000000357627869:\\
\;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(u2 \cdot \left(\pi \cdot \left(2 + \left(-1.3333333333333333 \cdot \left(u2 \cdot u2\right)\right) \cdot \left(\pi \cdot \pi\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin t\_0 \cdot \sqrt{u1 \cdot \left(1 + u1 \cdot \left(0.5 + u1 \cdot \left(0.3333333333333333 + u1 \cdot 0.25\right)\right)\right)}\\
\end{array}
\end{array}
if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 0.0900000036Initial program 59.6%
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
neg-lowering-neg.f32N/A
sub-negN/A
log1p-defineN/A
log1p-lowering-log1p.f32N/A
neg-lowering-neg.f32N/A
sin-lowering-sin.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.4%
Simplified98.4%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f32N/A
associate-*r*N/A
*-commutativeN/A
unpow3N/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
Simplified98.3%
if 0.0900000036 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) Initial program 57.3%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3295.5%
Simplified95.5%
Final simplification97.7%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (* u2 (* 2.0 PI))))
(if (<= t_0 0.09000000357627869)
(*
(sqrt (- (log1p (- u1))))
(* u2 (* PI (+ 2.0 (* (* -1.3333333333333333 (* u2 u2)) (* PI PI))))))
(*
(sin t_0)
(sqrt (* u1 (+ 1.0 (* u1 (+ 0.5 (* u1 0.3333333333333333))))))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = u2 * (2.0f * ((float) M_PI));
float tmp;
if (t_0 <= 0.09000000357627869f) {
tmp = sqrtf(-log1pf(-u1)) * (u2 * (((float) M_PI) * (2.0f + ((-1.3333333333333333f * (u2 * u2)) * (((float) M_PI) * ((float) M_PI))))));
} else {
tmp = sinf(t_0) * sqrtf((u1 * (1.0f + (u1 * (0.5f + (u1 * 0.3333333333333333f))))));
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi))) tmp = Float32(0.0) if (t_0 <= Float32(0.09000000357627869)) tmp = Float32(sqrt(Float32(-log1p(Float32(-u1)))) * Float32(u2 * Float32(Float32(pi) * Float32(Float32(2.0) + Float32(Float32(Float32(-1.3333333333333333) * Float32(u2 * u2)) * Float32(Float32(pi) * Float32(pi))))))); else tmp = Float32(sin(t_0) * sqrt(Float32(u1 * Float32(Float32(1.0) + Float32(u1 * Float32(Float32(0.5) + Float32(u1 * Float32(0.3333333333333333)))))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := u2 \cdot \left(2 \cdot \pi\right)\\
\mathbf{if}\;t\_0 \leq 0.09000000357627869:\\
\;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(u2 \cdot \left(\pi \cdot \left(2 + \left(-1.3333333333333333 \cdot \left(u2 \cdot u2\right)\right) \cdot \left(\pi \cdot \pi\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin t\_0 \cdot \sqrt{u1 \cdot \left(1 + u1 \cdot \left(0.5 + u1 \cdot 0.3333333333333333\right)\right)}\\
\end{array}
\end{array}
if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 0.0900000036Initial program 59.6%
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
neg-lowering-neg.f32N/A
sub-negN/A
log1p-defineN/A
log1p-lowering-log1p.f32N/A
neg-lowering-neg.f32N/A
sin-lowering-sin.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.4%
Simplified98.4%
Taylor expanded in u2 around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f32N/A
associate-*r*N/A
*-commutativeN/A
unpow3N/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
Simplified98.3%
if 0.0900000036 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) Initial program 57.3%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3293.9%
Simplified93.9%
Final simplification97.3%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (* u2 (* 2.0 PI))))
(if (<= t_0 0.0005000000237487257)
(* (sqrt (- (log1p (- u1)))) (* 2.0 (* PI u2)))
(*
(sin t_0)
(sqrt (* u1 (+ 1.0 (* u1 (+ 0.5 (* u1 0.3333333333333333))))))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = u2 * (2.0f * ((float) M_PI));
float tmp;
if (t_0 <= 0.0005000000237487257f) {
tmp = sqrtf(-log1pf(-u1)) * (2.0f * (((float) M_PI) * u2));
} else {
tmp = sinf(t_0) * sqrtf((u1 * (1.0f + (u1 * (0.5f + (u1 * 0.3333333333333333f))))));
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi))) tmp = Float32(0.0) if (t_0 <= Float32(0.0005000000237487257)) tmp = Float32(sqrt(Float32(-log1p(Float32(-u1)))) * Float32(Float32(2.0) * Float32(Float32(pi) * u2))); else tmp = Float32(sin(t_0) * sqrt(Float32(u1 * Float32(Float32(1.0) + Float32(u1 * Float32(Float32(0.5) + Float32(u1 * Float32(0.3333333333333333)))))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := u2 \cdot \left(2 \cdot \pi\right)\\
\mathbf{if}\;t\_0 \leq 0.0005000000237487257:\\
\;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(2 \cdot \left(\pi \cdot u2\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin t\_0 \cdot \sqrt{u1 \cdot \left(1 + u1 \cdot \left(0.5 + u1 \cdot 0.3333333333333333\right)\right)}\\
\end{array}
\end{array}
if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 5.00000024e-4Initial program 58.8%
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
neg-lowering-neg.f32N/A
sub-negN/A
log1p-defineN/A
log1p-lowering-log1p.f32N/A
neg-lowering-neg.f32N/A
sin-lowering-sin.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.4%
Simplified98.4%
Taylor expanded in u2 around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.4%
Simplified98.4%
if 5.00000024e-4 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) Initial program 59.5%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3293.1%
Simplified93.1%
Final simplification96.0%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (* u2 (* 2.0 PI))))
(if (<= t_0 0.0005000000237487257)
(* (sqrt (- (log1p (- u1)))) (* 2.0 (* PI u2)))
(* (sin t_0) (sqrt (* (* u1 u1) (+ 0.5 (/ 1.0 u1))))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = u2 * (2.0f * ((float) M_PI));
float tmp;
if (t_0 <= 0.0005000000237487257f) {
tmp = sqrtf(-log1pf(-u1)) * (2.0f * (((float) M_PI) * u2));
} else {
tmp = sinf(t_0) * sqrtf(((u1 * u1) * (0.5f + (1.0f / u1))));
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi))) tmp = Float32(0.0) if (t_0 <= Float32(0.0005000000237487257)) tmp = Float32(sqrt(Float32(-log1p(Float32(-u1)))) * Float32(Float32(2.0) * Float32(Float32(pi) * u2))); else tmp = Float32(sin(t_0) * sqrt(Float32(Float32(u1 * u1) * Float32(Float32(0.5) + Float32(Float32(1.0) / u1))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := u2 \cdot \left(2 \cdot \pi\right)\\
\mathbf{if}\;t\_0 \leq 0.0005000000237487257:\\
\;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(2 \cdot \left(\pi \cdot u2\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin t\_0 \cdot \sqrt{\left(u1 \cdot u1\right) \cdot \left(0.5 + \frac{1}{u1}\right)}\\
\end{array}
\end{array}
if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 5.00000024e-4Initial program 58.8%
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
neg-lowering-neg.f32N/A
sub-negN/A
log1p-defineN/A
log1p-lowering-log1p.f32N/A
neg-lowering-neg.f32N/A
sin-lowering-sin.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.4%
Simplified98.4%
Taylor expanded in u2 around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.4%
Simplified98.4%
if 5.00000024e-4 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) Initial program 59.5%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3289.6%
Simplified89.6%
Taylor expanded in u1 around inf
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f3289.7%
Simplified89.7%
Final simplification94.5%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (* u2 (* 2.0 PI))))
(if (<= t_0 0.0005000000237487257)
(* (sqrt (- (log1p (- u1)))) (* 2.0 (* PI u2)))
(* (sin t_0) (sqrt (+ u1 (* u1 (* u1 0.5))))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = u2 * (2.0f * ((float) M_PI));
float tmp;
if (t_0 <= 0.0005000000237487257f) {
tmp = sqrtf(-log1pf(-u1)) * (2.0f * (((float) M_PI) * u2));
} else {
tmp = sinf(t_0) * sqrtf((u1 + (u1 * (u1 * 0.5f))));
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi))) tmp = Float32(0.0) if (t_0 <= Float32(0.0005000000237487257)) tmp = Float32(sqrt(Float32(-log1p(Float32(-u1)))) * Float32(Float32(2.0) * Float32(Float32(pi) * u2))); else tmp = Float32(sin(t_0) * sqrt(Float32(u1 + Float32(u1 * Float32(u1 * Float32(0.5)))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := u2 \cdot \left(2 \cdot \pi\right)\\
\mathbf{if}\;t\_0 \leq 0.0005000000237487257:\\
\;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(2 \cdot \left(\pi \cdot u2\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin t\_0 \cdot \sqrt{u1 + u1 \cdot \left(u1 \cdot 0.5\right)}\\
\end{array}
\end{array}
if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 5.00000024e-4Initial program 58.8%
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
neg-lowering-neg.f32N/A
sub-negN/A
log1p-defineN/A
log1p-lowering-log1p.f32N/A
neg-lowering-neg.f32N/A
sin-lowering-sin.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.4%
Simplified98.4%
Taylor expanded in u2 around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.4%
Simplified98.4%
if 5.00000024e-4 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) Initial program 59.5%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3289.6%
Simplified89.6%
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f3289.7%
Applied egg-rr89.7%
Final simplification94.5%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (* u2 (* 2.0 PI))))
(if (<= t_0 0.0005000000237487257)
(* (sqrt (- (log1p (- u1)))) (* 2.0 (* PI u2)))
(* (sin t_0) (sqrt (* u1 (+ 1.0 (* u1 0.5))))))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = u2 * (2.0f * ((float) M_PI));
float tmp;
if (t_0 <= 0.0005000000237487257f) {
tmp = sqrtf(-log1pf(-u1)) * (2.0f * (((float) M_PI) * u2));
} else {
tmp = sinf(t_0) * sqrtf((u1 * (1.0f + (u1 * 0.5f))));
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi))) tmp = Float32(0.0) if (t_0 <= Float32(0.0005000000237487257)) tmp = Float32(sqrt(Float32(-log1p(Float32(-u1)))) * Float32(Float32(2.0) * Float32(Float32(pi) * u2))); else tmp = Float32(sin(t_0) * sqrt(Float32(u1 * Float32(Float32(1.0) + Float32(u1 * Float32(0.5)))))); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := u2 \cdot \left(2 \cdot \pi\right)\\
\mathbf{if}\;t\_0 \leq 0.0005000000237487257:\\
\;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(2 \cdot \left(\pi \cdot u2\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin t\_0 \cdot \sqrt{u1 \cdot \left(1 + u1 \cdot 0.5\right)}\\
\end{array}
\end{array}
if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 5.00000024e-4Initial program 58.8%
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
neg-lowering-neg.f32N/A
sub-negN/A
log1p-defineN/A
log1p-lowering-log1p.f32N/A
neg-lowering-neg.f32N/A
sin-lowering-sin.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.4%
Simplified98.4%
Taylor expanded in u2 around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.4%
Simplified98.4%
if 5.00000024e-4 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) Initial program 59.5%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3289.6%
Simplified89.6%
Final simplification94.4%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (* u2 (* 2.0 PI))))
(if (<= t_0 0.15000000596046448)
(*
(sqrt
(* u1 (+ 1.0 (* u1 (+ 0.5 (* u1 (+ 0.3333333333333333 (* u1 0.25))))))))
(*
u2
(+ (* 2.0 PI) (* (* -1.3333333333333333 (* u2 u2)) (* PI (* PI PI))))))
(* (sin t_0) (pow (* u1 u1) 0.25)))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = u2 * (2.0f * ((float) M_PI));
float tmp;
if (t_0 <= 0.15000000596046448f) {
tmp = sqrtf((u1 * (1.0f + (u1 * (0.5f + (u1 * (0.3333333333333333f + (u1 * 0.25f)))))))) * (u2 * ((2.0f * ((float) M_PI)) + ((-1.3333333333333333f * (u2 * u2)) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI))))));
} else {
tmp = sinf(t_0) * powf((u1 * u1), 0.25f);
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi))) tmp = Float32(0.0) if (t_0 <= Float32(0.15000000596046448)) tmp = Float32(sqrt(Float32(u1 * Float32(Float32(1.0) + Float32(u1 * Float32(Float32(0.5) + Float32(u1 * Float32(Float32(0.3333333333333333) + Float32(u1 * Float32(0.25))))))))) * Float32(u2 * Float32(Float32(Float32(2.0) * Float32(pi)) + Float32(Float32(Float32(-1.3333333333333333) * Float32(u2 * u2)) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))))))); else tmp = Float32(sin(t_0) * (Float32(u1 * u1) ^ Float32(0.25))); end return tmp end
function tmp_2 = code(cosTheta_i, u1, u2) t_0 = u2 * (single(2.0) * single(pi)); tmp = single(0.0); if (t_0 <= single(0.15000000596046448)) tmp = sqrt((u1 * (single(1.0) + (u1 * (single(0.5) + (u1 * (single(0.3333333333333333) + (u1 * single(0.25))))))))) * (u2 * ((single(2.0) * single(pi)) + ((single(-1.3333333333333333) * (u2 * u2)) * (single(pi) * (single(pi) * single(pi)))))); else tmp = sin(t_0) * ((u1 * u1) ^ single(0.25)); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := u2 \cdot \left(2 \cdot \pi\right)\\
\mathbf{if}\;t\_0 \leq 0.15000000596046448:\\
\;\;\;\;\sqrt{u1 \cdot \left(1 + u1 \cdot \left(0.5 + u1 \cdot \left(0.3333333333333333 + u1 \cdot 0.25\right)\right)\right)} \cdot \left(u2 \cdot \left(2 \cdot \pi + \left(-1.3333333333333333 \cdot \left(u2 \cdot u2\right)\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin t\_0 \cdot {\left(u1 \cdot u1\right)}^{0.25}\\
\end{array}
\end{array}
if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 0.150000006Initial program 60.0%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3293.7%
Simplified93.7%
Taylor expanded in u2 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
cube-multN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3293.3%
Simplified93.3%
if 0.150000006 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) Initial program 55.2%
Taylor expanded in u1 around 0
Simplified78.7%
pow1/2N/A
metadata-evalN/A
pow-prod-upN/A
pow-prod-downN/A
pow-lowering-pow.f32N/A
*-lowering-*.f3278.8%
Applied egg-rr78.8%
Final simplification90.6%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(let* ((t_0 (* u2 (* 2.0 PI))))
(if (<= t_0 0.15000000596046448)
(*
(sqrt
(* u1 (+ 1.0 (* u1 (+ 0.5 (* u1 (+ 0.3333333333333333 (* u1 0.25))))))))
(*
u2
(+ (* 2.0 PI) (* (* -1.3333333333333333 (* u2 u2)) (* PI (* PI PI))))))
(* (sin t_0) (sqrt u1)))))
float code(float cosTheta_i, float u1, float u2) {
float t_0 = u2 * (2.0f * ((float) M_PI));
float tmp;
if (t_0 <= 0.15000000596046448f) {
tmp = sqrtf((u1 * (1.0f + (u1 * (0.5f + (u1 * (0.3333333333333333f + (u1 * 0.25f)))))))) * (u2 * ((2.0f * ((float) M_PI)) + ((-1.3333333333333333f * (u2 * u2)) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI))))));
} else {
tmp = sinf(t_0) * sqrtf(u1);
}
return tmp;
}
function code(cosTheta_i, u1, u2) t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi))) tmp = Float32(0.0) if (t_0 <= Float32(0.15000000596046448)) tmp = Float32(sqrt(Float32(u1 * Float32(Float32(1.0) + Float32(u1 * Float32(Float32(0.5) + Float32(u1 * Float32(Float32(0.3333333333333333) + Float32(u1 * Float32(0.25))))))))) * Float32(u2 * Float32(Float32(Float32(2.0) * Float32(pi)) + Float32(Float32(Float32(-1.3333333333333333) * Float32(u2 * u2)) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))))))); else tmp = Float32(sin(t_0) * sqrt(u1)); end return tmp end
function tmp_2 = code(cosTheta_i, u1, u2) t_0 = u2 * (single(2.0) * single(pi)); tmp = single(0.0); if (t_0 <= single(0.15000000596046448)) tmp = sqrt((u1 * (single(1.0) + (u1 * (single(0.5) + (u1 * (single(0.3333333333333333) + (u1 * single(0.25))))))))) * (u2 * ((single(2.0) * single(pi)) + ((single(-1.3333333333333333) * (u2 * u2)) * (single(pi) * (single(pi) * single(pi)))))); else tmp = sin(t_0) * sqrt(u1); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := u2 \cdot \left(2 \cdot \pi\right)\\
\mathbf{if}\;t\_0 \leq 0.15000000596046448:\\
\;\;\;\;\sqrt{u1 \cdot \left(1 + u1 \cdot \left(0.5 + u1 \cdot \left(0.3333333333333333 + u1 \cdot 0.25\right)\right)\right)} \cdot \left(u2 \cdot \left(2 \cdot \pi + \left(-1.3333333333333333 \cdot \left(u2 \cdot u2\right)\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sin t\_0 \cdot \sqrt{u1}\\
\end{array}
\end{array}
if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 0.150000006Initial program 60.0%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3293.7%
Simplified93.7%
Taylor expanded in u2 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
cube-multN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3293.3%
Simplified93.3%
if 0.150000006 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) Initial program 55.2%
Taylor expanded in u1 around 0
Simplified78.7%
Final simplification90.5%
(FPCore (cosTheta_i u1 u2) :precision binary32 (* (sqrt (* u1 (+ 1.0 (* u1 (+ 0.5 (* u1 (+ 0.3333333333333333 (* u1 0.25)))))))) (* u2 (+ (* 2.0 PI) (* (* -1.3333333333333333 (* u2 u2)) (* PI (* PI PI)))))))
float code(float cosTheta_i, float u1, float u2) {
return sqrtf((u1 * (1.0f + (u1 * (0.5f + (u1 * (0.3333333333333333f + (u1 * 0.25f)))))))) * (u2 * ((2.0f * ((float) M_PI)) + ((-1.3333333333333333f * (u2 * u2)) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI))))));
}
function code(cosTheta_i, u1, u2) return Float32(sqrt(Float32(u1 * Float32(Float32(1.0) + Float32(u1 * Float32(Float32(0.5) + Float32(u1 * Float32(Float32(0.3333333333333333) + Float32(u1 * Float32(0.25))))))))) * Float32(u2 * Float32(Float32(Float32(2.0) * Float32(pi)) + Float32(Float32(Float32(-1.3333333333333333) * Float32(u2 * u2)) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))))))) end
function tmp = code(cosTheta_i, u1, u2) tmp = sqrt((u1 * (single(1.0) + (u1 * (single(0.5) + (u1 * (single(0.3333333333333333) + (u1 * single(0.25))))))))) * (u2 * ((single(2.0) * single(pi)) + ((single(-1.3333333333333333) * (u2 * u2)) * (single(pi) * (single(pi) * single(pi)))))); end
\begin{array}{l}
\\
\sqrt{u1 \cdot \left(1 + u1 \cdot \left(0.5 + u1 \cdot \left(0.3333333333333333 + u1 \cdot 0.25\right)\right)\right)} \cdot \left(u2 \cdot \left(2 \cdot \pi + \left(-1.3333333333333333 \cdot \left(u2 \cdot u2\right)\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right)\right)
\end{array}
Initial program 59.1%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3294.0%
Simplified94.0%
Taylor expanded in u2 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
cube-multN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3285.1%
Simplified85.1%
Final simplification85.1%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(if (<= (* u2 (* 2.0 PI)) 0.017999999225139618)
(*
(sqrt
(* u1 (+ 1.0 (* u1 (+ 0.5 (* u1 (+ 0.3333333333333333 (* u1 0.25))))))))
(* 2.0 (* PI u2)))
(*
(*
u2
(+ (* 2.0 PI) (* (* -1.3333333333333333 (* u2 u2)) (* PI (* PI PI)))))
(sqrt u1))))
float code(float cosTheta_i, float u1, float u2) {
float tmp;
if ((u2 * (2.0f * ((float) M_PI))) <= 0.017999999225139618f) {
tmp = sqrtf((u1 * (1.0f + (u1 * (0.5f + (u1 * (0.3333333333333333f + (u1 * 0.25f)))))))) * (2.0f * (((float) M_PI) * u2));
} else {
tmp = (u2 * ((2.0f * ((float) M_PI)) + ((-1.3333333333333333f * (u2 * u2)) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))))) * sqrtf(u1);
}
return tmp;
}
function code(cosTheta_i, u1, u2) tmp = Float32(0.0) if (Float32(u2 * Float32(Float32(2.0) * Float32(pi))) <= Float32(0.017999999225139618)) tmp = Float32(sqrt(Float32(u1 * Float32(Float32(1.0) + Float32(u1 * Float32(Float32(0.5) + Float32(u1 * Float32(Float32(0.3333333333333333) + Float32(u1 * Float32(0.25))))))))) * Float32(Float32(2.0) * Float32(Float32(pi) * u2))); else tmp = Float32(Float32(u2 * Float32(Float32(Float32(2.0) * Float32(pi)) + Float32(Float32(Float32(-1.3333333333333333) * Float32(u2 * u2)) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))))) * sqrt(u1)); end return tmp end
function tmp_2 = code(cosTheta_i, u1, u2) tmp = single(0.0); if ((u2 * (single(2.0) * single(pi))) <= single(0.017999999225139618)) tmp = sqrt((u1 * (single(1.0) + (u1 * (single(0.5) + (u1 * (single(0.3333333333333333) + (u1 * single(0.25))))))))) * (single(2.0) * (single(pi) * u2)); else tmp = (u2 * ((single(2.0) * single(pi)) + ((single(-1.3333333333333333) * (u2 * u2)) * (single(pi) * (single(pi) * single(pi)))))) * sqrt(u1); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;u2 \cdot \left(2 \cdot \pi\right) \leq 0.017999999225139618:\\
\;\;\;\;\sqrt{u1 \cdot \left(1 + u1 \cdot \left(0.5 + u1 \cdot \left(0.3333333333333333 + u1 \cdot 0.25\right)\right)\right)} \cdot \left(2 \cdot \left(\pi \cdot u2\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(u2 \cdot \left(2 \cdot \pi + \left(-1.3333333333333333 \cdot \left(u2 \cdot u2\right)\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right)\right) \cdot \sqrt{u1}\\
\end{array}
\end{array}
if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 0.0179999992Initial program 58.6%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3293.6%
Simplified93.6%
Taylor expanded in u2 around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3291.6%
Simplified91.6%
if 0.0179999992 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) Initial program 60.2%
Taylor expanded in u1 around 0
Simplified76.3%
Taylor expanded in u2 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
cube-multN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3256.1%
Simplified56.1%
Final simplification81.0%
(FPCore (cosTheta_i u1 u2)
:precision binary32
(if (<= (* u2 (* 2.0 PI)) 0.009999999776482582)
(*
(sqrt (* u1 (+ 1.0 (* u1 (+ 0.5 (* u1 0.3333333333333333))))))
(* 2.0 (* PI u2)))
(*
(*
u2
(+ (* 2.0 PI) (* (* -1.3333333333333333 (* u2 u2)) (* PI (* PI PI)))))
(sqrt u1))))
float code(float cosTheta_i, float u1, float u2) {
float tmp;
if ((u2 * (2.0f * ((float) M_PI))) <= 0.009999999776482582f) {
tmp = sqrtf((u1 * (1.0f + (u1 * (0.5f + (u1 * 0.3333333333333333f)))))) * (2.0f * (((float) M_PI) * u2));
} else {
tmp = (u2 * ((2.0f * ((float) M_PI)) + ((-1.3333333333333333f * (u2 * u2)) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI)))))) * sqrtf(u1);
}
return tmp;
}
function code(cosTheta_i, u1, u2) tmp = Float32(0.0) if (Float32(u2 * Float32(Float32(2.0) * Float32(pi))) <= Float32(0.009999999776482582)) tmp = Float32(sqrt(Float32(u1 * Float32(Float32(1.0) + Float32(u1 * Float32(Float32(0.5) + Float32(u1 * Float32(0.3333333333333333))))))) * Float32(Float32(2.0) * Float32(Float32(pi) * u2))); else tmp = Float32(Float32(u2 * Float32(Float32(Float32(2.0) * Float32(pi)) + Float32(Float32(Float32(-1.3333333333333333) * Float32(u2 * u2)) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi)))))) * sqrt(u1)); end return tmp end
function tmp_2 = code(cosTheta_i, u1, u2) tmp = single(0.0); if ((u2 * (single(2.0) * single(pi))) <= single(0.009999999776482582)) tmp = sqrt((u1 * (single(1.0) + (u1 * (single(0.5) + (u1 * single(0.3333333333333333))))))) * (single(2.0) * (single(pi) * u2)); else tmp = (u2 * ((single(2.0) * single(pi)) + ((single(-1.3333333333333333) * (u2 * u2)) * (single(pi) * (single(pi) * single(pi)))))) * sqrt(u1); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;u2 \cdot \left(2 \cdot \pi\right) \leq 0.009999999776482582:\\
\;\;\;\;\sqrt{u1 \cdot \left(1 + u1 \cdot \left(0.5 + u1 \cdot 0.3333333333333333\right)\right)} \cdot \left(2 \cdot \left(\pi \cdot u2\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(u2 \cdot \left(2 \cdot \pi + \left(-1.3333333333333333 \cdot \left(u2 \cdot u2\right)\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right)\right) \cdot \sqrt{u1}\\
\end{array}
\end{array}
if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 0.00999999978Initial program 58.5%
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
neg-lowering-neg.f32N/A
sub-negN/A
log1p-defineN/A
log1p-lowering-log1p.f32N/A
neg-lowering-neg.f32N/A
sin-lowering-sin.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.4%
Simplified98.4%
Taylor expanded in u2 around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3296.4%
Simplified96.4%
Taylor expanded in u1 around 0
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-+r+N/A
*-commutativeN/A
distribute-lft-inN/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3290.6%
Simplified90.6%
if 0.00999999978 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) Initial program 60.5%
Taylor expanded in u1 around 0
Simplified76.0%
Taylor expanded in u2 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
cube-multN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3257.1%
Simplified57.1%
Final simplification80.0%
(FPCore (cosTheta_i u1 u2) :precision binary32 (* (sqrt (* u1 (+ 1.0 (* u1 0.5)))) (* u2 (+ (* 2.0 PI) (* (* -1.3333333333333333 (* u2 u2)) (* PI (* PI PI)))))))
float code(float cosTheta_i, float u1, float u2) {
return sqrtf((u1 * (1.0f + (u1 * 0.5f)))) * (u2 * ((2.0f * ((float) M_PI)) + ((-1.3333333333333333f * (u2 * u2)) * (((float) M_PI) * (((float) M_PI) * ((float) M_PI))))));
}
function code(cosTheta_i, u1, u2) return Float32(sqrt(Float32(u1 * Float32(Float32(1.0) + Float32(u1 * Float32(0.5))))) * Float32(u2 * Float32(Float32(Float32(2.0) * Float32(pi)) + Float32(Float32(Float32(-1.3333333333333333) * Float32(u2 * u2)) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(pi))))))) end
function tmp = code(cosTheta_i, u1, u2) tmp = sqrt((u1 * (single(1.0) + (u1 * single(0.5))))) * (u2 * ((single(2.0) * single(pi)) + ((single(-1.3333333333333333) * (u2 * u2)) * (single(pi) * (single(pi) * single(pi)))))); end
\begin{array}{l}
\\
\sqrt{u1 \cdot \left(1 + u1 \cdot 0.5\right)} \cdot \left(u2 \cdot \left(2 \cdot \pi + \left(-1.3333333333333333 \cdot \left(u2 \cdot u2\right)\right) \cdot \left(\pi \cdot \left(\pi \cdot \pi\right)\right)\right)\right)
\end{array}
Initial program 59.1%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3288.7%
Simplified88.7%
Taylor expanded in u2 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
cube-multN/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3280.8%
Simplified80.8%
Final simplification80.8%
(FPCore (cosTheta_i u1 u2) :precision binary32 (* (sqrt (* u1 (+ 1.0 (* u1 (+ 0.5 (* u1 0.3333333333333333)))))) (* 2.0 (* PI u2))))
float code(float cosTheta_i, float u1, float u2) {
return sqrtf((u1 * (1.0f + (u1 * (0.5f + (u1 * 0.3333333333333333f)))))) * (2.0f * (((float) M_PI) * u2));
}
function code(cosTheta_i, u1, u2) return Float32(sqrt(Float32(u1 * Float32(Float32(1.0) + Float32(u1 * Float32(Float32(0.5) + Float32(u1 * Float32(0.3333333333333333))))))) * Float32(Float32(2.0) * Float32(Float32(pi) * u2))) end
function tmp = code(cosTheta_i, u1, u2) tmp = sqrt((u1 * (single(1.0) + (u1 * (single(0.5) + (u1 * single(0.3333333333333333))))))) * (single(2.0) * (single(pi) * u2)); end
\begin{array}{l}
\\
\sqrt{u1 \cdot \left(1 + u1 \cdot \left(0.5 + u1 \cdot 0.3333333333333333\right)\right)} \cdot \left(2 \cdot \left(\pi \cdot u2\right)\right)
\end{array}
Initial program 59.1%
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
neg-lowering-neg.f32N/A
sub-negN/A
log1p-defineN/A
log1p-lowering-log1p.f32N/A
neg-lowering-neg.f32N/A
sin-lowering-sin.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3298.3%
Simplified98.3%
Taylor expanded in u2 around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3280.2%
Simplified80.2%
Taylor expanded in u1 around 0
distribute-lft-inN/A
*-commutativeN/A
associate-+r+N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
associate-+r+N/A
*-commutativeN/A
distribute-lft-inN/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3276.1%
Simplified76.1%
Final simplification76.1%
(FPCore (cosTheta_i u1 u2) :precision binary32 (* (* 2.0 (* PI u2)) (sqrt (* u1 (+ 1.0 (* u1 0.5))))))
float code(float cosTheta_i, float u1, float u2) {
return (2.0f * (((float) M_PI) * u2)) * sqrtf((u1 * (1.0f + (u1 * 0.5f))));
}
function code(cosTheta_i, u1, u2) return Float32(Float32(Float32(2.0) * Float32(Float32(pi) * u2)) * sqrt(Float32(u1 * Float32(Float32(1.0) + Float32(u1 * Float32(0.5)))))) end
function tmp = code(cosTheta_i, u1, u2) tmp = (single(2.0) * (single(pi) * u2)) * sqrt((u1 * (single(1.0) + (u1 * single(0.5))))); end
\begin{array}{l}
\\
\left(2 \cdot \left(\pi \cdot u2\right)\right) \cdot \sqrt{u1 \cdot \left(1 + u1 \cdot 0.5\right)}
\end{array}
Initial program 59.1%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f3288.7%
Simplified88.7%
Taylor expanded in u2 around 0
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3273.7%
Simplified73.7%
Final simplification73.7%
(FPCore (cosTheta_i u1 u2) :precision binary32 (* u2 (* PI (* 2.0 (pow u1 0.5)))))
float code(float cosTheta_i, float u1, float u2) {
return u2 * (((float) M_PI) * (2.0f * powf(u1, 0.5f)));
}
function code(cosTheta_i, u1, u2) return Float32(u2 * Float32(Float32(pi) * Float32(Float32(2.0) * (u1 ^ Float32(0.5))))) end
function tmp = code(cosTheta_i, u1, u2) tmp = u2 * (single(pi) * (single(2.0) * (u1 ^ single(0.5)))); end
\begin{array}{l}
\\
u2 \cdot \left(\pi \cdot \left(2 \cdot {u1}^{0.5}\right)\right)
\end{array}
Initial program 59.1%
Taylor expanded in u1 around 0
Simplified76.4%
Taylor expanded in u2 around 0
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3265.9%
Simplified65.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
pow1/2N/A
pow-lowering-pow.f32N/A
PI-lowering-PI.f3266.0%
Applied egg-rr66.0%
Final simplification66.0%
(FPCore (cosTheta_i u1 u2) :precision binary32 (* 2.0 (* PI (* u2 (sqrt u1)))))
float code(float cosTheta_i, float u1, float u2) {
return 2.0f * (((float) M_PI) * (u2 * sqrtf(u1)));
}
function code(cosTheta_i, u1, u2) return Float32(Float32(2.0) * Float32(Float32(pi) * Float32(u2 * sqrt(u1)))) end
function tmp = code(cosTheta_i, u1, u2) tmp = single(2.0) * (single(pi) * (u2 * sqrt(u1))); end
\begin{array}{l}
\\
2 \cdot \left(\pi \cdot \left(u2 \cdot \sqrt{u1}\right)\right)
\end{array}
Initial program 59.1%
Taylor expanded in u1 around 0
Simplified76.4%
Taylor expanded in u2 around 0
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3265.9%
Simplified65.9%
Taylor expanded in u1 around 0
*-lowering-*.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
PI-lowering-PI.f3266.0%
Simplified66.0%
Final simplification66.0%
herbie shell --seed 2024152
(FPCore (cosTheta_i u1 u2)
:name "Beckmann Sample, near normal, slope_y"
: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)))) (sin (* (* 2.0 PI) u2))))