
(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 8 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)
(/
(/ (sqrt (- 1.0 (+ cosTheta cosTheta))) (* cosTheta (sqrt PI)))
(exp (* cosTheta cosTheta))))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + ((sqrtf((1.0f - (cosTheta + cosTheta))) / (cosTheta * sqrtf(((float) M_PI)))) / expf((cosTheta * cosTheta))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(sqrt(Float32(Float32(1.0) - Float32(cosTheta + cosTheta))) / Float32(cosTheta * sqrt(Float32(pi)))) / exp(Float32(cosTheta * cosTheta))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / ((single(1.0) + c) + ((sqrt((single(1.0) - (cosTheta + cosTheta))) / (cosTheta * sqrt(single(pi)))) / exp((cosTheta * cosTheta)))); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \frac{\frac{\sqrt{1 - \left(cosTheta + cosTheta\right)}}{cosTheta \cdot \sqrt{\pi}}}{e^{cosTheta \cdot cosTheta}}}
\end{array}
Initial program 97.3%
associate-*l*97.3%
*-lft-identity97.3%
associate-*r*97.3%
*-lft-identity97.3%
associate--l-97.3%
*-commutative97.3%
Simplified97.3%
associate--r+97.3%
*-commutative97.3%
associate-*l*97.3%
frac-times98.6%
*-un-lft-identity98.6%
associate--r+98.5%
associate-/l/97.8%
*-commutative97.8%
pow-exp97.8%
pow-neg97.8%
div-inv97.8%
associate-/l/98.5%
pow-exp98.5%
Applied egg-rr98.5%
Final simplification98.5%
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(+ 1.0 c)
(/
(/ (sqrt (/ (+ 1.0 (* cosTheta -2.0)) PI)) cosTheta)
(exp (* cosTheta cosTheta))))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + ((sqrtf(((1.0f + (cosTheta * -2.0f)) / ((float) M_PI))) / cosTheta) / expf((cosTheta * cosTheta))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(sqrt(Float32(Float32(Float32(1.0) + Float32(cosTheta * Float32(-2.0))) / Float32(pi))) / cosTheta) / exp(Float32(cosTheta * cosTheta))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / ((single(1.0) + c) + ((sqrt(((single(1.0) + (cosTheta * single(-2.0))) / single(pi))) / cosTheta) / exp((cosTheta * cosTheta)))); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \frac{\frac{\sqrt{\frac{1 + cosTheta \cdot -2}{\pi}}}{cosTheta}}{e^{cosTheta \cdot cosTheta}}}
\end{array}
Initial program 97.3%
associate-*l*97.3%
*-lft-identity97.3%
associate-*r*97.3%
*-lft-identity97.3%
associate--l-97.3%
*-commutative97.3%
Simplified97.3%
associate--r+97.3%
*-commutative97.3%
associate-*l*97.3%
frac-times98.6%
*-un-lft-identity98.6%
associate--r+98.5%
associate-/l/97.8%
*-commutative97.8%
pow-exp97.8%
pow-neg97.8%
div-inv97.8%
associate-/l/98.5%
pow-exp98.5%
Applied egg-rr98.5%
*-un-lft-identity98.5%
associate-/r*98.1%
count-298.1%
sqrt-div98.0%
cancel-sign-sub-inv98.0%
metadata-eval98.0%
Applied egg-rr98.0%
*-lft-identity98.0%
metadata-eval98.0%
cancel-sign-sub-inv98.0%
cancel-sign-sub-inv98.0%
metadata-eval98.0%
*-commutative98.0%
Simplified98.0%
Final simplification98.0%
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
1.0
(/
(sqrt (/ (+ 1.0 (* cosTheta -2.0)) PI))
(* cosTheta (exp (* cosTheta cosTheta)))))))
float code(float cosTheta, float c) {
return 1.0f / (1.0f + (sqrtf(((1.0f + (cosTheta * -2.0f)) / ((float) M_PI))) / (cosTheta * expf((cosTheta * cosTheta)))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(1.0) + Float32(sqrt(Float32(Float32(Float32(1.0) + Float32(cosTheta * Float32(-2.0))) / Float32(pi))) / Float32(cosTheta * exp(Float32(cosTheta * cosTheta)))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / (single(1.0) + (sqrt(((single(1.0) + (cosTheta * single(-2.0))) / single(pi))) / (cosTheta * exp((cosTheta * cosTheta))))); end
\begin{array}{l}
\\
\frac{1}{1 + \frac{\sqrt{\frac{1 + cosTheta \cdot -2}{\pi}}}{cosTheta \cdot e^{cosTheta \cdot cosTheta}}}
\end{array}
Initial program 97.3%
associate-*l*97.3%
*-lft-identity97.3%
associate-*r*97.3%
*-lft-identity97.3%
associate--l-97.3%
*-commutative97.3%
Simplified97.3%
associate--r+97.3%
*-commutative97.3%
associate-*l*97.3%
frac-times98.6%
*-un-lft-identity98.6%
associate--r+98.5%
associate-/l/97.8%
*-commutative97.8%
pow-exp97.8%
pow-neg97.8%
div-inv97.8%
associate-/l/98.5%
pow-exp98.5%
Applied egg-rr98.5%
Taylor expanded in c around 0 97.2%
associate-*r/97.8%
cancel-sign-sub-inv97.8%
metadata-eval97.8%
*-rgt-identity97.8%
metadata-eval97.8%
cancel-sign-sub-inv97.8%
cancel-sign-sub-inv97.8%
metadata-eval97.8%
*-commutative97.8%
*-commutative97.8%
unpow297.8%
Simplified97.8%
Final simplification97.8%
(FPCore (cosTheta c) :precision binary32 (/ 1.0 (+ (+ 1.0 c) (* (/ 1.0 (sqrt PI)) (+ (+ (* cosTheta -1.5) (/ 1.0 cosTheta)) -1.0)))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + ((1.0f / sqrtf(((float) M_PI))) * (((cosTheta * -1.5f) + (1.0f / cosTheta)) + -1.0f)));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(Float32(1.0) / sqrt(Float32(pi))) * Float32(Float32(Float32(cosTheta * Float32(-1.5)) + Float32(Float32(1.0) / cosTheta)) + Float32(-1.0))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / ((single(1.0) + c) + ((single(1.0) / sqrt(single(pi))) * (((cosTheta * single(-1.5)) + (single(1.0) / cosTheta)) + single(-1.0)))); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \frac{1}{\sqrt{\pi}} \cdot \left(\left(cosTheta \cdot -1.5 + \frac{1}{cosTheta}\right) + -1\right)}
\end{array}
Initial program 97.3%
associate-*l*97.3%
*-lft-identity97.3%
associate-*r*97.3%
*-lft-identity97.3%
associate--l-97.3%
*-commutative97.3%
Simplified97.3%
Taylor expanded in cosTheta around 0 95.2%
Final simplification95.2%
(FPCore (cosTheta c) :precision binary32 (/ 1.0 (+ (+ 1.0 c) (* (/ 1.0 (sqrt PI)) (+ (/ 1.0 cosTheta) -1.0)))))
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + ((1.0f / sqrtf(((float) M_PI))) * ((1.0f / cosTheta) + -1.0f)));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(Float32(1.0) / sqrt(Float32(pi))) * Float32(Float32(Float32(1.0) / cosTheta) + Float32(-1.0))))) 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(-1.0)))); end
\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \frac{1}{\sqrt{\pi}} \cdot \left(\frac{1}{cosTheta} + -1\right)}
\end{array}
Initial program 97.3%
associate-*l*97.3%
*-lft-identity97.3%
associate-*r*97.3%
*-lft-identity97.3%
associate--l-97.3%
*-commutative97.3%
Simplified97.3%
Taylor expanded in cosTheta around 0 93.5%
Final simplification93.5%
(FPCore (cosTheta c) :precision binary32 (/ 1.0 (+ 1.0 (* (+ (/ 1.0 cosTheta) -1.0) (sqrt (/ 1.0 PI))))))
float code(float cosTheta, float c) {
return 1.0f / (1.0f + (((1.0f / cosTheta) + -1.0f) * sqrtf((1.0f / ((float) M_PI)))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(1.0) + Float32(Float32(Float32(Float32(1.0) / cosTheta) + Float32(-1.0)) * sqrt(Float32(Float32(1.0) / Float32(pi)))))) end
function tmp = code(cosTheta, c) tmp = single(1.0) / (single(1.0) + (((single(1.0) / cosTheta) + single(-1.0)) * sqrt((single(1.0) / single(pi))))); end
\begin{array}{l}
\\
\frac{1}{1 + \left(\frac{1}{cosTheta} + -1\right) \cdot \sqrt{\frac{1}{\pi}}}
\end{array}
Initial program 97.3%
associate-*l*97.3%
*-lft-identity97.3%
associate-*r*97.3%
*-lft-identity97.3%
associate--l-97.3%
*-commutative97.3%
Simplified97.3%
associate--r+97.3%
*-commutative97.3%
associate-*l*97.3%
frac-times98.6%
*-un-lft-identity98.6%
associate--r+98.5%
associate-/l/97.8%
*-commutative97.8%
pow-exp97.8%
pow-neg97.8%
div-inv97.8%
associate-/l/98.5%
pow-exp98.5%
Applied egg-rr98.5%
Taylor expanded in c around 0 97.2%
associate-*r/97.8%
cancel-sign-sub-inv97.8%
metadata-eval97.8%
*-rgt-identity97.8%
metadata-eval97.8%
cancel-sign-sub-inv97.8%
cancel-sign-sub-inv97.8%
metadata-eval97.8%
*-commutative97.8%
*-commutative97.8%
unpow297.8%
Simplified97.8%
*-commutative97.8%
add-log-exp97.9%
+-commutative97.9%
*-commutative97.9%
fma-def97.9%
Applied egg-rr97.9%
Taylor expanded in cosTheta around 0 93.5%
distribute-rgt-out93.4%
Simplified93.4%
Final simplification93.4%
(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 97.3%
associate-*l*97.3%
*-lft-identity97.3%
associate-*r*97.3%
*-lft-identity97.3%
associate--l-97.3%
*-commutative97.3%
Simplified97.3%
associate--r+97.3%
*-commutative97.3%
associate-*l*97.3%
frac-times98.6%
*-un-lft-identity98.6%
associate--r+98.5%
associate-/l/97.8%
*-commutative97.8%
pow-exp97.8%
pow-neg97.8%
div-inv97.8%
associate-/l/98.5%
pow-exp98.5%
Applied egg-rr98.5%
Taylor expanded in cosTheta around 0 91.3%
Final simplification91.3%
(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 97.3%
Taylor expanded in c around 0 97.3%
Taylor expanded in cosTheta around inf 11.4%
Final simplification11.4%
herbie shell --seed 2023188
(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))))))