
(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 (/ (fabs (- (/ b (/ a b)) a)) a)))
double code(double a, double b) {
return sqrt((fabs(((b / (a / b)) - a)) / a));
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = sqrt((abs(((b / (a / b)) - a)) / a))
end function
public static double code(double a, double b) {
return Math.sqrt((Math.abs(((b / (a / b)) - a)) / a));
}
def code(a, b): return math.sqrt((math.fabs(((b / (a / b)) - a)) / a))
function code(a, b) return sqrt(Float64(abs(Float64(Float64(b / Float64(a / b)) - a)) / a)) end
function tmp = code(a, b) tmp = sqrt((abs(((b / (a / b)) - a)) / a)); end
code[a_, b_] := N[Sqrt[N[(N[Abs[N[(N[(b / N[(a / b), $MachinePrecision]), $MachinePrecision] - a), $MachinePrecision]], $MachinePrecision] / a), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\sqrt{\frac{\left|\frac{b}{\frac{a}{b}} - a\right|}{a}}
\end{array}
Initial program 74.2%
sqrt-lowering-sqrt.f64N/A
div-subN/A
fabs-subN/A
fabs-lowering-fabs.f64N/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64100.0%
Simplified100.0%
neg-fabsN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
metadata-evalN/A
associate-*r/N/A
associate-/r*N/A
sub-negN/A
metadata-evalN/A
*-inversesN/A
div-subN/A
associate-/r*N/A
fabs-divN/A
rem-sqrt-squareN/A
pow2N/A
sqrt-pow1N/A
metadata-evalN/A
unpow1N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in b around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6498.3%
Simplified98.3%
fabs-subN/A
fabs-lowering-fabs.f64N/A
--lowering--.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64100.0%
Applied egg-rr100.0%
(FPCore (a b) :precision binary64 (sqrt (fabs (+ (/ (* b (/ b a)) a) -1.0))))
double code(double a, double b) {
return sqrt(fabs((((b * (b / a)) / a) + -1.0)));
}
real(8) function code(a, b)
real(8), intent (in) :: a
real(8), intent (in) :: b
code = sqrt(abs((((b * (b / a)) / a) + (-1.0d0))))
end function
public static double code(double a, double b) {
return Math.sqrt(Math.abs((((b * (b / a)) / a) + -1.0)));
}
def code(a, b): return math.sqrt(math.fabs((((b * (b / a)) / a) + -1.0)))
function code(a, b) return sqrt(abs(Float64(Float64(Float64(b * Float64(b / a)) / a) + -1.0))) end
function tmp = code(a, b) tmp = sqrt(abs((((b * (b / a)) / a) + -1.0))); end
code[a_, b_] := N[Sqrt[N[Abs[N[(N[(N[(b * N[(b / a), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision] + -1.0), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\sqrt{\left|\frac{b \cdot \frac{b}{a}}{a} + -1\right|}
\end{array}
Initial program 74.2%
sqrt-lowering-sqrt.f64N/A
div-subN/A
fabs-subN/A
fabs-lowering-fabs.f64N/A
sub-negN/A
*-inversesN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-/r*N/A
/-lowering-/.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64100.0%
Simplified100.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.2%
Taylor expanded in a around inf
Simplified97.7%
metadata-evalN/A
metadata-eval97.7%
Applied egg-rr97.7%
herbie shell --seed 2024161
(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)))))