
(FPCore (x) :precision binary64 (/ x (+ 1.0 (sqrt (+ x 1.0)))))
double code(double x) {
return x / (1.0 + sqrt((x + 1.0)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x / (1.0d0 + sqrt((x + 1.0d0)))
end function
public static double code(double x) {
return x / (1.0 + Math.sqrt((x + 1.0)));
}
def code(x): return x / (1.0 + math.sqrt((x + 1.0)))
function code(x) return Float64(x / Float64(1.0 + sqrt(Float64(x + 1.0)))) end
function tmp = code(x) tmp = x / (1.0 + sqrt((x + 1.0))); end
code[x_] := N[(x / N[(1.0 + N[Sqrt[N[(x + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{1 + \sqrt{x + 1}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 5 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (/ x (+ 1.0 (sqrt (+ x 1.0)))))
double code(double x) {
return x / (1.0 + sqrt((x + 1.0)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x / (1.0d0 + sqrt((x + 1.0d0)))
end function
public static double code(double x) {
return x / (1.0 + Math.sqrt((x + 1.0)));
}
def code(x): return x / (1.0 + math.sqrt((x + 1.0)))
function code(x) return Float64(x / Float64(1.0 + sqrt(Float64(x + 1.0)))) end
function tmp = code(x) tmp = x / (1.0 + sqrt((x + 1.0))); end
code[x_] := N[(x / N[(1.0 + N[Sqrt[N[(x + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{1 + \sqrt{x + 1}}
\end{array}
(FPCore (x) :precision binary64 (/ x (+ 1.0 (sqrt (+ x 1.0)))))
double code(double x) {
return x / (1.0 + sqrt((x + 1.0)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x / (1.0d0 + sqrt((x + 1.0d0)))
end function
public static double code(double x) {
return x / (1.0 + Math.sqrt((x + 1.0)));
}
def code(x): return x / (1.0 + math.sqrt((x + 1.0)))
function code(x) return Float64(x / Float64(1.0 + sqrt(Float64(x + 1.0)))) end
function tmp = code(x) tmp = x / (1.0 + sqrt((x + 1.0))); end
code[x_] := N[(x / N[(1.0 + N[Sqrt[N[(x + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{1 + \sqrt{x + 1}}
\end{array}
Initial program 99.7%
Final simplification99.7%
(FPCore (x) :precision binary64 (if (<= x 7.2e-6) (/ x (+ 2.0 (* x 0.5))) (+ (sqrt (+ x 1.0)) -1.0)))
double code(double x) {
double tmp;
if (x <= 7.2e-6) {
tmp = x / (2.0 + (x * 0.5));
} else {
tmp = sqrt((x + 1.0)) + -1.0;
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 7.2d-6) then
tmp = x / (2.0d0 + (x * 0.5d0))
else
tmp = sqrt((x + 1.0d0)) + (-1.0d0)
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 7.2e-6) {
tmp = x / (2.0 + (x * 0.5));
} else {
tmp = Math.sqrt((x + 1.0)) + -1.0;
}
return tmp;
}
def code(x): tmp = 0 if x <= 7.2e-6: tmp = x / (2.0 + (x * 0.5)) else: tmp = math.sqrt((x + 1.0)) + -1.0 return tmp
function code(x) tmp = 0.0 if (x <= 7.2e-6) tmp = Float64(x / Float64(2.0 + Float64(x * 0.5))); else tmp = Float64(sqrt(Float64(x + 1.0)) + -1.0); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 7.2e-6) tmp = x / (2.0 + (x * 0.5)); else tmp = sqrt((x + 1.0)) + -1.0; end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 7.2e-6], N[(x / N[(2.0 + N[(x * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Sqrt[N[(x + 1.0), $MachinePrecision]], $MachinePrecision] + -1.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 7.2 \cdot 10^{-6}:\\
\;\;\;\;\frac{x}{2 + x \cdot 0.5}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{x + 1} + -1\\
\end{array}
\end{array}
if x < 7.19999999999999967e-6Initial program 100.0%
Taylor expanded in x around 0 99.6%
+-commutative99.6%
unpow299.6%
associate-*r*99.6%
distribute-rgt-out99.5%
*-commutative99.5%
Simplified99.5%
flip-+99.5%
associate-*r/99.6%
metadata-eval99.6%
swap-sqr99.6%
metadata-eval99.6%
*-commutative99.6%
cancel-sign-sub-inv99.6%
metadata-eval99.6%
Applied egg-rr99.6%
associate-/l*99.6%
*-commutative99.6%
associate-*l*99.6%
Simplified99.6%
Taylor expanded in x around 0 99.6%
if 7.19999999999999967e-6 < x Initial program 99.3%
flip-+99.1%
metadata-eval99.1%
add-sqr-sqrt99.8%
+-commutative99.8%
associate--r+99.8%
metadata-eval99.8%
neg-sub099.8%
associate-/r/99.8%
Applied egg-rr99.8%
sub-neg99.8%
+-commutative99.8%
distribute-rgt-in99.8%
distribute-lft-neg-out99.8%
distribute-rgt-neg-in99.8%
neg-mul-199.8%
*-commutative99.8%
associate-/r*99.8%
*-inverses99.8%
metadata-eval99.8%
metadata-eval99.8%
*-rgt-identity99.8%
associate-*r/99.8%
neg-mul-199.8%
times-frac99.8%
metadata-eval99.8%
mul-1-neg99.8%
*-inverses99.8%
metadata-eval99.8%
Simplified99.8%
Final simplification99.7%
(FPCore (x) :precision binary64 (/ x (+ 2.0 (* x 0.5))))
double code(double x) {
return x / (2.0 + (x * 0.5));
}
real(8) function code(x)
real(8), intent (in) :: x
code = x / (2.0d0 + (x * 0.5d0))
end function
public static double code(double x) {
return x / (2.0 + (x * 0.5));
}
def code(x): return x / (2.0 + (x * 0.5))
function code(x) return Float64(x / Float64(2.0 + Float64(x * 0.5))) end
function tmp = code(x) tmp = x / (2.0 + (x * 0.5)); end
code[x_] := N[(x / N[(2.0 + N[(x * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{2 + x \cdot 0.5}
\end{array}
Initial program 99.7%
Taylor expanded in x around 0 67.3%
+-commutative67.3%
unpow267.3%
associate-*r*67.3%
distribute-rgt-out67.3%
*-commutative67.3%
Simplified67.3%
flip-+67.3%
associate-*r/67.3%
metadata-eval67.3%
swap-sqr67.3%
metadata-eval67.3%
*-commutative67.3%
cancel-sign-sub-inv67.3%
metadata-eval67.3%
Applied egg-rr67.3%
associate-/l*67.3%
*-commutative67.3%
associate-*l*67.3%
Simplified67.3%
Taylor expanded in x around 0 69.3%
Final simplification69.3%
(FPCore (x) :precision binary64 (* x 0.5))
double code(double x) {
return x * 0.5;
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * 0.5d0
end function
public static double code(double x) {
return x * 0.5;
}
def code(x): return x * 0.5
function code(x) return Float64(x * 0.5) end
function tmp = code(x) tmp = x * 0.5; end
code[x_] := N[(x * 0.5), $MachinePrecision]
\begin{array}{l}
\\
x \cdot 0.5
\end{array}
Initial program 99.7%
Taylor expanded in x around 0 68.5%
Final simplification68.5%
(FPCore (x) :precision binary64 2.0)
double code(double x) {
return 2.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 2.0d0
end function
public static double code(double x) {
return 2.0;
}
def code(x): return 2.0
function code(x) return 2.0 end
function tmp = code(x) tmp = 2.0; end
code[x_] := 2.0
\begin{array}{l}
\\
2
\end{array}
Initial program 99.7%
Taylor expanded in x around 0 67.3%
+-commutative67.3%
unpow267.3%
associate-*r*67.3%
distribute-rgt-out67.3%
*-commutative67.3%
Simplified67.3%
flip-+67.3%
associate-*r/67.3%
metadata-eval67.3%
swap-sqr67.3%
metadata-eval67.3%
*-commutative67.3%
cancel-sign-sub-inv67.3%
metadata-eval67.3%
Applied egg-rr67.3%
associate-/l*67.3%
*-commutative67.3%
associate-*l*67.3%
Simplified67.3%
Taylor expanded in x around 0 69.3%
Taylor expanded in x around inf 4.8%
Final simplification4.8%
herbie shell --seed 2023287
(FPCore (x)
:name "Numeric.Log:$clog1p from log-domain-0.10.2.1, B"
:precision binary64
(/ x (+ 1.0 (sqrt (+ x 1.0)))))