
(FPCore (cosTheta alpha)
:precision binary32
(let* ((t_0 (- (* alpha alpha) 1.0)))
(/
t_0
(* (* (PI) (log (* alpha alpha))) (+ 1.0 (* (* t_0 cosTheta) cosTheta))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \alpha \cdot \alpha - 1\\
\frac{t\_0}{\left(\mathsf{PI}\left(\right) \cdot \log \left(\alpha \cdot \alpha\right)\right) \cdot \left(1 + \left(t\_0 \cdot cosTheta\right) \cdot cosTheta\right)}
\end{array}
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 21 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (cosTheta alpha)
:precision binary32
(let* ((t_0 (- (* alpha alpha) 1.0)))
(/
t_0
(* (* (PI) (log (* alpha alpha))) (+ 1.0 (* (* t_0 cosTheta) cosTheta))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \alpha \cdot \alpha - 1\\
\frac{t\_0}{\left(\mathsf{PI}\left(\right) \cdot \log \left(\alpha \cdot \alpha\right)\right) \cdot \left(1 + \left(t\_0 \cdot cosTheta\right) \cdot cosTheta\right)}
\end{array}
\end{array}
(FPCore (cosTheta alpha)
:precision binary32
(let* ((t_0 (- (* alpha alpha) 1.0)))
(/
t_0
(*
(- (* (* cosTheta t_0) cosTheta) -1.0)
(log (pow (* alpha alpha) (PI)))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \alpha \cdot \alpha - 1\\
\frac{t\_0}{\left(\left(cosTheta \cdot t\_0\right) \cdot cosTheta - -1\right) \cdot \log \left({\left(\alpha \cdot \alpha\right)}^{\mathsf{PI}\left(\right)}\right)}
\end{array}
\end{array}
Initial program 98.4%
lift-*.f32N/A
lift-log.f32N/A
log-pow-revN/A
lower-log.f32N/A
lower-pow.f3298.8
Applied rewrites98.8%
Final simplification98.8%
(FPCore (cosTheta alpha)
:precision binary32
(let* ((t_0 (- (* alpha alpha) 1.0)))
(/
t_0
(*
(log (pow alpha (+ (PI) (PI))))
(- (* (* cosTheta t_0) cosTheta) -1.0)))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \alpha \cdot \alpha - 1\\
\frac{t\_0}{\log \left({\alpha}^{\left(\mathsf{PI}\left(\right) + \mathsf{PI}\left(\right)\right)}\right) \cdot \left(\left(cosTheta \cdot t\_0\right) \cdot cosTheta - -1\right)}
\end{array}
\end{array}
Initial program 98.4%
lift-*.f32N/A
lift-log.f32N/A
log-pow-revN/A
lower-log.f32N/A
lower-pow.f3298.8
Applied rewrites98.8%
lift-pow.f32N/A
lift-*.f32N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f32N/A
lower-+.f3298.8
Applied rewrites98.8%
Final simplification98.8%
(FPCore (cosTheta alpha)
:precision binary32
(let* ((t_0 (- (* alpha alpha) 1.0)))
(/
t_0
(*
(* (log (* alpha alpha)) (PI))
(- (* (* cosTheta t_0) cosTheta) -1.0)))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \alpha \cdot \alpha - 1\\
\frac{t\_0}{\left(\log \left(\alpha \cdot \alpha\right) \cdot \mathsf{PI}\left(\right)\right) \cdot \left(\left(cosTheta \cdot t\_0\right) \cdot cosTheta - -1\right)}
\end{array}
\end{array}
Initial program 98.4%
Final simplification98.4%
(FPCore (cosTheta alpha)
:precision binary32
(let* ((t_0 (- (* alpha alpha) 1.0)))
(/
t_0
(* (* (log alpha) (+ (PI) (PI))) (- (* (* cosTheta t_0) cosTheta) -1.0)))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \alpha \cdot \alpha - 1\\
\frac{t\_0}{\left(\log \alpha \cdot \left(\mathsf{PI}\left(\right) + \mathsf{PI}\left(\right)\right)\right) \cdot \left(\left(cosTheta \cdot t\_0\right) \cdot cosTheta - -1\right)}
\end{array}
\end{array}
Initial program 98.4%
lift-*.f32N/A
*-commutativeN/A
lift-PI.f32N/A
add-cube-cbrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
pow2N/A
lower-pow.f32N/A
lift-PI.f32N/A
lower-cbrt.f32N/A
lift-PI.f32N/A
lower-cbrt.f3298.0
Applied rewrites98.0%
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
lift-pow.f32N/A
pow-plusN/A
lift-cbrt.f32N/A
metadata-evalN/A
rem-cube-cbrtN/A
*-commutativeN/A
lift-log.f32N/A
log-pow-revN/A
lift-*.f32N/A
unpow-prod-downN/A
pow-prod-upN/A
log-powN/A
lift-log.f32N/A
lower-*.f32N/A
lower-+.f3298.4
Applied rewrites98.4%
Final simplification98.4%
(FPCore (cosTheta alpha) :precision binary32 (/ (- (* alpha alpha) 1.0) (* (- (* (- cosTheta) cosTheta) -1.0) (* (log (* alpha alpha)) (PI)))))
\begin{array}{l}
\\
\frac{\alpha \cdot \alpha - 1}{\left(\left(-cosTheta\right) \cdot cosTheta - -1\right) \cdot \left(\log \left(\alpha \cdot \alpha\right) \cdot \mathsf{PI}\left(\right)\right)}
\end{array}
Initial program 98.4%
Taylor expanded in alpha around 0
mul-1-negN/A
lower-neg.f3297.3
Applied rewrites97.3%
Final simplification97.3%
(FPCore (cosTheta alpha) :precision binary32 (/ (- (* alpha alpha) 1.0) (* (- (* (- cosTheta) cosTheta) -1.0) (* (log alpha) (+ (PI) (PI))))))
\begin{array}{l}
\\
\frac{\alpha \cdot \alpha - 1}{\left(\left(-cosTheta\right) \cdot cosTheta - -1\right) \cdot \left(\log \alpha \cdot \left(\mathsf{PI}\left(\right) + \mathsf{PI}\left(\right)\right)\right)}
\end{array}
Initial program 98.4%
lift-*.f32N/A
*-commutativeN/A
lift-PI.f32N/A
add-cube-cbrtN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f32N/A
pow2N/A
lower-pow.f32N/A
lift-PI.f32N/A
lower-cbrt.f32N/A
lift-PI.f32N/A
lower-cbrt.f3298.0
Applied rewrites98.0%
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
lift-pow.f32N/A
pow-plusN/A
lift-cbrt.f32N/A
metadata-evalN/A
rem-cube-cbrtN/A
*-commutativeN/A
lift-log.f32N/A
log-pow-revN/A
lift-*.f32N/A
unpow-prod-downN/A
pow-prod-upN/A
log-powN/A
lift-log.f32N/A
lower-*.f32N/A
lower-+.f3298.4
Applied rewrites98.4%
Taylor expanded in alpha around 0
mul-1-negN/A
lower-neg.f3297.2
Applied rewrites97.2%
Final simplification97.2%
(FPCore (cosTheta alpha) :precision binary32 (* (/ (- alpha 1.0) (PI)) (/ (- alpha -1.0) (* 2.0 (log alpha)))))
\begin{array}{l}
\\
\frac{\alpha - 1}{\mathsf{PI}\left(\right)} \cdot \frac{\alpha - -1}{2 \cdot \log \alpha}
\end{array}
Initial program 98.4%
lift-*.f32N/A
lift-log.f32N/A
log-pow-revN/A
lower-log.f32N/A
lower-pow.f3298.8
Applied rewrites98.8%
lift-pow.f32N/A
lift-*.f32N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f32N/A
lower-+.f3298.8
Applied rewrites98.8%
Applied rewrites3.3%
Applied rewrites95.3%
Final simplification95.3%
(FPCore (cosTheta alpha) :precision binary32 (* (/ (- alpha 1.0) (* 2.0 (log alpha))) (/ (- alpha -1.0) (PI))))
\begin{array}{l}
\\
\frac{\alpha - 1}{2 \cdot \log \alpha} \cdot \frac{\alpha - -1}{\mathsf{PI}\left(\right)}
\end{array}
Initial program 98.4%
lift-*.f32N/A
lift-log.f32N/A
log-pow-revN/A
lower-log.f32N/A
lower-pow.f3298.8
Applied rewrites98.8%
lift-pow.f32N/A
lift-*.f32N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f32N/A
lower-+.f3298.8
Applied rewrites98.8%
Applied rewrites3.3%
Applied rewrites95.2%
Final simplification95.2%
(FPCore (cosTheta alpha) :precision binary32 (/ (/ -0.5 (PI)) (* (- 1.0 (* cosTheta cosTheta)) (log alpha))))
\begin{array}{l}
\\
\frac{\frac{-0.5}{\mathsf{PI}\left(\right)}}{\left(1 - cosTheta \cdot cosTheta\right) \cdot \log \alpha}
\end{array}
Initial program 98.4%
lift-*.f32N/A
lift-log.f32N/A
log-pow-revN/A
lower-log.f32N/A
lower-pow.f3298.8
Applied rewrites98.8%
Taylor expanded in alpha around 0
associate-/r*N/A
lower-/.f32N/A
lower-/.f32N/A
lower-PI.f32N/A
*-commutativeN/A
lower-*.f32N/A
mul-1-negN/A
unsub-negN/A
lower--.f32N/A
unpow2N/A
lower-*.f32N/A
lower-log.f3268.3
Applied rewrites68.3%
(FPCore (cosTheta alpha) :precision binary32 (/ 0.5 (* (* (log alpha) (PI)) (- (* cosTheta cosTheta) 1.0))))
\begin{array}{l}
\\
\frac{0.5}{\left(\log \alpha \cdot \mathsf{PI}\left(\right)\right) \cdot \left(cosTheta \cdot cosTheta - 1\right)}
\end{array}
Initial program 98.4%
Taylor expanded in alpha around 0
metadata-evalN/A
distribute-neg-fracN/A
distribute-frac-neg2N/A
distribute-lft-neg-inN/A
lower-/.f32N/A
distribute-lft-neg-inN/A
associate-*r*N/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-*.f32N/A
+-commutativeN/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
unpow2N/A
metadata-evalN/A
lower-fma.f32N/A
*-commutativeN/A
Applied rewrites10.3%
Applied rewrites68.2%
Final simplification68.2%
(FPCore (cosTheta alpha) :precision binary32 (/ (/ -0.5 (PI)) (log alpha)))
\begin{array}{l}
\\
\frac{\frac{-0.5}{\mathsf{PI}\left(\right)}}{\log \alpha}
\end{array}
Initial program 98.4%
Taylor expanded in alpha around 0
metadata-evalN/A
distribute-neg-fracN/A
distribute-frac-neg2N/A
distribute-lft-neg-inN/A
lower-/.f32N/A
distribute-lft-neg-inN/A
associate-*r*N/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-*.f32N/A
+-commutativeN/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
unpow2N/A
metadata-evalN/A
lower-fma.f32N/A
*-commutativeN/A
Applied rewrites10.3%
Taylor expanded in cosTheta around 0
Applied rewrites67.3%
Applied rewrites67.3%
(FPCore (cosTheta alpha) :precision binary32 (/ -0.5 (* (log alpha) (PI))))
\begin{array}{l}
\\
\frac{-0.5}{\log \alpha \cdot \mathsf{PI}\left(\right)}
\end{array}
Initial program 98.4%
Taylor expanded in alpha around 0
metadata-evalN/A
distribute-neg-fracN/A
distribute-frac-neg2N/A
distribute-lft-neg-inN/A
lower-/.f32N/A
distribute-lft-neg-inN/A
associate-*r*N/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-*.f32N/A
+-commutativeN/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
unpow2N/A
metadata-evalN/A
lower-fma.f32N/A
*-commutativeN/A
Applied rewrites10.7%
Taylor expanded in cosTheta around 0
Applied rewrites67.3%
(FPCore (cosTheta alpha) :precision binary32 (/ (- (/ (* alpha alpha) 2.0) -0.5) (PI)))
\begin{array}{l}
\\
\frac{\frac{\alpha \cdot \alpha}{2} - -0.5}{\mathsf{PI}\left(\right)}
\end{array}
Initial program 98.4%
lift-*.f32N/A
lift-log.f32N/A
log-pow-revN/A
lower-log.f32N/A
lower-pow.f3298.8
Applied rewrites98.8%
lift-pow.f32N/A
lift-*.f32N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f32N/A
lower-+.f3298.8
Applied rewrites98.8%
Applied rewrites3.3%
Applied rewrites26.6%
(FPCore (cosTheta alpha) :precision binary32 (/ (- 0.5 (/ (* alpha alpha) 2.0)) (PI)))
\begin{array}{l}
\\
\frac{0.5 - \frac{\alpha \cdot \alpha}{2}}{\mathsf{PI}\left(\right)}
\end{array}
Initial program 98.4%
lift-*.f32N/A
lift-log.f32N/A
log-pow-revN/A
lower-log.f32N/A
lower-pow.f3298.8
Applied rewrites98.8%
lift-pow.f32N/A
lift-*.f32N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f32N/A
lower-+.f3298.8
Applied rewrites98.8%
Applied rewrites3.3%
Applied rewrites25.8%
(FPCore (cosTheta alpha) :precision binary32 (/ (- (* alpha alpha) -0.5) (PI)))
\begin{array}{l}
\\
\frac{\alpha \cdot \alpha - -0.5}{\mathsf{PI}\left(\right)}
\end{array}
Initial program 98.4%
lift-*.f32N/A
lift-log.f32N/A
log-pow-revN/A
lower-log.f32N/A
lower-pow.f3298.8
Applied rewrites98.8%
lift-pow.f32N/A
lift-*.f32N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f32N/A
lower-+.f3298.8
Applied rewrites98.8%
Applied rewrites3.3%
Applied rewrites25.6%
(FPCore (cosTheta alpha) :precision binary32 (/ (- 0.5 (* alpha alpha)) (PI)))
\begin{array}{l}
\\
\frac{0.5 - \alpha \cdot \alpha}{\mathsf{PI}\left(\right)}
\end{array}
Initial program 98.4%
lift-*.f32N/A
lift-log.f32N/A
log-pow-revN/A
lower-log.f32N/A
lower-pow.f3298.8
Applied rewrites98.8%
lift-pow.f32N/A
lift-*.f32N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f32N/A
lower-+.f3298.8
Applied rewrites98.8%
Applied rewrites3.3%
Applied rewrites25.1%
(FPCore (cosTheta alpha) :precision binary32 (/ (fma alpha alpha 1.0) (PI)))
\begin{array}{l}
\\
\frac{\mathsf{fma}\left(\alpha, \alpha, 1\right)}{\mathsf{PI}\left(\right)}
\end{array}
Initial program 98.4%
lift-*.f32N/A
lift-log.f32N/A
log-pow-revN/A
lower-log.f32N/A
lower-pow.f3298.8
Applied rewrites98.8%
lift-pow.f32N/A
lift-*.f32N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f32N/A
lower-+.f3298.8
Applied rewrites98.8%
Applied rewrites3.3%
Applied rewrites23.9%
(FPCore (cosTheta alpha) :precision binary32 (let* ((t_0 (* (PI) alpha))) (* t_0 t_0)))
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{PI}\left(\right) \cdot \alpha\\
t\_0 \cdot t\_0
\end{array}
\end{array}
Initial program 98.4%
lift-*.f32N/A
lift-log.f32N/A
log-pow-revN/A
lower-log.f32N/A
lower-pow.f3298.8
Applied rewrites98.8%
lift-pow.f32N/A
lift-*.f32N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f32N/A
lower-+.f3298.8
Applied rewrites98.8%
Applied rewrites3.3%
Applied rewrites21.3%
(FPCore (cosTheta alpha) :precision binary32 (* (* (PI) alpha) alpha))
\begin{array}{l}
\\
\left(\mathsf{PI}\left(\right) \cdot \alpha\right) \cdot \alpha
\end{array}
Initial program 98.4%
lift-*.f32N/A
lift-log.f32N/A
log-pow-revN/A
lower-log.f32N/A
lower-pow.f3298.8
Applied rewrites98.8%
lift-pow.f32N/A
lift-*.f32N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f32N/A
lower-+.f3298.8
Applied rewrites98.8%
Applied rewrites3.3%
Applied rewrites21.0%
(FPCore (cosTheta alpha) :precision binary32 (* (PI) (* alpha alpha)))
\begin{array}{l}
\\
\mathsf{PI}\left(\right) \cdot \left(\alpha \cdot \alpha\right)
\end{array}
Initial program 98.4%
lift-*.f32N/A
lift-log.f32N/A
log-pow-revN/A
lower-log.f32N/A
lower-pow.f3298.8
Applied rewrites98.8%
lift-pow.f32N/A
lift-*.f32N/A
unpow-prod-downN/A
pow-prod-upN/A
lower-pow.f32N/A
lower-+.f3298.8
Applied rewrites98.8%
Applied rewrites3.3%
Applied rewrites21.0%
Final simplification21.0%
(FPCore (cosTheta alpha) :precision binary32 (PI))
\begin{array}{l}
\\
\mathsf{PI}\left(\right)
\end{array}
Initial program 98.4%
Taylor expanded in alpha around 0
metadata-evalN/A
distribute-neg-fracN/A
distribute-frac-neg2N/A
distribute-lft-neg-inN/A
lower-/.f32N/A
distribute-lft-neg-inN/A
associate-*r*N/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-*.f32N/A
+-commutativeN/A
distribute-neg-inN/A
mul-1-negN/A
remove-double-negN/A
unpow2N/A
metadata-evalN/A
lower-fma.f32N/A
*-commutativeN/A
Applied rewrites10.4%
Taylor expanded in cosTheta around 0
Applied rewrites67.3%
Applied rewrites3.2%
Applied rewrites20.0%
herbie shell --seed 2024312
(FPCore (cosTheta alpha)
:name "GTR1 distribution"
:precision binary32
:pre (and (and (<= 0.0 cosTheta) (<= cosTheta 1.0)) (and (<= 0.0001 alpha) (<= alpha 1.0)))
(/ (- (* alpha alpha) 1.0) (* (* (PI) (log (* alpha alpha))) (+ 1.0 (* (* (- (* alpha alpha) 1.0) cosTheta) cosTheta)))))