
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(+ 1.0 c)
(*
(* (/ 1.0 (sqrt PI)) (/ (sqrt (- (- 1.0 cosTheta) cosTheta)) cosTheta))
(exp (* (- cosTheta) cosTheta))))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + (((1.0f / sqrtf(((float) M_PI))) * (sqrtf(((1.0f - cosTheta) - cosTheta)) / cosTheta)) * expf((-cosTheta * cosTheta))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(Float32(Float32(1.0) / sqrt(Float32(pi))) * Float32(sqrt(Float32(Float32(Float32(1.0) - cosTheta) - cosTheta)) / cosTheta)) * exp(Float32(Float32(-cosTheta) * cosTheta))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / ((single(1.0) + c) + (((single(1.0) / sqrt(single(pi))) * (sqrt(((single(1.0) - cosTheta) - cosTheta)) / cosTheta)) * exp((-cosTheta * cosTheta)))); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \left(\frac{1}{\sqrt{\pi}} \cdot \frac{\sqrt{\left(1 - cosTheta\right) - cosTheta}}{cosTheta}\right) \cdot e^{\left(-cosTheta\right) \cdot cosTheta}}
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 18 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(+ 1.0 c)
(*
(* (/ 1.0 (sqrt PI)) (/ (sqrt (- (- 1.0 cosTheta) cosTheta)) cosTheta))
(exp (* (- cosTheta) cosTheta))))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + (((1.0f / sqrtf(((float) M_PI))) * (sqrtf(((1.0f - cosTheta) - cosTheta)) / cosTheta)) * expf((-cosTheta * cosTheta))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(Float32(Float32(1.0) / sqrt(Float32(pi))) * Float32(sqrt(Float32(Float32(Float32(1.0) - cosTheta) - cosTheta)) / cosTheta)) * exp(Float32(Float32(-cosTheta) * cosTheta))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / ((single(1.0) + c) + (((single(1.0) / sqrt(single(pi))) * (sqrt(((single(1.0) - cosTheta) - cosTheta)) / cosTheta)) * exp((-cosTheta * cosTheta)))); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \left(\frac{1}{\sqrt{\pi}} \cdot \frac{\sqrt{\left(1 - cosTheta\right) - cosTheta}}{cosTheta}\right) \cdot e^{\left(-cosTheta\right) \cdot cosTheta}}
\end{array}
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(+ 1.0 c)
(/
(pow (+ 1.0 (* cosTheta -2.0)) 0.5)
(* (exp (* cosTheta cosTheta)) (* cosTheta (pow PI 0.5)))))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + (powf((1.0f + (cosTheta * -2.0f)), 0.5f) / (expf((cosTheta * cosTheta)) * (cosTheta * powf(((float) M_PI), 0.5f)))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32((Float32(Float32(1.0) + Float32(cosTheta * Float32(-2.0))) ^ Float32(0.5)) / Float32(exp(Float32(cosTheta * cosTheta)) * Float32(cosTheta * (Float32(pi) ^ Float32(0.5))))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / ((single(1.0) + c) + (((single(1.0) + (cosTheta * single(-2.0))) ^ single(0.5)) / (exp((cosTheta * cosTheta)) * (cosTheta * (single(pi) ^ single(0.5)))))); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \frac{{\left(1 + cosTheta \cdot -2\right)}^{0.5}}{e^{cosTheta \cdot cosTheta} \cdot \left(cosTheta \cdot {\pi}^{0.5}\right)}}
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
associate-/l/N/A
/-lowering-/.f32N/A
pow1/2N/A
pow-lowering-pow.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
pow1/2N/A
pow-lowering-pow.f32N/A
PI-lowering-PI.f3298.6%
Applied egg-rr98.6%
Final simplification98.6%
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(+ 1.0 c)
(/
1.0
(/
(pow (/ (+ 1.0 (* cosTheta -2.0)) PI) -0.5)
(/ (/ 1.0 cosTheta) (exp (* cosTheta cosTheta))))))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + (1.0f / (powf(((1.0f + (cosTheta * -2.0f)) / ((float) M_PI)), -0.5f) / ((1.0f / cosTheta) / expf((cosTheta * cosTheta))))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(1.0) / Float32((Float32(Float32(Float32(1.0) + Float32(cosTheta * Float32(-2.0))) / Float32(pi)) ^ Float32(-0.5)) / Float32(Float32(Float32(1.0) / cosTheta) / exp(Float32(cosTheta * cosTheta))))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / ((single(1.0) + c) + (single(1.0) / ((((single(1.0) + (cosTheta * single(-2.0))) / single(pi)) ^ single(-0.5)) / ((single(1.0) / cosTheta) / exp((cosTheta * cosTheta)))))); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \frac{1}{\frac{{\left(\frac{1 + cosTheta \cdot -2}{\pi}\right)}^{-0.5}}{\frac{\frac{1}{cosTheta}}{e^{cosTheta \cdot cosTheta}}}}}
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
clear-numN/A
/-lowering-/.f32N/A
div-invN/A
times-fracN/A
*-lowering-*.f32N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f3298.5%
Applied egg-rr98.5%
clear-numN/A
un-div-invN/A
associate-/r*N/A
/-lowering-/.f32N/A
pow1/2N/A
clear-numN/A
inv-powN/A
pow-powN/A
pow-lowering-pow.f32N/A
/-lowering-/.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
metadata-evalN/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
*-lowering-*.f3298.6%
Applied egg-rr98.6%
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(+ 1.0 c)
(/
1.0
(*
(sqrt (/ PI (+ 1.0 (* cosTheta -2.0))))
(*
cosTheta
(+
1.0
(*
(* cosTheta cosTheta)
(+
1.0
(*
cosTheta
(*
cosTheta
(+ 0.5 (* (* cosTheta cosTheta) 0.16666666666666666)))))))))))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + (1.0f / (sqrtf((((float) M_PI) / (1.0f + (cosTheta * -2.0f)))) * (cosTheta * (1.0f + ((cosTheta * cosTheta) * (1.0f + (cosTheta * (cosTheta * (0.5f + ((cosTheta * cosTheta) * 0.16666666666666666f)))))))))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(1.0) / Float32(sqrt(Float32(Float32(pi) / Float32(Float32(1.0) + Float32(cosTheta * Float32(-2.0))))) * Float32(cosTheta * Float32(Float32(1.0) + Float32(Float32(cosTheta * cosTheta) * Float32(Float32(1.0) + Float32(cosTheta * Float32(cosTheta * Float32(Float32(0.5) + Float32(Float32(cosTheta * cosTheta) * Float32(0.16666666666666666))))))))))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / ((single(1.0) + c) + (single(1.0) / (sqrt((single(pi) / (single(1.0) + (cosTheta * single(-2.0))))) * (cosTheta * (single(1.0) + ((cosTheta * cosTheta) * (single(1.0) + (cosTheta * (cosTheta * (single(0.5) + ((cosTheta * cosTheta) * single(0.16666666666666666)))))))))))); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \frac{1}{\sqrt{\frac{\pi}{1 + cosTheta \cdot -2}} \cdot \left(cosTheta \cdot \left(1 + \left(cosTheta \cdot cosTheta\right) \cdot \left(1 + cosTheta \cdot \left(cosTheta \cdot \left(0.5 + \left(cosTheta \cdot cosTheta\right) \cdot 0.16666666666666666\right)\right)\right)\right)\right)}}
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
clear-numN/A
/-lowering-/.f32N/A
div-invN/A
times-fracN/A
*-lowering-*.f32N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f3298.5%
Applied egg-rr98.5%
Taylor expanded in cosTheta around 0
*-lowering-*.f32N/A
distribute-lft-inN/A
*-rgt-identityN/A
+-lowering-+.f32N/A
*-rgt-identityN/A
distribute-lft-inN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
Simplified98.5%
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
1.0
(/
1.0
(*
(sqrt (/ PI (+ 1.0 (* cosTheta -2.0))))
(*
cosTheta
(+
1.0
(*
(* cosTheta cosTheta)
(+
1.0
(*
cosTheta
(*
cosTheta
(+ 0.5 (* (* cosTheta cosTheta) 0.16666666666666666)))))))))))))
float code(float cosTheta, float c) {
return 1.0f / (1.0f + (1.0f / (sqrtf((((float) M_PI) / (1.0f + (cosTheta * -2.0f)))) * (cosTheta * (1.0f + ((cosTheta * cosTheta) * (1.0f + (cosTheta * (cosTheta * (0.5f + ((cosTheta * cosTheta) * 0.16666666666666666f)))))))))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(1.0) + Float32(Float32(1.0) / Float32(sqrt(Float32(Float32(pi) / Float32(Float32(1.0) + Float32(cosTheta * Float32(-2.0))))) * Float32(cosTheta * Float32(Float32(1.0) + Float32(Float32(cosTheta * cosTheta) * Float32(Float32(1.0) + Float32(cosTheta * Float32(cosTheta * Float32(Float32(0.5) + Float32(Float32(cosTheta * cosTheta) * Float32(0.16666666666666666))))))))))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / (single(1.0) + (single(1.0) / (sqrt((single(pi) / (single(1.0) + (cosTheta * single(-2.0))))) * (cosTheta * (single(1.0) + ((cosTheta * cosTheta) * (single(1.0) + (cosTheta * (cosTheta * (single(0.5) + ((cosTheta * cosTheta) * single(0.16666666666666666)))))))))))); end
\begin{array}{l}
\\
\frac{1}{1 + \frac{1}{\sqrt{\frac{\pi}{1 + cosTheta \cdot -2}} \cdot \left(cosTheta \cdot \left(1 + \left(cosTheta \cdot cosTheta\right) \cdot \left(1 + cosTheta \cdot \left(cosTheta \cdot \left(0.5 + \left(cosTheta \cdot cosTheta\right) \cdot 0.16666666666666666\right)\right)\right)\right)\right)}}
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
clear-numN/A
/-lowering-/.f32N/A
div-invN/A
times-fracN/A
*-lowering-*.f32N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f3298.5%
Applied egg-rr98.5%
Taylor expanded in cosTheta around 0
*-lowering-*.f32N/A
distribute-lft-inN/A
*-rgt-identityN/A
+-lowering-+.f32N/A
*-rgt-identityN/A
distribute-lft-inN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
Simplified98.5%
Taylor expanded in c around 0
Simplified98.3%
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(+ 1.0 c)
(/
1.0
(*
(sqrt (/ PI (+ 1.0 (* cosTheta -2.0))))
(*
cosTheta
(+
1.0
(* (* cosTheta cosTheta) (+ 1.0 (* cosTheta (* cosTheta 0.5)))))))))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + (1.0f / (sqrtf((((float) M_PI) / (1.0f + (cosTheta * -2.0f)))) * (cosTheta * (1.0f + ((cosTheta * cosTheta) * (1.0f + (cosTheta * (cosTheta * 0.5f)))))))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(1.0) / Float32(sqrt(Float32(Float32(pi) / Float32(Float32(1.0) + Float32(cosTheta * Float32(-2.0))))) * Float32(cosTheta * Float32(Float32(1.0) + Float32(Float32(cosTheta * cosTheta) * Float32(Float32(1.0) + Float32(cosTheta * Float32(cosTheta * Float32(0.5))))))))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / ((single(1.0) + c) + (single(1.0) / (sqrt((single(pi) / (single(1.0) + (cosTheta * single(-2.0))))) * (cosTheta * (single(1.0) + ((cosTheta * cosTheta) * (single(1.0) + (cosTheta * (cosTheta * single(0.5)))))))))); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \frac{1}{\sqrt{\frac{\pi}{1 + cosTheta \cdot -2}} \cdot \left(cosTheta \cdot \left(1 + \left(cosTheta \cdot cosTheta\right) \cdot \left(1 + cosTheta \cdot \left(cosTheta \cdot 0.5\right)\right)\right)\right)}}
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
clear-numN/A
/-lowering-/.f32N/A
div-invN/A
times-fracN/A
*-lowering-*.f32N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f3298.5%
Applied egg-rr98.5%
Taylor expanded in cosTheta around 0
*-lowering-*.f32N/A
distribute-rgt-inN/A
*-lft-identityN/A
*-commutativeN/A
associate-*r*N/A
+-lowering-+.f32N/A
*-rgt-identityN/A
distribute-lft-inN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f3298.3%
Simplified98.3%
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(+ 1.0 c)
(/
1.0
(*
(sqrt (/ PI (+ 1.0 (* cosTheta -2.0))))
(* cosTheta (+ 1.0 (* cosTheta cosTheta))))))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + (1.0f / (sqrtf((((float) M_PI) / (1.0f + (cosTheta * -2.0f)))) * (cosTheta * (1.0f + (cosTheta * cosTheta))))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(1.0) / Float32(sqrt(Float32(Float32(pi) / Float32(Float32(1.0) + Float32(cosTheta * Float32(-2.0))))) * Float32(cosTheta * Float32(Float32(1.0) + Float32(cosTheta * cosTheta))))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / ((single(1.0) + c) + (single(1.0) / (sqrt((single(pi) / (single(1.0) + (cosTheta * single(-2.0))))) * (cosTheta * (single(1.0) + (cosTheta * cosTheta)))))); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \frac{1}{\sqrt{\frac{\pi}{1 + cosTheta \cdot -2}} \cdot \left(cosTheta \cdot \left(1 + cosTheta \cdot cosTheta\right)\right)}}
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
clear-numN/A
/-lowering-/.f32N/A
div-invN/A
times-fracN/A
*-lowering-*.f32N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f3298.5%
Applied egg-rr98.5%
Taylor expanded in cosTheta around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
unpow2N/A
*-lowering-*.f3297.7%
Simplified97.7%
(FPCore (cosTheta c) :precision binary32 (/ 1.0 (+ (+ 1.0 c) (/ 1.0 (* cosTheta (sqrt (/ PI (+ 1.0 (* cosTheta -2.0)))))))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + (1.0f / (cosTheta * sqrtf((((float) M_PI) / (1.0f + (cosTheta * -2.0f)))))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(1.0) / Float32(cosTheta * sqrt(Float32(Float32(pi) / Float32(Float32(1.0) + Float32(cosTheta * Float32(-2.0))))))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / ((single(1.0) + c) + (single(1.0) / (cosTheta * sqrt((single(pi) / (single(1.0) + (cosTheta * single(-2.0)))))))); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \frac{1}{cosTheta \cdot \sqrt{\frac{\pi}{1 + cosTheta \cdot -2}}}}
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
clear-numN/A
/-lowering-/.f32N/A
div-invN/A
times-fracN/A
*-lowering-*.f32N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f3298.5%
Applied egg-rr98.5%
Taylor expanded in cosTheta around 0
Simplified96.2%
Final simplification96.2%
(FPCore (cosTheta c) :precision binary32 (/ 1.0 (+ (+ 1.0 c) (/ 1.0 (* cosTheta (* (+ 1.0 cosTheta) (sqrt PI)))))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + (1.0f / (cosTheta * ((1.0f + cosTheta) * sqrtf(((float) M_PI))))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(1.0) / Float32(cosTheta * Float32(Float32(Float32(1.0) + cosTheta) * sqrt(Float32(pi))))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / ((single(1.0) + c) + (single(1.0) / (cosTheta * ((single(1.0) + cosTheta) * sqrt(single(pi)))))); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \frac{1}{cosTheta \cdot \left(\left(1 + cosTheta\right) \cdot \sqrt{\pi}\right)}}
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
clear-numN/A
/-lowering-/.f32N/A
div-invN/A
times-fracN/A
*-lowering-*.f32N/A
sqrt-undivN/A
sqrt-lowering-sqrt.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f3298.5%
Applied egg-rr98.5%
Taylor expanded in cosTheta around 0
*-lowering-*.f32N/A
distribute-rgt1-inN/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
sqrt-lowering-sqrt.f32N/A
PI-lowering-PI.f3295.9%
Simplified95.9%
Final simplification95.9%
(FPCore (cosTheta c) :precision binary32 (/ 1.0 (+ (+ 1.0 c) (/ (/ (- 1.0 cosTheta) (pow PI 0.5)) cosTheta))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + (((1.0f - cosTheta) / powf(((float) M_PI), 0.5f)) / cosTheta));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(Float32(Float32(1.0) - cosTheta) / (Float32(pi) ^ Float32(0.5))) / cosTheta))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / ((single(1.0) + c) + (((single(1.0) - cosTheta) / (single(pi) ^ single(0.5))) / cosTheta)); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \frac{\frac{1 - cosTheta}{{\pi}^{0.5}}}{cosTheta}}
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
Taylor expanded in c around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
+-lowering-+.f32N/A
Simplified98.1%
Taylor expanded in cosTheta around 0
associate-*r*N/A
/-lowering-/.f32N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f32N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
sqrt-lowering-sqrt.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f3295.4%
Simplified95.4%
/-lowering-/.f32N/A
sqrt-divN/A
metadata-evalN/A
un-div-invN/A
/-lowering-/.f32N/A
--lowering--.f32N/A
pow1/2N/A
pow-lowering-pow.f32N/A
PI-lowering-PI.f3295.4%
Applied egg-rr95.4%
Final simplification95.4%
(FPCore (cosTheta c) :precision binary32 (/ 1.0 (+ 1.0 (+ c (* (- 1.0 cosTheta) (/ (pow PI -0.5) cosTheta))))))
float code(float cosTheta, float c) {
return 1.0f / (1.0f + (c + ((1.0f - cosTheta) * (powf(((float) M_PI), -0.5f) / cosTheta))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(1.0) + Float32(c + Float32(Float32(Float32(1.0) - cosTheta) * Float32((Float32(pi) ^ Float32(-0.5)) / cosTheta))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / (single(1.0) + (c + ((single(1.0) - cosTheta) * ((single(pi) ^ single(-0.5)) / cosTheta)))); end
\begin{array}{l}
\\
\frac{1}{1 + \left(c + \left(1 - cosTheta\right) \cdot \frac{{\pi}^{-0.5}}{cosTheta}\right)}
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
Taylor expanded in c around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
+-lowering-+.f32N/A
Simplified98.1%
Taylor expanded in cosTheta around 0
associate-*r*N/A
/-lowering-/.f32N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f32N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
sqrt-lowering-sqrt.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f3295.4%
Simplified95.4%
+-commutativeN/A
associate-+r+N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-/l*N/A
*-lowering-*.f32N/A
--lowering--.f32N/A
/-lowering-/.f32N/A
inv-powN/A
sqrt-pow1N/A
metadata-evalN/A
pow-lowering-pow.f32N/A
PI-lowering-PI.f3295.4%
Applied egg-rr95.4%
Final simplification95.4%
(FPCore (cosTheta c) :precision binary32 (/ 1.0 (+ (+ 1.0 c) (* (pow PI -0.5) (/ (- 1.0 cosTheta) cosTheta)))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + (powf(((float) M_PI), -0.5f) * ((1.0f - cosTheta) / cosTheta)));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32((Float32(pi) ^ Float32(-0.5)) * Float32(Float32(Float32(1.0) - cosTheta) / cosTheta)))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / ((single(1.0) + c) + ((single(pi) ^ single(-0.5)) * ((single(1.0) - cosTheta) / cosTheta))); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + {\pi}^{-0.5} \cdot \frac{1 - cosTheta}{cosTheta}}
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
Taylor expanded in c around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
+-lowering-+.f32N/A
Simplified98.1%
Taylor expanded in cosTheta around 0
associate-*r*N/A
/-lowering-/.f32N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f32N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
sqrt-lowering-sqrt.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f3295.4%
Simplified95.4%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f32N/A
inv-powN/A
sqrt-pow1N/A
metadata-evalN/A
pow-lowering-pow.f32N/A
PI-lowering-PI.f32N/A
/-lowering-/.f32N/A
--lowering--.f3295.4%
Applied egg-rr95.4%
Final simplification95.4%
(FPCore (cosTheta c) :precision binary32 (/ 1.0 (+ 1.0 (* (sqrt (/ 1.0 PI)) (+ (/ 1.0 cosTheta) -1.0)))))
float code(float cosTheta, float c) {
return 1.0f / (1.0f + (sqrtf((1.0f / ((float) M_PI))) * ((1.0f / cosTheta) + -1.0f)));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(1.0) + Float32(sqrt(Float32(Float32(1.0) / Float32(pi))) * Float32(Float32(Float32(1.0) / cosTheta) + Float32(-1.0))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / (single(1.0) + (sqrt((single(1.0) / single(pi))) * ((single(1.0) / cosTheta) + single(-1.0)))); end
\begin{array}{l}
\\
\frac{1}{1 + \sqrt{\frac{1}{\pi}} \cdot \left(\frac{1}{cosTheta} + -1\right)}
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
Taylor expanded in c around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
+-lowering-+.f32N/A
Simplified98.1%
Taylor expanded in cosTheta around 0
associate-*r*N/A
/-lowering-/.f32N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f32N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
sqrt-lowering-sqrt.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f3295.4%
Simplified95.4%
Taylor expanded in c around 0
/-lowering-/.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f32N/A
div-subN/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
+-lowering-+.f32N/A
/-lowering-/.f3295.3%
Simplified95.3%
(FPCore (cosTheta c) :precision binary32 (/ 1.0 (+ 1.0 (* (pow PI -0.5) (/ (- 1.0 cosTheta) cosTheta)))))
float code(float cosTheta, float c) {
return 1.0f / (1.0f + (powf(((float) M_PI), -0.5f) * ((1.0f - cosTheta) / cosTheta)));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(1.0) + Float32((Float32(pi) ^ Float32(-0.5)) * Float32(Float32(Float32(1.0) - cosTheta) / cosTheta)))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / (single(1.0) + ((single(pi) ^ single(-0.5)) * ((single(1.0) - cosTheta) / cosTheta))); end
\begin{array}{l}
\\
\frac{1}{1 + {\pi}^{-0.5} \cdot \frac{1 - cosTheta}{cosTheta}}
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
Taylor expanded in c around 0
+-commutativeN/A
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
+-lowering-+.f32N/A
Simplified98.1%
Taylor expanded in cosTheta around 0
associate-*r*N/A
/-lowering-/.f32N/A
distribute-rgt1-inN/A
+-commutativeN/A
*-lowering-*.f32N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
sqrt-lowering-sqrt.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f3295.4%
Simplified95.4%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f32N/A
inv-powN/A
sqrt-pow1N/A
metadata-evalN/A
pow-lowering-pow.f32N/A
PI-lowering-PI.f32N/A
/-lowering-/.f32N/A
--lowering--.f3295.4%
Applied egg-rr95.4%
Taylor expanded in c around 0
Simplified95.2%
Final simplification95.2%
(FPCore (cosTheta c) :precision binary32 (* cosTheta (- (sqrt PI) (* cosTheta (* c PI)))))
float code(float cosTheta, float c) {
return cosTheta * (sqrtf(((float) M_PI)) - (cosTheta * (c * ((float) M_PI))));
}
function code(cosTheta, c) return Float32(cosTheta * Float32(sqrt(Float32(pi)) - Float32(cosTheta * Float32(c * Float32(pi))))) end
function tmp = code(cosTheta, c) tmp = cosTheta * (sqrt(single(pi)) - (cosTheta * (c * single(pi)))); end
\begin{array}{l}
\\
cosTheta \cdot \left(\sqrt{\pi} - cosTheta \cdot \left(c \cdot \pi\right)\right)
\end{array}
Initial program 98.0%
Taylor expanded in cosTheta around 0
*-lowering-*.f32N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
sqrt-lowering-sqrt.f32N/A
PI-lowering-PI.f32N/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f32N/A
Simplified96.1%
Taylor expanded in c around inf
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3293.2%
Simplified93.2%
Final simplification93.2%
(FPCore (cosTheta c) :precision binary32 (* cosTheta (sqrt PI)))
float code(float cosTheta, float c) {
return cosTheta * sqrtf(((float) M_PI));
}
function code(cosTheta, c) return Float32(cosTheta * sqrt(Float32(pi))) end
function tmp = code(cosTheta, c) tmp = cosTheta * sqrt(single(pi)); end
\begin{array}{l}
\\
cosTheta \cdot \sqrt{\pi}
\end{array}
Initial program 98.0%
Taylor expanded in cosTheta around 0
*-lowering-*.f32N/A
sqrt-lowering-sqrt.f32N/A
PI-lowering-PI.f3293.2%
Simplified93.2%
(FPCore (cosTheta c) :precision binary32 (- 0.0 (* c (* cosTheta (* cosTheta PI)))))
float code(float cosTheta, float c) {
return 0.0f - (c * (cosTheta * (cosTheta * ((float) M_PI))));
}
function code(cosTheta, c) return Float32(Float32(0.0) - Float32(c * Float32(cosTheta * Float32(cosTheta * Float32(pi))))) end
function tmp = code(cosTheta, c) tmp = single(0.0) - (c * (cosTheta * (cosTheta * single(pi)))); end
\begin{array}{l}
\\
0 - c \cdot \left(cosTheta \cdot \left(cosTheta \cdot \pi\right)\right)
\end{array}
Initial program 98.0%
Taylor expanded in cosTheta around 0
*-lowering-*.f32N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
sqrt-lowering-sqrt.f32N/A
PI-lowering-PI.f32N/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f32N/A
Simplified96.1%
Taylor expanded in c around inf
mul-1-negN/A
neg-lowering-neg.f32N/A
*-lowering-*.f32N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3211.0%
Simplified11.0%
Final simplification11.0%
(FPCore (cosTheta c) :precision binary32 (- 1.0 c))
float code(float cosTheta, float c) {
return 1.0f - c;
}
real(4) function code(costheta, c)
real(4), intent (in) :: costheta
real(4), intent (in) :: c
code = 1.0e0 - c
end function
function code(cosTheta, c) return Float32(Float32(1.0) - c) end
function tmp = code(cosTheta, c) tmp = single(1.0) - c; end
\begin{array}{l}
\\
1 - c
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
Taylor expanded in cosTheta around 0
/-lowering-/.f3292.7%
Simplified92.7%
Taylor expanded in cosTheta around inf
/-lowering-/.f32N/A
+-lowering-+.f3210.8%
Simplified10.8%
Taylor expanded in c around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f3210.8%
Simplified10.8%
(FPCore (cosTheta c) :precision binary32 1.0)
float code(float cosTheta, float c) {
return 1.0f;
}
real(4) function code(costheta, c)
real(4), intent (in) :: costheta
real(4), intent (in) :: c
code = 1.0e0
end function
function code(cosTheta, c) return Float32(1.0) end
function tmp = code(cosTheta, c) tmp = single(1.0); end
\begin{array}{l}
\\
1
\end{array}
Initial program 98.0%
/-lowering-/.f32N/A
+-lowering-+.f32N/A
+-lowering-+.f32N/A
associate-*l*N/A
distribute-lft-neg-outN/A
exp-negN/A
associate-*r/N/A
*-rgt-identityN/A
times-fracN/A
*-lft-identityN/A
/-lowering-/.f32N/A
Simplified98.4%
Taylor expanded in cosTheta around 0
/-lowering-/.f3292.7%
Simplified92.7%
Taylor expanded in cosTheta around inf
/-lowering-/.f32N/A
+-lowering-+.f3210.8%
Simplified10.8%
Taylor expanded in c around 0
Simplified10.8%
herbie shell --seed 2024138
(FPCore (cosTheta c)
:name "Beckmann Sample, normalization factor"
:precision binary32
:pre (and (and (< 0.0 cosTheta) (< cosTheta 0.9999)) (and (< -1.0 c) (< c 1.0)))
(/ 1.0 (+ (+ 1.0 c) (* (* (/ 1.0 (sqrt PI)) (/ (sqrt (- (- 1.0 cosTheta) cosTheta)) cosTheta)) (exp (* (- cosTheta) cosTheta))))))