
(FPCore (a b) :precision binary64 (* (* (/ (PI) 2.0) (/ 1.0 (- (* b b) (* a a)))) (- (/ 1.0 a) (/ 1.0 b))))
\begin{array}{l}
\\
\left(\frac{\mathsf{PI}\left(\right)}{2} \cdot \frac{1}{b \cdot b - a \cdot a}\right) \cdot \left(\frac{1}{a} - \frac{1}{b}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 9 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a b) :precision binary64 (* (* (/ (PI) 2.0) (/ 1.0 (- (* b b) (* a a)))) (- (/ 1.0 a) (/ 1.0 b))))
\begin{array}{l}
\\
\left(\frac{\mathsf{PI}\left(\right)}{2} \cdot \frac{1}{b \cdot b - a \cdot a}\right) \cdot \left(\frac{1}{a} - \frac{1}{b}\right)
\end{array}
(FPCore (a b) :precision binary64 (/ (/ (* 0.5 (PI)) (* a b)) (+ a b)))
\begin{array}{l}
\\
\frac{\frac{0.5 \cdot \mathsf{PI}\left(\right)}{a \cdot b}}{a + b}
\end{array}
Initial program 76.6%
Applied rewrites88.7%
Applied rewrites99.7%
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
*-rgt-identityN/A
lift-/.f64N/A
associate-*l/N/A
metadata-evalN/A
div-invN/A
lift-PI.f64N/A
lower-/.f64N/A
lift-PI.f64N/A
div-invN/A
metadata-evalN/A
lower-*.f6499.7
lift-*.f64N/A
*-commutativeN/A
lower-*.f6499.7
Applied rewrites99.7%
Final simplification99.7%
(FPCore (a b) :precision binary64 (if (<= a -5e+121) (* (/ (PI) (* (* a b) a)) 0.5) (* (/ 0.5 (* (* (+ a b) a) b)) (PI))))
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -5 \cdot 10^{+121}:\\
\;\;\;\;\frac{\mathsf{PI}\left(\right)}{\left(a \cdot b\right) \cdot a} \cdot 0.5\\
\mathbf{else}:\\
\;\;\;\;\frac{0.5}{\left(\left(a + b\right) \cdot a\right) \cdot b} \cdot \mathsf{PI}\left(\right)\\
\end{array}
\end{array}
if a < -5.00000000000000007e121Initial program 62.3%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-PI.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6475.5
Applied rewrites75.5%
Applied rewrites99.8%
if -5.00000000000000007e121 < a Initial program 79.1%
Applied rewrites90.9%
Applied rewrites99.7%
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
*-rgt-identityN/A
lift-/.f64N/A
associate-*l/N/A
metadata-evalN/A
div-invN/A
lift-PI.f64N/A
lower-/.f64N/A
lift-PI.f64N/A
div-invN/A
metadata-evalN/A
lower-*.f6499.7
lift-*.f64N/A
*-commutativeN/A
lower-*.f6499.7
Applied rewrites99.7%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
lift-*.f64N/A
associate-/l*N/A
lower-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lift-+.f64N/A
*-commutativeN/A
lower-/.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-+.f64N/A
+-commutativeN/A
distribute-rgt-outN/A
lower-*.f64N/A
distribute-rgt-outN/A
+-commutativeN/A
lift-+.f64N/A
lower-*.f6495.4
lift-+.f64N/A
+-commutativeN/A
lift-+.f6495.4
Applied rewrites95.4%
Final simplification96.0%
(FPCore (a b) :precision binary64 (/ (/ (* 0.5 (PI)) (+ a b)) (* a b)))
\begin{array}{l}
\\
\frac{\frac{0.5 \cdot \mathsf{PI}\left(\right)}{a + b}}{a \cdot b}
\end{array}
Initial program 76.6%
Applied rewrites88.7%
Applied rewrites99.7%
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
*-rgt-identityN/A
lift-/.f64N/A
associate-*l/N/A
metadata-evalN/A
div-invN/A
lift-PI.f64N/A
lower-/.f64N/A
lift-PI.f64N/A
div-invN/A
metadata-evalN/A
lower-*.f6499.7
lift-*.f64N/A
*-commutativeN/A
lower-*.f6499.7
Applied rewrites99.7%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
lift-+.f64N/A
+-commutativeN/A
lift-+.f64N/A
associate-/r*N/A
lower-/.f64N/A
lower-/.f6499.6
lift-*.f64N/A
*-commutativeN/A
lower-*.f6499.6
lift-+.f64N/A
+-commutativeN/A
lift-+.f6499.6
lift-*.f64N/A
*-commutativeN/A
lower-*.f6499.6
Applied rewrites99.6%
(FPCore (a b) :precision binary64 (* (/ 0.5 (+ a b)) (/ (PI) (* a b))))
\begin{array}{l}
\\
\frac{0.5}{a + b} \cdot \frac{\mathsf{PI}\left(\right)}{a \cdot b}
\end{array}
Initial program 76.6%
Applied rewrites88.7%
Applied rewrites99.7%
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lift-*.f64N/A
*-rgt-identityN/A
lift-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6499.6
lift-+.f64N/A
+-commutativeN/A
lower-+.f6499.6
Applied rewrites99.6%
Final simplification99.6%
(FPCore (a b) :precision binary64 (if (<= a -2.9e-30) (* (/ (PI) (* (* a b) a)) 0.5) (* (/ (PI) (* (* b b) a)) 0.5)))
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -2.9 \cdot 10^{-30}:\\
\;\;\;\;\frac{\mathsf{PI}\left(\right)}{\left(a \cdot b\right) \cdot a} \cdot 0.5\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{PI}\left(\right)}{\left(b \cdot b\right) \cdot a} \cdot 0.5\\
\end{array}
\end{array}
if a < -2.89999999999999989e-30Initial program 76.3%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-PI.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6479.2
Applied rewrites79.2%
Applied rewrites94.3%
if -2.89999999999999989e-30 < a Initial program 76.7%
Taylor expanded in a around 0
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-PI.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6460.9
Applied rewrites60.9%
(FPCore (a b) :precision binary64 (/ (PI) (* (* 2.0 (+ a b)) (* a b))))
\begin{array}{l}
\\
\frac{\mathsf{PI}\left(\right)}{\left(2 \cdot \left(a + b\right)\right) \cdot \left(a \cdot b\right)}
\end{array}
Initial program 76.6%
Applied rewrites88.7%
Applied rewrites99.7%
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
*-rgt-identityN/A
lift-/.f64N/A
associate-*l/N/A
metadata-evalN/A
div-invN/A
lift-PI.f64N/A
lower-/.f64N/A
lift-PI.f64N/A
div-invN/A
metadata-evalN/A
lower-*.f6499.7
lift-*.f64N/A
*-commutativeN/A
lower-*.f6499.7
Applied rewrites99.7%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
lift-*.f64N/A
lift-+.f64N/A
+-commutativeN/A
lift-+.f64N/A
metadata-evalN/A
div-invN/A
*-commutativeN/A
associate-/r*N/A
associate-*l*N/A
lift-*.f64N/A
lift-*.f64N/A
lower-/.f6499.4
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
Applied rewrites99.4%
Final simplification99.4%
(FPCore (a b) :precision binary64 (* (/ (PI) (* (* a b) a)) 0.5))
\begin{array}{l}
\\
\frac{\mathsf{PI}\left(\right)}{\left(a \cdot b\right) \cdot a} \cdot 0.5
\end{array}
Initial program 76.6%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-PI.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6458.8
Applied rewrites58.8%
Applied rewrites65.7%
(FPCore (a b) :precision binary64 (* (/ 0.5 (* (* a b) a)) (PI)))
\begin{array}{l}
\\
\frac{0.5}{\left(a \cdot b\right) \cdot a} \cdot \mathsf{PI}\left(\right)
\end{array}
Initial program 76.6%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-PI.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6458.8
Applied rewrites58.8%
Applied rewrites58.8%
Applied rewrites65.7%
Final simplification65.7%
(FPCore (a b) :precision binary64 (* (/ 0.5 (* (* a a) b)) (PI)))
\begin{array}{l}
\\
\frac{0.5}{\left(a \cdot a\right) \cdot b} \cdot \mathsf{PI}\left(\right)
\end{array}
Initial program 76.6%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-PI.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6458.8
Applied rewrites58.8%
Applied rewrites58.8%
Final simplification58.8%
herbie shell --seed 2024295
(FPCore (a b)
:name "NMSE Section 6.1 mentioned, B"
:precision binary64
(* (* (/ (PI) 2.0) (/ 1.0 (- (* b b) (* a a)))) (- (/ 1.0 a) (/ 1.0 b))))