
(FPCore (a b) :precision binary64 (sqrt (fabs (/ (- (* a a) (* b b)) (* a a)))))
double code(double a, double b) {
return sqrt(fabs((((a * a) - (b * b)) / (a * a))));
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = sqrt(abs((((a * a) - (b * b)) / (a * a))))
end function
public static double code(double a, double b) {
return Math.sqrt(Math.abs((((a * a) - (b * b)) / (a * a))));
}
def code(a, b): return math.sqrt(math.fabs((((a * a) - (b * b)) / (a * a))))
function code(a, b) return sqrt(abs(Float64(Float64(Float64(a * a) - Float64(b * b)) / Float64(a * a)))) end
function tmp = code(a, b) tmp = sqrt(abs((((a * a) - (b * b)) / (a * a)))); end
code[a_, b_] := N[Sqrt[N[Abs[N[(N[(N[(a * a), $MachinePrecision] - N[(b * b), $MachinePrecision]), $MachinePrecision] / N[(a * a), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\sqrt{\left|\frac{a \cdot a - b \cdot b}{a \cdot a}\right|}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 3 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a b) :precision binary64 (sqrt (fabs (/ (- (* a a) (* b b)) (* a a)))))
double code(double a, double b) {
return sqrt(fabs((((a * a) - (b * b)) / (a * a))));
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = sqrt(abs((((a * a) - (b * b)) / (a * a))))
end function
public static double code(double a, double b) {
return Math.sqrt(Math.abs((((a * a) - (b * b)) / (a * a))));
}
def code(a, b): return math.sqrt(math.fabs((((a * a) - (b * b)) / (a * a))))
function code(a, b) return sqrt(abs(Float64(Float64(Float64(a * a) - Float64(b * b)) / Float64(a * a)))) end
function tmp = code(a, b) tmp = sqrt(abs((((a * a) - (b * b)) / (a * a)))); end
code[a_, b_] := N[Sqrt[N[Abs[N[(N[(N[(a * a), $MachinePrecision] - N[(b * b), $MachinePrecision]), $MachinePrecision] / N[(a * a), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\sqrt{\left|\frac{a \cdot a - b \cdot b}{a \cdot a}\right|}
\end{array}
(FPCore (a b) :precision binary64 (sqrt (/ (/ (+ b a) a) (/ a (- a b)))))
double code(double a, double b) {
return sqrt((((b + a) / a) / (a / (a - b))));
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = sqrt((((b + a) / a) / (a / (a - b))))
end function
public static double code(double a, double b) {
return Math.sqrt((((b + a) / a) / (a / (a - b))));
}
def code(a, b): return math.sqrt((((b + a) / a) / (a / (a - b))))
function code(a, b) return sqrt(Float64(Float64(Float64(b + a) / a) / Float64(a / Float64(a - b)))) end
function tmp = code(a, b) tmp = sqrt((((b + a) / a) / (a / (a - b)))); end
code[a_, b_] := N[Sqrt[N[(N[(N[(b + a), $MachinePrecision] / a), $MachinePrecision] / N[(a / N[(a - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\sqrt{\frac{\frac{b + a}{a}}{\frac{a}{a - b}}}
\end{array}
Initial program 76.2%
sqr-neg76.2%
sqr-neg76.2%
div-sub76.2%
fabs-sub76.2%
times-frac76.2%
*-inverses100.0%
difference-of-sqr-199.9%
difference-of-sqr--1100.0%
fma-define100.0%
Simplified100.0%
add-sqr-sqrt100.0%
associate-/l*100.0%
add-sqr-sqrt100.0%
sqrt-prod76.9%
sqrt-div76.9%
sqrt-prod76.9%
add-sqr-sqrt76.9%
associate-/l*76.9%
add-sqr-sqrt76.9%
sqrt-prod76.9%
sqrt-div76.9%
sqrt-prod76.9%
metadata-eval76.9%
fmm-def76.9%
add-sqr-sqrt76.9%
fabs-sub76.9%
*-inverses76.2%
associate-*r/76.2%
Applied egg-rr100.0%
clear-num100.0%
un-div-inv100.0%
Applied egg-rr100.0%
(FPCore (a b) :precision binary64 (sqrt (* (/ (+ b a) a) (/ (- a b) a))))
double code(double a, double b) {
return sqrt((((b + a) / a) * ((a - b) / a)));
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = sqrt((((b + a) / a) * ((a - b) / a)))
end function
public static double code(double a, double b) {
return Math.sqrt((((b + a) / a) * ((a - b) / a)));
}
def code(a, b): return math.sqrt((((b + a) / a) * ((a - b) / a)))
function code(a, b) return sqrt(Float64(Float64(Float64(b + a) / a) * Float64(Float64(a - b) / a))) end
function tmp = code(a, b) tmp = sqrt((((b + a) / a) * ((a - b) / a))); end
code[a_, b_] := N[Sqrt[N[(N[(N[(b + a), $MachinePrecision] / a), $MachinePrecision] * N[(N[(a - b), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\sqrt{\frac{b + a}{a} \cdot \frac{a - b}{a}}
\end{array}
Initial program 76.2%
sqr-neg76.2%
sqr-neg76.2%
div-sub76.2%
fabs-sub76.2%
times-frac76.2%
*-inverses100.0%
difference-of-sqr-199.9%
difference-of-sqr--1100.0%
fma-define100.0%
Simplified100.0%
add-sqr-sqrt100.0%
associate-/l*100.0%
add-sqr-sqrt100.0%
sqrt-prod76.9%
sqrt-div76.9%
sqrt-prod76.9%
add-sqr-sqrt76.9%
associate-/l*76.9%
add-sqr-sqrt76.9%
sqrt-prod76.9%
sqrt-div76.9%
sqrt-prod76.9%
metadata-eval76.9%
fmm-def76.9%
add-sqr-sqrt76.9%
fabs-sub76.9%
*-inverses76.2%
associate-*r/76.2%
Applied egg-rr100.0%
(FPCore (a b) :precision binary64 1.0)
double code(double a, double b) {
return 1.0;
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = 1.0d0
end function
public static double code(double a, double b) {
return 1.0;
}
def code(a, b): return 1.0
function code(a, b) return 1.0 end
function tmp = code(a, b) tmp = 1.0; end
code[a_, b_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 76.2%
sqr-neg76.2%
fabs-div76.2%
sqr-neg76.2%
fabs-sub76.2%
sqr-neg76.2%
distribute-rgt-neg-out76.2%
fabs-neg76.2%
fabs-div76.2%
cancel-sign-sub-inv76.2%
+-commutative76.2%
sqr-neg76.2%
cancel-sign-sub-inv76.2%
div-sub76.2%
Simplified76.9%
pow1/276.9%
pow-to-exp76.9%
add-sqr-sqrt76.2%
fabs-sqr76.2%
add-sqr-sqrt76.2%
sub-neg76.2%
log1p-define76.2%
associate-*r/76.2%
frac-times100.0%
pow2100.0%
Applied egg-rr100.0%
Taylor expanded in b around 0 98.4%
herbie shell --seed 2024156
(FPCore (a b)
:name "Eccentricity of an ellipse"
:precision binary64
:pre (and (and (<= 0.0 b) (<= b a)) (<= a 1.0))
(sqrt (fabs (/ (- (* a a) (* b b)) (* a a)))))