
(FPCore (a b) :precision binary64 (* (* (/ PI 2.0) (/ 1.0 (- (* b b) (* a a)))) (- (/ 1.0 a) (/ 1.0 b))))
double code(double a, double b) {
return ((((double) M_PI) / 2.0) * (1.0 / ((b * b) - (a * a)))) * ((1.0 / a) - (1.0 / b));
}
public static double code(double a, double b) {
return ((Math.PI / 2.0) * (1.0 / ((b * b) - (a * a)))) * ((1.0 / a) - (1.0 / b));
}
def code(a, b): return ((math.pi / 2.0) * (1.0 / ((b * b) - (a * a)))) * ((1.0 / a) - (1.0 / b))
function code(a, b) return Float64(Float64(Float64(pi / 2.0) * Float64(1.0 / Float64(Float64(b * b) - Float64(a * a)))) * Float64(Float64(1.0 / a) - Float64(1.0 / b))) end
function tmp = code(a, b) tmp = ((pi / 2.0) * (1.0 / ((b * b) - (a * a)))) * ((1.0 / a) - (1.0 / b)); end
code[a_, b_] := N[(N[(N[(Pi / 2.0), $MachinePrecision] * N[(1.0 / N[(N[(b * b), $MachinePrecision] - N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 / a), $MachinePrecision] - N[(1.0 / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\frac{\pi}{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 7 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a b) :precision binary64 (* (* (/ PI 2.0) (/ 1.0 (- (* b b) (* a a)))) (- (/ 1.0 a) (/ 1.0 b))))
double code(double a, double b) {
return ((((double) M_PI) / 2.0) * (1.0 / ((b * b) - (a * a)))) * ((1.0 / a) - (1.0 / b));
}
public static double code(double a, double b) {
return ((Math.PI / 2.0) * (1.0 / ((b * b) - (a * a)))) * ((1.0 / a) - (1.0 / b));
}
def code(a, b): return ((math.pi / 2.0) * (1.0 / ((b * b) - (a * a)))) * ((1.0 / a) - (1.0 / b))
function code(a, b) return Float64(Float64(Float64(pi / 2.0) * Float64(1.0 / Float64(Float64(b * b) - Float64(a * a)))) * Float64(Float64(1.0 / a) - Float64(1.0 / b))) end
function tmp = code(a, b) tmp = ((pi / 2.0) * (1.0 / ((b * b) - (a * a)))) * ((1.0 / a) - (1.0 / b)); end
code[a_, b_] := N[(N[(N[(Pi / 2.0), $MachinePrecision] * N[(1.0 / N[(N[(b * b), $MachinePrecision] - N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 / a), $MachinePrecision] - N[(1.0 / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\frac{\pi}{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 (* (/ PI (+ a b)) (/ 0.5 (* a b))))
double code(double a, double b) {
return (((double) M_PI) / (a + b)) * (0.5 / (a * b));
}
public static double code(double a, double b) {
return (Math.PI / (a + b)) * (0.5 / (a * b));
}
def code(a, b): return (math.pi / (a + b)) * (0.5 / (a * b))
function code(a, b) return Float64(Float64(pi / Float64(a + b)) * Float64(0.5 / Float64(a * b))) end
function tmp = code(a, b) tmp = (pi / (a + b)) * (0.5 / (a * b)); end
code[a_, b_] := N[(N[(Pi / N[(a + b), $MachinePrecision]), $MachinePrecision] * N[(0.5 / N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\pi}{a + b} \cdot \frac{0.5}{a \cdot b}
\end{array}
Initial program 77.2%
associate-*l*77.2%
*-rgt-identity77.2%
associate-/l*77.2%
metadata-eval77.2%
associate-*l/77.2%
*-lft-identity77.2%
sub-neg77.2%
distribute-neg-frac77.2%
metadata-eval77.2%
Simplified77.2%
metadata-eval77.2%
div-inv77.2%
*-commutative77.2%
clear-num77.0%
frac-times77.0%
*-un-lft-identity77.0%
frac-add76.9%
associate-/r/77.0%
*-un-lft-identity77.0%
*-commutative77.0%
neg-mul-177.0%
sub-neg77.0%
flip-+99.3%
+-commutative99.3%
Applied egg-rr99.3%
*-commutative99.3%
associate-/r*99.3%
Simplified99.3%
div-inv99.3%
metadata-eval99.3%
times-frac99.7%
Applied egg-rr99.7%
(FPCore (a b) :precision binary64 (if (<= a -1.3e-42) (/ (/ PI 2.0) (* a (* a b))) (/ (/ (* PI 0.5) b) (* a b))))
double code(double a, double b) {
double tmp;
if (a <= -1.3e-42) {
tmp = (((double) M_PI) / 2.0) / (a * (a * b));
} else {
tmp = ((((double) M_PI) * 0.5) / b) / (a * b);
}
return tmp;
}
public static double code(double a, double b) {
double tmp;
if (a <= -1.3e-42) {
tmp = (Math.PI / 2.0) / (a * (a * b));
} else {
tmp = ((Math.PI * 0.5) / b) / (a * b);
}
return tmp;
}
def code(a, b): tmp = 0 if a <= -1.3e-42: tmp = (math.pi / 2.0) / (a * (a * b)) else: tmp = ((math.pi * 0.5) / b) / (a * b) return tmp
function code(a, b) tmp = 0.0 if (a <= -1.3e-42) tmp = Float64(Float64(pi / 2.0) / Float64(a * Float64(a * b))); else tmp = Float64(Float64(Float64(pi * 0.5) / b) / Float64(a * b)); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (a <= -1.3e-42) tmp = (pi / 2.0) / (a * (a * b)); else tmp = ((pi * 0.5) / b) / (a * b); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[a, -1.3e-42], N[(N[(Pi / 2.0), $MachinePrecision] / N[(a * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(Pi * 0.5), $MachinePrecision] / b), $MachinePrecision] / N[(a * b), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.3 \cdot 10^{-42}:\\
\;\;\;\;\frac{\frac{\pi}{2}}{a \cdot \left(a \cdot b\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\pi \cdot 0.5}{b}}{a \cdot b}\\
\end{array}
\end{array}
if a < -1.3e-42Initial program 85.7%
associate-*l*85.8%
*-rgt-identity85.8%
associate-/l*85.8%
metadata-eval85.8%
associate-*l/85.7%
*-lft-identity85.7%
sub-neg85.7%
distribute-neg-frac85.7%
metadata-eval85.7%
Simplified85.7%
metadata-eval85.7%
div-inv85.7%
*-commutative85.7%
clear-num85.8%
frac-times85.8%
*-un-lft-identity85.8%
frac-add85.7%
associate-/r/85.8%
*-un-lft-identity85.8%
*-commutative85.8%
neg-mul-185.8%
sub-neg85.8%
flip-+99.6%
+-commutative99.6%
Applied egg-rr99.6%
*-commutative99.6%
associate-/r*99.6%
Simplified99.6%
Taylor expanded in a around inf 79.8%
if -1.3e-42 < a Initial program 73.8%
associate-*l*73.8%
*-rgt-identity73.8%
associate-/l*73.8%
metadata-eval73.8%
associate-*l/73.8%
*-lft-identity73.8%
sub-neg73.8%
distribute-neg-frac73.8%
metadata-eval73.8%
Simplified73.8%
metadata-eval73.8%
div-inv73.8%
*-commutative73.8%
clear-num73.4%
frac-times73.4%
*-un-lft-identity73.4%
frac-add73.4%
associate-/r/73.5%
*-un-lft-identity73.5%
*-commutative73.5%
neg-mul-173.5%
sub-neg73.5%
flip-+99.2%
+-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-/r*99.2%
Simplified99.2%
Taylor expanded in a around 0 72.1%
*-un-lft-identity72.1%
*-commutative72.1%
times-frac72.0%
div-inv72.0%
metadata-eval72.0%
associate-*r/72.0%
Applied egg-rr72.0%
associate-*l/72.0%
*-lft-identity72.0%
associate-*r/72.1%
*-commutative72.1%
Simplified72.1%
Final simplification74.3%
(FPCore (a b) :precision binary64 (if (<= a -1.3e-42) (/ (/ PI 2.0) (* a (* a b))) (* (/ 0.5 b) (/ (/ PI b) a))))
double code(double a, double b) {
double tmp;
if (a <= -1.3e-42) {
tmp = (((double) M_PI) / 2.0) / (a * (a * b));
} else {
tmp = (0.5 / b) * ((((double) M_PI) / b) / a);
}
return tmp;
}
public static double code(double a, double b) {
double tmp;
if (a <= -1.3e-42) {
tmp = (Math.PI / 2.0) / (a * (a * b));
} else {
tmp = (0.5 / b) * ((Math.PI / b) / a);
}
return tmp;
}
def code(a, b): tmp = 0 if a <= -1.3e-42: tmp = (math.pi / 2.0) / (a * (a * b)) else: tmp = (0.5 / b) * ((math.pi / b) / a) return tmp
function code(a, b) tmp = 0.0 if (a <= -1.3e-42) tmp = Float64(Float64(pi / 2.0) / Float64(a * Float64(a * b))); else tmp = Float64(Float64(0.5 / b) * Float64(Float64(pi / b) / a)); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (a <= -1.3e-42) tmp = (pi / 2.0) / (a * (a * b)); else tmp = (0.5 / b) * ((pi / b) / a); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[a, -1.3e-42], N[(N[(Pi / 2.0), $MachinePrecision] / N[(a * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.5 / b), $MachinePrecision] * N[(N[(Pi / b), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.3 \cdot 10^{-42}:\\
\;\;\;\;\frac{\frac{\pi}{2}}{a \cdot \left(a \cdot b\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.5}{b} \cdot \frac{\frac{\pi}{b}}{a}\\
\end{array}
\end{array}
if a < -1.3e-42Initial program 85.7%
associate-*l*85.8%
*-rgt-identity85.8%
associate-/l*85.8%
metadata-eval85.8%
associate-*l/85.7%
*-lft-identity85.7%
sub-neg85.7%
distribute-neg-frac85.7%
metadata-eval85.7%
Simplified85.7%
metadata-eval85.7%
div-inv85.7%
*-commutative85.7%
clear-num85.8%
frac-times85.8%
*-un-lft-identity85.8%
frac-add85.7%
associate-/r/85.8%
*-un-lft-identity85.8%
*-commutative85.8%
neg-mul-185.8%
sub-neg85.8%
flip-+99.6%
+-commutative99.6%
Applied egg-rr99.6%
*-commutative99.6%
associate-/r*99.6%
Simplified99.6%
Taylor expanded in a around inf 79.8%
if -1.3e-42 < a Initial program 73.8%
associate-*l*73.8%
*-rgt-identity73.8%
associate-/l*73.8%
metadata-eval73.8%
associate-*l/73.8%
*-lft-identity73.8%
sub-neg73.8%
distribute-neg-frac73.8%
metadata-eval73.8%
Simplified73.8%
metadata-eval73.8%
div-inv73.8%
*-commutative73.8%
clear-num73.4%
frac-times73.4%
*-un-lft-identity73.4%
frac-add73.4%
associate-/r/73.5%
*-un-lft-identity73.5%
*-commutative73.5%
neg-mul-173.5%
sub-neg73.5%
flip-+99.2%
+-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-/r*99.2%
Simplified99.2%
Taylor expanded in a around 0 72.1%
*-un-lft-identity72.1%
*-commutative72.1%
times-frac72.0%
div-inv72.0%
metadata-eval72.0%
associate-*r/72.0%
Applied egg-rr72.0%
associate-*l/72.0%
*-lft-identity72.0%
associate-*r/72.1%
*-commutative72.1%
Simplified72.1%
associate-/l*72.1%
*-commutative72.1%
times-frac72.1%
Applied egg-rr72.1%
(FPCore (a b) :precision binary64 (if (<= a -4.2e-43) (/ 0.5 (* (* a b) (/ a PI))) (* (/ 0.5 b) (/ (/ PI b) a))))
double code(double a, double b) {
double tmp;
if (a <= -4.2e-43) {
tmp = 0.5 / ((a * b) * (a / ((double) M_PI)));
} else {
tmp = (0.5 / b) * ((((double) M_PI) / b) / a);
}
return tmp;
}
public static double code(double a, double b) {
double tmp;
if (a <= -4.2e-43) {
tmp = 0.5 / ((a * b) * (a / Math.PI));
} else {
tmp = (0.5 / b) * ((Math.PI / b) / a);
}
return tmp;
}
def code(a, b): tmp = 0 if a <= -4.2e-43: tmp = 0.5 / ((a * b) * (a / math.pi)) else: tmp = (0.5 / b) * ((math.pi / b) / a) return tmp
function code(a, b) tmp = 0.0 if (a <= -4.2e-43) tmp = Float64(0.5 / Float64(Float64(a * b) * Float64(a / pi))); else tmp = Float64(Float64(0.5 / b) * Float64(Float64(pi / b) / a)); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (a <= -4.2e-43) tmp = 0.5 / ((a * b) * (a / pi)); else tmp = (0.5 / b) * ((pi / b) / a); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[a, -4.2e-43], N[(0.5 / N[(N[(a * b), $MachinePrecision] * N[(a / Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.5 / b), $MachinePrecision] * N[(N[(Pi / b), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -4.2 \cdot 10^{-43}:\\
\;\;\;\;\frac{0.5}{\left(a \cdot b\right) \cdot \frac{a}{\pi}}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.5}{b} \cdot \frac{\frac{\pi}{b}}{a}\\
\end{array}
\end{array}
if a < -4.2000000000000001e-43Initial program 85.7%
associate-*l*85.8%
*-rgt-identity85.8%
associate-/l*85.8%
metadata-eval85.8%
associate-*l/85.7%
*-lft-identity85.7%
sub-neg85.7%
distribute-neg-frac85.7%
metadata-eval85.7%
Simplified85.7%
metadata-eval85.7%
div-inv85.7%
*-commutative85.7%
clear-num85.8%
frac-times85.8%
*-un-lft-identity85.8%
frac-add85.7%
associate-/r/85.8%
*-un-lft-identity85.8%
*-commutative85.8%
neg-mul-185.8%
sub-neg85.8%
flip-+99.6%
+-commutative99.6%
Applied egg-rr99.6%
*-commutative99.6%
associate-/r*99.6%
Simplified99.6%
div-inv99.6%
metadata-eval99.6%
times-frac99.6%
Applied egg-rr99.6%
Taylor expanded in a around inf 78.7%
clear-num78.6%
frac-times79.7%
metadata-eval79.7%
Applied egg-rr79.7%
if -4.2000000000000001e-43 < a Initial program 73.8%
associate-*l*73.8%
*-rgt-identity73.8%
associate-/l*73.8%
metadata-eval73.8%
associate-*l/73.8%
*-lft-identity73.8%
sub-neg73.8%
distribute-neg-frac73.8%
metadata-eval73.8%
Simplified73.8%
metadata-eval73.8%
div-inv73.8%
*-commutative73.8%
clear-num73.4%
frac-times73.4%
*-un-lft-identity73.4%
frac-add73.4%
associate-/r/73.5%
*-un-lft-identity73.5%
*-commutative73.5%
neg-mul-173.5%
sub-neg73.5%
flip-+99.2%
+-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-/r*99.2%
Simplified99.2%
Taylor expanded in a around 0 72.1%
*-un-lft-identity72.1%
*-commutative72.1%
times-frac72.0%
div-inv72.0%
metadata-eval72.0%
associate-*r/72.0%
Applied egg-rr72.0%
associate-*l/72.0%
*-lft-identity72.0%
associate-*r/72.1%
*-commutative72.1%
Simplified72.1%
associate-/l*72.1%
*-commutative72.1%
times-frac72.1%
Applied egg-rr72.1%
Final simplification74.3%
(FPCore (a b) :precision binary64 (if (<= a -8.5e-43) (* PI (/ 0.5 (* a (* a b)))) (* (/ 0.5 b) (/ (/ PI b) a))))
double code(double a, double b) {
double tmp;
if (a <= -8.5e-43) {
tmp = ((double) M_PI) * (0.5 / (a * (a * b)));
} else {
tmp = (0.5 / b) * ((((double) M_PI) / b) / a);
}
return tmp;
}
public static double code(double a, double b) {
double tmp;
if (a <= -8.5e-43) {
tmp = Math.PI * (0.5 / (a * (a * b)));
} else {
tmp = (0.5 / b) * ((Math.PI / b) / a);
}
return tmp;
}
def code(a, b): tmp = 0 if a <= -8.5e-43: tmp = math.pi * (0.5 / (a * (a * b))) else: tmp = (0.5 / b) * ((math.pi / b) / a) return tmp
function code(a, b) tmp = 0.0 if (a <= -8.5e-43) tmp = Float64(pi * Float64(0.5 / Float64(a * Float64(a * b)))); else tmp = Float64(Float64(0.5 / b) * Float64(Float64(pi / b) / a)); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (a <= -8.5e-43) tmp = pi * (0.5 / (a * (a * b))); else tmp = (0.5 / b) * ((pi / b) / a); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[a, -8.5e-43], N[(Pi * N[(0.5 / N[(a * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.5 / b), $MachinePrecision] * N[(N[(Pi / b), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -8.5 \cdot 10^{-43}:\\
\;\;\;\;\pi \cdot \frac{0.5}{a \cdot \left(a \cdot b\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.5}{b} \cdot \frac{\frac{\pi}{b}}{a}\\
\end{array}
\end{array}
if a < -8.50000000000000056e-43Initial program 85.7%
*-commutative85.7%
associate-*r*85.8%
associate-*r/85.8%
associate-*r*85.8%
*-rgt-identity85.8%
sub-neg85.8%
distribute-neg-frac85.8%
metadata-eval85.8%
Simplified85.8%
*-commutative85.8%
associate-*r/85.7%
div-inv85.7%
metadata-eval85.7%
associate-*l*85.7%
*-commutative85.7%
Applied egg-rr99.5%
Taylor expanded in a around inf 79.6%
if -8.50000000000000056e-43 < a Initial program 73.8%
associate-*l*73.8%
*-rgt-identity73.8%
associate-/l*73.8%
metadata-eval73.8%
associate-*l/73.8%
*-lft-identity73.8%
sub-neg73.8%
distribute-neg-frac73.8%
metadata-eval73.8%
Simplified73.8%
metadata-eval73.8%
div-inv73.8%
*-commutative73.8%
clear-num73.4%
frac-times73.4%
*-un-lft-identity73.4%
frac-add73.4%
associate-/r/73.5%
*-un-lft-identity73.5%
*-commutative73.5%
neg-mul-173.5%
sub-neg73.5%
flip-+99.2%
+-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-/r*99.2%
Simplified99.2%
Taylor expanded in a around 0 72.1%
*-un-lft-identity72.1%
*-commutative72.1%
times-frac72.0%
div-inv72.0%
metadata-eval72.0%
associate-*r/72.0%
Applied egg-rr72.0%
associate-*l/72.0%
*-lft-identity72.0%
associate-*r/72.1%
*-commutative72.1%
Simplified72.1%
associate-/l*72.1%
*-commutative72.1%
times-frac72.1%
Applied egg-rr72.1%
Final simplification74.2%
(FPCore (a b) :precision binary64 (* PI (/ 0.5 (* (+ a b) (* a b)))))
double code(double a, double b) {
return ((double) M_PI) * (0.5 / ((a + b) * (a * b)));
}
public static double code(double a, double b) {
return Math.PI * (0.5 / ((a + b) * (a * b)));
}
def code(a, b): return math.pi * (0.5 / ((a + b) * (a * b)))
function code(a, b) return Float64(pi * Float64(0.5 / Float64(Float64(a + b) * Float64(a * b)))) end
function tmp = code(a, b) tmp = pi * (0.5 / ((a + b) * (a * b))); end
code[a_, b_] := N[(Pi * N[(0.5 / N[(N[(a + b), $MachinePrecision] * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\pi \cdot \frac{0.5}{\left(a + b\right) \cdot \left(a \cdot b\right)}
\end{array}
Initial program 77.2%
*-commutative77.2%
associate-*r*77.2%
associate-*r/77.3%
associate-*r*77.3%
*-rgt-identity77.3%
sub-neg77.3%
distribute-neg-frac77.3%
metadata-eval77.3%
Simplified77.3%
*-commutative77.3%
associate-*r/77.2%
div-inv77.2%
metadata-eval77.2%
associate-*l*77.2%
*-commutative77.2%
Applied egg-rr99.3%
Final simplification99.3%
(FPCore (a b) :precision binary64 (* (/ 0.5 b) (/ (/ PI b) a)))
double code(double a, double b) {
return (0.5 / b) * ((((double) M_PI) / b) / a);
}
public static double code(double a, double b) {
return (0.5 / b) * ((Math.PI / b) / a);
}
def code(a, b): return (0.5 / b) * ((math.pi / b) / a)
function code(a, b) return Float64(Float64(0.5 / b) * Float64(Float64(pi / b) / a)) end
function tmp = code(a, b) tmp = (0.5 / b) * ((pi / b) / a); end
code[a_, b_] := N[(N[(0.5 / b), $MachinePrecision] * N[(N[(Pi / b), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.5}{b} \cdot \frac{\frac{\pi}{b}}{a}
\end{array}
Initial program 77.2%
associate-*l*77.2%
*-rgt-identity77.2%
associate-/l*77.2%
metadata-eval77.2%
associate-*l/77.2%
*-lft-identity77.2%
sub-neg77.2%
distribute-neg-frac77.2%
metadata-eval77.2%
Simplified77.2%
metadata-eval77.2%
div-inv77.2%
*-commutative77.2%
clear-num77.0%
frac-times77.0%
*-un-lft-identity77.0%
frac-add76.9%
associate-/r/77.0%
*-un-lft-identity77.0%
*-commutative77.0%
neg-mul-177.0%
sub-neg77.0%
flip-+99.3%
+-commutative99.3%
Applied egg-rr99.3%
*-commutative99.3%
associate-/r*99.3%
Simplified99.3%
Taylor expanded in a around 0 61.9%
*-un-lft-identity61.9%
*-commutative61.9%
times-frac61.7%
div-inv61.7%
metadata-eval61.7%
associate-*r/61.7%
Applied egg-rr61.7%
associate-*l/61.7%
*-lft-identity61.7%
associate-*r/61.7%
*-commutative61.7%
Simplified61.7%
associate-/l*61.7%
*-commutative61.7%
times-frac61.7%
Applied egg-rr61.7%
herbie shell --seed 2024150
(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))))