
(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 79.3%
sqr-neg79.3%
fabs-div79.3%
sqr-neg79.3%
fabs-sub79.3%
sqr-neg79.3%
distribute-rgt-neg-out79.3%
fabs-neg79.3%
fabs-div79.3%
cancel-sign-sub-inv79.3%
+-commutative79.3%
sqr-neg79.3%
cancel-sign-sub-inv79.3%
div-sub79.3%
Simplified79.9%
associate-*r/79.3%
frac-times100.0%
clear-num100.0%
un-div-inv100.0%
Applied egg-rr100.0%
Taylor expanded in b around 0 79.3%
fabs-sub79.3%
unpow279.3%
unpow279.3%
times-frac100.0%
unpow2100.0%
fabs-sub100.0%
rem-square-sqrt100.0%
fabs-sqr100.0%
rem-square-sqrt100.0%
Simplified100.0%
unpow2100.0%
clear-num100.0%
un-div-inv100.0%
Applied egg-rr100.0%
Final simplification100.0%
(FPCore (a b) :precision binary64 (+ 1.0 (/ -0.5 (pow (/ a b) 2.0))))
double code(double a, double b) {
return 1.0 + (-0.5 / pow((a / b), 2.0));
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = 1.0d0 + ((-0.5d0) / ((a / b) ** 2.0d0))
end function
public static double code(double a, double b) {
return 1.0 + (-0.5 / Math.pow((a / b), 2.0));
}
def code(a, b): return 1.0 + (-0.5 / math.pow((a / b), 2.0))
function code(a, b) return Float64(1.0 + Float64(-0.5 / (Float64(a / b) ^ 2.0))) end
function tmp = code(a, b) tmp = 1.0 + (-0.5 / ((a / b) ^ 2.0)); end
code[a_, b_] := N[(1.0 + N[(-0.5 / N[Power[N[(a / b), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 + \frac{-0.5}{{\left(\frac{a}{b}\right)}^{2}}
\end{array}
Initial program 79.3%
sqr-neg79.3%
fabs-div79.3%
sqr-neg79.3%
fabs-sub79.3%
sqr-neg79.3%
distribute-rgt-neg-out79.3%
fabs-neg79.3%
fabs-div79.3%
cancel-sign-sub-inv79.3%
+-commutative79.3%
sqr-neg79.3%
cancel-sign-sub-inv79.3%
div-sub79.3%
Simplified79.9%
associate-*r/79.3%
frac-times100.0%
clear-num100.0%
un-div-inv100.0%
Applied egg-rr100.0%
Taylor expanded in b around 0 79.3%
fabs-sub79.3%
unpow279.3%
unpow279.3%
times-frac100.0%
unpow2100.0%
fabs-sub100.0%
rem-square-sqrt100.0%
fabs-sqr100.0%
rem-square-sqrt100.0%
Simplified100.0%
Taylor expanded in b around 0 78.5%
+-commutative78.5%
fma-define78.5%
unpow278.5%
unpow278.5%
times-frac98.9%
unpow298.9%
Simplified98.9%
fma-undefine98.9%
unpow298.9%
clear-num98.9%
associate-/r/98.9%
un-div-inv98.9%
div-inv98.9%
clear-num98.9%
pow298.9%
Applied egg-rr98.9%
Final simplification98.9%
(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 79.3%
sqr-neg79.3%
fabs-div79.3%
sqr-neg79.3%
fabs-sub79.3%
sqr-neg79.3%
distribute-rgt-neg-out79.3%
fabs-neg79.3%
fabs-div79.3%
cancel-sign-sub-inv79.3%
+-commutative79.3%
sqr-neg79.3%
cancel-sign-sub-inv79.3%
div-sub79.3%
Simplified79.9%
associate-*r/79.3%
frac-times100.0%
clear-num100.0%
un-div-inv100.0%
Applied egg-rr100.0%
Taylor expanded in b around 0 79.3%
fabs-sub79.3%
unpow279.3%
unpow279.3%
times-frac100.0%
unpow2100.0%
fabs-sub100.0%
rem-square-sqrt100.0%
fabs-sqr100.0%
rem-square-sqrt100.0%
Simplified100.0%
unpow2100.0%
clear-num100.0%
un-div-inv100.0%
Applied egg-rr100.0%
Taylor expanded in b around 0 98.0%
Final simplification98.0%
herbie shell --seed 2024130
(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)))))