
(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 7 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 (if (<= x 0.0013) (/ x (+ 1.0 (+ 1.0 (* x (+ 0.5 (* x (- (* x 0.0625) 0.125))))))) (+ (sqrt (+ x 1.0)) -1.0)))
double code(double x) {
double tmp;
if (x <= 0.0013) {
tmp = x / (1.0 + (1.0 + (x * (0.5 + (x * ((x * 0.0625) - 0.125))))));
} 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 <= 0.0013d0) then
tmp = x / (1.0d0 + (1.0d0 + (x * (0.5d0 + (x * ((x * 0.0625d0) - 0.125d0))))))
else
tmp = sqrt((x + 1.0d0)) + (-1.0d0)
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 0.0013) {
tmp = x / (1.0 + (1.0 + (x * (0.5 + (x * ((x * 0.0625) - 0.125))))));
} else {
tmp = Math.sqrt((x + 1.0)) + -1.0;
}
return tmp;
}
def code(x): tmp = 0 if x <= 0.0013: tmp = x / (1.0 + (1.0 + (x * (0.5 + (x * ((x * 0.0625) - 0.125)))))) else: tmp = math.sqrt((x + 1.0)) + -1.0 return tmp
function code(x) tmp = 0.0 if (x <= 0.0013) tmp = Float64(x / Float64(1.0 + Float64(1.0 + Float64(x * Float64(0.5 + Float64(x * Float64(Float64(x * 0.0625) - 0.125))))))); else tmp = Float64(sqrt(Float64(x + 1.0)) + -1.0); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 0.0013) tmp = x / (1.0 + (1.0 + (x * (0.5 + (x * ((x * 0.0625) - 0.125)))))); else tmp = sqrt((x + 1.0)) + -1.0; end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 0.0013], N[(x / N[(1.0 + N[(1.0 + N[(x * N[(0.5 + N[(x * N[(N[(x * 0.0625), $MachinePrecision] - 0.125), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $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 0.0013:\\
\;\;\;\;\frac{x}{1 + \left(1 + x \cdot \left(0.5 + x \cdot \left(x \cdot 0.0625 - 0.125\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{x + 1} + -1\\
\end{array}
\end{array}
if x < 0.0012999999999999999Initial program 100.0%
Taylor expanded in x around 0 100.0%
if 0.0012999999999999999 < x Initial program 99.4%
add-log-exp5.0%
*-un-lft-identity5.0%
log-prod5.0%
metadata-eval5.0%
add-log-exp99.4%
frac-2neg99.4%
distribute-frac-neg299.4%
neg-sub099.4%
metadata-eval99.4%
associate--r+99.3%
metadata-eval99.3%
+-commutative99.3%
add-sqr-sqrt99.9%
flip--100.0%
Applied egg-rr100.0%
unsub-neg100.0%
associate--r-100.0%
metadata-eval100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (x) :precision binary64 (if (<= x 3.1) (/ x (+ 1.0 (+ 1.0 (* x (+ 0.5 (* x (- (* x 0.0625) 0.125))))))) (+ -1.0 (sqrt x))))
double code(double x) {
double tmp;
if (x <= 3.1) {
tmp = x / (1.0 + (1.0 + (x * (0.5 + (x * ((x * 0.0625) - 0.125))))));
} else {
tmp = -1.0 + sqrt(x);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 3.1d0) then
tmp = x / (1.0d0 + (1.0d0 + (x * (0.5d0 + (x * ((x * 0.0625d0) - 0.125d0))))))
else
tmp = (-1.0d0) + sqrt(x)
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 3.1) {
tmp = x / (1.0 + (1.0 + (x * (0.5 + (x * ((x * 0.0625) - 0.125))))));
} else {
tmp = -1.0 + Math.sqrt(x);
}
return tmp;
}
def code(x): tmp = 0 if x <= 3.1: tmp = x / (1.0 + (1.0 + (x * (0.5 + (x * ((x * 0.0625) - 0.125)))))) else: tmp = -1.0 + math.sqrt(x) return tmp
function code(x) tmp = 0.0 if (x <= 3.1) tmp = Float64(x / Float64(1.0 + Float64(1.0 + Float64(x * Float64(0.5 + Float64(x * Float64(Float64(x * 0.0625) - 0.125))))))); else tmp = Float64(-1.0 + sqrt(x)); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 3.1) tmp = x / (1.0 + (1.0 + (x * (0.5 + (x * ((x * 0.0625) - 0.125)))))); else tmp = -1.0 + sqrt(x); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 3.1], N[(x / N[(1.0 + N[(1.0 + N[(x * N[(0.5 + N[(x * N[(N[(x * 0.0625), $MachinePrecision] - 0.125), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-1.0 + N[Sqrt[x], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 3.1:\\
\;\;\;\;\frac{x}{1 + \left(1 + x \cdot \left(0.5 + x \cdot \left(x \cdot 0.0625 - 0.125\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;-1 + \sqrt{x}\\
\end{array}
\end{array}
if x < 3.10000000000000009Initial program 100.0%
Taylor expanded in x around 0 100.0%
if 3.10000000000000009 < x Initial program 99.4%
Taylor expanded in x around inf 98.3%
Final simplification99.4%
(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.8%
Final simplification99.8%
(FPCore (x) :precision binary64 (if (<= x 3.7) (/ x (+ 1.0 (+ 1.0 (* x (+ 0.5 (* x (- (* x 0.0625) 0.125))))))) (sqrt x)))
double code(double x) {
double tmp;
if (x <= 3.7) {
tmp = x / (1.0 + (1.0 + (x * (0.5 + (x * ((x * 0.0625) - 0.125))))));
} else {
tmp = sqrt(x);
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 3.7d0) then
tmp = x / (1.0d0 + (1.0d0 + (x * (0.5d0 + (x * ((x * 0.0625d0) - 0.125d0))))))
else
tmp = sqrt(x)
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 3.7) {
tmp = x / (1.0 + (1.0 + (x * (0.5 + (x * ((x * 0.0625) - 0.125))))));
} else {
tmp = Math.sqrt(x);
}
return tmp;
}
def code(x): tmp = 0 if x <= 3.7: tmp = x / (1.0 + (1.0 + (x * (0.5 + (x * ((x * 0.0625) - 0.125)))))) else: tmp = math.sqrt(x) return tmp
function code(x) tmp = 0.0 if (x <= 3.7) tmp = Float64(x / Float64(1.0 + Float64(1.0 + Float64(x * Float64(0.5 + Float64(x * Float64(Float64(x * 0.0625) - 0.125))))))); else tmp = sqrt(x); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 3.7) tmp = x / (1.0 + (1.0 + (x * (0.5 + (x * ((x * 0.0625) - 0.125)))))); else tmp = sqrt(x); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 3.7], N[(x / N[(1.0 + N[(1.0 + N[(x * N[(0.5 + N[(x * N[(N[(x * 0.0625), $MachinePrecision] - 0.125), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Sqrt[x], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 3.7:\\
\;\;\;\;\frac{x}{1 + \left(1 + x \cdot \left(0.5 + x \cdot \left(x \cdot 0.0625 - 0.125\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{x}\\
\end{array}
\end{array}
if x < 3.7000000000000002Initial program 100.0%
Taylor expanded in x around 0 100.0%
if 3.7000000000000002 < x Initial program 99.4%
Taylor expanded in x around inf 94.9%
Final simplification98.3%
(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.8%
Taylor expanded in x around 0 67.8%
+-commutative67.8%
*-commutative67.8%
Simplified67.8%
Taylor expanded in x around 0 69.4%
*-commutative69.4%
Simplified69.4%
Final simplification69.4%
(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.8%
Taylor expanded in x around 0 69.1%
Final simplification69.1%
(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.8%
Taylor expanded in x around 0 69.4%
+-commutative69.4%
Simplified69.4%
Taylor expanded in x around inf 4.9%
Final simplification4.9%
herbie shell --seed 2024078
(FPCore (x)
:name "Numeric.Log:$clog1p from log-domain-0.10.2.1, B"
:precision binary64
(/ x (+ 1.0 (sqrt (+ x 1.0)))))