
(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 4 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 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(Float64(Float64(pi * 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[(N[(N[(Pi * 0.5), $MachinePrecision] / N[(a + b), $MachinePrecision]), $MachinePrecision] / N[(a * b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{\pi \cdot 0.5}{a + b}}{a \cdot b}
\end{array}
Initial program 78.0%
un-div-inv78.0%
difference-of-squares90.1%
associate-/r*90.3%
div-inv90.3%
metadata-eval90.3%
Applied egg-rr90.3%
associate-*l/99.6%
associate-/l*99.6%
Applied egg-rr99.6%
associate-/l*99.5%
+-commutative99.5%
sub-neg99.5%
distribute-neg-frac99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in a around 0 99.5%
un-div-inv99.6%
associate-*r/99.6%
Applied egg-rr99.6%
(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 \cdot b} \cdot \frac{0.5}{a + b}
\end{array}
Initial program 78.0%
un-div-inv78.0%
difference-of-squares90.1%
associate-/r*90.3%
div-inv90.3%
metadata-eval90.3%
Applied egg-rr90.3%
associate-*l/99.6%
associate-/l*99.6%
Applied egg-rr99.6%
associate-/l*99.5%
+-commutative99.5%
sub-neg99.5%
distribute-neg-frac99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in a around 0 99.5%
un-div-inv99.6%
associate-*r/99.6%
Applied egg-rr99.6%
associate-/l/99.5%
times-frac99.6%
Simplified99.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(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 78.0%
un-div-inv78.0%
difference-of-squares90.1%
associate-/r*90.3%
div-inv90.3%
metadata-eval90.3%
Applied egg-rr90.3%
associate-*l/99.6%
associate-/l*99.6%
Applied egg-rr99.6%
associate-/l*99.5%
+-commutative99.5%
sub-neg99.5%
distribute-neg-frac99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in a around 0 99.5%
un-div-inv99.6%
associate-*r/99.6%
Applied egg-rr99.6%
associate-/l/99.5%
associate-/l*99.5%
Simplified99.5%
Final simplification99.5%
(FPCore (a b) :precision binary64 (* (/ PI b) (/ 0.5 (* a b))))
double code(double a, double b) {
return (((double) M_PI) / b) * (0.5 / (a * b));
}
public static double code(double a, double b) {
return (Math.PI / b) * (0.5 / (a * b));
}
def code(a, b): return (math.pi / b) * (0.5 / (a * b))
function code(a, b) return Float64(Float64(pi / b) * Float64(0.5 / Float64(a * b))) end
function tmp = code(a, b) tmp = (pi / b) * (0.5 / (a * b)); end
code[a_, b_] := N[(N[(Pi / b), $MachinePrecision] * N[(0.5 / N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\pi}{b} \cdot \frac{0.5}{a \cdot b}
\end{array}
Initial program 78.0%
associate-*l*78.0%
*-rgt-identity78.0%
associate-/l*78.0%
metadata-eval78.0%
associate-*l/78.0%
*-lft-identity78.0%
sub-neg78.0%
distribute-neg-frac78.0%
metadata-eval78.0%
Simplified78.0%
metadata-eval78.0%
div-inv78.0%
associate-*r/78.1%
*-commutative78.1%
difference-of-squares90.2%
associate-/r*99.6%
Applied egg-rr73.2%
Taylor expanded in a around 0 73.3%
associate-*r/73.3%
*-commutative73.3%
*-commutative73.3%
times-frac73.2%
Simplified73.2%
*-un-lft-identity73.2%
frac-times73.3%
Applied egg-rr73.3%
*-lft-identity73.3%
times-frac73.2%
associate-*r/73.3%
associate-/l/73.3%
Simplified73.3%
Taylor expanded in b around inf 71.9%
herbie shell --seed 2024170
(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))))