
(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 8 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 (/ (* 0.5 (/ PI (+ a b))) (* a b)))
double code(double a, double b) {
return (0.5 * (((double) M_PI) / (a + b))) / (a * b);
}
public static double code(double a, double b) {
return (0.5 * (Math.PI / (a + b))) / (a * b);
}
def code(a, b): return (0.5 * (math.pi / (a + b))) / (a * b)
function code(a, b) return Float64(Float64(0.5 * Float64(pi / Float64(a + b))) / Float64(a * b)) end
function tmp = code(a, b) tmp = (0.5 * (pi / (a + b))) / (a * b); end
code[a_, b_] := N[(N[(0.5 * N[(Pi / N[(a + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(a * b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.5 \cdot \frac{\pi}{a + b}}{a \cdot b}
\end{array}
Initial program 84.5%
*-commutative84.5%
associate-*r*84.5%
associate-*r/84.5%
associate-*r*84.5%
*-rgt-identity84.5%
sub-neg84.5%
distribute-neg-frac84.5%
metadata-eval84.5%
Simplified84.5%
associate-*l/84.5%
div-inv84.5%
metadata-eval84.5%
*-commutative84.5%
associate-*r*84.5%
*-commutative84.5%
Applied egg-rr99.2%
associate-*l/99.2%
*-commutative99.2%
metadata-eval99.2%
div-inv99.2%
associate-/r*99.7%
div-inv99.7%
metadata-eval99.7%
*-commutative99.7%
*-un-lft-identity99.7%
times-frac99.7%
metadata-eval99.7%
Applied egg-rr99.7%
(FPCore (a b) :precision binary64 (if (<= a -9.2e-47) (* (/ 0.5 a) (/ PI (* a b))) (* PI (/ (/ 0.5 b) (* a b)))))
double code(double a, double b) {
double tmp;
if (a <= -9.2e-47) {
tmp = (0.5 / a) * (((double) M_PI) / (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 <= -9.2e-47) {
tmp = (0.5 / a) * (Math.PI / (a * b));
} else {
tmp = Math.PI * ((0.5 / b) / (a * b));
}
return tmp;
}
def code(a, b): tmp = 0 if a <= -9.2e-47: tmp = (0.5 / a) * (math.pi / (a * b)) else: tmp = math.pi * ((0.5 / b) / (a * b)) return tmp
function code(a, b) tmp = 0.0 if (a <= -9.2e-47) tmp = Float64(Float64(0.5 / a) * Float64(pi / Float64(a * b))); else tmp = Float64(pi * Float64(Float64(0.5 / b) / Float64(a * b))); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (a <= -9.2e-47) tmp = (0.5 / a) * (pi / (a * b)); else tmp = pi * ((0.5 / b) / (a * b)); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[a, -9.2e-47], N[(N[(0.5 / a), $MachinePrecision] * N[(Pi / N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(Pi * N[(N[(0.5 / b), $MachinePrecision] / N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -9.2 \cdot 10^{-47}:\\
\;\;\;\;\frac{0.5}{a} \cdot \frac{\pi}{a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;\pi \cdot \frac{\frac{0.5}{b}}{a \cdot b}\\
\end{array}
\end{array}
if a < -9.19999999999999928e-47Initial program 80.0%
associate-*l*80.0%
*-rgt-identity80.0%
associate-/l*80.0%
metadata-eval80.0%
associate-*l/80.1%
*-lft-identity80.1%
sub-neg80.1%
distribute-neg-frac80.1%
metadata-eval80.1%
Simplified80.1%
metadata-eval80.1%
div-inv80.1%
*-commutative80.1%
clear-num80.0%
frac-times80.0%
*-un-lft-identity80.0%
frac-add80.0%
associate-/r/80.1%
*-un-lft-identity80.1%
*-commutative80.1%
neg-mul-180.1%
sub-neg80.1%
flip-+98.9%
+-commutative98.9%
Applied egg-rr98.9%
*-commutative98.9%
associate-/r*98.9%
Simplified98.9%
associate-/l/98.9%
div-inv99.0%
metadata-eval99.0%
div-inv99.0%
clear-num99.0%
associate-/r*99.5%
Applied egg-rr99.5%
associate-*r/99.6%
associate-*l/99.7%
*-commutative99.7%
associate-/r*99.7%
Simplified99.7%
Taylor expanded in a around inf 95.0%
Taylor expanded in a around 0 95.0%
if -9.19999999999999928e-47 < a Initial program 86.0%
associate-*l*86.0%
*-rgt-identity86.0%
associate-/l*86.0%
metadata-eval86.0%
associate-*l/86.1%
*-lft-identity86.1%
sub-neg86.1%
distribute-neg-frac86.1%
metadata-eval86.1%
Simplified86.1%
metadata-eval86.1%
div-inv86.1%
*-commutative86.1%
clear-num85.9%
frac-times85.9%
*-un-lft-identity85.9%
frac-add85.8%
associate-/r/85.9%
*-un-lft-identity85.9%
*-commutative85.9%
neg-mul-185.9%
sub-neg85.9%
flip-+99.2%
+-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-/r*99.2%
Simplified99.2%
div-inv99.2%
metadata-eval99.2%
*-commutative99.2%
associate-*l/99.2%
associate-/r*99.8%
Applied egg-rr99.8%
Taylor expanded in a around 0 63.5%
Final simplification71.6%
(FPCore (a b) :precision binary64 (if (<= a -1.2e-45) (* (/ 0.5 a) (/ PI (* a b))) (* PI (/ 0.5 (* b (* a b))))))
double code(double a, double b) {
double tmp;
if (a <= -1.2e-45) {
tmp = (0.5 / a) * (((double) M_PI) / (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.2e-45) {
tmp = (0.5 / a) * (Math.PI / (a * b));
} else {
tmp = Math.PI * (0.5 / (b * (a * b)));
}
return tmp;
}
def code(a, b): tmp = 0 if a <= -1.2e-45: tmp = (0.5 / a) * (math.pi / (a * b)) else: tmp = math.pi * (0.5 / (b * (a * b))) return tmp
function code(a, b) tmp = 0.0 if (a <= -1.2e-45) tmp = Float64(Float64(0.5 / a) * Float64(pi / Float64(a * b))); else tmp = Float64(pi * Float64(0.5 / Float64(b * Float64(a * b)))); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (a <= -1.2e-45) tmp = (0.5 / a) * (pi / (a * b)); else tmp = pi * (0.5 / (b * (a * b))); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[a, -1.2e-45], N[(N[(0.5 / a), $MachinePrecision] * N[(Pi / N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(Pi * N[(0.5 / N[(b * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.2 \cdot 10^{-45}:\\
\;\;\;\;\frac{0.5}{a} \cdot \frac{\pi}{a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;\pi \cdot \frac{0.5}{b \cdot \left(a \cdot b\right)}\\
\end{array}
\end{array}
if a < -1.19999999999999995e-45Initial program 80.0%
associate-*l*80.0%
*-rgt-identity80.0%
associate-/l*80.0%
metadata-eval80.0%
associate-*l/80.1%
*-lft-identity80.1%
sub-neg80.1%
distribute-neg-frac80.1%
metadata-eval80.1%
Simplified80.1%
metadata-eval80.1%
div-inv80.1%
*-commutative80.1%
clear-num80.0%
frac-times80.0%
*-un-lft-identity80.0%
frac-add80.0%
associate-/r/80.1%
*-un-lft-identity80.1%
*-commutative80.1%
neg-mul-180.1%
sub-neg80.1%
flip-+98.9%
+-commutative98.9%
Applied egg-rr98.9%
*-commutative98.9%
associate-/r*98.9%
Simplified98.9%
associate-/l/98.9%
div-inv99.0%
metadata-eval99.0%
div-inv99.0%
clear-num99.0%
associate-/r*99.5%
Applied egg-rr99.5%
associate-*r/99.6%
associate-*l/99.7%
*-commutative99.7%
associate-/r*99.7%
Simplified99.7%
Taylor expanded in a around inf 95.0%
Taylor expanded in a around 0 95.0%
if -1.19999999999999995e-45 < a Initial program 86.0%
*-commutative86.0%
associate-*r*86.1%
associate-*r/86.1%
associate-*r*86.1%
*-rgt-identity86.1%
sub-neg86.1%
distribute-neg-frac86.1%
metadata-eval86.1%
Simplified86.1%
associate-*l/86.1%
div-inv86.1%
metadata-eval86.1%
*-commutative86.1%
associate-*r*86.1%
*-commutative86.1%
Applied egg-rr99.2%
Taylor expanded in a around 0 62.9%
Final simplification71.2%
(FPCore (a b) :precision binary64 (if (<= a -5.1e-46) (* (/ 0.5 a) (/ PI (* a b))) (* (/ (/ PI a) b) (/ 0.5 b))))
double code(double a, double b) {
double tmp;
if (a <= -5.1e-46) {
tmp = (0.5 / a) * (((double) M_PI) / (a * b));
} else {
tmp = ((((double) M_PI) / a) / b) * (0.5 / b);
}
return tmp;
}
public static double code(double a, double b) {
double tmp;
if (a <= -5.1e-46) {
tmp = (0.5 / a) * (Math.PI / (a * b));
} else {
tmp = ((Math.PI / a) / b) * (0.5 / b);
}
return tmp;
}
def code(a, b): tmp = 0 if a <= -5.1e-46: tmp = (0.5 / a) * (math.pi / (a * b)) else: tmp = ((math.pi / a) / b) * (0.5 / b) return tmp
function code(a, b) tmp = 0.0 if (a <= -5.1e-46) tmp = Float64(Float64(0.5 / a) * Float64(pi / Float64(a * b))); else tmp = Float64(Float64(Float64(pi / a) / b) * Float64(0.5 / b)); end return tmp end
function tmp_2 = code(a, b) tmp = 0.0; if (a <= -5.1e-46) tmp = (0.5 / a) * (pi / (a * b)); else tmp = ((pi / a) / b) * (0.5 / b); end tmp_2 = tmp; end
code[a_, b_] := If[LessEqual[a, -5.1e-46], N[(N[(0.5 / a), $MachinePrecision] * N[(Pi / N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(Pi / a), $MachinePrecision] / b), $MachinePrecision] * N[(0.5 / b), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -5.1 \cdot 10^{-46}:\\
\;\;\;\;\frac{0.5}{a} \cdot \frac{\pi}{a \cdot b}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\pi}{a}}{b} \cdot \frac{0.5}{b}\\
\end{array}
\end{array}
if a < -5.0999999999999997e-46Initial program 80.0%
associate-*l*80.0%
*-rgt-identity80.0%
associate-/l*80.0%
metadata-eval80.0%
associate-*l/80.1%
*-lft-identity80.1%
sub-neg80.1%
distribute-neg-frac80.1%
metadata-eval80.1%
Simplified80.1%
metadata-eval80.1%
div-inv80.1%
*-commutative80.1%
clear-num80.0%
frac-times80.0%
*-un-lft-identity80.0%
frac-add80.0%
associate-/r/80.1%
*-un-lft-identity80.1%
*-commutative80.1%
neg-mul-180.1%
sub-neg80.1%
flip-+98.9%
+-commutative98.9%
Applied egg-rr98.9%
*-commutative98.9%
associate-/r*98.9%
Simplified98.9%
associate-/l/98.9%
div-inv99.0%
metadata-eval99.0%
div-inv99.0%
clear-num99.0%
associate-/r*99.5%
Applied egg-rr99.5%
associate-*r/99.6%
associate-*l/99.7%
*-commutative99.7%
associate-/r*99.7%
Simplified99.7%
Taylor expanded in a around inf 95.0%
Taylor expanded in a around 0 95.0%
if -5.0999999999999997e-46 < a Initial program 86.0%
associate-*l*86.0%
*-rgt-identity86.0%
associate-/l*86.0%
metadata-eval86.0%
associate-*l/86.1%
*-lft-identity86.1%
sub-neg86.1%
distribute-neg-frac86.1%
metadata-eval86.1%
Simplified86.1%
metadata-eval86.1%
div-inv86.1%
*-commutative86.1%
clear-num85.9%
frac-times85.9%
*-un-lft-identity85.9%
frac-add85.8%
associate-/r/85.9%
*-un-lft-identity85.9%
*-commutative85.9%
neg-mul-185.9%
sub-neg85.9%
flip-+99.2%
+-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-/r*99.2%
Simplified99.2%
associate-/l/99.2%
div-inv99.2%
metadata-eval99.2%
div-inv99.2%
clear-num99.2%
associate-/r*99.8%
Applied egg-rr99.8%
associate-*r/99.8%
associate-*l/99.8%
*-commutative99.8%
associate-/r*99.8%
Simplified99.8%
Taylor expanded in a around 0 63.5%
Final simplification71.6%
(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(Float64(0.5 / 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[(N[(0.5 / N[(a + b), $MachinePrecision]), $MachinePrecision] / N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\pi \cdot \frac{\frac{0.5}{a + b}}{a \cdot b}
\end{array}
Initial program 84.5%
associate-*l*84.5%
*-rgt-identity84.5%
associate-/l*84.5%
metadata-eval84.5%
associate-*l/84.5%
*-lft-identity84.5%
sub-neg84.5%
distribute-neg-frac84.5%
metadata-eval84.5%
Simplified84.5%
metadata-eval84.5%
div-inv84.5%
*-commutative84.5%
clear-num84.4%
frac-times84.4%
*-un-lft-identity84.4%
frac-add84.3%
associate-/r/84.4%
*-un-lft-identity84.4%
*-commutative84.4%
neg-mul-184.4%
sub-neg84.4%
flip-+99.2%
+-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-/r*99.2%
Simplified99.2%
div-inv99.2%
metadata-eval99.2%
*-commutative99.2%
associate-*l/99.2%
associate-/r*99.7%
Applied egg-rr99.7%
Final simplification99.7%
(FPCore (a b) :precision binary64 (* (/ 0.5 (+ a b)) (/ (/ PI a) b)))
double code(double a, double b) {
return (0.5 / (a + b)) * ((((double) M_PI) / a) / b);
}
public static double code(double a, double b) {
return (0.5 / (a + b)) * ((Math.PI / a) / b);
}
def code(a, b): return (0.5 / (a + b)) * ((math.pi / a) / b)
function code(a, b) return Float64(Float64(0.5 / Float64(a + b)) * Float64(Float64(pi / a) / b)) end
function tmp = code(a, b) tmp = (0.5 / (a + b)) * ((pi / a) / b); end
code[a_, b_] := N[(N[(0.5 / N[(a + b), $MachinePrecision]), $MachinePrecision] * N[(N[(Pi / a), $MachinePrecision] / b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.5}{a + b} \cdot \frac{\frac{\pi}{a}}{b}
\end{array}
Initial program 84.5%
associate-*l*84.5%
*-rgt-identity84.5%
associate-/l*84.5%
metadata-eval84.5%
associate-*l/84.5%
*-lft-identity84.5%
sub-neg84.5%
distribute-neg-frac84.5%
metadata-eval84.5%
Simplified84.5%
metadata-eval84.5%
div-inv84.5%
*-commutative84.5%
clear-num84.4%
frac-times84.4%
*-un-lft-identity84.4%
frac-add84.3%
associate-/r/84.4%
*-un-lft-identity84.4%
*-commutative84.4%
neg-mul-184.4%
sub-neg84.4%
flip-+99.2%
+-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-/r*99.2%
Simplified99.2%
associate-/l/99.2%
div-inv99.2%
metadata-eval99.2%
div-inv99.2%
clear-num99.2%
associate-/r*99.7%
Applied egg-rr99.7%
associate-*r/99.7%
associate-*l/99.7%
*-commutative99.7%
associate-/r*99.7%
Simplified99.7%
(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 84.5%
*-commutative84.5%
associate-*r*84.5%
associate-*r/84.5%
associate-*r*84.5%
*-rgt-identity84.5%
sub-neg84.5%
distribute-neg-frac84.5%
metadata-eval84.5%
Simplified84.5%
associate-*l/84.5%
div-inv84.5%
metadata-eval84.5%
*-commutative84.5%
associate-*r*84.5%
*-commutative84.5%
Applied egg-rr99.2%
Final simplification99.2%
(FPCore (a b) :precision binary64 (* (/ 0.5 a) (/ PI (* a b))))
double code(double a, double b) {
return (0.5 / a) * (((double) M_PI) / (a * b));
}
public static double code(double a, double b) {
return (0.5 / a) * (Math.PI / (a * b));
}
def code(a, b): return (0.5 / a) * (math.pi / (a * b))
function code(a, b) return Float64(Float64(0.5 / a) * Float64(pi / Float64(a * b))) end
function tmp = code(a, b) tmp = (0.5 / a) * (pi / (a * b)); end
code[a_, b_] := N[(N[(0.5 / a), $MachinePrecision] * N[(Pi / N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.5}{a} \cdot \frac{\pi}{a \cdot b}
\end{array}
Initial program 84.5%
associate-*l*84.5%
*-rgt-identity84.5%
associate-/l*84.5%
metadata-eval84.5%
associate-*l/84.5%
*-lft-identity84.5%
sub-neg84.5%
distribute-neg-frac84.5%
metadata-eval84.5%
Simplified84.5%
metadata-eval84.5%
div-inv84.5%
*-commutative84.5%
clear-num84.4%
frac-times84.4%
*-un-lft-identity84.4%
frac-add84.3%
associate-/r/84.4%
*-un-lft-identity84.4%
*-commutative84.4%
neg-mul-184.4%
sub-neg84.4%
flip-+99.2%
+-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-/r*99.2%
Simplified99.2%
associate-/l/99.2%
div-inv99.2%
metadata-eval99.2%
div-inv99.2%
clear-num99.2%
associate-/r*99.7%
Applied egg-rr99.7%
associate-*r/99.7%
associate-*l/99.7%
*-commutative99.7%
associate-/r*99.7%
Simplified99.7%
Taylor expanded in a around inf 72.2%
Taylor expanded in a around 0 72.2%
herbie shell --seed 2024146
(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))))