
(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 (- 1.0 (/ (/ b a) (/ a b)))))
double code(double a, double b) {
return sqrt((1.0 - ((b / a) / (a / b))));
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = sqrt((1.0d0 - ((b / a) / (a / b))))
end function
public static double code(double a, double b) {
return Math.sqrt((1.0 - ((b / a) / (a / b))));
}
def code(a, b): return math.sqrt((1.0 - ((b / a) / (a / b))))
function code(a, b) return sqrt(Float64(1.0 - Float64(Float64(b / a) / Float64(a / b)))) end
function tmp = code(a, b) tmp = sqrt((1.0 - ((b / a) / (a / b)))); end
code[a_, b_] := N[Sqrt[N[(1.0 - N[(N[(b / a), $MachinePrecision] / N[(a / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\sqrt{1 - \frac{\frac{b}{a}}{\frac{a}{b}}}
\end{array}
Initial program 80.1%
sqr-neg80.1%
associate-/r*79.6%
sqr-neg79.6%
associate-/r*80.1%
div-sub80.1%
fabs-sub80.1%
times-frac80.1%
*-inverses100.0%
difference-of-sqr-1100.0%
difference-of-sqr--1100.0%
fma-define100.0%
Simplified100.0%
add-sqr-sqrt100.0%
associate-/l*100.0%
add-sqr-sqrt100.0%
sqrt-prod80.7%
sqrt-div80.7%
sqrt-prod80.7%
add-sqr-sqrt80.7%
associate-/l*80.7%
add-sqr-sqrt80.7%
sqrt-prod80.7%
sqrt-div80.7%
sqrt-prod80.7%
metadata-eval80.7%
fma-neg80.7%
add-sqr-sqrt80.7%
fabs-sub80.7%
add-sqr-sqrt80.1%
Applied egg-rr100.0%
unpow2100.0%
clear-num100.0%
un-div-inv100.0%
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 80.1%
sqr-neg80.1%
associate-/r*79.6%
sqr-neg79.6%
associate-/r*80.1%
div-sub80.1%
fabs-sub80.1%
times-frac80.1%
*-inverses100.0%
difference-of-sqr-1100.0%
difference-of-sqr--1100.0%
fma-define100.0%
Simplified100.0%
add-sqr-sqrt100.0%
associate-/l*100.0%
add-sqr-sqrt100.0%
sqrt-prod80.7%
sqrt-div80.7%
sqrt-prod80.7%
add-sqr-sqrt80.7%
associate-/l*80.7%
add-sqr-sqrt80.7%
sqrt-prod80.7%
sqrt-div80.7%
sqrt-prod80.7%
metadata-eval80.7%
fma-neg80.7%
add-sqr-sqrt80.7%
fabs-sub80.7%
add-sqr-sqrt80.1%
Applied egg-rr100.0%
unpow2100.0%
clear-num100.0%
un-div-inv100.0%
Applied egg-rr100.0%
Taylor expanded in b around 0 98.6%
(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 80.1%
sqr-neg80.1%
associate-/r*79.6%
sqr-neg79.6%
associate-/r*80.1%
div-sub80.1%
fabs-sub80.1%
times-frac80.1%
*-inverses100.0%
difference-of-sqr-1100.0%
difference-of-sqr--1100.0%
fma-define100.0%
Simplified100.0%
add-sqr-sqrt100.0%
associate-/l*100.0%
add-sqr-sqrt100.0%
sqrt-prod80.7%
sqrt-div80.7%
sqrt-prod80.7%
add-sqr-sqrt80.7%
associate-/l*80.7%
add-sqr-sqrt80.7%
sqrt-prod80.7%
sqrt-div80.7%
sqrt-prod80.7%
metadata-eval80.7%
fma-neg80.7%
add-sqr-sqrt80.7%
fabs-sub80.7%
*-inverses80.1%
associate-*r/80.1%
Applied egg-rr100.0%
Taylor expanded in a around -inf 0.0%
unpow20.0%
rem-square-sqrt1.6%
Simplified1.6%
herbie shell --seed 2024135
(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)))))