
(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 5 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)) (* (+ b a) (* b a))))
\begin{array}{l}
\\
\frac{0.5 \cdot \mathsf{PI}\left(\right)}{\left(b + a\right) \cdot \left(b \cdot a\right)}
\end{array}
Initial program 83.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
difference-of-squaresN/A
times-fracN/A
*-rgt-identityN/A
lower-*.f64N/A
lower-/.f64N/A
*-rgt-identityN/A
+-commutativeN/A
lower-+.f64N/A
*-lft-identityN/A
*-rgt-identityN/A
lower-/.f64N/A
*-lft-identityN/A
*-rgt-identityN/A
lower--.f6492.7
Applied rewrites92.7%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
lift--.f64N/A
lift-/.f64N/A
lift-/.f64N/A
frac-subN/A
frac-timesN/A
lower-/.f64N/A
Applied rewrites99.3%
Taylor expanded in a around 0
Applied rewrites99.3%
(FPCore (a b) :precision binary64 (if (or (<= a -1.06e-64) (not (<= a 56000000.0))) (/ (* (PI) 0.5) (* (* a a) b)) (* (/ (PI) (* (* a b) b)) 0.5)))
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.06 \cdot 10^{-64} \lor \neg \left(a \leq 56000000\right):\\
\;\;\;\;\frac{\mathsf{PI}\left(\right) \cdot 0.5}{\left(a \cdot a\right) \cdot b}\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{PI}\left(\right)}{\left(a \cdot b\right) \cdot b} \cdot 0.5\\
\end{array}
\end{array}
if a < -1.06000000000000007e-64 or 5.6e7 < a Initial program 81.0%
Taylor expanded in a around inf
Applied rewrites81.4%
Applied rewrites81.8%
if -1.06000000000000007e-64 < a < 5.6e7Initial program 86.2%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
difference-of-squaresN/A
times-fracN/A
*-rgt-identityN/A
lower-*.f64N/A
lower-/.f64N/A
*-rgt-identityN/A
+-commutativeN/A
lower-+.f64N/A
*-lft-identityN/A
*-rgt-identityN/A
lower-/.f64N/A
*-lft-identityN/A
*-rgt-identityN/A
lower--.f6492.6
Applied rewrites92.6%
Taylor expanded in a around 0
Applied rewrites80.6%
Applied rewrites87.8%
Final simplification84.5%
(FPCore (a b) :precision binary64 (if (<= b -1.2e+132) (* (PI) (/ 0.5 (* (* a b) b))) (/ (* 0.5 (PI)) (* (* (+ a b) b) a))))
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.2 \cdot 10^{+132}:\\
\;\;\;\;\mathsf{PI}\left(\right) \cdot \frac{0.5}{\left(a \cdot b\right) \cdot b}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.5 \cdot \mathsf{PI}\left(\right)}{\left(\left(a + b\right) \cdot b\right) \cdot a}\\
\end{array}
\end{array}
if b < -1.2000000000000001e132Initial program 59.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
difference-of-squaresN/A
times-fracN/A
*-rgt-identityN/A
lower-*.f64N/A
lower-/.f64N/A
*-rgt-identityN/A
+-commutativeN/A
lower-+.f64N/A
*-lft-identityN/A
*-rgt-identityN/A
lower-/.f64N/A
*-lft-identityN/A
*-rgt-identityN/A
lower--.f6483.8
Applied rewrites83.8%
Taylor expanded in a around 0
Applied rewrites81.5%
Applied rewrites81.5%
Applied rewrites99.6%
if -1.2000000000000001e132 < b Initial program 88.5%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
difference-of-squaresN/A
times-fracN/A
*-rgt-identityN/A
lower-*.f64N/A
lower-/.f64N/A
*-rgt-identityN/A
+-commutativeN/A
lower-+.f64N/A
*-lft-identityN/A
*-rgt-identityN/A
lower-/.f64N/A
*-lft-identityN/A
*-rgt-identityN/A
lower--.f6494.6
Applied rewrites94.6%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
lift--.f64N/A
lift-/.f64N/A
lift-/.f64N/A
frac-subN/A
frac-timesN/A
lower-/.f64N/A
Applied rewrites99.2%
Taylor expanded in a around 0
Applied rewrites99.2%
lift-*.f64N/A
lift-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6498.8
lift-+.f64N/A
+-commutativeN/A
lift-+.f6498.8
Applied rewrites98.8%
(FPCore (a b) :precision binary64 (* (/ (PI) (* (* a b) b)) 0.5))
\begin{array}{l}
\\
\frac{\mathsf{PI}\left(\right)}{\left(a \cdot b\right) \cdot b} \cdot 0.5
\end{array}
Initial program 83.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
difference-of-squaresN/A
times-fracN/A
*-rgt-identityN/A
lower-*.f64N/A
lower-/.f64N/A
*-rgt-identityN/A
+-commutativeN/A
lower-+.f64N/A
*-lft-identityN/A
*-rgt-identityN/A
lower-/.f64N/A
*-lft-identityN/A
*-rgt-identityN/A
lower--.f6492.7
Applied rewrites92.7%
Taylor expanded in a around 0
Applied rewrites62.9%
Applied rewrites66.6%
(FPCore (a b) :precision binary64 (* (PI) (/ 0.5 (* (* a b) b))))
\begin{array}{l}
\\
\mathsf{PI}\left(\right) \cdot \frac{0.5}{\left(a \cdot b\right) \cdot b}
\end{array}
Initial program 83.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
difference-of-squaresN/A
times-fracN/A
*-rgt-identityN/A
lower-*.f64N/A
lower-/.f64N/A
*-rgt-identityN/A
+-commutativeN/A
lower-+.f64N/A
*-lft-identityN/A
*-rgt-identityN/A
lower-/.f64N/A
*-lft-identityN/A
*-rgt-identityN/A
lower--.f6492.7
Applied rewrites92.7%
Taylor expanded in a around 0
Applied rewrites62.9%
Applied rewrites62.9%
Applied rewrites66.5%
herbie shell --seed 2025019
(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))))