
(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 (+ (+ 1.0 (/ (/ b a) (/ a b))) -1.0))))
double code(double a, double b) {
return sqrt((1.0 - ((1.0 + ((b / a) / (a / b))) + -1.0)));
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = sqrt((1.0d0 - ((1.0d0 + ((b / a) / (a / b))) + (-1.0d0))))
end function
public static double code(double a, double b) {
return Math.sqrt((1.0 - ((1.0 + ((b / a) / (a / b))) + -1.0)));
}
def code(a, b): return math.sqrt((1.0 - ((1.0 + ((b / a) / (a / b))) + -1.0)))
function code(a, b) return sqrt(Float64(1.0 - Float64(Float64(1.0 + Float64(Float64(b / a) / Float64(a / b))) + -1.0))) end
function tmp = code(a, b) tmp = sqrt((1.0 - ((1.0 + ((b / a) / (a / b))) + -1.0))); end
code[a_, b_] := N[Sqrt[N[(1.0 - N[(N[(1.0 + N[(N[(b / a), $MachinePrecision] / N[(a / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\sqrt{1 - \left(\left(1 + \frac{\frac{b}{a}}{\frac{a}{b}}\right) + -1\right)}
\end{array}
Initial program 77.3%
sqr-neg77.3%
sqr-neg77.3%
div-sub77.3%
*-inverses77.3%
times-frac100.0%
Simplified100.0%
Taylor expanded in b around 0 77.3%
fabs-sub77.3%
*-inverses77.3%
fabs-sub77.3%
div-sub77.3%
fabs-div77.3%
rem-square-sqrt77.3%
fabs-sqr77.3%
unpow277.3%
fabs-sqr77.3%
unpow277.3%
rem-square-sqrt77.3%
div-sub77.3%
*-inverses77.3%
unpow277.3%
unpow277.3%
times-frac100.0%
unpow2100.0%
Simplified100.0%
expm1-log1p-u100.0%
expm1-udef100.0%
log1p-udef100.0%
add-exp-log100.0%
+-commutative100.0%
Applied egg-rr100.0%
unpow2100.0%
clear-num100.0%
un-div-inv100.0%
Applied egg-rr100.0%
Final simplification100.0%
(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 77.3%
sqr-neg77.3%
sqr-neg77.3%
div-sub77.3%
*-inverses77.3%
times-frac100.0%
Simplified100.0%
Taylor expanded in b around 0 77.3%
fabs-sub77.3%
*-inverses77.3%
fabs-sub77.3%
div-sub77.3%
fabs-div77.3%
rem-square-sqrt77.3%
fabs-sqr77.3%
unpow277.3%
fabs-sqr77.3%
unpow277.3%
rem-square-sqrt77.3%
div-sub77.3%
*-inverses77.3%
unpow277.3%
unpow277.3%
times-frac100.0%
unpow2100.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)
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 77.3%
sqr-neg77.3%
sqr-neg77.3%
div-sub77.3%
*-inverses77.3%
times-frac100.0%
Simplified100.0%
Taylor expanded in b around 0 77.3%
fabs-sub77.3%
*-inverses77.3%
fabs-sub77.3%
div-sub77.3%
fabs-div77.3%
rem-square-sqrt77.3%
fabs-sqr77.3%
unpow277.3%
fabs-sqr77.3%
unpow277.3%
rem-square-sqrt77.3%
div-sub77.3%
*-inverses77.3%
unpow277.3%
unpow277.3%
times-frac100.0%
unpow2100.0%
Simplified100.0%
Taylor expanded in b around 0 96.4%
Final simplification96.4%
herbie shell --seed 2024013
(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)))))