
(FPCore (x) :precision binary64 (sqrt (* 2.0 (pow x 2.0))))
double code(double x) {
return sqrt((2.0 * pow(x, 2.0)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = sqrt((2.0d0 * (x ** 2.0d0)))
end function
public static double code(double x) {
return Math.sqrt((2.0 * Math.pow(x, 2.0)));
}
def code(x): return math.sqrt((2.0 * math.pow(x, 2.0)))
function code(x) return sqrt(Float64(2.0 * (x ^ 2.0))) end
function tmp = code(x) tmp = sqrt((2.0 * (x ^ 2.0))); end
code[x_] := N[Sqrt[N[(2.0 * N[Power[x, 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\sqrt{2 \cdot {x}^{2}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 4 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (sqrt (* 2.0 (pow x 2.0))))
double code(double x) {
return sqrt((2.0 * pow(x, 2.0)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = sqrt((2.0d0 * (x ** 2.0d0)))
end function
public static double code(double x) {
return Math.sqrt((2.0 * Math.pow(x, 2.0)));
}
def code(x): return math.sqrt((2.0 * math.pow(x, 2.0)))
function code(x) return sqrt(Float64(2.0 * (x ^ 2.0))) end
function tmp = code(x) tmp = sqrt((2.0 * (x ^ 2.0))); end
code[x_] := N[Sqrt[N[(2.0 * N[Power[x, 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\sqrt{2 \cdot {x}^{2}}
\end{array}
x_m = (fabs.f64 x) (FPCore (x_m) :precision binary64 (* (pow 2.0 0.375) (* (pow 2.0 0.125) x_m)))
x_m = fabs(x);
double code(double x_m) {
return pow(2.0, 0.375) * (pow(2.0, 0.125) * x_m);
}
x_m = abs(x)
real(8) function code(x_m)
real(8), intent (in) :: x_m
code = (2.0d0 ** 0.375d0) * ((2.0d0 ** 0.125d0) * x_m)
end function
x_m = Math.abs(x);
public static double code(double x_m) {
return Math.pow(2.0, 0.375) * (Math.pow(2.0, 0.125) * x_m);
}
x_m = math.fabs(x) def code(x_m): return math.pow(2.0, 0.375) * (math.pow(2.0, 0.125) * x_m)
x_m = abs(x) function code(x_m) return Float64((2.0 ^ 0.375) * Float64((2.0 ^ 0.125) * x_m)) end
x_m = abs(x); function tmp = code(x_m) tmp = (2.0 ^ 0.375) * ((2.0 ^ 0.125) * x_m); end
x_m = N[Abs[x], $MachinePrecision] code[x$95$m_] := N[(N[Power[2.0, 0.375], $MachinePrecision] * N[(N[Power[2.0, 0.125], $MachinePrecision] * x$95$m), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
x_m = \left|x\right|
\\
{2}^{0.375} \cdot \left({2}^{0.125} \cdot x\_m\right)
\end{array}
Initial program 52.7%
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.7%
Simplified52.7%
pow1/2N/A
sqr-powN/A
associate-*r*N/A
unpow-prod-downN/A
unpow-prod-downN/A
associate-*r*N/A
sqr-powN/A
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
sqr-powN/A
associate-*l*N/A
*-lowering-*.f64N/A
Applied egg-rr48.9%
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
pow-prod-upN/A
pow-lowering-pow.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
pow-lowering-pow.f6449.0%
Applied egg-rr49.0%
x_m = (fabs.f64 x) (FPCore (x_m) :precision binary64 (* (sqrt (* 2.0 x_m)) (sqrt x_m)))
x_m = fabs(x);
double code(double x_m) {
return sqrt((2.0 * x_m)) * sqrt(x_m);
}
x_m = abs(x)
real(8) function code(x_m)
real(8), intent (in) :: x_m
code = sqrt((2.0d0 * x_m)) * sqrt(x_m)
end function
x_m = Math.abs(x);
public static double code(double x_m) {
return Math.sqrt((2.0 * x_m)) * Math.sqrt(x_m);
}
x_m = math.fabs(x) def code(x_m): return math.sqrt((2.0 * x_m)) * math.sqrt(x_m)
x_m = abs(x) function code(x_m) return Float64(sqrt(Float64(2.0 * x_m)) * sqrt(x_m)) end
x_m = abs(x); function tmp = code(x_m) tmp = sqrt((2.0 * x_m)) * sqrt(x_m); end
x_m = N[Abs[x], $MachinePrecision] code[x$95$m_] := N[(N[Sqrt[N[(2.0 * x$95$m), $MachinePrecision]], $MachinePrecision] * N[Sqrt[x$95$m], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
x_m = \left|x\right|
\\
\sqrt{2 \cdot x\_m} \cdot \sqrt{x\_m}
\end{array}
Initial program 52.7%
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.7%
Simplified52.7%
pow1/2N/A
associate-*r*N/A
unpow-prod-downN/A
*-lowering-*.f64N/A
pow1/2N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
pow1/2N/A
sqrt-lowering-sqrt.f6447.8%
Applied egg-rr47.8%
x_m = (fabs.f64 x) (FPCore (x_m) :precision binary64 (* x_m (sqrt 2.0)))
x_m = fabs(x);
double code(double x_m) {
return x_m * sqrt(2.0);
}
x_m = abs(x)
real(8) function code(x_m)
real(8), intent (in) :: x_m
code = x_m * sqrt(2.0d0)
end function
x_m = Math.abs(x);
public static double code(double x_m) {
return x_m * Math.sqrt(2.0);
}
x_m = math.fabs(x) def code(x_m): return x_m * math.sqrt(2.0)
x_m = abs(x) function code(x_m) return Float64(x_m * sqrt(2.0)) end
x_m = abs(x); function tmp = code(x_m) tmp = x_m * sqrt(2.0); end
x_m = N[Abs[x], $MachinePrecision] code[x$95$m_] := N[(x$95$m * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
x_m = \left|x\right|
\\
x\_m \cdot \sqrt{2}
\end{array}
Initial program 52.7%
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.7%
Simplified52.7%
sqrt-prodN/A
pow1/2N/A
sqrt-prodN/A
rem-square-sqrtN/A
*-lowering-*.f64N/A
pow1/2N/A
sqrt-lowering-sqrt.f6448.9%
Applied egg-rr48.9%
Final simplification48.9%
x_m = (fabs.f64 x) (FPCore (x_m) :precision binary64 x_m)
x_m = fabs(x);
double code(double x_m) {
return x_m;
}
x_m = abs(x)
real(8) function code(x_m)
real(8), intent (in) :: x_m
code = x_m
end function
x_m = Math.abs(x);
public static double code(double x_m) {
return x_m;
}
x_m = math.fabs(x) def code(x_m): return x_m
x_m = abs(x) function code(x_m) return x_m end
x_m = abs(x); function tmp = code(x_m) tmp = x_m; end
x_m = N[Abs[x], $MachinePrecision] code[x$95$m_] := x$95$m
\begin{array}{l}
x_m = \left|x\right|
\\
x\_m
\end{array}
Initial program 52.7%
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6452.7%
Simplified52.7%
sqrt-prodN/A
pow1/2N/A
sqrt-prodN/A
rem-square-sqrtN/A
*-lowering-*.f64N/A
pow1/2N/A
sqrt-lowering-sqrt.f6448.9%
Applied egg-rr48.9%
pow1/2N/A
metadata-evalN/A
pow-prod-upN/A
pow-prod-downN/A
metadata-evalN/A
metadata-evalN/A
pow-sqrN/A
pow-prod-downN/A
pow-to-expN/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
metadata-eval48.7%
Applied egg-rr48.7%
Applied egg-rr10.8%
*-rgt-identity10.8%
Applied egg-rr10.8%
herbie shell --seed 2024148
(FPCore (x)
:name "sqrt D (should all be same)"
:precision binary64
(sqrt (* 2.0 (pow x 2.0))))