
(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 9 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 (/ (/ 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}
Initial program 100.0%
Final simplification100.0%
(FPCore (t)
:precision binary64
(let* ((t_1 (* (* 2.0 t) (* 2.0 t))) (t_2 (/ (- 8.0 (/ 4.0 t)) (- -1.0 t))))
(if (or (<= t -0.52) (not (<= t 0.76)))
(/ (+ 5.0 t_2) (+ 6.0 t_2))
(/ (+ 1.0 t_1) (+ 2.0 t_1)))))
double code(double t) {
double t_1 = (2.0 * t) * (2.0 * t);
double t_2 = (8.0 - (4.0 / t)) / (-1.0 - t);
double tmp;
if ((t <= -0.52) || !(t <= 0.76)) {
tmp = (5.0 + t_2) / (6.0 + t_2);
} else {
tmp = (1.0 + t_1) / (2.0 + t_1);
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (2.0d0 * t) * (2.0d0 * t)
t_2 = (8.0d0 - (4.0d0 / t)) / ((-1.0d0) - t)
if ((t <= (-0.52d0)) .or. (.not. (t <= 0.76d0))) then
tmp = (5.0d0 + t_2) / (6.0d0 + t_2)
else
tmp = (1.0d0 + t_1) / (2.0d0 + t_1)
end if
code = tmp
end function
public static double code(double t) {
double t_1 = (2.0 * t) * (2.0 * t);
double t_2 = (8.0 - (4.0 / t)) / (-1.0 - t);
double tmp;
if ((t <= -0.52) || !(t <= 0.76)) {
tmp = (5.0 + t_2) / (6.0 + t_2);
} else {
tmp = (1.0 + t_1) / (2.0 + t_1);
}
return tmp;
}
def code(t): t_1 = (2.0 * t) * (2.0 * t) t_2 = (8.0 - (4.0 / t)) / (-1.0 - t) tmp = 0 if (t <= -0.52) or not (t <= 0.76): tmp = (5.0 + t_2) / (6.0 + t_2) else: tmp = (1.0 + t_1) / (2.0 + t_1) return tmp
function code(t) t_1 = Float64(Float64(2.0 * t) * Float64(2.0 * t)) t_2 = Float64(Float64(8.0 - Float64(4.0 / t)) / Float64(-1.0 - t)) tmp = 0.0 if ((t <= -0.52) || !(t <= 0.76)) tmp = Float64(Float64(5.0 + t_2) / Float64(6.0 + t_2)); else tmp = Float64(Float64(1.0 + t_1) / Float64(2.0 + t_1)); end return tmp end
function tmp_2 = code(t) t_1 = (2.0 * t) * (2.0 * t); t_2 = (8.0 - (4.0 / t)) / (-1.0 - t); tmp = 0.0; if ((t <= -0.52) || ~((t <= 0.76))) tmp = (5.0 + t_2) / (6.0 + t_2); else tmp = (1.0 + t_1) / (2.0 + t_1); end tmp_2 = tmp; end
code[t_] := Block[{t$95$1 = N[(N[(2.0 * t), $MachinePrecision] * N[(2.0 * t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(8.0 - N[(4.0 / t), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - t), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[t, -0.52], N[Not[LessEqual[t, 0.76]], $MachinePrecision]], N[(N[(5.0 + t$95$2), $MachinePrecision] / N[(6.0 + t$95$2), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 + t$95$1), $MachinePrecision] / N[(2.0 + t$95$1), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(2 \cdot t\right) \cdot \left(2 \cdot t\right)\\
t_2 := \frac{8 - \frac{4}{t}}{-1 - t}\\
\mathbf{if}\;t \leq -0.52 \lor \neg \left(t \leq 0.76\right):\\
\;\;\;\;\frac{5 + t_2}{6 + t_2}\\
\mathbf{else}:\\
\;\;\;\;\frac{1 + t_1}{2 + t_1}\\
\end{array}
\end{array}
if t < -0.52000000000000002 or 0.76000000000000001 < t Initial program 100.0%
Simplified100.0%
associate-*l/100.0%
frac-2neg100.0%
sub-neg100.0%
distribute-rgt-in100.0%
+-commutative100.0%
metadata-eval100.0%
metadata-eval100.0%
distribute-neg-in100.0%
metadata-eval100.0%
Applied egg-rr100.0%
+-commutative100.0%
distribute-neg-in100.0%
metadata-eval100.0%
associate-*l/100.0%
metadata-eval100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
unsub-neg100.0%
Simplified100.0%
associate-*l/100.0%
frac-2neg100.0%
sub-neg100.0%
distribute-rgt-in100.0%
+-commutative100.0%
metadata-eval100.0%
metadata-eval100.0%
distribute-neg-in100.0%
metadata-eval100.0%
Applied egg-rr100.0%
+-commutative100.0%
distribute-neg-in100.0%
metadata-eval100.0%
associate-*l/100.0%
metadata-eval100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
unsub-neg100.0%
Simplified100.0%
Taylor expanded in t around inf 98.6%
associate-*r/98.6%
metadata-eval98.6%
Simplified98.6%
Taylor expanded in t around inf 98.7%
associate-*r/98.6%
metadata-eval98.6%
Simplified98.7%
if -0.52000000000000002 < t < 0.76000000000000001Initial program 100.0%
Taylor expanded in t around 0 99.4%
Taylor expanded in t around 0 99.4%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.4%
Final simplification99.0%
(FPCore (t) :precision binary64 (let* ((t_1 (/ 2.0 (+ 1.0 t))) (t_2 (* t_1 (- t_1 4.0)))) (/ (+ 5.0 t_2) (+ t_2 6.0))))
double code(double t) {
double t_1 = 2.0 / (1.0 + t);
double t_2 = t_1 * (t_1 - 4.0);
return (5.0 + t_2) / (t_2 + 6.0);
}
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)
t_2 = t_1 * (t_1 - 4.0d0)
code = (5.0d0 + t_2) / (t_2 + 6.0d0)
end function
public static double code(double t) {
double t_1 = 2.0 / (1.0 + t);
double t_2 = t_1 * (t_1 - 4.0);
return (5.0 + t_2) / (t_2 + 6.0);
}
def code(t): t_1 = 2.0 / (1.0 + t) t_2 = t_1 * (t_1 - 4.0) return (5.0 + t_2) / (t_2 + 6.0)
function code(t) t_1 = Float64(2.0 / Float64(1.0 + t)) t_2 = Float64(t_1 * Float64(t_1 - 4.0)) return Float64(Float64(5.0 + t_2) / Float64(t_2 + 6.0)) end
function tmp = code(t) t_1 = 2.0 / (1.0 + t); t_2 = t_1 * (t_1 - 4.0); tmp = (5.0 + t_2) / (t_2 + 6.0); end
code[t_] := Block[{t$95$1 = N[(2.0 / N[(1.0 + t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * N[(t$95$1 - 4.0), $MachinePrecision]), $MachinePrecision]}, N[(N[(5.0 + t$95$2), $MachinePrecision] / N[(t$95$2 + 6.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{2}{1 + t}\\
t_2 := t_1 \cdot \left(t_1 - 4\right)\\
\frac{5 + t_2}{t_2 + 6}
\end{array}
\end{array}
Initial program 100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (t)
:precision binary64
(let* ((t_1 (* (* 2.0 t) (* 2.0 t))))
(if (or (<= t -0.59) (not (<= t 0.68)))
(- 0.8333333333333334 (/ 0.2222222222222222 t))
(/ (+ 1.0 t_1) (+ 2.0 t_1)))))
double code(double t) {
double t_1 = (2.0 * t) * (2.0 * t);
double tmp;
if ((t <= -0.59) || !(t <= 0.68)) {
tmp = 0.8333333333333334 - (0.2222222222222222 / t);
} else {
tmp = (1.0 + t_1) / (2.0 + t_1);
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (2.0d0 * t) * (2.0d0 * t)
if ((t <= (-0.59d0)) .or. (.not. (t <= 0.68d0))) then
tmp = 0.8333333333333334d0 - (0.2222222222222222d0 / t)
else
tmp = (1.0d0 + t_1) / (2.0d0 + t_1)
end if
code = tmp
end function
public static double code(double t) {
double t_1 = (2.0 * t) * (2.0 * t);
double tmp;
if ((t <= -0.59) || !(t <= 0.68)) {
tmp = 0.8333333333333334 - (0.2222222222222222 / t);
} else {
tmp = (1.0 + t_1) / (2.0 + t_1);
}
return tmp;
}
def code(t): t_1 = (2.0 * t) * (2.0 * t) tmp = 0 if (t <= -0.59) or not (t <= 0.68): tmp = 0.8333333333333334 - (0.2222222222222222 / t) else: tmp = (1.0 + t_1) / (2.0 + t_1) return tmp
function code(t) t_1 = Float64(Float64(2.0 * t) * Float64(2.0 * t)) tmp = 0.0 if ((t <= -0.59) || !(t <= 0.68)) tmp = Float64(0.8333333333333334 - Float64(0.2222222222222222 / t)); else tmp = Float64(Float64(1.0 + t_1) / Float64(2.0 + t_1)); end return tmp end
function tmp_2 = code(t) t_1 = (2.0 * t) * (2.0 * t); tmp = 0.0; if ((t <= -0.59) || ~((t <= 0.68))) tmp = 0.8333333333333334 - (0.2222222222222222 / t); else tmp = (1.0 + t_1) / (2.0 + t_1); end tmp_2 = tmp; end
code[t_] := Block[{t$95$1 = N[(N[(2.0 * t), $MachinePrecision] * N[(2.0 * t), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[t, -0.59], N[Not[LessEqual[t, 0.68]], $MachinePrecision]], N[(0.8333333333333334 - N[(0.2222222222222222 / t), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 + t$95$1), $MachinePrecision] / N[(2.0 + t$95$1), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(2 \cdot t\right) \cdot \left(2 \cdot t\right)\\
\mathbf{if}\;t \leq -0.59 \lor \neg \left(t \leq 0.68\right):\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222}{t}\\
\mathbf{else}:\\
\;\;\;\;\frac{1 + t_1}{2 + t_1}\\
\end{array}
\end{array}
if t < -0.589999999999999969 or 0.680000000000000049 < t Initial program 100.0%
Simplified100.0%
Taylor expanded in t around inf 97.9%
Taylor expanded in t around inf 98.2%
associate-*r/98.2%
metadata-eval98.2%
Simplified98.2%
if -0.589999999999999969 < t < 0.680000000000000049Initial program 100.0%
Taylor expanded in t around 0 99.4%
Taylor expanded in t around 0 99.4%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.4%
Final simplification98.8%
(FPCore (t) :precision binary64 (let* ((t_1 (/ (+ 8.0 (/ -4.0 (+ 1.0 t))) (- -1.0 t)))) (/ (+ 5.0 t_1) (+ 6.0 t_1))))
double code(double t) {
double t_1 = (8.0 + (-4.0 / (1.0 + t))) / (-1.0 - t);
return (5.0 + t_1) / (6.0 + t_1);
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
t_1 = (8.0d0 + ((-4.0d0) / (1.0d0 + t))) / ((-1.0d0) - t)
code = (5.0d0 + t_1) / (6.0d0 + t_1)
end function
public static double code(double t) {
double t_1 = (8.0 + (-4.0 / (1.0 + t))) / (-1.0 - t);
return (5.0 + t_1) / (6.0 + t_1);
}
def code(t): t_1 = (8.0 + (-4.0 / (1.0 + t))) / (-1.0 - t) return (5.0 + t_1) / (6.0 + t_1)
function code(t) t_1 = Float64(Float64(8.0 + Float64(-4.0 / Float64(1.0 + t))) / Float64(-1.0 - t)) return Float64(Float64(5.0 + t_1) / Float64(6.0 + t_1)) end
function tmp = code(t) t_1 = (8.0 + (-4.0 / (1.0 + t))) / (-1.0 - t); tmp = (5.0 + t_1) / (6.0 + t_1); end
code[t_] := Block[{t$95$1 = N[(N[(8.0 + N[(-4.0 / N[(1.0 + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(-1.0 - t), $MachinePrecision]), $MachinePrecision]}, N[(N[(5.0 + t$95$1), $MachinePrecision] / N[(6.0 + t$95$1), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{8 + \frac{-4}{1 + t}}{-1 - t}\\
\frac{5 + t_1}{6 + t_1}
\end{array}
\end{array}
Initial program 100.0%
Simplified100.0%
associate-*l/100.0%
frac-2neg100.0%
sub-neg100.0%
distribute-rgt-in100.0%
+-commutative100.0%
metadata-eval100.0%
metadata-eval100.0%
distribute-neg-in100.0%
metadata-eval100.0%
Applied egg-rr100.0%
+-commutative100.0%
distribute-neg-in100.0%
metadata-eval100.0%
associate-*l/100.0%
metadata-eval100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
unsub-neg100.0%
Simplified100.0%
associate-*l/100.0%
frac-2neg100.0%
sub-neg100.0%
distribute-rgt-in100.0%
+-commutative100.0%
metadata-eval100.0%
metadata-eval100.0%
distribute-neg-in100.0%
metadata-eval100.0%
Applied egg-rr100.0%
+-commutative100.0%
distribute-neg-in100.0%
metadata-eval100.0%
associate-*l/100.0%
metadata-eval100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
unsub-neg100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (t) :precision binary64 (let* ((t_1 (/ 2.0 (+ 1.0 t)))) (/ (+ 5.0 (* t_1 (- t_1 4.0))) (+ 6.0 (* t_1 -2.0)))))
double code(double t) {
double t_1 = 2.0 / (1.0 + t);
return (5.0 + (t_1 * (t_1 - 4.0))) / (6.0 + (t_1 * -2.0));
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
t_1 = 2.0d0 / (1.0d0 + t)
code = (5.0d0 + (t_1 * (t_1 - 4.0d0))) / (6.0d0 + (t_1 * (-2.0d0)))
end function
public static double code(double t) {
double t_1 = 2.0 / (1.0 + t);
return (5.0 + (t_1 * (t_1 - 4.0))) / (6.0 + (t_1 * -2.0));
}
def code(t): t_1 = 2.0 / (1.0 + t) return (5.0 + (t_1 * (t_1 - 4.0))) / (6.0 + (t_1 * -2.0))
function code(t) t_1 = Float64(2.0 / Float64(1.0 + t)) return Float64(Float64(5.0 + Float64(t_1 * Float64(t_1 - 4.0))) / Float64(6.0 + Float64(t_1 * -2.0))) end
function tmp = code(t) t_1 = 2.0 / (1.0 + t); tmp = (5.0 + (t_1 * (t_1 - 4.0))) / (6.0 + (t_1 * -2.0)); end
code[t_] := Block[{t$95$1 = N[(2.0 / N[(1.0 + t), $MachinePrecision]), $MachinePrecision]}, N[(N[(5.0 + N[(t$95$1 * N[(t$95$1 - 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(6.0 + N[(t$95$1 * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{2}{1 + t}\\
\frac{5 + t_1 \cdot \left(t_1 - 4\right)}{6 + t_1 \cdot -2}
\end{array}
\end{array}
Initial program 100.0%
Simplified100.0%
Taylor expanded in t around 0 97.1%
Final simplification97.1%
(FPCore (t) :precision binary64 (if (or (<= t -0.48) (not (<= t 0.65))) (- 0.8333333333333334 (/ 0.2222222222222222 t)) 0.5))
double code(double t) {
double tmp;
if ((t <= -0.48) || !(t <= 0.65)) {
tmp = 0.8333333333333334 - (0.2222222222222222 / t);
} else {
tmp = 0.5;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if ((t <= (-0.48d0)) .or. (.not. (t <= 0.65d0))) then
tmp = 0.8333333333333334d0 - (0.2222222222222222d0 / t)
else
tmp = 0.5d0
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if ((t <= -0.48) || !(t <= 0.65)) {
tmp = 0.8333333333333334 - (0.2222222222222222 / t);
} else {
tmp = 0.5;
}
return tmp;
}
def code(t): tmp = 0 if (t <= -0.48) or not (t <= 0.65): tmp = 0.8333333333333334 - (0.2222222222222222 / t) else: tmp = 0.5 return tmp
function code(t) tmp = 0.0 if ((t <= -0.48) || !(t <= 0.65)) tmp = Float64(0.8333333333333334 - Float64(0.2222222222222222 / t)); else tmp = 0.5; end return tmp end
function tmp_2 = code(t) tmp = 0.0; if ((t <= -0.48) || ~((t <= 0.65))) tmp = 0.8333333333333334 - (0.2222222222222222 / t); else tmp = 0.5; end tmp_2 = tmp; end
code[t_] := If[Or[LessEqual[t, -0.48], N[Not[LessEqual[t, 0.65]], $MachinePrecision]], N[(0.8333333333333334 - N[(0.2222222222222222 / t), $MachinePrecision]), $MachinePrecision], 0.5]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.48 \lor \neg \left(t \leq 0.65\right):\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222}{t}\\
\mathbf{else}:\\
\;\;\;\;0.5\\
\end{array}
\end{array}
if t < -0.47999999999999998 or 0.650000000000000022 < t Initial program 100.0%
Simplified100.0%
Taylor expanded in t around inf 97.9%
Taylor expanded in t around inf 98.2%
associate-*r/98.2%
metadata-eval98.2%
Simplified98.2%
if -0.47999999999999998 < t < 0.650000000000000022Initial program 100.0%
Simplified100.0%
Taylor expanded in t around 0 97.5%
Taylor expanded in t around 0 99.1%
Final simplification98.7%
(FPCore (t) :precision binary64 (if (<= t -0.335) 0.8333333333333334 (if (<= t 1.0) 0.5 0.8333333333333334)))
double code(double t) {
double tmp;
if (t <= -0.335) {
tmp = 0.8333333333333334;
} else if (t <= 1.0) {
tmp = 0.5;
} else {
tmp = 0.8333333333333334;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-0.335d0)) then
tmp = 0.8333333333333334d0
else if (t <= 1.0d0) then
tmp = 0.5d0
else
tmp = 0.8333333333333334d0
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if (t <= -0.335) {
tmp = 0.8333333333333334;
} else if (t <= 1.0) {
tmp = 0.5;
} else {
tmp = 0.8333333333333334;
}
return tmp;
}
def code(t): tmp = 0 if t <= -0.335: tmp = 0.8333333333333334 elif t <= 1.0: tmp = 0.5 else: tmp = 0.8333333333333334 return tmp
function code(t) tmp = 0.0 if (t <= -0.335) tmp = 0.8333333333333334; elseif (t <= 1.0) tmp = 0.5; else tmp = 0.8333333333333334; end return tmp end
function tmp_2 = code(t) tmp = 0.0; if (t <= -0.335) tmp = 0.8333333333333334; elseif (t <= 1.0) tmp = 0.5; else tmp = 0.8333333333333334; end tmp_2 = tmp; end
code[t_] := If[LessEqual[t, -0.335], 0.8333333333333334, If[LessEqual[t, 1.0], 0.5, 0.8333333333333334]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.335:\\
\;\;\;\;0.8333333333333334\\
\mathbf{elif}\;t \leq 1:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;0.8333333333333334\\
\end{array}
\end{array}
if t < -0.33500000000000002 or 1 < t Initial program 100.0%
Simplified100.0%
Taylor expanded in t around inf 97.9%
Taylor expanded in t around inf 96.9%
if -0.33500000000000002 < t < 1Initial program 100.0%
Simplified100.0%
Taylor expanded in t around 0 97.5%
Taylor expanded in t around 0 99.1%
Final simplification98.0%
(FPCore (t) :precision binary64 0.5)
double code(double t) {
return 0.5;
}
real(8) function code(t)
real(8), intent (in) :: t
code = 0.5d0
end function
public static double code(double t) {
return 0.5;
}
def code(t): return 0.5
function code(t) return 0.5 end
function tmp = code(t) tmp = 0.5; end
code[t_] := 0.5
\begin{array}{l}
\\
0.5
\end{array}
Initial program 100.0%
Simplified100.0%
Taylor expanded in t around 0 97.1%
Taylor expanded in t around 0 58.1%
Final simplification58.1%
herbie shell --seed 2024021
(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))))))))