
(FPCore (x) :precision binary64 (sqrt (* (* 2.0 x) x)))
double code(double x) {
return sqrt(((2.0 * x) * x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = sqrt(((2.0d0 * x) * x))
end function
public static double code(double x) {
return Math.sqrt(((2.0 * x) * x));
}
def code(x): return math.sqrt(((2.0 * x) * x))
function code(x) return sqrt(Float64(Float64(2.0 * x) * x)) end
function tmp = code(x) tmp = sqrt(((2.0 * x) * x)); end
code[x_] := N[Sqrt[N[(N[(2.0 * x), $MachinePrecision] * x), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\sqrt{\left(2 \cdot x\right) \cdot x}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 4 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (sqrt (* (* 2.0 x) x)))
double code(double x) {
return sqrt(((2.0 * x) * x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = sqrt(((2.0d0 * x) * x))
end function
public static double code(double x) {
return Math.sqrt(((2.0 * x) * x));
}
def code(x): return math.sqrt(((2.0 * x) * x))
function code(x) return sqrt(Float64(Float64(2.0 * x) * x)) end
function tmp = code(x) tmp = sqrt(((2.0 * x) * x)); end
code[x_] := N[Sqrt[N[(N[(2.0 * x), $MachinePrecision] * x), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\sqrt{\left(2 \cdot x\right) \cdot x}
\end{array}
(FPCore (x) :precision binary64 (hypot x x))
double code(double x) {
return hypot(x, x);
}
public static double code(double x) {
return Math.hypot(x, x);
}
def code(x): return math.hypot(x, x)
function code(x) return hypot(x, x) end
function tmp = code(x) tmp = hypot(x, x); end
code[x_] := N[Sqrt[x ^ 2 + x ^ 2], $MachinePrecision]
\begin{array}{l}
\\
\mathsf{hypot}\left(x, x\right)
\end{array}
Initial program 50.7%
add-sqr-sqrt50.3%
pow250.3%
pow1/250.3%
associate-*l*50.3%
*-commutative50.3%
unpow-prod-down50.3%
pow1/250.3%
sqrt-unprod50.9%
add-sqr-sqrt51.0%
pow1/251.0%
Applied egg-rr51.0%
unpow251.0%
add-sqr-sqrt52.2%
*-commutative52.2%
add-sqr-sqrt52.1%
associate-*l*52.3%
pow1/252.3%
sqrt-pow152.3%
metadata-eval52.3%
pow1/252.3%
sqrt-pow152.3%
metadata-eval52.3%
Applied egg-rr52.3%
Taylor expanded in x around 0 52.2%
rem-square-sqrt51.0%
fabs-sqr51.0%
rem-square-sqrt99.2%
rem-sqrt-square50.5%
unpow1/250.5%
*-commutative50.5%
unpow1/250.5%
*-commutative50.5%
swap-sqr50.3%
unpow1/250.3%
unpow1/250.3%
rem-square-sqrt50.7%
unpow250.7%
count-250.7%
unpow250.7%
unpow250.7%
hypot-define100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (x) :precision binary64 (* x (* x 2.0)))
double code(double x) {
return x * (x * 2.0);
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * (x * 2.0d0)
end function
public static double code(double x) {
return x * (x * 2.0);
}
def code(x): return x * (x * 2.0)
function code(x) return Float64(x * Float64(x * 2.0)) end
function tmp = code(x) tmp = x * (x * 2.0); end
code[x_] := N[(x * N[(x * 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(x \cdot 2\right)
\end{array}
Initial program 50.7%
add-sqr-sqrt50.3%
pow250.3%
pow1/250.3%
associate-*l*50.3%
*-commutative50.3%
unpow-prod-down50.3%
pow1/250.3%
sqrt-unprod50.9%
add-sqr-sqrt51.0%
pow1/251.0%
Applied egg-rr51.0%
unpow251.0%
add-sqr-sqrt52.2%
*-commutative52.2%
add-sqr-sqrt52.1%
associate-*l*52.3%
pow1/252.3%
sqrt-pow152.3%
metadata-eval52.3%
pow1/252.3%
sqrt-pow152.3%
metadata-eval52.3%
Applied egg-rr52.3%
Taylor expanded in x around 0 52.2%
rem-square-sqrt51.0%
fabs-sqr51.0%
rem-square-sqrt99.2%
rem-sqrt-square50.5%
unpow1/250.5%
*-commutative50.5%
unpow1/250.5%
*-commutative50.5%
swap-sqr50.3%
unpow1/250.3%
unpow1/250.3%
rem-square-sqrt50.7%
unpow250.7%
count-250.7%
unpow250.7%
unpow250.7%
hypot-define100.0%
Simplified100.0%
hypot-undefine50.7%
flip-+0.0%
difference-of-squares0.0%
+-inverses0.0%
metadata-eval0.0%
+-inverses0.0%
metadata-eval0.0%
associate-*r/0.0%
metadata-eval0.0%
+-inverses0.0%
metadata-eval0.0%
+-inverses0.0%
flip-+6.8%
sqrt-unprod7.1%
add-sqr-sqrt7.1%
distribute-lft-out7.1%
*-commutative7.1%
*-un-lft-identity7.1%
distribute-rgt-out7.1%
distribute-lft-out7.1%
metadata-eval7.1%
Applied egg-rr7.1%
Final simplification7.1%
(FPCore (x) :precision binary64 32.0)
double code(double x) {
return 32.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 32.0d0
end function
public static double code(double x) {
return 32.0;
}
def code(x): return 32.0
function code(x) return 32.0 end
function tmp = code(x) tmp = 32.0; end
code[x_] := 32.0
\begin{array}{l}
\\
32
\end{array}
Initial program 50.7%
add-sqr-sqrt50.3%
pow250.3%
pow1/250.3%
associate-*l*50.3%
*-commutative50.3%
unpow-prod-down50.3%
pow1/250.3%
sqrt-unprod50.9%
add-sqr-sqrt51.0%
pow1/251.0%
Applied egg-rr51.0%
unpow251.0%
add-sqr-sqrt52.2%
*-commutative52.2%
add-sqr-sqrt52.1%
associate-*l*52.3%
pow1/252.3%
sqrt-pow152.3%
metadata-eval52.3%
pow1/252.3%
sqrt-pow152.3%
metadata-eval52.3%
Applied egg-rr52.3%
Taylor expanded in x around 0 52.2%
rem-square-sqrt51.0%
fabs-sqr51.0%
rem-square-sqrt99.2%
rem-sqrt-square50.5%
unpow1/250.5%
*-commutative50.5%
unpow1/250.5%
*-commutative50.5%
swap-sqr50.3%
unpow1/250.3%
unpow1/250.3%
rem-square-sqrt50.7%
unpow250.7%
count-250.7%
unpow250.7%
unpow250.7%
hypot-define100.0%
Simplified100.0%
Applied egg-rr0.0%
Simplified5.4%
Final simplification5.4%
(FPCore (x) :precision binary64 64.0)
double code(double x) {
return 64.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 64.0d0
end function
public static double code(double x) {
return 64.0;
}
def code(x): return 64.0
function code(x) return 64.0 end
function tmp = code(x) tmp = 64.0; end
code[x_] := 64.0
\begin{array}{l}
\\
64
\end{array}
Initial program 50.7%
add-sqr-sqrt50.3%
pow250.3%
pow1/250.3%
associate-*l*50.3%
*-commutative50.3%
unpow-prod-down50.3%
pow1/250.3%
sqrt-unprod50.9%
add-sqr-sqrt51.0%
pow1/251.0%
Applied egg-rr51.0%
unpow251.0%
add-sqr-sqrt52.2%
*-commutative52.2%
add-sqr-sqrt52.1%
associate-*l*52.3%
pow1/252.3%
sqrt-pow152.3%
metadata-eval52.3%
pow1/252.3%
sqrt-pow152.3%
metadata-eval52.3%
Applied egg-rr52.3%
Taylor expanded in x around 0 52.2%
rem-square-sqrt51.0%
fabs-sqr51.0%
rem-square-sqrt99.2%
rem-sqrt-square50.5%
unpow1/250.5%
*-commutative50.5%
unpow1/250.5%
*-commutative50.5%
swap-sqr50.3%
unpow1/250.3%
unpow1/250.3%
rem-square-sqrt50.7%
unpow250.7%
count-250.7%
unpow250.7%
unpow250.7%
hypot-define100.0%
Simplified100.0%
Applied egg-rr0.0%
Simplified5.4%
Final simplification5.4%
herbie shell --seed 2024053
(FPCore (x)
:name "sqrt B (should all be same)"
:precision binary64
(sqrt (* (* 2.0 x) x)))