
(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 74.6%
sqr-neg74.6%
fabs-div74.6%
sqr-neg74.6%
fabs-sub74.6%
sqr-neg74.6%
distribute-rgt-neg-out74.6%
fabs-neg74.6%
fabs-div74.6%
cancel-sign-sub-inv74.6%
+-commutative74.6%
sqr-neg74.6%
cancel-sign-sub-inv74.6%
Simplified75.4%
*-inverses74.6%
associate-*r/74.6%
div-sub74.6%
difference-of-squares74.6%
times-frac100.0%
+-commutative100.0%
Applied egg-rr100.0%
Taylor expanded in b around 0 74.6%
fabs-div74.6%
difference-of-squares74.6%
unpow274.6%
unpow274.6%
*-rgt-identity74.6%
fabs-div74.6%
*-rgt-identity74.6%
div-sub74.6%
fabs-sub74.6%
unpow274.6%
unpow274.6%
times-frac74.6%
unpow274.6%
*-inverses100.0%
fabs-sub100.0%
rem-square-sqrt100.0%
fabs-sqr100.0%
Simplified100.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 74.6%
sqr-neg74.6%
fabs-div74.6%
sqr-neg74.6%
fabs-sub74.6%
sqr-neg74.6%
distribute-rgt-neg-out74.6%
fabs-neg74.6%
fabs-div74.6%
cancel-sign-sub-inv74.6%
+-commutative74.6%
sqr-neg74.6%
cancel-sign-sub-inv74.6%
Simplified75.4%
*-inverses74.6%
associate-*r/74.6%
div-sub74.6%
difference-of-squares74.6%
times-frac100.0%
+-commutative100.0%
Applied egg-rr100.0%
add-sqr-sqrt100.0%
fabs-sqr100.0%
add-sqr-sqrt100.0%
associate-*l/100.0%
associate-/l*99.7%
div-sub99.7%
*-inverses99.7%
Applied egg-rr99.7%
Taylor expanded in a around inf 96.5%
associate-*r/96.5%
distribute-rgt1-in96.5%
metadata-eval96.5%
mul0-lft96.5%
metadata-eval96.5%
Simplified96.5%
Taylor expanded in a around 0 96.5%
Final simplification96.5%
(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 74.6%
sqr-neg74.6%
fabs-div74.6%
sqr-neg74.6%
fabs-sub74.6%
sqr-neg74.6%
distribute-rgt-neg-out74.6%
fabs-neg74.6%
fabs-div74.6%
cancel-sign-sub-inv74.6%
+-commutative74.6%
sqr-neg74.6%
cancel-sign-sub-inv74.6%
Simplified75.4%
*-inverses74.6%
associate-*r/74.6%
div-sub74.6%
difference-of-squares74.6%
times-frac100.0%
+-commutative100.0%
Applied egg-rr100.0%
add-sqr-sqrt100.0%
fabs-sqr100.0%
add-sqr-sqrt100.0%
associate-*l/100.0%
associate-/l*99.7%
div-sub99.7%
*-inverses99.7%
Applied egg-rr99.7%
Taylor expanded in a around -inf 0.0%
unpow20.0%
rem-square-sqrt1.6%
Simplified1.6%
herbie shell --seed 2024112
(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)))))