
(FPCore (t) :precision binary64 (let* ((t_1 (- 2.0 (/ (/ 2.0 t) (+ 1.0 (/ 1.0 t))))) (t_2 (* t_1 t_1))) (/ (+ 1.0 t_2) (+ 2.0 t_2))))
double code(double t) {
double t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t)));
double t_2 = t_1 * t_1;
return (1.0 + t_2) / (2.0 + t_2);
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: t_2
t_1 = 2.0d0 - ((2.0d0 / t) / (1.0d0 + (1.0d0 / t)))
t_2 = t_1 * t_1
code = (1.0d0 + t_2) / (2.0d0 + t_2)
end function
public static double code(double t) {
double t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t)));
double t_2 = t_1 * t_1;
return (1.0 + t_2) / (2.0 + t_2);
}
def code(t): t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t))) t_2 = t_1 * t_1 return (1.0 + t_2) / (2.0 + t_2)
function code(t) t_1 = Float64(2.0 - Float64(Float64(2.0 / t) / Float64(1.0 + Float64(1.0 / t)))) t_2 = Float64(t_1 * t_1) return Float64(Float64(1.0 + t_2) / Float64(2.0 + t_2)) end
function tmp = code(t) t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t))); t_2 = t_1 * t_1; tmp = (1.0 + t_2) / (2.0 + t_2); end
code[t_] := Block[{t$95$1 = N[(2.0 - N[(N[(2.0 / t), $MachinePrecision] / N[(1.0 + N[(1.0 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * t$95$1), $MachinePrecision]}, N[(N[(1.0 + t$95$2), $MachinePrecision] / N[(2.0 + t$95$2), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 2 - \frac{\frac{2}{t}}{1 + \frac{1}{t}}\\
t_2 := t\_1 \cdot t\_1\\
\frac{1 + t\_2}{2 + t\_2}
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 6 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (t) :precision binary64 (let* ((t_1 (- 2.0 (/ (/ 2.0 t) (+ 1.0 (/ 1.0 t))))) (t_2 (* t_1 t_1))) (/ (+ 1.0 t_2) (+ 2.0 t_2))))
double code(double t) {
double t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t)));
double t_2 = t_1 * t_1;
return (1.0 + t_2) / (2.0 + t_2);
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: t_2
t_1 = 2.0d0 - ((2.0d0 / t) / (1.0d0 + (1.0d0 / t)))
t_2 = t_1 * t_1
code = (1.0d0 + t_2) / (2.0d0 + t_2)
end function
public static double code(double t) {
double t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t)));
double t_2 = t_1 * t_1;
return (1.0 + t_2) / (2.0 + t_2);
}
def code(t): t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t))) t_2 = t_1 * t_1 return (1.0 + t_2) / (2.0 + t_2)
function code(t) t_1 = Float64(2.0 - Float64(Float64(2.0 / t) / Float64(1.0 + Float64(1.0 / t)))) t_2 = Float64(t_1 * t_1) return Float64(Float64(1.0 + t_2) / Float64(2.0 + t_2)) end
function tmp = code(t) t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t))); t_2 = t_1 * t_1; tmp = (1.0 + t_2) / (2.0 + t_2); end
code[t_] := Block[{t$95$1 = N[(2.0 - N[(N[(2.0 / t), $MachinePrecision] / N[(1.0 + N[(1.0 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * t$95$1), $MachinePrecision]}, N[(N[(1.0 + t$95$2), $MachinePrecision] / N[(2.0 + t$95$2), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 2 - \frac{\frac{2}{t}}{1 + \frac{1}{t}}\\
t_2 := t\_1 \cdot t\_1\\
\frac{1 + t\_2}{2 + t\_2}
\end{array}
\end{array}
(FPCore (t) :precision binary64 (let* ((t_1 (- (/ -2.0 (- -1.0 t)) 2.0)) (t_2 (+ 2.0 (/ 2.0 (- -1.0 t))))) (/ (+ 1.0 (* t_1 t_1)) (+ 2.0 (* t_2 t_2)))))
double code(double t) {
double t_1 = (-2.0 / (-1.0 - t)) - 2.0;
double t_2 = 2.0 + (2.0 / (-1.0 - t));
return (1.0 + (t_1 * t_1)) / (2.0 + (t_2 * t_2));
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: t_2
t_1 = ((-2.0d0) / ((-1.0d0) - t)) - 2.0d0
t_2 = 2.0d0 + (2.0d0 / ((-1.0d0) - t))
code = (1.0d0 + (t_1 * t_1)) / (2.0d0 + (t_2 * t_2))
end function
public static double code(double t) {
double t_1 = (-2.0 / (-1.0 - t)) - 2.0;
double t_2 = 2.0 + (2.0 / (-1.0 - t));
return (1.0 + (t_1 * t_1)) / (2.0 + (t_2 * t_2));
}
def code(t): t_1 = (-2.0 / (-1.0 - t)) - 2.0 t_2 = 2.0 + (2.0 / (-1.0 - t)) return (1.0 + (t_1 * t_1)) / (2.0 + (t_2 * t_2))
function code(t) t_1 = Float64(Float64(-2.0 / Float64(-1.0 - t)) - 2.0) t_2 = Float64(2.0 + Float64(2.0 / Float64(-1.0 - t))) return Float64(Float64(1.0 + Float64(t_1 * t_1)) / Float64(2.0 + Float64(t_2 * t_2))) end
function tmp = code(t) t_1 = (-2.0 / (-1.0 - t)) - 2.0; t_2 = 2.0 + (2.0 / (-1.0 - t)); tmp = (1.0 + (t_1 * t_1)) / (2.0 + (t_2 * t_2)); end
code[t_] := Block[{t$95$1 = N[(N[(-2.0 / N[(-1.0 - t), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision]}, Block[{t$95$2 = N[(2.0 + N[(2.0 / N[(-1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(N[(1.0 + N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[(t$95$2 * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{-2}{-1 - t} - 2\\
t_2 := 2 + \frac{2}{-1 - t}\\
\frac{1 + t\_1 \cdot t\_1}{2 + t\_2 \cdot t\_2}
\end{array}
\end{array}
Initial program 100.0%
add-log-exp100.0%
*-un-lft-identity100.0%
log-prod100.0%
metadata-eval100.0%
add-log-exp100.0%
associate-/l/100.0%
*-commutative100.0%
Applied egg-rr100.0%
+-lft-identity100.0%
distribute-rgt-in100.0%
lft-mult-inverse100.0%
*-lft-identity100.0%
Simplified100.0%
add-log-exp100.0%
*-un-lft-identity100.0%
log-prod100.0%
metadata-eval100.0%
add-log-exp100.0%
associate-/l/100.0%
*-commutative100.0%
Applied egg-rr100.0%
+-lft-identity100.0%
distribute-rgt-in100.0%
lft-mult-inverse100.0%
*-lft-identity100.0%
Simplified100.0%
sub-neg100.0%
distribute-neg-frac100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Applied egg-rr100.0%
associate-/r*100.0%
distribute-rgt-in100.0%
lft-mult-inverse100.0%
*-lft-identity100.0%
Simplified100.0%
sub-neg100.0%
distribute-neg-frac100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Applied egg-rr100.0%
associate-/r*100.0%
distribute-rgt-in100.0%
lft-mult-inverse100.0%
*-lft-identity100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (t)
:precision binary64
(let* ((t_1 (+ 2.0 (/ 2.0 (- -1.0 t)))))
(if (<= (+ 2.0 (/ (/ 2.0 t) (- -1.0 (/ 1.0 t)))) 0.2)
(/ (- 1.0 (* (* 2.0 t) (- (/ -2.0 (- -1.0 t)) 2.0))) (+ 2.0 (* t_1 t_1)))
(+
0.8333333333333334
(/
(-
(/ (+ 0.037037037037037035 (/ 0.04938271604938271 t)) t)
0.2222222222222222)
t)))))
double code(double t) {
double t_1 = 2.0 + (2.0 / (-1.0 - t));
double tmp;
if ((2.0 + ((2.0 / t) / (-1.0 - (1.0 / t)))) <= 0.2) {
tmp = (1.0 - ((2.0 * t) * ((-2.0 / (-1.0 - t)) - 2.0))) / (2.0 + (t_1 * t_1));
} else {
tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t);
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = 2.0d0 + (2.0d0 / ((-1.0d0) - t))
if ((2.0d0 + ((2.0d0 / t) / ((-1.0d0) - (1.0d0 / t)))) <= 0.2d0) then
tmp = (1.0d0 - ((2.0d0 * t) * (((-2.0d0) / ((-1.0d0) - t)) - 2.0d0))) / (2.0d0 + (t_1 * t_1))
else
tmp = 0.8333333333333334d0 + ((((0.037037037037037035d0 + (0.04938271604938271d0 / t)) / t) - 0.2222222222222222d0) / t)
end if
code = tmp
end function
public static double code(double t) {
double t_1 = 2.0 + (2.0 / (-1.0 - t));
double tmp;
if ((2.0 + ((2.0 / t) / (-1.0 - (1.0 / t)))) <= 0.2) {
tmp = (1.0 - ((2.0 * t) * ((-2.0 / (-1.0 - t)) - 2.0))) / (2.0 + (t_1 * t_1));
} else {
tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t);
}
return tmp;
}
def code(t): t_1 = 2.0 + (2.0 / (-1.0 - t)) tmp = 0 if (2.0 + ((2.0 / t) / (-1.0 - (1.0 / t)))) <= 0.2: tmp = (1.0 - ((2.0 * t) * ((-2.0 / (-1.0 - t)) - 2.0))) / (2.0 + (t_1 * t_1)) else: tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t) return tmp
function code(t) t_1 = Float64(2.0 + Float64(2.0 / Float64(-1.0 - t))) tmp = 0.0 if (Float64(2.0 + Float64(Float64(2.0 / t) / Float64(-1.0 - Float64(1.0 / t)))) <= 0.2) tmp = Float64(Float64(1.0 - Float64(Float64(2.0 * t) * Float64(Float64(-2.0 / Float64(-1.0 - t)) - 2.0))) / Float64(2.0 + Float64(t_1 * t_1))); else tmp = Float64(0.8333333333333334 + Float64(Float64(Float64(Float64(0.037037037037037035 + Float64(0.04938271604938271 / t)) / t) - 0.2222222222222222) / t)); end return tmp end
function tmp_2 = code(t) t_1 = 2.0 + (2.0 / (-1.0 - t)); tmp = 0.0; if ((2.0 + ((2.0 / t) / (-1.0 - (1.0 / t)))) <= 0.2) tmp = (1.0 - ((2.0 * t) * ((-2.0 / (-1.0 - t)) - 2.0))) / (2.0 + (t_1 * t_1)); else tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t); end tmp_2 = tmp; end
code[t_] := Block[{t$95$1 = N[(2.0 + N[(2.0 / N[(-1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(2.0 + N[(N[(2.0 / t), $MachinePrecision] / N[(-1.0 - N[(1.0 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.2], N[(N[(1.0 - N[(N[(2.0 * t), $MachinePrecision] * N[(N[(-2.0 / N[(-1.0 - t), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.8333333333333334 + N[(N[(N[(N[(0.037037037037037035 + N[(0.04938271604938271 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] - 0.2222222222222222), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 2 + \frac{2}{-1 - t}\\
\mathbf{if}\;2 + \frac{\frac{2}{t}}{-1 - \frac{1}{t}} \leq 0.2:\\
\;\;\;\;\frac{1 - \left(2 \cdot t\right) \cdot \left(\frac{-2}{-1 - t} - 2\right)}{2 + t\_1 \cdot t\_1}\\
\mathbf{else}:\\
\;\;\;\;0.8333333333333334 + \frac{\frac{0.037037037037037035 + \frac{0.04938271604938271}{t}}{t} - 0.2222222222222222}{t}\\
\end{array}
\end{array}
if (-.f64 #s(literal 2 binary64) (/.f64 (/.f64 #s(literal 2 binary64) t) (+.f64 #s(literal 1 binary64) (/.f64 #s(literal 1 binary64) t)))) < 0.20000000000000001Initial program 100.0%
add-log-exp100.0%
*-un-lft-identity100.0%
log-prod100.0%
metadata-eval100.0%
add-log-exp100.0%
associate-/l/100.0%
*-commutative100.0%
Applied egg-rr100.0%
+-lft-identity100.0%
distribute-rgt-in100.0%
lft-mult-inverse100.0%
*-lft-identity100.0%
Simplified100.0%
add-log-exp100.0%
*-un-lft-identity100.0%
log-prod100.0%
metadata-eval100.0%
add-log-exp100.0%
associate-/l/100.0%
*-commutative100.0%
Applied egg-rr100.0%
+-lft-identity100.0%
distribute-rgt-in100.0%
lft-mult-inverse100.0%
*-lft-identity100.0%
Simplified100.0%
sub-neg100.0%
distribute-neg-frac100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Applied egg-rr100.0%
associate-/r*100.0%
distribute-rgt-in100.0%
lft-mult-inverse100.0%
*-lft-identity100.0%
Simplified100.0%
Taylor expanded in t around 0 99.1%
if 0.20000000000000001 < (-.f64 #s(literal 2 binary64) (/.f64 (/.f64 #s(literal 2 binary64) t) (+.f64 #s(literal 1 binary64) (/.f64 #s(literal 1 binary64) t)))) Initial program 100.0%
Taylor expanded in t around -inf 99.1%
mul-1-neg99.1%
unsub-neg99.1%
mul-1-neg99.1%
unsub-neg99.1%
sub-neg99.1%
associate-*r/99.1%
metadata-eval99.1%
distribute-neg-frac99.1%
metadata-eval99.1%
Simplified99.1%
Taylor expanded in t around -inf 99.4%
mul-1-neg99.4%
unsub-neg99.4%
mul-1-neg99.4%
unsub-neg99.4%
associate-*r/99.4%
metadata-eval99.4%
Simplified99.4%
Final simplification99.3%
(FPCore (t)
:precision binary64
(let* ((t_1 (- 4.0 (/ (+ 8.0 (/ -12.0 t)) t))))
(if (or (<= t -5.5e-17) (not (<= t 1.12e-16)))
(/ (+ 1.0 t_1) (+ 2.0 t_1))
(/ 1.0 (+ 2.0 (* (- (/ -2.0 (- -1.0 t)) 2.0) (- (/ 2.0 t) 2.0)))))))
double code(double t) {
double t_1 = 4.0 - ((8.0 + (-12.0 / t)) / t);
double tmp;
if ((t <= -5.5e-17) || !(t <= 1.12e-16)) {
tmp = (1.0 + t_1) / (2.0 + t_1);
} else {
tmp = 1.0 / (2.0 + (((-2.0 / (-1.0 - t)) - 2.0) * ((2.0 / t) - 2.0)));
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = 4.0d0 - ((8.0d0 + ((-12.0d0) / t)) / t)
if ((t <= (-5.5d-17)) .or. (.not. (t <= 1.12d-16))) then
tmp = (1.0d0 + t_1) / (2.0d0 + t_1)
else
tmp = 1.0d0 / (2.0d0 + ((((-2.0d0) / ((-1.0d0) - t)) - 2.0d0) * ((2.0d0 / t) - 2.0d0)))
end if
code = tmp
end function
public static double code(double t) {
double t_1 = 4.0 - ((8.0 + (-12.0 / t)) / t);
double tmp;
if ((t <= -5.5e-17) || !(t <= 1.12e-16)) {
tmp = (1.0 + t_1) / (2.0 + t_1);
} else {
tmp = 1.0 / (2.0 + (((-2.0 / (-1.0 - t)) - 2.0) * ((2.0 / t) - 2.0)));
}
return tmp;
}
def code(t): t_1 = 4.0 - ((8.0 + (-12.0 / t)) / t) tmp = 0 if (t <= -5.5e-17) or not (t <= 1.12e-16): tmp = (1.0 + t_1) / (2.0 + t_1) else: tmp = 1.0 / (2.0 + (((-2.0 / (-1.0 - t)) - 2.0) * ((2.0 / t) - 2.0))) return tmp
function code(t) t_1 = Float64(4.0 - Float64(Float64(8.0 + Float64(-12.0 / t)) / t)) tmp = 0.0 if ((t <= -5.5e-17) || !(t <= 1.12e-16)) tmp = Float64(Float64(1.0 + t_1) / Float64(2.0 + t_1)); else tmp = Float64(1.0 / Float64(2.0 + Float64(Float64(Float64(-2.0 / Float64(-1.0 - t)) - 2.0) * Float64(Float64(2.0 / t) - 2.0)))); end return tmp end
function tmp_2 = code(t) t_1 = 4.0 - ((8.0 + (-12.0 / t)) / t); tmp = 0.0; if ((t <= -5.5e-17) || ~((t <= 1.12e-16))) tmp = (1.0 + t_1) / (2.0 + t_1); else tmp = 1.0 / (2.0 + (((-2.0 / (-1.0 - t)) - 2.0) * ((2.0 / t) - 2.0))); end tmp_2 = tmp; end
code[t_] := Block[{t$95$1 = N[(4.0 - N[(N[(8.0 + N[(-12.0 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[t, -5.5e-17], N[Not[LessEqual[t, 1.12e-16]], $MachinePrecision]], N[(N[(1.0 + t$95$1), $MachinePrecision] / N[(2.0 + t$95$1), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(2.0 + N[(N[(N[(-2.0 / N[(-1.0 - t), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision] * N[(N[(2.0 / t), $MachinePrecision] - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 4 - \frac{8 + \frac{-12}{t}}{t}\\
\mathbf{if}\;t \leq -5.5 \cdot 10^{-17} \lor \neg \left(t \leq 1.12 \cdot 10^{-16}\right):\\
\;\;\;\;\frac{1 + t\_1}{2 + t\_1}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{2 + \left(\frac{-2}{-1 - t} - 2\right) \cdot \left(\frac{2}{t} - 2\right)}\\
\end{array}
\end{array}
if t < -5.50000000000000001e-17 or 1.12e-16 < t Initial program 100.0%
Taylor expanded in t around inf 93.4%
+-commutative93.4%
associate--l+93.4%
+-commutative93.4%
associate--r-93.4%
associate-*r/93.4%
metadata-eval93.4%
unpow293.4%
associate-/r*93.4%
metadata-eval93.4%
associate-*r/93.4%
div-sub93.4%
sub-neg93.4%
associate-*r/93.4%
metadata-eval93.4%
distribute-neg-frac93.4%
metadata-eval93.4%
Simplified93.4%
Taylor expanded in t around inf 94.1%
+-commutative93.4%
associate--l+93.4%
+-commutative93.4%
associate--r-93.4%
associate-*r/93.4%
metadata-eval93.4%
unpow293.4%
associate-/r*93.4%
metadata-eval93.4%
associate-*r/93.4%
div-sub93.4%
sub-neg93.4%
associate-*r/93.4%
metadata-eval93.4%
distribute-neg-frac93.4%
metadata-eval93.4%
Simplified94.1%
if -5.50000000000000001e-17 < t < 1.12e-16Initial program 100.0%
Taylor expanded in t around inf 100.0%
associate-*r/100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
Taylor expanded in t around 0 100.0%
sub-neg100.0%
distribute-neg-frac100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Applied egg-rr100.0%
associate-/r*100.0%
distribute-rgt-in100.0%
lft-mult-inverse100.0%
*-lft-identity100.0%
Simplified100.0%
Final simplification96.6%
(FPCore (t) :precision binary64 (let* ((t_1 (+ 2.0 (/ 2.0 (- -1.0 t))))) (/ (- 1.0 (* 2.0 (- (/ -2.0 (- -1.0 t)) 2.0))) (+ 2.0 (* t_1 t_1)))))
double code(double t) {
double t_1 = 2.0 + (2.0 / (-1.0 - t));
return (1.0 - (2.0 * ((-2.0 / (-1.0 - t)) - 2.0))) / (2.0 + (t_1 * t_1));
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
t_1 = 2.0d0 + (2.0d0 / ((-1.0d0) - t))
code = (1.0d0 - (2.0d0 * (((-2.0d0) / ((-1.0d0) - t)) - 2.0d0))) / (2.0d0 + (t_1 * t_1))
end function
public static double code(double t) {
double t_1 = 2.0 + (2.0 / (-1.0 - t));
return (1.0 - (2.0 * ((-2.0 / (-1.0 - t)) - 2.0))) / (2.0 + (t_1 * t_1));
}
def code(t): t_1 = 2.0 + (2.0 / (-1.0 - t)) return (1.0 - (2.0 * ((-2.0 / (-1.0 - t)) - 2.0))) / (2.0 + (t_1 * t_1))
function code(t) t_1 = Float64(2.0 + Float64(2.0 / Float64(-1.0 - t))) return Float64(Float64(1.0 - Float64(2.0 * Float64(Float64(-2.0 / Float64(-1.0 - t)) - 2.0))) / Float64(2.0 + Float64(t_1 * t_1))) end
function tmp = code(t) t_1 = 2.0 + (2.0 / (-1.0 - t)); tmp = (1.0 - (2.0 * ((-2.0 / (-1.0 - t)) - 2.0))) / (2.0 + (t_1 * t_1)); end
code[t_] := Block[{t$95$1 = N[(2.0 + N[(2.0 / N[(-1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(N[(1.0 - N[(2.0 * N[(N[(-2.0 / N[(-1.0 - t), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 2 + \frac{2}{-1 - t}\\
\frac{1 - 2 \cdot \left(\frac{-2}{-1 - t} - 2\right)}{2 + t\_1 \cdot t\_1}
\end{array}
\end{array}
Initial program 100.0%
add-log-exp100.0%
*-un-lft-identity100.0%
log-prod100.0%
metadata-eval100.0%
add-log-exp100.0%
associate-/l/100.0%
*-commutative100.0%
Applied egg-rr100.0%
+-lft-identity100.0%
distribute-rgt-in100.0%
lft-mult-inverse100.0%
*-lft-identity100.0%
Simplified100.0%
add-log-exp100.0%
*-un-lft-identity100.0%
log-prod100.0%
metadata-eval100.0%
add-log-exp100.0%
associate-/l/100.0%
*-commutative100.0%
Applied egg-rr100.0%
+-lft-identity100.0%
distribute-rgt-in100.0%
lft-mult-inverse100.0%
*-lft-identity100.0%
Simplified100.0%
sub-neg100.0%
distribute-neg-frac100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Applied egg-rr100.0%
associate-/r*100.0%
distribute-rgt-in100.0%
lft-mult-inverse100.0%
*-lft-identity100.0%
Simplified100.0%
Taylor expanded in t around inf 97.3%
Final simplification97.3%
(FPCore (t)
:precision binary64
(if (<= t -5.5e-17)
(-
0.8333333333333334
(/ (+ 0.2222222222222222 (/ -0.037037037037037035 t)) t))
(if (<= t 1.12e-16)
(/ 1.0 (+ 2.0 (* (- (/ -2.0 (- -1.0 t)) 2.0) (- (/ 2.0 t) 2.0))))
(+
0.8333333333333334
(/
(-
(/ (+ 0.037037037037037035 (/ 0.04938271604938271 t)) t)
0.2222222222222222)
t)))))
double code(double t) {
double tmp;
if (t <= -5.5e-17) {
tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t);
} else if (t <= 1.12e-16) {
tmp = 1.0 / (2.0 + (((-2.0 / (-1.0 - t)) - 2.0) * ((2.0 / t) - 2.0)));
} else {
tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t);
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-5.5d-17)) then
tmp = 0.8333333333333334d0 - ((0.2222222222222222d0 + ((-0.037037037037037035d0) / t)) / t)
else if (t <= 1.12d-16) then
tmp = 1.0d0 / (2.0d0 + ((((-2.0d0) / ((-1.0d0) - t)) - 2.0d0) * ((2.0d0 / t) - 2.0d0)))
else
tmp = 0.8333333333333334d0 + ((((0.037037037037037035d0 + (0.04938271604938271d0 / t)) / t) - 0.2222222222222222d0) / t)
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if (t <= -5.5e-17) {
tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t);
} else if (t <= 1.12e-16) {
tmp = 1.0 / (2.0 + (((-2.0 / (-1.0 - t)) - 2.0) * ((2.0 / t) - 2.0)));
} else {
tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t);
}
return tmp;
}
def code(t): tmp = 0 if t <= -5.5e-17: tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t) elif t <= 1.12e-16: tmp = 1.0 / (2.0 + (((-2.0 / (-1.0 - t)) - 2.0) * ((2.0 / t) - 2.0))) else: tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t) return tmp
function code(t) tmp = 0.0 if (t <= -5.5e-17) tmp = Float64(0.8333333333333334 - Float64(Float64(0.2222222222222222 + Float64(-0.037037037037037035 / t)) / t)); elseif (t <= 1.12e-16) tmp = Float64(1.0 / Float64(2.0 + Float64(Float64(Float64(-2.0 / Float64(-1.0 - t)) - 2.0) * Float64(Float64(2.0 / t) - 2.0)))); else tmp = Float64(0.8333333333333334 + Float64(Float64(Float64(Float64(0.037037037037037035 + Float64(0.04938271604938271 / t)) / t) - 0.2222222222222222) / t)); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if (t <= -5.5e-17) tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t); elseif (t <= 1.12e-16) tmp = 1.0 / (2.0 + (((-2.0 / (-1.0 - t)) - 2.0) * ((2.0 / t) - 2.0))); else tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t); end tmp_2 = tmp; end
code[t_] := If[LessEqual[t, -5.5e-17], N[(0.8333333333333334 - N[(N[(0.2222222222222222 + N[(-0.037037037037037035 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.12e-16], N[(1.0 / N[(2.0 + N[(N[(N[(-2.0 / N[(-1.0 - t), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision] * N[(N[(2.0 / t), $MachinePrecision] - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.8333333333333334 + N[(N[(N[(N[(0.037037037037037035 + N[(0.04938271604938271 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] - 0.2222222222222222), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -5.5 \cdot 10^{-17}:\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222 + \frac{-0.037037037037037035}{t}}{t}\\
\mathbf{elif}\;t \leq 1.12 \cdot 10^{-16}:\\
\;\;\;\;\frac{1}{2 + \left(\frac{-2}{-1 - t} - 2\right) \cdot \left(\frac{2}{t} - 2\right)}\\
\mathbf{else}:\\
\;\;\;\;0.8333333333333334 + \frac{\frac{0.037037037037037035 + \frac{0.04938271604938271}{t}}{t} - 0.2222222222222222}{t}\\
\end{array}
\end{array}
if t < -5.50000000000000001e-17Initial program 100.0%
Taylor expanded in t around inf 92.7%
+-commutative92.7%
associate--l+92.7%
+-commutative92.7%
associate--r-92.8%
associate-*r/92.8%
metadata-eval92.8%
unpow292.8%
associate-/r*92.8%
metadata-eval92.8%
associate-*r/92.8%
div-sub92.8%
sub-neg92.8%
associate-*r/92.8%
metadata-eval92.8%
distribute-neg-frac92.8%
metadata-eval92.8%
Simplified92.8%
Taylor expanded in t around inf 92.9%
sub-neg92.9%
+-commutative92.9%
associate-+r+92.9%
+-commutative92.9%
sub-neg92.9%
associate--r-92.9%
associate-*r/92.9%
metadata-eval92.9%
unpow292.9%
associate-/r*92.9%
metadata-eval92.9%
associate-*r/92.9%
div-sub92.9%
sub-neg92.9%
associate-*r/92.9%
metadata-eval92.9%
distribute-neg-frac92.9%
metadata-eval92.9%
Simplified92.9%
if -5.50000000000000001e-17 < t < 1.12e-16Initial program 100.0%
Taylor expanded in t around inf 100.0%
associate-*r/100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
Taylor expanded in t around 0 100.0%
sub-neg100.0%
distribute-neg-frac100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Applied egg-rr100.0%
associate-/r*100.0%
distribute-rgt-in100.0%
lft-mult-inverse100.0%
*-lft-identity100.0%
Simplified100.0%
if 1.12e-16 < t Initial program 100.0%
Taylor expanded in t around -inf 93.6%
mul-1-neg93.6%
unsub-neg93.6%
mul-1-neg93.6%
unsub-neg93.6%
sub-neg93.6%
associate-*r/93.6%
metadata-eval93.6%
distribute-neg-frac93.6%
metadata-eval93.6%
Simplified93.6%
Taylor expanded in t around -inf 94.4%
mul-1-neg94.4%
unsub-neg94.4%
mul-1-neg94.4%
unsub-neg94.4%
associate-*r/94.4%
metadata-eval94.4%
Simplified94.4%
Final simplification96.4%
(FPCore (t) :precision binary64 0.8333333333333334)
double code(double t) {
return 0.8333333333333334;
}
real(8) function code(t)
real(8), intent (in) :: t
code = 0.8333333333333334d0
end function
public static double code(double t) {
return 0.8333333333333334;
}
def code(t): return 0.8333333333333334
function code(t) return 0.8333333333333334 end
function tmp = code(t) tmp = 0.8333333333333334; end
code[t_] := 0.8333333333333334
\begin{array}{l}
\\
0.8333333333333334
\end{array}
Initial program 100.0%
Taylor expanded in t around inf 54.7%
+-commutative54.7%
associate--l+54.7%
+-commutative54.7%
associate--r-54.7%
associate-*r/54.7%
metadata-eval54.7%
unpow254.7%
associate-/r*54.7%
metadata-eval54.7%
associate-*r/54.7%
div-sub54.7%
sub-neg54.7%
associate-*r/54.7%
metadata-eval54.7%
distribute-neg-frac54.7%
metadata-eval54.7%
Simplified54.7%
Taylor expanded in t around inf 61.1%
Final simplification61.1%
herbie shell --seed 2024071
(FPCore (t)
:name "Kahan p13 Example 2"
:precision binary64
(/ (+ 1.0 (* (- 2.0 (/ (/ 2.0 t) (+ 1.0 (/ 1.0 t)))) (- 2.0 (/ (/ 2.0 t) (+ 1.0 (/ 1.0 t)))))) (+ 2.0 (* (- 2.0 (/ (/ 2.0 t) (+ 1.0 (/ 1.0 t)))) (- 2.0 (/ (/ 2.0 t) (+ 1.0 (/ 1.0 t))))))))