
(FPCore (t) :precision binary64 (let* ((t_1 (- 2.0 (/ (/ 2.0 t) (+ 1.0 (/ 1.0 t)))))) (- 1.0 (/ 1.0 (+ 2.0 (* t_1 t_1))))))
double code(double t) {
double t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t)));
return 1.0 - (1.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 / t) / (1.0d0 + (1.0d0 / t)))
code = 1.0d0 - (1.0d0 / (2.0d0 + (t_1 * t_1)))
end function
public static double code(double t) {
double t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t)));
return 1.0 - (1.0 / (2.0 + (t_1 * t_1)));
}
def code(t): t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t))) return 1.0 - (1.0 / (2.0 + (t_1 * t_1)))
function code(t) t_1 = Float64(2.0 - Float64(Float64(2.0 / t) / Float64(1.0 + Float64(1.0 / t)))) return Float64(1.0 - Float64(1.0 / Float64(2.0 + Float64(t_1 * t_1)))) end
function tmp = code(t) t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t))); tmp = 1.0 - (1.0 / (2.0 + (t_1 * t_1))); 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]}, N[(1.0 - N[(1.0 / N[(2.0 + N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 2 - \frac{\frac{2}{t}}{1 + \frac{1}{t}}\\
1 - \frac{1}{2 + t\_1 \cdot t\_1}
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 12 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (t) :precision binary64 (let* ((t_1 (- 2.0 (/ (/ 2.0 t) (+ 1.0 (/ 1.0 t)))))) (- 1.0 (/ 1.0 (+ 2.0 (* t_1 t_1))))))
double code(double t) {
double t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t)));
return 1.0 - (1.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 / t) / (1.0d0 + (1.0d0 / t)))
code = 1.0d0 - (1.0d0 / (2.0d0 + (t_1 * t_1)))
end function
public static double code(double t) {
double t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t)));
return 1.0 - (1.0 / (2.0 + (t_1 * t_1)));
}
def code(t): t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t))) return 1.0 - (1.0 / (2.0 + (t_1 * t_1)))
function code(t) t_1 = Float64(2.0 - Float64(Float64(2.0 / t) / Float64(1.0 + Float64(1.0 / t)))) return Float64(1.0 - Float64(1.0 / Float64(2.0 + Float64(t_1 * t_1)))) end
function tmp = code(t) t_1 = 2.0 - ((2.0 / t) / (1.0 + (1.0 / t))); tmp = 1.0 - (1.0 / (2.0 + (t_1 * t_1))); 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]}, N[(1.0 - N[(1.0 / N[(2.0 + N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 2 - \frac{\frac{2}{t}}{1 + \frac{1}{t}}\\
1 - \frac{1}{2 + t\_1 \cdot t\_1}
\end{array}
\end{array}
(FPCore (t) :precision binary64 (let* ((t_1 (+ 2.0 (/ -2.0 (+ 1.0 t))))) (+ 1.0 (/ -1.0 (+ 2.0 (* t_1 t_1))))))
double code(double t) {
double t_1 = 2.0 + (-2.0 / (1.0 + t));
return 1.0 + (-1.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 + ((-1.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 + (-1.0 / (2.0 + (t_1 * t_1)));
}
def code(t): t_1 = 2.0 + (-2.0 / (1.0 + t)) return 1.0 + (-1.0 / (2.0 + (t_1 * t_1)))
function code(t) t_1 = Float64(2.0 + Float64(-2.0 / Float64(1.0 + t))) return Float64(1.0 + Float64(-1.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 + (-1.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[(1.0 + N[(-1.0 / N[(2.0 + N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 2 + \frac{-2}{1 + t}\\
1 + \frac{-1}{2 + t\_1 \cdot t\_1}
\end{array}
\end{array}
Initial program 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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (t)
:precision binary64
(let* ((t_1
(/
(-
0.2222222222222222
(/ (+ 0.037037037037037035 (/ 0.04938271604938271 t)) t))
t)))
(if (<= t -0.58)
(- 1.0 (+ 0.16666666666666666 t_1))
(if (<= t 0.84)
(+ 1.0 (/ -1.0 (+ 2.0 (* (* t (+ 2.0 (* -2.0 t))) (* 2.0 t)))))
(- 0.8333333333333334 t_1)))))
double code(double t) {
double t_1 = (0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t;
double tmp;
if (t <= -0.58) {
tmp = 1.0 - (0.16666666666666666 + t_1);
} else if (t <= 0.84) {
tmp = 1.0 + (-1.0 / (2.0 + ((t * (2.0 + (-2.0 * t))) * (2.0 * t))));
} else {
tmp = 0.8333333333333334 - t_1;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (0.2222222222222222d0 - ((0.037037037037037035d0 + (0.04938271604938271d0 / t)) / t)) / t
if (t <= (-0.58d0)) then
tmp = 1.0d0 - (0.16666666666666666d0 + t_1)
else if (t <= 0.84d0) then
tmp = 1.0d0 + ((-1.0d0) / (2.0d0 + ((t * (2.0d0 + ((-2.0d0) * t))) * (2.0d0 * t))))
else
tmp = 0.8333333333333334d0 - t_1
end if
code = tmp
end function
public static double code(double t) {
double t_1 = (0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t;
double tmp;
if (t <= -0.58) {
tmp = 1.0 - (0.16666666666666666 + t_1);
} else if (t <= 0.84) {
tmp = 1.0 + (-1.0 / (2.0 + ((t * (2.0 + (-2.0 * t))) * (2.0 * t))));
} else {
tmp = 0.8333333333333334 - t_1;
}
return tmp;
}
def code(t): t_1 = (0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t tmp = 0 if t <= -0.58: tmp = 1.0 - (0.16666666666666666 + t_1) elif t <= 0.84: tmp = 1.0 + (-1.0 / (2.0 + ((t * (2.0 + (-2.0 * t))) * (2.0 * t)))) else: tmp = 0.8333333333333334 - t_1 return tmp
function code(t) t_1 = Float64(Float64(0.2222222222222222 - Float64(Float64(0.037037037037037035 + Float64(0.04938271604938271 / t)) / t)) / t) tmp = 0.0 if (t <= -0.58) tmp = Float64(1.0 - Float64(0.16666666666666666 + t_1)); elseif (t <= 0.84) tmp = Float64(1.0 + Float64(-1.0 / Float64(2.0 + Float64(Float64(t * Float64(2.0 + Float64(-2.0 * t))) * Float64(2.0 * t))))); else tmp = Float64(0.8333333333333334 - t_1); end return tmp end
function tmp_2 = code(t) t_1 = (0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t; tmp = 0.0; if (t <= -0.58) tmp = 1.0 - (0.16666666666666666 + t_1); elseif (t <= 0.84) tmp = 1.0 + (-1.0 / (2.0 + ((t * (2.0 + (-2.0 * t))) * (2.0 * t)))); else tmp = 0.8333333333333334 - t_1; end tmp_2 = tmp; end
code[t_] := Block[{t$95$1 = N[(N[(0.2222222222222222 - N[(N[(0.037037037037037035 + N[(0.04938271604938271 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]}, If[LessEqual[t, -0.58], N[(1.0 - N[(0.16666666666666666 + t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.84], N[(1.0 + N[(-1.0 / N[(2.0 + N[(N[(t * N[(2.0 + N[(-2.0 * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(2.0 * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.8333333333333334 - t$95$1), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.2222222222222222 - \frac{0.037037037037037035 + \frac{0.04938271604938271}{t}}{t}}{t}\\
\mathbf{if}\;t \leq -0.58:\\
\;\;\;\;1 - \left(0.16666666666666666 + t\_1\right)\\
\mathbf{elif}\;t \leq 0.84:\\
\;\;\;\;1 + \frac{-1}{2 + \left(t \cdot \left(2 + -2 \cdot t\right)\right) \cdot \left(2 \cdot t\right)}\\
\mathbf{else}:\\
\;\;\;\;0.8333333333333334 - t\_1\\
\end{array}
\end{array}
if t < -0.57999999999999996Initial program 100.0%
Taylor expanded in t around inf 99.5%
+-commutative99.5%
mul-1-neg99.5%
unsub-neg99.5%
associate-*r/99.5%
metadata-eval99.5%
unpow299.5%
associate-/r*99.5%
div-sub99.5%
associate-*r/99.5%
metadata-eval99.5%
Simplified99.5%
if -0.57999999999999996 < t < 0.839999999999999969Initial program 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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.5%
*-commutative99.5%
Simplified99.5%
if 0.839999999999999969 < t Initial program 100.0%
Taylor expanded in t around -inf 99.8%
mul-1-neg99.8%
unsub-neg99.8%
mul-1-neg99.8%
unsub-neg99.8%
associate-*r/99.8%
metadata-eval99.8%
Simplified99.8%
Final simplification99.6%
(FPCore (t)
:precision binary64
(let* ((t_1
(/
(-
0.2222222222222222
(/ (+ 0.037037037037037035 (/ 0.04938271604938271 t)) t))
t)))
(if (<= t -1.12)
(- 1.0 (+ 0.16666666666666666 t_1))
(if (<= t 0.8)
(+ 1.0 (/ -1.0 (+ 2.0 (* t (+ 4.0 (/ -4.0 (+ 1.0 t)))))))
(- 0.8333333333333334 t_1)))))
double code(double t) {
double t_1 = (0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t;
double tmp;
if (t <= -1.12) {
tmp = 1.0 - (0.16666666666666666 + t_1);
} else if (t <= 0.8) {
tmp = 1.0 + (-1.0 / (2.0 + (t * (4.0 + (-4.0 / (1.0 + t))))));
} else {
tmp = 0.8333333333333334 - t_1;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (0.2222222222222222d0 - ((0.037037037037037035d0 + (0.04938271604938271d0 / t)) / t)) / t
if (t <= (-1.12d0)) then
tmp = 1.0d0 - (0.16666666666666666d0 + t_1)
else if (t <= 0.8d0) then
tmp = 1.0d0 + ((-1.0d0) / (2.0d0 + (t * (4.0d0 + ((-4.0d0) / (1.0d0 + t))))))
else
tmp = 0.8333333333333334d0 - t_1
end if
code = tmp
end function
public static double code(double t) {
double t_1 = (0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t;
double tmp;
if (t <= -1.12) {
tmp = 1.0 - (0.16666666666666666 + t_1);
} else if (t <= 0.8) {
tmp = 1.0 + (-1.0 / (2.0 + (t * (4.0 + (-4.0 / (1.0 + t))))));
} else {
tmp = 0.8333333333333334 - t_1;
}
return tmp;
}
def code(t): t_1 = (0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t tmp = 0 if t <= -1.12: tmp = 1.0 - (0.16666666666666666 + t_1) elif t <= 0.8: tmp = 1.0 + (-1.0 / (2.0 + (t * (4.0 + (-4.0 / (1.0 + t)))))) else: tmp = 0.8333333333333334 - t_1 return tmp
function code(t) t_1 = Float64(Float64(0.2222222222222222 - Float64(Float64(0.037037037037037035 + Float64(0.04938271604938271 / t)) / t)) / t) tmp = 0.0 if (t <= -1.12) tmp = Float64(1.0 - Float64(0.16666666666666666 + t_1)); elseif (t <= 0.8) tmp = Float64(1.0 + Float64(-1.0 / Float64(2.0 + Float64(t * Float64(4.0 + Float64(-4.0 / Float64(1.0 + t))))))); else tmp = Float64(0.8333333333333334 - t_1); end return tmp end
function tmp_2 = code(t) t_1 = (0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t; tmp = 0.0; if (t <= -1.12) tmp = 1.0 - (0.16666666666666666 + t_1); elseif (t <= 0.8) tmp = 1.0 + (-1.0 / (2.0 + (t * (4.0 + (-4.0 / (1.0 + t)))))); else tmp = 0.8333333333333334 - t_1; end tmp_2 = tmp; end
code[t_] := Block[{t$95$1 = N[(N[(0.2222222222222222 - N[(N[(0.037037037037037035 + N[(0.04938271604938271 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]}, If[LessEqual[t, -1.12], N[(1.0 - N[(0.16666666666666666 + t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.8], N[(1.0 + N[(-1.0 / N[(2.0 + N[(t * N[(4.0 + N[(-4.0 / N[(1.0 + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.8333333333333334 - t$95$1), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.2222222222222222 - \frac{0.037037037037037035 + \frac{0.04938271604938271}{t}}{t}}{t}\\
\mathbf{if}\;t \leq -1.12:\\
\;\;\;\;1 - \left(0.16666666666666666 + t\_1\right)\\
\mathbf{elif}\;t \leq 0.8:\\
\;\;\;\;1 + \frac{-1}{2 + t \cdot \left(4 + \frac{-4}{1 + t}\right)}\\
\mathbf{else}:\\
\;\;\;\;0.8333333333333334 - t\_1\\
\end{array}
\end{array}
if t < -1.1200000000000001Initial program 100.0%
Taylor expanded in t around inf 99.5%
+-commutative99.5%
mul-1-neg99.5%
unsub-neg99.5%
associate-*r/99.5%
metadata-eval99.5%
unpow299.5%
associate-/r*99.5%
div-sub99.5%
associate-*r/99.5%
metadata-eval99.5%
Simplified99.5%
if -1.1200000000000001 < t < 0.80000000000000004Initial program 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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
*-commutative99.5%
distribute-lft-in99.5%
Applied egg-rr99.5%
distribute-lft-out99.5%
distribute-rgt-out99.5%
associate-*r*99.5%
metadata-eval99.5%
associate-*r*99.5%
distribute-rgt-out99.5%
associate-*l/99.5%
metadata-eval99.5%
Simplified99.5%
if 0.80000000000000004 < t Initial program 100.0%
Taylor expanded in t around -inf 99.8%
mul-1-neg99.8%
unsub-neg99.8%
mul-1-neg99.8%
unsub-neg99.8%
associate-*r/99.8%
metadata-eval99.8%
Simplified99.8%
Final simplification99.6%
(FPCore (t)
:precision binary64
(if (or (<= t -0.34) (not (<= t 0.69)))
(-
0.8333333333333334
(/
(-
0.2222222222222222
(/ (+ 0.037037037037037035 (/ 0.04938271604938271 t)) t))
t))
(+ 1.0 (/ -1.0 (+ 2.0 (* (* 2.0 t) (* 2.0 t)))))))
double code(double t) {
double tmp;
if ((t <= -0.34) || !(t <= 0.69)) {
tmp = 0.8333333333333334 - ((0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t);
} else {
tmp = 1.0 + (-1.0 / (2.0 + ((2.0 * t) * (2.0 * t))));
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if ((t <= (-0.34d0)) .or. (.not. (t <= 0.69d0))) then
tmp = 0.8333333333333334d0 - ((0.2222222222222222d0 - ((0.037037037037037035d0 + (0.04938271604938271d0 / t)) / t)) / t)
else
tmp = 1.0d0 + ((-1.0d0) / (2.0d0 + ((2.0d0 * t) * (2.0d0 * t))))
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if ((t <= -0.34) || !(t <= 0.69)) {
tmp = 0.8333333333333334 - ((0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t);
} else {
tmp = 1.0 + (-1.0 / (2.0 + ((2.0 * t) * (2.0 * t))));
}
return tmp;
}
def code(t): tmp = 0 if (t <= -0.34) or not (t <= 0.69): tmp = 0.8333333333333334 - ((0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t) else: tmp = 1.0 + (-1.0 / (2.0 + ((2.0 * t) * (2.0 * t)))) return tmp
function code(t) tmp = 0.0 if ((t <= -0.34) || !(t <= 0.69)) tmp = Float64(0.8333333333333334 - Float64(Float64(0.2222222222222222 - Float64(Float64(0.037037037037037035 + Float64(0.04938271604938271 / t)) / t)) / t)); else tmp = Float64(1.0 + Float64(-1.0 / Float64(2.0 + Float64(Float64(2.0 * t) * Float64(2.0 * t))))); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if ((t <= -0.34) || ~((t <= 0.69))) tmp = 0.8333333333333334 - ((0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t); else tmp = 1.0 + (-1.0 / (2.0 + ((2.0 * t) * (2.0 * t)))); end tmp_2 = tmp; end
code[t_] := If[Or[LessEqual[t, -0.34], N[Not[LessEqual[t, 0.69]], $MachinePrecision]], N[(0.8333333333333334 - N[(N[(0.2222222222222222 - N[(N[(0.037037037037037035 + N[(0.04938271604938271 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-1.0 / N[(2.0 + N[(N[(2.0 * t), $MachinePrecision] * N[(2.0 * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.34 \lor \neg \left(t \leq 0.69\right):\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222 - \frac{0.037037037037037035 + \frac{0.04938271604938271}{t}}{t}}{t}\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{-1}{2 + \left(2 \cdot t\right) \cdot \left(2 \cdot t\right)}\\
\end{array}
\end{array}
if t < -0.340000000000000024 or 0.68999999999999995 < t Initial program 100.0%
Taylor expanded in t around -inf 99.6%
mul-1-neg99.6%
unsub-neg99.6%
mul-1-neg99.6%
unsub-neg99.6%
associate-*r/99.6%
metadata-eval99.6%
Simplified99.6%
if -0.340000000000000024 < t < 0.68999999999999995Initial program 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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.4%
Final simplification99.6%
(FPCore (t)
:precision binary64
(if (or (<= t -0.33) (not (<= t 0.68)))
(-
0.8333333333333334
(/
(-
0.2222222222222222
(/ (+ 0.037037037037037035 (/ 0.04938271604938271 t)) t))
t))
0.5))
double code(double t) {
double tmp;
if ((t <= -0.33) || !(t <= 0.68)) {
tmp = 0.8333333333333334 - ((0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t);
} else {
tmp = 0.5;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if ((t <= (-0.33d0)) .or. (.not. (t <= 0.68d0))) then
tmp = 0.8333333333333334d0 - ((0.2222222222222222d0 - ((0.037037037037037035d0 + (0.04938271604938271d0 / t)) / t)) / t)
else
tmp = 0.5d0
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if ((t <= -0.33) || !(t <= 0.68)) {
tmp = 0.8333333333333334 - ((0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t);
} else {
tmp = 0.5;
}
return tmp;
}
def code(t): tmp = 0 if (t <= -0.33) or not (t <= 0.68): tmp = 0.8333333333333334 - ((0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t) else: tmp = 0.5 return tmp
function code(t) tmp = 0.0 if ((t <= -0.33) || !(t <= 0.68)) tmp = Float64(0.8333333333333334 - Float64(Float64(0.2222222222222222 - Float64(Float64(0.037037037037037035 + Float64(0.04938271604938271 / t)) / t)) / t)); else tmp = 0.5; end return tmp end
function tmp_2 = code(t) tmp = 0.0; if ((t <= -0.33) || ~((t <= 0.68))) tmp = 0.8333333333333334 - ((0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t); else tmp = 0.5; end tmp_2 = tmp; end
code[t_] := If[Or[LessEqual[t, -0.33], N[Not[LessEqual[t, 0.68]], $MachinePrecision]], N[(0.8333333333333334 - N[(N[(0.2222222222222222 - N[(N[(0.037037037037037035 + N[(0.04938271604938271 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], 0.5]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.33 \lor \neg \left(t \leq 0.68\right):\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222 - \frac{0.037037037037037035 + \frac{0.04938271604938271}{t}}{t}}{t}\\
\mathbf{else}:\\
\;\;\;\;0.5\\
\end{array}
\end{array}
if t < -0.330000000000000016 or 0.680000000000000049 < t Initial program 100.0%
Taylor expanded in t around -inf 99.6%
mul-1-neg99.6%
unsub-neg99.6%
mul-1-neg99.6%
unsub-neg99.6%
associate-*r/99.6%
metadata-eval99.6%
Simplified99.6%
if -0.330000000000000016 < t < 0.680000000000000049Initial program 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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.1%
Final simplification99.4%
(FPCore (t)
:precision binary64
(let* ((t_1
(/
(-
0.2222222222222222
(/ (+ 0.037037037037037035 (/ 0.04938271604938271 t)) t))
t)))
(if (<= t -0.34)
(- 1.0 (+ 0.16666666666666666 t_1))
(if (<= t 0.69)
(+ 1.0 (/ -1.0 (+ 2.0 (* (* 2.0 t) (* 2.0 t)))))
(- 0.8333333333333334 t_1)))))
double code(double t) {
double t_1 = (0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t;
double tmp;
if (t <= -0.34) {
tmp = 1.0 - (0.16666666666666666 + t_1);
} else if (t <= 0.69) {
tmp = 1.0 + (-1.0 / (2.0 + ((2.0 * t) * (2.0 * t))));
} else {
tmp = 0.8333333333333334 - t_1;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (0.2222222222222222d0 - ((0.037037037037037035d0 + (0.04938271604938271d0 / t)) / t)) / t
if (t <= (-0.34d0)) then
tmp = 1.0d0 - (0.16666666666666666d0 + t_1)
else if (t <= 0.69d0) then
tmp = 1.0d0 + ((-1.0d0) / (2.0d0 + ((2.0d0 * t) * (2.0d0 * t))))
else
tmp = 0.8333333333333334d0 - t_1
end if
code = tmp
end function
public static double code(double t) {
double t_1 = (0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t;
double tmp;
if (t <= -0.34) {
tmp = 1.0 - (0.16666666666666666 + t_1);
} else if (t <= 0.69) {
tmp = 1.0 + (-1.0 / (2.0 + ((2.0 * t) * (2.0 * t))));
} else {
tmp = 0.8333333333333334 - t_1;
}
return tmp;
}
def code(t): t_1 = (0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t tmp = 0 if t <= -0.34: tmp = 1.0 - (0.16666666666666666 + t_1) elif t <= 0.69: tmp = 1.0 + (-1.0 / (2.0 + ((2.0 * t) * (2.0 * t)))) else: tmp = 0.8333333333333334 - t_1 return tmp
function code(t) t_1 = Float64(Float64(0.2222222222222222 - Float64(Float64(0.037037037037037035 + Float64(0.04938271604938271 / t)) / t)) / t) tmp = 0.0 if (t <= -0.34) tmp = Float64(1.0 - Float64(0.16666666666666666 + t_1)); elseif (t <= 0.69) tmp = Float64(1.0 + Float64(-1.0 / Float64(2.0 + Float64(Float64(2.0 * t) * Float64(2.0 * t))))); else tmp = Float64(0.8333333333333334 - t_1); end return tmp end
function tmp_2 = code(t) t_1 = (0.2222222222222222 - ((0.037037037037037035 + (0.04938271604938271 / t)) / t)) / t; tmp = 0.0; if (t <= -0.34) tmp = 1.0 - (0.16666666666666666 + t_1); elseif (t <= 0.69) tmp = 1.0 + (-1.0 / (2.0 + ((2.0 * t) * (2.0 * t)))); else tmp = 0.8333333333333334 - t_1; end tmp_2 = tmp; end
code[t_] := Block[{t$95$1 = N[(N[(0.2222222222222222 - N[(N[(0.037037037037037035 + N[(0.04938271604938271 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]}, If[LessEqual[t, -0.34], N[(1.0 - N[(0.16666666666666666 + t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.69], N[(1.0 + N[(-1.0 / N[(2.0 + N[(N[(2.0 * t), $MachinePrecision] * N[(2.0 * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.8333333333333334 - t$95$1), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.2222222222222222 - \frac{0.037037037037037035 + \frac{0.04938271604938271}{t}}{t}}{t}\\
\mathbf{if}\;t \leq -0.34:\\
\;\;\;\;1 - \left(0.16666666666666666 + t\_1\right)\\
\mathbf{elif}\;t \leq 0.69:\\
\;\;\;\;1 + \frac{-1}{2 + \left(2 \cdot t\right) \cdot \left(2 \cdot t\right)}\\
\mathbf{else}:\\
\;\;\;\;0.8333333333333334 - t\_1\\
\end{array}
\end{array}
if t < -0.340000000000000024Initial program 100.0%
Taylor expanded in t around inf 99.5%
+-commutative99.5%
mul-1-neg99.5%
unsub-neg99.5%
associate-*r/99.5%
metadata-eval99.5%
unpow299.5%
associate-/r*99.5%
div-sub99.5%
associate-*r/99.5%
metadata-eval99.5%
Simplified99.5%
if -0.340000000000000024 < t < 0.68999999999999995Initial program 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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.4%
if 0.68999999999999995 < t Initial program 100.0%
Taylor expanded in t around -inf 99.8%
mul-1-neg99.8%
unsub-neg99.8%
mul-1-neg99.8%
unsub-neg99.8%
associate-*r/99.8%
metadata-eval99.8%
Simplified99.8%
Final simplification99.6%
(FPCore (t)
:precision binary64
(if (or (<= t -0.52) (not (<= t 0.23)))
(-
0.8333333333333334
(/ (+ 0.2222222222222222 (/ -0.037037037037037035 t)) t))
0.5))
double code(double t) {
double tmp;
if ((t <= -0.52) || !(t <= 0.23)) {
tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t);
} else {
tmp = 0.5;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if ((t <= (-0.52d0)) .or. (.not. (t <= 0.23d0))) then
tmp = 0.8333333333333334d0 - ((0.2222222222222222d0 + ((-0.037037037037037035d0) / t)) / t)
else
tmp = 0.5d0
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if ((t <= -0.52) || !(t <= 0.23)) {
tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t);
} else {
tmp = 0.5;
}
return tmp;
}
def code(t): tmp = 0 if (t <= -0.52) or not (t <= 0.23): tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t) else: tmp = 0.5 return tmp
function code(t) tmp = 0.0 if ((t <= -0.52) || !(t <= 0.23)) tmp = Float64(0.8333333333333334 - Float64(Float64(0.2222222222222222 + Float64(-0.037037037037037035 / t)) / t)); else tmp = 0.5; end return tmp end
function tmp_2 = code(t) tmp = 0.0; if ((t <= -0.52) || ~((t <= 0.23))) tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t); else tmp = 0.5; end tmp_2 = tmp; end
code[t_] := If[Or[LessEqual[t, -0.52], N[Not[LessEqual[t, 0.23]], $MachinePrecision]], N[(0.8333333333333334 - N[(N[(0.2222222222222222 + N[(-0.037037037037037035 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], 0.5]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.52 \lor \neg \left(t \leq 0.23\right):\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222 + \frac{-0.037037037037037035}{t}}{t}\\
\mathbf{else}:\\
\;\;\;\;0.5\\
\end{array}
\end{array}
if t < -0.52000000000000002 or 0.23000000000000001 < t Initial program 100.0%
Taylor expanded in t around inf 99.5%
associate--l+99.5%
unpow299.5%
associate-/r*99.5%
metadata-eval99.5%
associate-*r/99.5%
associate-*r/99.5%
metadata-eval99.5%
div-sub99.5%
remove-double-neg99.5%
sub-neg99.5%
+-commutative99.5%
distribute-neg-in99.5%
metadata-eval99.5%
metadata-eval99.5%
sub-neg99.5%
distribute-neg-frac99.5%
unsub-neg99.5%
Simplified99.5%
if -0.52000000000000002 < t < 0.23000000000000001Initial program 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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.1%
Final simplification99.3%
(FPCore (t)
:precision binary64
(let* ((t_1 (/ (+ 0.2222222222222222 (/ -0.037037037037037035 t)) t)))
(if (<= t -0.52)
(- 1.0 (+ 0.16666666666666666 t_1))
(if (<= t 0.23) 0.5 (- 0.8333333333333334 t_1)))))
double code(double t) {
double t_1 = (0.2222222222222222 + (-0.037037037037037035 / t)) / t;
double tmp;
if (t <= -0.52) {
tmp = 1.0 - (0.16666666666666666 + t_1);
} else if (t <= 0.23) {
tmp = 0.5;
} else {
tmp = 0.8333333333333334 - t_1;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (0.2222222222222222d0 + ((-0.037037037037037035d0) / t)) / t
if (t <= (-0.52d0)) then
tmp = 1.0d0 - (0.16666666666666666d0 + t_1)
else if (t <= 0.23d0) then
tmp = 0.5d0
else
tmp = 0.8333333333333334d0 - t_1
end if
code = tmp
end function
public static double code(double t) {
double t_1 = (0.2222222222222222 + (-0.037037037037037035 / t)) / t;
double tmp;
if (t <= -0.52) {
tmp = 1.0 - (0.16666666666666666 + t_1);
} else if (t <= 0.23) {
tmp = 0.5;
} else {
tmp = 0.8333333333333334 - t_1;
}
return tmp;
}
def code(t): t_1 = (0.2222222222222222 + (-0.037037037037037035 / t)) / t tmp = 0 if t <= -0.52: tmp = 1.0 - (0.16666666666666666 + t_1) elif t <= 0.23: tmp = 0.5 else: tmp = 0.8333333333333334 - t_1 return tmp
function code(t) t_1 = Float64(Float64(0.2222222222222222 + Float64(-0.037037037037037035 / t)) / t) tmp = 0.0 if (t <= -0.52) tmp = Float64(1.0 - Float64(0.16666666666666666 + t_1)); elseif (t <= 0.23) tmp = 0.5; else tmp = Float64(0.8333333333333334 - t_1); end return tmp end
function tmp_2 = code(t) t_1 = (0.2222222222222222 + (-0.037037037037037035 / t)) / t; tmp = 0.0; if (t <= -0.52) tmp = 1.0 - (0.16666666666666666 + t_1); elseif (t <= 0.23) tmp = 0.5; else tmp = 0.8333333333333334 - t_1; end tmp_2 = tmp; end
code[t_] := Block[{t$95$1 = N[(N[(0.2222222222222222 + N[(-0.037037037037037035 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]}, If[LessEqual[t, -0.52], N[(1.0 - N[(0.16666666666666666 + t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.23], 0.5, N[(0.8333333333333334 - t$95$1), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.2222222222222222 + \frac{-0.037037037037037035}{t}}{t}\\
\mathbf{if}\;t \leq -0.52:\\
\;\;\;\;1 - \left(0.16666666666666666 + t\_1\right)\\
\mathbf{elif}\;t \leq 0.23:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;0.8333333333333334 - t\_1\\
\end{array}
\end{array}
if t < -0.52000000000000002Initial program 100.0%
Taylor expanded in t around inf 99.3%
associate--l+99.3%
unpow299.3%
associate-/r*99.3%
metadata-eval99.3%
associate-*r/99.3%
associate-*r/99.3%
metadata-eval99.3%
div-sub99.3%
sub-neg99.3%
associate-*r/99.3%
metadata-eval99.3%
distribute-neg-frac99.3%
metadata-eval99.3%
Simplified99.3%
if -0.52000000000000002 < t < 0.23000000000000001Initial program 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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.1%
if 0.23000000000000001 < t Initial program 100.0%
Taylor expanded in t around inf 99.6%
associate--l+99.6%
unpow299.6%
associate-/r*99.6%
metadata-eval99.6%
associate-*r/99.6%
associate-*r/99.6%
metadata-eval99.6%
div-sub99.6%
remove-double-neg99.6%
sub-neg99.6%
+-commutative99.6%
distribute-neg-in99.6%
metadata-eval99.6%
metadata-eval99.6%
sub-neg99.6%
distribute-neg-frac99.6%
unsub-neg99.6%
Simplified99.6%
(FPCore (t) :precision binary64 (if (or (<= t -0.48) (not (<= t 0.68))) (- 0.8333333333333334 (/ 0.2222222222222222 t)) 0.5))
double code(double t) {
double tmp;
if ((t <= -0.48) || !(t <= 0.68)) {
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.68d0))) 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.68)) {
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.68): tmp = 0.8333333333333334 - (0.2222222222222222 / t) else: tmp = 0.5 return tmp
function code(t) tmp = 0.0 if ((t <= -0.48) || !(t <= 0.68)) 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.68))) 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.68]], $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.68\right):\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222}{t}\\
\mathbf{else}:\\
\;\;\;\;0.5\\
\end{array}
\end{array}
if t < -0.47999999999999998 or 0.680000000000000049 < t Initial program 100.0%
Taylor expanded in t around inf 99.2%
associate-*r/99.2%
metadata-eval99.2%
Simplified99.2%
if -0.47999999999999998 < t < 0.680000000000000049Initial program 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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.1%
Final simplification99.2%
(FPCore (t) :precision binary64 (if (<= t -0.48) (- 1.0 (+ 0.16666666666666666 (/ 0.2222222222222222 t))) (if (<= t 0.68) 0.5 (- 0.8333333333333334 (/ 0.2222222222222222 t)))))
double code(double t) {
double tmp;
if (t <= -0.48) {
tmp = 1.0 - (0.16666666666666666 + (0.2222222222222222 / t));
} else if (t <= 0.68) {
tmp = 0.5;
} else {
tmp = 0.8333333333333334 - (0.2222222222222222 / t);
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-0.48d0)) then
tmp = 1.0d0 - (0.16666666666666666d0 + (0.2222222222222222d0 / t))
else if (t <= 0.68d0) then
tmp = 0.5d0
else
tmp = 0.8333333333333334d0 - (0.2222222222222222d0 / t)
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if (t <= -0.48) {
tmp = 1.0 - (0.16666666666666666 + (0.2222222222222222 / t));
} else if (t <= 0.68) {
tmp = 0.5;
} else {
tmp = 0.8333333333333334 - (0.2222222222222222 / t);
}
return tmp;
}
def code(t): tmp = 0 if t <= -0.48: tmp = 1.0 - (0.16666666666666666 + (0.2222222222222222 / t)) elif t <= 0.68: tmp = 0.5 else: tmp = 0.8333333333333334 - (0.2222222222222222 / t) return tmp
function code(t) tmp = 0.0 if (t <= -0.48) tmp = Float64(1.0 - Float64(0.16666666666666666 + Float64(0.2222222222222222 / t))); elseif (t <= 0.68) tmp = 0.5; else tmp = Float64(0.8333333333333334 - Float64(0.2222222222222222 / t)); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if (t <= -0.48) tmp = 1.0 - (0.16666666666666666 + (0.2222222222222222 / t)); elseif (t <= 0.68) tmp = 0.5; else tmp = 0.8333333333333334 - (0.2222222222222222 / t); end tmp_2 = tmp; end
code[t_] := If[LessEqual[t, -0.48], N[(1.0 - N[(0.16666666666666666 + N[(0.2222222222222222 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.68], 0.5, N[(0.8333333333333334 - N[(0.2222222222222222 / t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.48:\\
\;\;\;\;1 - \left(0.16666666666666666 + \frac{0.2222222222222222}{t}\right)\\
\mathbf{elif}\;t \leq 0.68:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222}{t}\\
\end{array}
\end{array}
if t < -0.47999999999999998Initial program 100.0%
Taylor expanded in t around inf 99.2%
associate-*r/99.2%
metadata-eval99.2%
Simplified99.2%
if -0.47999999999999998 < t < 0.680000000000000049Initial program 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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.1%
if 0.680000000000000049 < t Initial program 100.0%
Taylor expanded in t around inf 99.3%
associate-*r/99.3%
metadata-eval99.3%
Simplified99.3%
(FPCore (t) :precision binary64 (if (<= t -0.33) 0.8333333333333334 (if (<= t 1.0) 0.5 0.8333333333333334)))
double code(double t) {
double tmp;
if (t <= -0.33) {
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.33d0)) 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.33) {
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.33: 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.33) 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.33) 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.33], 0.8333333333333334, If[LessEqual[t, 1.0], 0.5, 0.8333333333333334]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.33:\\
\;\;\;\;0.8333333333333334\\
\mathbf{elif}\;t \leq 1:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;0.8333333333333334\\
\end{array}
\end{array}
if t < -0.330000000000000016 or 1 < t Initial program 100.0%
Taylor expanded in t around inf 98.2%
if -0.330000000000000016 < t < 1Initial program 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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.1%
(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%
sub-neg100.0%
distribute-neg-frac100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Applied egg-rr100.0%
associate-/r*100.0%
distribute-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.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-lft-in100.0%
*-rgt-identity100.0%
rgt-mult-inverse100.0%
Simplified100.0%
Taylor expanded in t around 0 56.4%
Taylor expanded in t around 0 55.3%
herbie shell --seed 2024106
(FPCore (t)
:name "Kahan p13 Example 3"
:precision binary64
(- 1.0 (/ 1.0 (+ 2.0 (* (- 2.0 (/ (/ 2.0 t) (+ 1.0 (/ 1.0 t)))) (- 2.0 (/ (/ 2.0 t) (+ 1.0 (/ 1.0 t)))))))))