
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(+ 1.0 c)
(*
(* (/ 1.0 (sqrt (PI))) (/ (sqrt (- (- 1.0 cosTheta) cosTheta)) cosTheta))
(exp (* (- cosTheta) cosTheta))))))\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \left(\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \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 12 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))))))\begin{array}{l}
\\
\frac{1}{\left(1 + c\right) + \left(\frac{1}{\sqrt{\mathsf{PI}\left(\right)}} \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
(+
(*
(exp (* (- cosTheta) cosTheta))
(/ (/ (sqrt (* (- 1.0 (* 2.0 cosTheta)) (PI))) cosTheta) (PI)))
(+ c 1.0))))\begin{array}{l}
\\
\frac{1}{e^{\left(-cosTheta\right) \cdot cosTheta} \cdot \frac{\frac{\sqrt{\left(1 - 2 \cdot cosTheta\right) \cdot \mathsf{PI}\left(\right)}}{cosTheta}}{\mathsf{PI}\left(\right)} + \left(c + 1\right)}
\end{array}
Initial program 97.9%
lift-*.f32N/A
lift-/.f32N/A
associate-*r/N/A
lower-/.f32N/A
lift-/.f32N/A
associate-*l/N/A
*-lft-identityN/A
lift-sqrt.f32N/A
lift-sqrt.f32N/A
sqrt-undivN/A
lower-sqrt.f32N/A
lower-/.f3298.0
Applied rewrites98.0%
lift-sqrt.f32N/A
lift-/.f32N/A
lift--.f32N/A
div-subN/A
frac-subN/A
lift-PI.f32N/A
lift-PI.f32N/A
sqrt-divN/A
lift-PI.f32N/A
lift-PI.f32N/A
sqrt-unprodN/A
rem-square-sqrtN/A
lower-/.f32N/A
Applied rewrites98.2%
lift-/.f32N/A
lift-/.f32N/A
associate-/l/N/A
associate-/r*N/A
lower-/.f32N/A
Applied rewrites98.9%
Final simplification98.9%
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(*
(/ (/ (sqrt (* (- (- 1.0 cosTheta) cosTheta) (PI))) cosTheta) (PI))
(exp (* (- cosTheta) cosTheta)))
(+ c 1.0))))\begin{array}{l}
\\
\frac{1}{\frac{\frac{\sqrt{\left(\left(1 - cosTheta\right) - cosTheta\right) \cdot \mathsf{PI}\left(\right)}}{cosTheta}}{\mathsf{PI}\left(\right)} \cdot e^{\left(-cosTheta\right) \cdot cosTheta} + \left(c + 1\right)}
\end{array}
Initial program 97.9%
lift-*.f32N/A
lift-/.f32N/A
associate-*r/N/A
lower-/.f32N/A
lift-/.f32N/A
associate-*l/N/A
*-lft-identityN/A
lift-sqrt.f32N/A
lift-sqrt.f32N/A
sqrt-undivN/A
lower-sqrt.f32N/A
lower-/.f3298.0
Applied rewrites98.0%
lift-sqrt.f32N/A
lift-/.f32N/A
lift--.f32N/A
div-subN/A
frac-subN/A
lift-PI.f32N/A
lift-PI.f32N/A
sqrt-divN/A
lift-PI.f32N/A
lift-PI.f32N/A
sqrt-unprodN/A
rem-square-sqrtN/A
lower-/.f32N/A
Applied rewrites98.2%
lift-/.f32N/A
lift-/.f32N/A
associate-/l/N/A
associate-/r*N/A
lower-/.f32N/A
Applied rewrites98.9%
lift--.f32N/A
lift-*.f32N/A
*-commutativeN/A
count-2N/A
associate--l-N/A
lift--.f32N/A
lift--.f3298.8
Applied rewrites98.8%
Final simplification98.8%
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(*
(/ (sqrt (* (- 1.0 (* 2.0 cosTheta)) (PI))) (* cosTheta (PI)))
(exp (* (- cosTheta) cosTheta)))
(+ c 1.0))))\begin{array}{l}
\\
\frac{1}{\frac{\sqrt{\left(1 - 2 \cdot cosTheta\right) \cdot \mathsf{PI}\left(\right)}}{cosTheta \cdot \mathsf{PI}\left(\right)} \cdot e^{\left(-cosTheta\right) \cdot cosTheta} + \left(c + 1\right)}
\end{array}
Initial program 97.9%
lift-*.f32N/A
lift-/.f32N/A
associate-*r/N/A
lower-/.f32N/A
lift-/.f32N/A
associate-*l/N/A
*-lft-identityN/A
lift-sqrt.f32N/A
lift-sqrt.f32N/A
sqrt-undivN/A
lower-sqrt.f32N/A
lower-/.f3298.0
Applied rewrites98.0%
lift-sqrt.f32N/A
lift-/.f32N/A
lift--.f32N/A
div-subN/A
frac-subN/A
lift-PI.f32N/A
lift-PI.f32N/A
sqrt-divN/A
lift-PI.f32N/A
lift-PI.f32N/A
sqrt-unprodN/A
rem-square-sqrtN/A
lower-/.f32N/A
Applied rewrites98.2%
lift-/.f32N/A
lift-/.f32N/A
associate-/r*N/A
lift-*.f32N/A
lower-/.f3298.8
Applied rewrites98.8%
Final simplification98.8%
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(*
(/ (sqrt (- (- 1.0 cosTheta) cosTheta)) (* (sqrt (PI)) cosTheta))
(exp (* (- cosTheta) cosTheta)))
1.0)))\begin{array}{l}
\\
\frac{1}{\frac{\sqrt{\left(1 - cosTheta\right) - cosTheta}}{\sqrt{\mathsf{PI}\left(\right)} \cdot cosTheta} \cdot e^{\left(-cosTheta\right) \cdot cosTheta} + 1}
\end{array}
Initial program 97.9%
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
frac-timesN/A
*-lft-identityN/A
lower-/.f32N/A
*-commutativeN/A
lower-*.f3298.7
Applied rewrites98.7%
Taylor expanded in c around 0
Applied rewrites98.3%
Final simplification98.3%
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(*
(/ (sqrt (/ (- (- 1.0 cosTheta) cosTheta) (PI))) cosTheta)
(exp (* (- cosTheta) cosTheta)))
(+ c 1.0))))\begin{array}{l}
\\
\frac{1}{\frac{\sqrt{\frac{\left(1 - cosTheta\right) - cosTheta}{\mathsf{PI}\left(\right)}}}{cosTheta} \cdot e^{\left(-cosTheta\right) \cdot cosTheta} + \left(c + 1\right)}
\end{array}
Initial program 97.9%
lift-*.f32N/A
lift-/.f32N/A
associate-*r/N/A
lower-/.f32N/A
lift-/.f32N/A
associate-*l/N/A
*-lft-identityN/A
lift-sqrt.f32N/A
lift-sqrt.f32N/A
sqrt-undivN/A
lower-sqrt.f32N/A
lower-/.f3298.0
Applied rewrites98.0%
Final simplification98.0%
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(*
(/ (sqrt (/ (- (- 1.0 cosTheta) cosTheta) (PI))) cosTheta)
(exp (* (- cosTheta) cosTheta)))
1.0)))\begin{array}{l}
\\
\frac{1}{\frac{\sqrt{\frac{\left(1 - cosTheta\right) - cosTheta}{\mathsf{PI}\left(\right)}}}{cosTheta} \cdot e^{\left(-cosTheta\right) \cdot cosTheta} + 1}
\end{array}
Initial program 97.9%
lift-*.f32N/A
lift-/.f32N/A
associate-*r/N/A
lower-/.f32N/A
lift-/.f32N/A
associate-*l/N/A
*-lft-identityN/A
lift-sqrt.f32N/A
lift-sqrt.f32N/A
sqrt-undivN/A
lower-sqrt.f32N/A
lower-/.f3298.0
Applied rewrites98.0%
Taylor expanded in c around 0
Applied rewrites97.7%
Final simplification97.7%
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(*
(/ (- 1.0 cosTheta) (* (sqrt (PI)) cosTheta))
(exp (* (- cosTheta) cosTheta)))
(+ c 1.0))))\begin{array}{l}
\\
\frac{1}{\frac{1 - cosTheta}{\sqrt{\mathsf{PI}\left(\right)} \cdot cosTheta} \cdot e^{\left(-cosTheta\right) \cdot cosTheta} + \left(c + 1\right)}
\end{array}
Initial program 97.9%
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
frac-timesN/A
*-lft-identityN/A
lower-/.f32N/A
*-commutativeN/A
lower-*.f3298.7
Applied rewrites98.7%
Taylor expanded in cosTheta around 0
mul-1-negN/A
unsub-negN/A
lower--.f3295.9
Applied rewrites95.9%
Final simplification95.9%
(FPCore (cosTheta c) :precision binary32 (* (/ 1.0 (* (+ (/ (/ (- 1.0 cosTheta) (sqrt (PI))) (+ c 1.0)) cosTheta) (+ c 1.0))) cosTheta))
\begin{array}{l}
\\
\frac{1}{\left(\frac{\frac{1 - cosTheta}{\sqrt{\mathsf{PI}\left(\right)}}}{c + 1} + cosTheta\right) \cdot \left(c + 1\right)} \cdot cosTheta
\end{array}
Initial program 97.9%
lift-+.f32N/A
lift-+.f32N/A
flip3-+N/A
clear-numN/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
lift-/.f32N/A
associate-*r/N/A
associate-*r/N/A
Applied rewrites98.1%
Taylor expanded in cosTheta around 0
associate-*r/N/A
*-rgt-identityN/A
lower-/.f32N/A
lower-sqrt.f32N/A
lower-/.f32N/A
lower-PI.f32N/A
+-commutativeN/A
lower-+.f3291.3
Applied rewrites91.3%
Taylor expanded in cosTheta around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f32N/A
associate-*r/N/A
*-rgt-identityN/A
lower-/.f32N/A
lower-sqrt.f32N/A
lower-/.f32N/A
lower-PI.f32N/A
+-commutativeN/A
lower-+.f32N/A
lower-*.f32N/A
lower-/.f32N/A
+-commutativeN/A
lower-+.f32N/A
lower-sqrt.f32N/A
Applied rewrites94.8%
lift-/.f32N/A
lift-/.f32N/A
lift-*.f32N/A
associate-/r*N/A
Applied rewrites95.6%
(FPCore (cosTheta c) :precision binary32 (/ cosTheta (* (+ (/ (/ (- 1.0 cosTheta) (sqrt (PI))) (+ c 1.0)) cosTheta) (+ c 1.0))))
\begin{array}{l}
\\
\frac{cosTheta}{\left(\frac{\frac{1 - cosTheta}{\sqrt{\mathsf{PI}\left(\right)}}}{c + 1} + cosTheta\right) \cdot \left(c + 1\right)}
\end{array}
Initial program 97.9%
lift-+.f32N/A
lift-+.f32N/A
flip3-+N/A
clear-numN/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
lift-/.f32N/A
associate-*r/N/A
associate-*r/N/A
Applied rewrites98.1%
Taylor expanded in cosTheta around 0
associate-*r/N/A
*-rgt-identityN/A
lower-/.f32N/A
lower-sqrt.f32N/A
lower-/.f32N/A
lower-PI.f32N/A
+-commutativeN/A
lower-+.f3291.3
Applied rewrites91.3%
Taylor expanded in cosTheta around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f32N/A
associate-*r/N/A
*-rgt-identityN/A
lower-/.f32N/A
lower-sqrt.f32N/A
lower-/.f32N/A
lower-PI.f32N/A
+-commutativeN/A
lower-+.f32N/A
lower-*.f32N/A
lower-/.f32N/A
+-commutativeN/A
lower-+.f32N/A
lower-sqrt.f32N/A
Applied rewrites94.8%
Applied rewrites94.9%
(FPCore (cosTheta c) :precision binary32 (/ 1.0 (+ (/ (* (sqrt (/ 1.0 (PI))) (- 1.0 cosTheta)) cosTheta) (+ c 1.0))))
\begin{array}{l}
\\
\frac{1}{\frac{\sqrt{\frac{1}{\mathsf{PI}\left(\right)}} \cdot \left(1 - cosTheta\right)}{cosTheta} + \left(c + 1\right)}
\end{array}
Initial program 97.9%
lift-*.f32N/A
lift-/.f32N/A
associate-*r/N/A
lower-/.f32N/A
lift-/.f32N/A
associate-*l/N/A
*-lft-identityN/A
lift-sqrt.f32N/A
lift-sqrt.f32N/A
sqrt-undivN/A
lower-sqrt.f32N/A
lower-/.f3298.0
Applied rewrites98.0%
Taylor expanded in cosTheta around 0
associate-*r*N/A
lower-/.f32N/A
distribute-rgt1-inN/A
+-commutativeN/A
lower-*.f32N/A
mul-1-negN/A
unsub-negN/A
lower--.f32N/A
lower-sqrt.f32N/A
lower-/.f32N/A
lower-PI.f3294.7
Applied rewrites94.7%
Final simplification94.7%
(FPCore (cosTheta c) :precision binary32 (* (sqrt (PI)) cosTheta))
\begin{array}{l}
\\
\sqrt{\mathsf{PI}\left(\right)} \cdot cosTheta
\end{array}
Initial program 97.9%
Taylor expanded in cosTheta around 0
*-commutativeN/A
lower-*.f32N/A
lower-sqrt.f32N/A
lower-PI.f3292.4
Applied rewrites92.4%
(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}
\\
\frac{1}{c}
\end{array}
Initial program 97.9%
Taylor expanded in c around inf
lower-/.f324.9
Applied rewrites4.9%
herbie shell --seed 2024255
(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))))))