
(FPCore (x) :precision binary64 (- 1.0 (sqrt (* 0.5 (+ 1.0 (/ 1.0 (hypot 1.0 x)))))))
double code(double x) {
return 1.0 - sqrt((0.5 * (1.0 + (1.0 / hypot(1.0, x)))));
}
public static double code(double x) {
return 1.0 - Math.sqrt((0.5 * (1.0 + (1.0 / Math.hypot(1.0, x)))));
}
def code(x): return 1.0 - math.sqrt((0.5 * (1.0 + (1.0 / math.hypot(1.0, x)))))
function code(x) return Float64(1.0 - sqrt(Float64(0.5 * Float64(1.0 + Float64(1.0 / hypot(1.0, x)))))) end
function tmp = code(x) tmp = 1.0 - sqrt((0.5 * (1.0 + (1.0 / hypot(1.0, x))))); end
code[x_] := N[(1.0 - N[Sqrt[N[(0.5 * N[(1.0 + N[(1.0 / N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 - \sqrt{0.5 \cdot \left(1 + \frac{1}{\mathsf{hypot}\left(1, x\right)}\right)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 13 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (- 1.0 (sqrt (* 0.5 (+ 1.0 (/ 1.0 (hypot 1.0 x)))))))
double code(double x) {
return 1.0 - sqrt((0.5 * (1.0 + (1.0 / hypot(1.0, x)))));
}
public static double code(double x) {
return 1.0 - Math.sqrt((0.5 * (1.0 + (1.0 / Math.hypot(1.0, x)))));
}
def code(x): return 1.0 - math.sqrt((0.5 * (1.0 + (1.0 / math.hypot(1.0, x)))))
function code(x) return Float64(1.0 - sqrt(Float64(0.5 * Float64(1.0 + Float64(1.0 / hypot(1.0, x)))))) end
function tmp = code(x) tmp = 1.0 - sqrt((0.5 * (1.0 + (1.0 / hypot(1.0, x))))); end
code[x_] := N[(1.0 - N[Sqrt[N[(0.5 * N[(1.0 + N[(1.0 / N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 - \sqrt{0.5 \cdot \left(1 + \frac{1}{\mathsf{hypot}\left(1, x\right)}\right)}
\end{array}
(FPCore (x)
:precision binary64
(let* ((t_0 (/ 0.25 (fma x x 1.0))) (t_1 (/ 0.5 (hypot 1.0 x))))
(/
(exp (log (/ (+ t_0 -0.25) (- -0.5 t_1))))
(+ 1.0 (/ (sqrt (- 0.25 t_0)) (sqrt (- 0.5 t_1)))))))
double code(double x) {
double t_0 = 0.25 / fma(x, x, 1.0);
double t_1 = 0.5 / hypot(1.0, x);
return exp(log(((t_0 + -0.25) / (-0.5 - t_1)))) / (1.0 + (sqrt((0.25 - t_0)) / sqrt((0.5 - t_1))));
}
function code(x) t_0 = Float64(0.25 / fma(x, x, 1.0)) t_1 = Float64(0.5 / hypot(1.0, x)) return Float64(exp(log(Float64(Float64(t_0 + -0.25) / Float64(-0.5 - t_1)))) / Float64(1.0 + Float64(sqrt(Float64(0.25 - t_0)) / sqrt(Float64(0.5 - t_1))))) end
code[x_] := Block[{t$95$0 = N[(0.25 / N[(x * x + 1.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(0.5 / N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]}, N[(N[Exp[N[Log[N[(N[(t$95$0 + -0.25), $MachinePrecision] / N[(-0.5 - t$95$1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] / N[(1.0 + N[(N[Sqrt[N[(0.25 - t$95$0), $MachinePrecision]], $MachinePrecision] / N[Sqrt[N[(0.5 - t$95$1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{0.25}{\mathsf{fma}\left(x, x, 1\right)}\\
t_1 := \frac{0.5}{\mathsf{hypot}\left(1, x\right)}\\
\frac{e^{\log \left(\frac{t\_0 + -0.25}{-0.5 - t\_1}\right)}}{1 + \frac{\sqrt{0.25 - t\_0}}{\sqrt{0.5 - t\_1}}}
\end{array}
\end{array}
Initial program 98.4%
distribute-lft-in98.4%
metadata-eval98.4%
associate-*r/98.4%
metadata-eval98.4%
Simplified98.4%
flip--98.4%
metadata-eval98.4%
add-sqr-sqrt99.9%
associate--r+99.9%
metadata-eval99.9%
Applied egg-rr99.9%
flip--99.9%
div-inv99.9%
metadata-eval99.9%
frac-times99.9%
metadata-eval99.9%
pow299.9%
Applied egg-rr99.9%
associate-*r/99.9%
*-rgt-identity99.9%
times-frac99.9%
metadata-eval99.9%
metadata-eval99.9%
times-frac99.9%
*-commutative99.9%
*-commutative99.9%
neg-mul-199.9%
neg-sub099.9%
neg-mul-199.9%
Simplified99.9%
flip-+99.9%
sqrt-div99.9%
metadata-eval99.9%
frac-times99.9%
metadata-eval99.9%
hypot-undefine99.9%
hypot-undefine99.9%
rem-square-sqrt99.9%
metadata-eval99.9%
+-commutative99.9%
fma-define99.9%
Applied egg-rr99.9%
add-exp-log100.0%
+-commutative100.0%
+-commutative100.0%
unpow2100.0%
fma-undefine100.0%
Applied egg-rr100.0%
(FPCore (x)
:precision binary64
(let* ((t_0 (/ 0.25 (fma x x 1.0))) (t_1 (/ 0.5 (hypot 1.0 x))))
(pow
(pow
(/
(+ t_0 -0.25)
(* (- -0.5 t_1) (+ 1.0 (sqrt (/ (- 0.25 t_0) (- 0.5 t_1))))))
3.0)
0.3333333333333333)))
double code(double x) {
double t_0 = 0.25 / fma(x, x, 1.0);
double t_1 = 0.5 / hypot(1.0, x);
return pow(pow(((t_0 + -0.25) / ((-0.5 - t_1) * (1.0 + sqrt(((0.25 - t_0) / (0.5 - t_1)))))), 3.0), 0.3333333333333333);
}
function code(x) t_0 = Float64(0.25 / fma(x, x, 1.0)) t_1 = Float64(0.5 / hypot(1.0, x)) return (Float64(Float64(t_0 + -0.25) / Float64(Float64(-0.5 - t_1) * Float64(1.0 + sqrt(Float64(Float64(0.25 - t_0) / Float64(0.5 - t_1)))))) ^ 3.0) ^ 0.3333333333333333 end
code[x_] := Block[{t$95$0 = N[(0.25 / N[(x * x + 1.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(0.5 / N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]}, N[Power[N[Power[N[(N[(t$95$0 + -0.25), $MachinePrecision] / N[(N[(-0.5 - t$95$1), $MachinePrecision] * N[(1.0 + N[Sqrt[N[(N[(0.25 - t$95$0), $MachinePrecision] / N[(0.5 - t$95$1), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision], 0.3333333333333333], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{0.25}{\mathsf{fma}\left(x, x, 1\right)}\\
t_1 := \frac{0.5}{\mathsf{hypot}\left(1, x\right)}\\
{\left({\left(\frac{t\_0 + -0.25}{\left(-0.5 - t\_1\right) \cdot \left(1 + \sqrt{\frac{0.25 - t\_0}{0.5 - t\_1}}\right)}\right)}^{3}\right)}^{0.3333333333333333}
\end{array}
\end{array}
Initial program 98.4%
distribute-lft-in98.4%
metadata-eval98.4%
associate-*r/98.4%
metadata-eval98.4%
Simplified98.4%
flip--98.4%
metadata-eval98.4%
add-sqr-sqrt99.9%
associate--r+99.9%
metadata-eval99.9%
Applied egg-rr99.9%
flip--99.9%
div-inv99.9%
metadata-eval99.9%
frac-times99.9%
metadata-eval99.9%
pow299.9%
Applied egg-rr99.9%
associate-*r/99.9%
*-rgt-identity99.9%
times-frac99.9%
metadata-eval99.9%
metadata-eval99.9%
times-frac99.9%
*-commutative99.9%
*-commutative99.9%
neg-mul-199.9%
neg-sub099.9%
neg-mul-199.9%
Simplified99.9%
flip-+99.9%
sqrt-div99.9%
metadata-eval99.9%
frac-times99.9%
metadata-eval99.9%
hypot-undefine99.9%
hypot-undefine99.9%
rem-square-sqrt99.9%
metadata-eval99.9%
+-commutative99.9%
fma-define99.9%
Applied egg-rr99.9%
add-cbrt-cube98.5%
pow1/399.9%
Applied egg-rr99.9%
Final simplification99.9%
(FPCore (x)
:precision binary64
(let* ((t_0 (/ -0.5 (hypot 1.0 x))))
(/
(/ (+ (/ 0.25 (fma x x 1.0)) -0.25) (+ -0.5 t_0))
(+ 1.0 (sqrt (/ (+ 0.25 (/ -0.25 (fma x x 1.0))) (+ 0.5 t_0)))))))
double code(double x) {
double t_0 = -0.5 / hypot(1.0, x);
return (((0.25 / fma(x, x, 1.0)) + -0.25) / (-0.5 + t_0)) / (1.0 + sqrt(((0.25 + (-0.25 / fma(x, x, 1.0))) / (0.5 + t_0))));
}
function code(x) t_0 = Float64(-0.5 / hypot(1.0, x)) return Float64(Float64(Float64(Float64(0.25 / fma(x, x, 1.0)) + -0.25) / Float64(-0.5 + t_0)) / Float64(1.0 + sqrt(Float64(Float64(0.25 + Float64(-0.25 / fma(x, x, 1.0))) / Float64(0.5 + t_0))))) end
code[x_] := Block[{t$95$0 = N[(-0.5 / N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]}, N[(N[(N[(N[(0.25 / N[(x * x + 1.0), $MachinePrecision]), $MachinePrecision] + -0.25), $MachinePrecision] / N[(-0.5 + t$95$0), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[Sqrt[N[(N[(0.25 + N[(-0.25 / N[(x * x + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(0.5 + t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{-0.5}{\mathsf{hypot}\left(1, x\right)}\\
\frac{\frac{\frac{0.25}{\mathsf{fma}\left(x, x, 1\right)} + -0.25}{-0.5 + t\_0}}{1 + \sqrt{\frac{0.25 + \frac{-0.25}{\mathsf{fma}\left(x, x, 1\right)}}{0.5 + t\_0}}}
\end{array}
\end{array}
Initial program 98.4%
distribute-lft-in98.4%
metadata-eval98.4%
associate-*r/98.4%
metadata-eval98.4%
Simplified98.4%
flip--98.4%
metadata-eval98.4%
add-sqr-sqrt99.9%
associate--r+99.9%
metadata-eval99.9%
Applied egg-rr99.9%
flip--99.9%
div-inv99.9%
metadata-eval99.9%
frac-times99.9%
metadata-eval99.9%
pow299.9%
Applied egg-rr99.9%
associate-*r/99.9%
*-rgt-identity99.9%
times-frac99.9%
metadata-eval99.9%
metadata-eval99.9%
times-frac99.9%
*-commutative99.9%
*-commutative99.9%
neg-mul-199.9%
neg-sub099.9%
neg-mul-199.9%
Simplified99.9%
flip-+99.9%
sqrt-div99.9%
metadata-eval99.9%
frac-times99.9%
metadata-eval99.9%
hypot-undefine99.9%
hypot-undefine99.9%
rem-square-sqrt99.9%
metadata-eval99.9%
+-commutative99.9%
fma-define99.9%
Applied egg-rr99.9%
*-un-lft-identity99.9%
associate-/l/99.9%
+-commutative99.9%
+-commutative99.9%
unpow299.9%
fma-undefine99.9%
sqrt-undiv99.9%
Applied egg-rr99.9%
associate-*r/99.9%
times-frac99.9%
*-commutative99.9%
associate-*r/99.9%
*-rgt-identity99.9%
+-commutative99.9%
sub-neg99.9%
distribute-neg-frac99.9%
metadata-eval99.9%
Simplified99.9%
Final simplification99.9%
(FPCore (x) :precision binary64 (/ (- 0.5 (sqrt (/ 0.25 (fma x x 1.0)))) (+ 1.0 (sqrt (+ 0.5 (/ 0.5 (hypot 1.0 x)))))))
double code(double x) {
return (0.5 - sqrt((0.25 / fma(x, x, 1.0)))) / (1.0 + sqrt((0.5 + (0.5 / hypot(1.0, x)))));
}
function code(x) return Float64(Float64(0.5 - sqrt(Float64(0.25 / fma(x, x, 1.0)))) / Float64(1.0 + sqrt(Float64(0.5 + Float64(0.5 / hypot(1.0, x)))))) end
code[x_] := N[(N[(0.5 - N[Sqrt[N[(0.25 / N[(x * x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[Sqrt[N[(0.5 + N[(0.5 / N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.5 - \sqrt{\frac{0.25}{\mathsf{fma}\left(x, x, 1\right)}}}{1 + \sqrt{0.5 + \frac{0.5}{\mathsf{hypot}\left(1, x\right)}}}
\end{array}
Initial program 98.4%
distribute-lft-in98.4%
metadata-eval98.4%
associate-*r/98.4%
metadata-eval98.4%
Simplified98.4%
flip--98.4%
metadata-eval98.4%
add-sqr-sqrt99.9%
associate--r+99.9%
metadata-eval99.9%
Applied egg-rr99.9%
add-sqr-sqrt99.9%
sqrt-unprod99.9%
frac-times99.9%
metadata-eval99.9%
hypot-undefine99.9%
hypot-undefine99.9%
rem-square-sqrt99.9%
metadata-eval99.9%
+-commutative99.9%
fma-define99.9%
Applied egg-rr99.9%
(FPCore (x) :precision binary64 (let* ((t_0 (/ 0.5 (hypot 1.0 x)))) (/ (- 0.5 t_0) (+ 1.0 (sqrt (+ 0.5 t_0))))))
double code(double x) {
double t_0 = 0.5 / hypot(1.0, x);
return (0.5 - t_0) / (1.0 + sqrt((0.5 + t_0)));
}
public static double code(double x) {
double t_0 = 0.5 / Math.hypot(1.0, x);
return (0.5 - t_0) / (1.0 + Math.sqrt((0.5 + t_0)));
}
def code(x): t_0 = 0.5 / math.hypot(1.0, x) return (0.5 - t_0) / (1.0 + math.sqrt((0.5 + t_0)))
function code(x) t_0 = Float64(0.5 / hypot(1.0, x)) return Float64(Float64(0.5 - t_0) / Float64(1.0 + sqrt(Float64(0.5 + t_0)))) end
function tmp = code(x) t_0 = 0.5 / hypot(1.0, x); tmp = (0.5 - t_0) / (1.0 + sqrt((0.5 + t_0))); end
code[x_] := Block[{t$95$0 = N[(0.5 / N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]}, N[(N[(0.5 - t$95$0), $MachinePrecision] / N[(1.0 + N[Sqrt[N[(0.5 + t$95$0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{0.5}{\mathsf{hypot}\left(1, x\right)}\\
\frac{0.5 - t\_0}{1 + \sqrt{0.5 + t\_0}}
\end{array}
\end{array}
Initial program 98.4%
distribute-lft-in98.4%
metadata-eval98.4%
associate-*r/98.4%
metadata-eval98.4%
Simplified98.4%
flip--98.4%
metadata-eval98.4%
add-sqr-sqrt99.9%
associate--r+99.9%
metadata-eval99.9%
Applied egg-rr99.9%
(FPCore (x) :precision binary64 (- 1.0 (sqrt (+ 0.5 (sqrt (/ 0.25 (fma x x 1.0)))))))
double code(double x) {
return 1.0 - sqrt((0.5 + sqrt((0.25 / fma(x, x, 1.0)))));
}
function code(x) return Float64(1.0 - sqrt(Float64(0.5 + sqrt(Float64(0.25 / fma(x, x, 1.0)))))) end
code[x_] := N[(1.0 - N[Sqrt[N[(0.5 + N[Sqrt[N[(0.25 / N[(x * x + 1.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 - \sqrt{0.5 + \sqrt{\frac{0.25}{\mathsf{fma}\left(x, x, 1\right)}}}
\end{array}
Initial program 98.4%
distribute-lft-in98.4%
metadata-eval98.4%
associate-*r/98.4%
metadata-eval98.4%
Simplified98.4%
add-sqr-sqrt99.9%
sqrt-unprod99.9%
frac-times99.9%
metadata-eval99.9%
hypot-undefine99.9%
hypot-undefine99.9%
rem-square-sqrt99.9%
metadata-eval99.9%
+-commutative99.9%
fma-define99.9%
Applied egg-rr98.4%
(FPCore (x) :precision binary64 (- 1.0 (sqrt (+ 0.5 (/ 0.5 (hypot 1.0 x))))))
double code(double x) {
return 1.0 - sqrt((0.5 + (0.5 / hypot(1.0, x))));
}
public static double code(double x) {
return 1.0 - Math.sqrt((0.5 + (0.5 / Math.hypot(1.0, x))));
}
def code(x): return 1.0 - math.sqrt((0.5 + (0.5 / math.hypot(1.0, x))))
function code(x) return Float64(1.0 - sqrt(Float64(0.5 + Float64(0.5 / hypot(1.0, x))))) end
function tmp = code(x) tmp = 1.0 - sqrt((0.5 + (0.5 / hypot(1.0, x)))); end
code[x_] := N[(1.0 - N[Sqrt[N[(0.5 + N[(0.5 / N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 - \sqrt{0.5 + \frac{0.5}{\mathsf{hypot}\left(1, x\right)}}
\end{array}
Initial program 98.4%
distribute-lft-in98.4%
metadata-eval98.4%
associate-*r/98.4%
metadata-eval98.4%
Simplified98.4%
(FPCore (x) :precision binary64 (/ (- 0.5 (/ 0.5 x)) (+ 1.0 (sqrt (+ 0.5 (/ 0.5 x))))))
double code(double x) {
return (0.5 - (0.5 / x)) / (1.0 + sqrt((0.5 + (0.5 / x))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = (0.5d0 - (0.5d0 / x)) / (1.0d0 + sqrt((0.5d0 + (0.5d0 / x))))
end function
public static double code(double x) {
return (0.5 - (0.5 / x)) / (1.0 + Math.sqrt((0.5 + (0.5 / x))));
}
def code(x): return (0.5 - (0.5 / x)) / (1.0 + math.sqrt((0.5 + (0.5 / x))))
function code(x) return Float64(Float64(0.5 - Float64(0.5 / x)) / Float64(1.0 + sqrt(Float64(0.5 + Float64(0.5 / x))))) end
function tmp = code(x) tmp = (0.5 - (0.5 / x)) / (1.0 + sqrt((0.5 + (0.5 / x)))); end
code[x_] := N[(N[(0.5 - N[(0.5 / x), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[Sqrt[N[(0.5 + N[(0.5 / x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.5 - \frac{0.5}{x}}{1 + \sqrt{0.5 + \frac{0.5}{x}}}
\end{array}
Initial program 98.4%
distribute-lft-in98.4%
metadata-eval98.4%
associate-*r/98.4%
metadata-eval98.4%
Simplified98.4%
flip--98.4%
metadata-eval98.4%
add-sqr-sqrt99.9%
associate--r+99.9%
metadata-eval99.9%
Applied egg-rr99.9%
Taylor expanded in x around inf 97.2%
Taylor expanded in x around inf 96.9%
associate-*r/96.9%
metadata-eval96.9%
Simplified96.9%
(FPCore (x) :precision binary64 (/ 0.5 (+ 1.0 (sqrt 0.5))))
double code(double x) {
return 0.5 / (1.0 + sqrt(0.5));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 0.5d0 / (1.0d0 + sqrt(0.5d0))
end function
public static double code(double x) {
return 0.5 / (1.0 + Math.sqrt(0.5));
}
def code(x): return 0.5 / (1.0 + math.sqrt(0.5))
function code(x) return Float64(0.5 / Float64(1.0 + sqrt(0.5))) end
function tmp = code(x) tmp = 0.5 / (1.0 + sqrt(0.5)); end
code[x_] := N[(0.5 / N[(1.0 + N[Sqrt[0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.5}{1 + \sqrt{0.5}}
\end{array}
Initial program 98.4%
distribute-lft-in98.4%
metadata-eval98.4%
associate-*r/98.4%
metadata-eval98.4%
Simplified98.4%
flip--98.4%
div-inv98.4%
metadata-eval98.4%
add-sqr-sqrt99.9%
associate--r+99.9%
metadata-eval99.9%
Applied egg-rr99.9%
Taylor expanded in x around inf 96.6%
(FPCore (x) :precision binary64 (- 1.0 (sqrt 0.5)))
double code(double x) {
return 1.0 - sqrt(0.5);
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0 - sqrt(0.5d0)
end function
public static double code(double x) {
return 1.0 - Math.sqrt(0.5);
}
def code(x): return 1.0 - math.sqrt(0.5)
function code(x) return Float64(1.0 - sqrt(0.5)) end
function tmp = code(x) tmp = 1.0 - sqrt(0.5); end
code[x_] := N[(1.0 - N[Sqrt[0.5], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 - \sqrt{0.5}
\end{array}
Initial program 98.4%
distribute-lft-in98.4%
metadata-eval98.4%
associate-*r/98.4%
metadata-eval98.4%
Simplified98.4%
Taylor expanded in x around inf 95.1%
(FPCore (x) :precision binary64 (+ 1.0 (- -1.0 (* (* x x) -0.125))))
double code(double x) {
return 1.0 + (-1.0 - ((x * x) * -0.125));
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.0d0 + ((-1.0d0) - ((x * x) * (-0.125d0)))
end function
public static double code(double x) {
return 1.0 + (-1.0 - ((x * x) * -0.125));
}
def code(x): return 1.0 + (-1.0 - ((x * x) * -0.125))
function code(x) return Float64(1.0 + Float64(-1.0 - Float64(Float64(x * x) * -0.125))) end
function tmp = code(x) tmp = 1.0 + (-1.0 - ((x * x) * -0.125)); end
code[x_] := N[(1.0 + N[(-1.0 - N[(N[(x * x), $MachinePrecision] * -0.125), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 + \left(-1 - \left(x \cdot x\right) \cdot -0.125\right)
\end{array}
Initial program 98.4%
distribute-lft-in98.4%
metadata-eval98.4%
associate-*r/98.4%
metadata-eval98.4%
Simplified98.4%
Taylor expanded in x around 0 4.5%
*-commutative4.5%
Simplified4.5%
unpow24.5%
Applied egg-rr4.5%
Final simplification4.5%
(FPCore (x) :precision binary64 (* x (* x 0.125)))
double code(double x) {
return x * (x * 0.125);
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * (x * 0.125d0)
end function
public static double code(double x) {
return x * (x * 0.125);
}
def code(x): return x * (x * 0.125)
function code(x) return Float64(x * Float64(x * 0.125)) end
function tmp = code(x) tmp = x * (x * 0.125); end
code[x_] := N[(x * N[(x * 0.125), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(x \cdot 0.125\right)
\end{array}
Initial program 98.4%
distribute-lft-in98.4%
metadata-eval98.4%
associate-*r/98.4%
metadata-eval98.4%
Simplified98.4%
Taylor expanded in x around 0 4.5%
*-commutative4.5%
Simplified4.5%
+-commutative4.5%
associate--r+4.5%
add-exp-log4.5%
*-commutative4.5%
cancel-sign-sub-inv4.5%
metadata-eval4.5%
log1p-undefine4.5%
expm1-undefine4.5%
expm1-log1p-u4.5%
unpow24.5%
associate-*r*4.5%
Applied egg-rr4.5%
Final simplification4.5%
(FPCore (x) :precision binary64 0.0)
double code(double x) {
return 0.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 0.0d0
end function
public static double code(double x) {
return 0.0;
}
def code(x): return 0.0
function code(x) return 0.0 end
function tmp = code(x) tmp = 0.0; end
code[x_] := 0.0
\begin{array}{l}
\\
0
\end{array}
Initial program 98.4%
distribute-lft-in98.4%
metadata-eval98.4%
associate-*r/98.4%
metadata-eval98.4%
Simplified98.4%
Taylor expanded in x around 0 3.1%
metadata-eval3.1%
Applied egg-rr3.1%
herbie shell --seed 2024139
(FPCore (x)
:name "Given's Rotation SVD example, simplified"
:precision binary64
(- 1.0 (sqrt (* 0.5 (+ 1.0 (/ 1.0 (hypot 1.0 x)))))))