
(FPCore (t) :precision binary64 (let* ((t_1 (/ (* 2.0 t) (+ 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 * t) / (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 * t) / (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 * t) / (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 * t) / (1.0 + t) t_2 = t_1 * t_1 return (1.0 + t_2) / (2.0 + t_2)
function code(t) t_1 = Float64(Float64(2.0 * t) / 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 * t) / (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[(N[(2.0 * t), $MachinePrecision] / N[(1.0 + t), $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 := \frac{2 \cdot t}{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 8 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (t) :precision binary64 (let* ((t_1 (/ (* 2.0 t) (+ 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 * t) / (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 * t) / (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 * t) / (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 * t) / (1.0 + t) t_2 = t_1 * t_1 return (1.0 + t_2) / (2.0 + t_2)
function code(t) t_1 = Float64(Float64(2.0 * t) / 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 * t) / (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[(N[(2.0 * t), $MachinePrecision] / N[(1.0 + t), $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 := \frac{2 \cdot t}{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 (/ t (+ t 1.0))) (t_2 (* 4.0 (* t_1 t_1)))) (/ (+ t_2 1.0) (+ t_2 2.0))))
double code(double t) {
double t_1 = t / (t + 1.0);
double t_2 = 4.0 * (t_1 * t_1);
return (t_2 + 1.0) / (t_2 + 2.0);
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: t_2
t_1 = t / (t + 1.0d0)
t_2 = 4.0d0 * (t_1 * t_1)
code = (t_2 + 1.0d0) / (t_2 + 2.0d0)
end function
public static double code(double t) {
double t_1 = t / (t + 1.0);
double t_2 = 4.0 * (t_1 * t_1);
return (t_2 + 1.0) / (t_2 + 2.0);
}
def code(t): t_1 = t / (t + 1.0) t_2 = 4.0 * (t_1 * t_1) return (t_2 + 1.0) / (t_2 + 2.0)
function code(t) t_1 = Float64(t / Float64(t + 1.0)) t_2 = Float64(4.0 * Float64(t_1 * t_1)) return Float64(Float64(t_2 + 1.0) / Float64(t_2 + 2.0)) end
function tmp = code(t) t_1 = t / (t + 1.0); t_2 = 4.0 * (t_1 * t_1); tmp = (t_2 + 1.0) / (t_2 + 2.0); end
code[t_] := Block[{t$95$1 = N[(t / N[(t + 1.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(4.0 * N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]}, N[(N[(t$95$2 + 1.0), $MachinePrecision] / N[(t$95$2 + 2.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{t + 1}\\
t_2 := 4 \cdot \left(t\_1 \cdot t\_1\right)\\
\frac{t\_2 + 1}{t\_2 + 2}
\end{array}
\end{array}
Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (t)
:precision binary64
(if (<= t -0.38)
0.8333333333333334
(if (<= t 0.86)
(/
(+ (* 4.0 (* t (* t (- 1.0 t)))) 1.0)
(- 2.0 (* 4.0 (* t (/ t (- -1.0 t))))))
(/
(+ 5.0 (/ (- (/ (+ 12.0 (/ -16.0 t)) t) 8.0) t))
(- 6.0 (/ (- (/ (- (/ 16.0 t) 12.0) t) -8.0) t))))))
double code(double t) {
double tmp;
if (t <= -0.38) {
tmp = 0.8333333333333334;
} else if (t <= 0.86) {
tmp = ((4.0 * (t * (t * (1.0 - t)))) + 1.0) / (2.0 - (4.0 * (t * (t / (-1.0 - t)))));
} else {
tmp = (5.0 + ((((12.0 + (-16.0 / t)) / t) - 8.0) / t)) / (6.0 - (((((16.0 / t) - 12.0) / t) - -8.0) / t));
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-0.38d0)) then
tmp = 0.8333333333333334d0
else if (t <= 0.86d0) then
tmp = ((4.0d0 * (t * (t * (1.0d0 - t)))) + 1.0d0) / (2.0d0 - (4.0d0 * (t * (t / ((-1.0d0) - t)))))
else
tmp = (5.0d0 + ((((12.0d0 + ((-16.0d0) / t)) / t) - 8.0d0) / t)) / (6.0d0 - (((((16.0d0 / t) - 12.0d0) / t) - (-8.0d0)) / t))
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if (t <= -0.38) {
tmp = 0.8333333333333334;
} else if (t <= 0.86) {
tmp = ((4.0 * (t * (t * (1.0 - t)))) + 1.0) / (2.0 - (4.0 * (t * (t / (-1.0 - t)))));
} else {
tmp = (5.0 + ((((12.0 + (-16.0 / t)) / t) - 8.0) / t)) / (6.0 - (((((16.0 / t) - 12.0) / t) - -8.0) / t));
}
return tmp;
}
def code(t): tmp = 0 if t <= -0.38: tmp = 0.8333333333333334 elif t <= 0.86: tmp = ((4.0 * (t * (t * (1.0 - t)))) + 1.0) / (2.0 - (4.0 * (t * (t / (-1.0 - t))))) else: tmp = (5.0 + ((((12.0 + (-16.0 / t)) / t) - 8.0) / t)) / (6.0 - (((((16.0 / t) - 12.0) / t) - -8.0) / t)) return tmp
function code(t) tmp = 0.0 if (t <= -0.38) tmp = 0.8333333333333334; elseif (t <= 0.86) tmp = Float64(Float64(Float64(4.0 * Float64(t * Float64(t * Float64(1.0 - t)))) + 1.0) / Float64(2.0 - Float64(4.0 * Float64(t * Float64(t / Float64(-1.0 - t)))))); else tmp = Float64(Float64(5.0 + Float64(Float64(Float64(Float64(12.0 + Float64(-16.0 / t)) / t) - 8.0) / t)) / Float64(6.0 - Float64(Float64(Float64(Float64(Float64(16.0 / t) - 12.0) / t) - -8.0) / t))); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if (t <= -0.38) tmp = 0.8333333333333334; elseif (t <= 0.86) tmp = ((4.0 * (t * (t * (1.0 - t)))) + 1.0) / (2.0 - (4.0 * (t * (t / (-1.0 - t))))); else tmp = (5.0 + ((((12.0 + (-16.0 / t)) / t) - 8.0) / t)) / (6.0 - (((((16.0 / t) - 12.0) / t) - -8.0) / t)); end tmp_2 = tmp; end
code[t_] := If[LessEqual[t, -0.38], 0.8333333333333334, If[LessEqual[t, 0.86], N[(N[(N[(4.0 * N[(t * N[(t * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / N[(2.0 - N[(4.0 * N[(t * N[(t / N[(-1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(5.0 + N[(N[(N[(N[(12.0 + N[(-16.0 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] - 8.0), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] / N[(6.0 - N[(N[(N[(N[(N[(16.0 / t), $MachinePrecision] - 12.0), $MachinePrecision] / t), $MachinePrecision] - -8.0), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.38:\\
\;\;\;\;0.8333333333333334\\
\mathbf{elif}\;t \leq 0.86:\\
\;\;\;\;\frac{4 \cdot \left(t \cdot \left(t \cdot \left(1 - t\right)\right)\right) + 1}{2 - 4 \cdot \left(t \cdot \frac{t}{-1 - t}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{5 + \frac{\frac{12 + \frac{-16}{t}}{t} - 8}{t}}{6 - \frac{\frac{\frac{16}{t} - 12}{t} - -8}{t}}\\
\end{array}
\end{array}
if t < -0.38Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
Taylor expanded in t around -inf 100.0%
mul-1-neg100.0%
unsub-neg100.0%
sub-neg100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
if -0.38 < t < 0.859999999999999987Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.5%
neg-mul-199.5%
sub-neg99.5%
Simplified99.5%
if 0.859999999999999987 < t Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around -inf 98.6%
mul-1-neg98.6%
unsub-neg98.6%
mul-1-neg98.6%
unsub-neg98.6%
sub-neg98.6%
associate-*r/98.6%
metadata-eval98.6%
distribute-neg-frac98.6%
metadata-eval98.6%
Simplified98.6%
Taylor expanded in t around inf 98.7%
Simplified98.7%
Final simplification99.4%
(FPCore (t)
:precision binary64
(if (<= t -0.38)
0.8333333333333334
(if (<= t 0.82)
(/
(+ (* 4.0 (* t (* t (- 1.0 t)))) 1.0)
(- 2.0 (* 4.0 (* t (/ t (- -1.0 t))))))
(/ (- 5.0 (/ (+ 8.0 (/ -12.0 t)) t)) (+ 6.0 (/ (- (/ 12.0 t) 8.0) t))))))
double code(double t) {
double tmp;
if (t <= -0.38) {
tmp = 0.8333333333333334;
} else if (t <= 0.82) {
tmp = ((4.0 * (t * (t * (1.0 - t)))) + 1.0) / (2.0 - (4.0 * (t * (t / (-1.0 - t)))));
} else {
tmp = (5.0 - ((8.0 + (-12.0 / t)) / t)) / (6.0 + (((12.0 / t) - 8.0) / t));
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-0.38d0)) then
tmp = 0.8333333333333334d0
else if (t <= 0.82d0) then
tmp = ((4.0d0 * (t * (t * (1.0d0 - t)))) + 1.0d0) / (2.0d0 - (4.0d0 * (t * (t / ((-1.0d0) - t)))))
else
tmp = (5.0d0 - ((8.0d0 + ((-12.0d0) / t)) / t)) / (6.0d0 + (((12.0d0 / t) - 8.0d0) / t))
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if (t <= -0.38) {
tmp = 0.8333333333333334;
} else if (t <= 0.82) {
tmp = ((4.0 * (t * (t * (1.0 - t)))) + 1.0) / (2.0 - (4.0 * (t * (t / (-1.0 - t)))));
} else {
tmp = (5.0 - ((8.0 + (-12.0 / t)) / t)) / (6.0 + (((12.0 / t) - 8.0) / t));
}
return tmp;
}
def code(t): tmp = 0 if t <= -0.38: tmp = 0.8333333333333334 elif t <= 0.82: tmp = ((4.0 * (t * (t * (1.0 - t)))) + 1.0) / (2.0 - (4.0 * (t * (t / (-1.0 - t))))) else: tmp = (5.0 - ((8.0 + (-12.0 / t)) / t)) / (6.0 + (((12.0 / t) - 8.0) / t)) return tmp
function code(t) tmp = 0.0 if (t <= -0.38) tmp = 0.8333333333333334; elseif (t <= 0.82) tmp = Float64(Float64(Float64(4.0 * Float64(t * Float64(t * Float64(1.0 - t)))) + 1.0) / Float64(2.0 - Float64(4.0 * Float64(t * Float64(t / Float64(-1.0 - t)))))); else tmp = Float64(Float64(5.0 - Float64(Float64(8.0 + Float64(-12.0 / t)) / t)) / Float64(6.0 + Float64(Float64(Float64(12.0 / t) - 8.0) / t))); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if (t <= -0.38) tmp = 0.8333333333333334; elseif (t <= 0.82) tmp = ((4.0 * (t * (t * (1.0 - t)))) + 1.0) / (2.0 - (4.0 * (t * (t / (-1.0 - t))))); else tmp = (5.0 - ((8.0 + (-12.0 / t)) / t)) / (6.0 + (((12.0 / t) - 8.0) / t)); end tmp_2 = tmp; end
code[t_] := If[LessEqual[t, -0.38], 0.8333333333333334, If[LessEqual[t, 0.82], N[(N[(N[(4.0 * N[(t * N[(t * N[(1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / N[(2.0 - N[(4.0 * N[(t * N[(t / N[(-1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(5.0 - N[(N[(8.0 + N[(-12.0 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] / N[(6.0 + N[(N[(N[(12.0 / t), $MachinePrecision] - 8.0), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.38:\\
\;\;\;\;0.8333333333333334\\
\mathbf{elif}\;t \leq 0.82:\\
\;\;\;\;\frac{4 \cdot \left(t \cdot \left(t \cdot \left(1 - t\right)\right)\right) + 1}{2 - 4 \cdot \left(t \cdot \frac{t}{-1 - t}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{5 - \frac{8 + \frac{-12}{t}}{t}}{6 + \frac{\frac{12}{t} - 8}{t}}\\
\end{array}
\end{array}
if t < -0.38Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
Taylor expanded in t around -inf 100.0%
mul-1-neg100.0%
unsub-neg100.0%
sub-neg100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
if -0.38 < t < 0.819999999999999951Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.5%
neg-mul-199.5%
sub-neg99.5%
Simplified99.5%
if 0.819999999999999951 < t Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around -inf 98.3%
mul-1-neg98.3%
unsub-neg98.3%
sub-neg98.3%
associate-*r/98.3%
metadata-eval98.3%
distribute-neg-frac98.3%
metadata-eval98.3%
Simplified98.3%
Taylor expanded in t around inf 98.4%
associate--l+98.4%
associate-*r/98.4%
metadata-eval98.4%
unpow298.4%
associate-/r*98.4%
metadata-eval98.4%
associate-*r/98.4%
div-sub98.4%
sub-neg98.4%
associate-*r/98.4%
metadata-eval98.4%
metadata-eval98.4%
distribute-neg-frac98.4%
distribute-neg-in98.4%
+-commutative98.4%
distribute-neg-frac98.4%
sub-neg98.4%
Simplified98.4%
Final simplification99.4%
(FPCore (t)
:precision binary64
(if (<= t -0.43)
0.8333333333333334
(if (<= t 1.85)
(/ (+ (* 4.0 (* t (/ t (+ t 1.0)))) 1.0) (+ 2.0 (* 4.0 (* t t))))
(/ (- 5.0 (/ 8.0 t)) (- 6.0 (/ 8.0 t))))))
double code(double t) {
double tmp;
if (t <= -0.43) {
tmp = 0.8333333333333334;
} else if (t <= 1.85) {
tmp = ((4.0 * (t * (t / (t + 1.0)))) + 1.0) / (2.0 + (4.0 * (t * t)));
} else {
tmp = (5.0 - (8.0 / t)) / (6.0 - (8.0 / t));
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-0.43d0)) then
tmp = 0.8333333333333334d0
else if (t <= 1.85d0) then
tmp = ((4.0d0 * (t * (t / (t + 1.0d0)))) + 1.0d0) / (2.0d0 + (4.0d0 * (t * t)))
else
tmp = (5.0d0 - (8.0d0 / t)) / (6.0d0 - (8.0d0 / t))
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if (t <= -0.43) {
tmp = 0.8333333333333334;
} else if (t <= 1.85) {
tmp = ((4.0 * (t * (t / (t + 1.0)))) + 1.0) / (2.0 + (4.0 * (t * t)));
} else {
tmp = (5.0 - (8.0 / t)) / (6.0 - (8.0 / t));
}
return tmp;
}
def code(t): tmp = 0 if t <= -0.43: tmp = 0.8333333333333334 elif t <= 1.85: tmp = ((4.0 * (t * (t / (t + 1.0)))) + 1.0) / (2.0 + (4.0 * (t * t))) else: tmp = (5.0 - (8.0 / t)) / (6.0 - (8.0 / t)) return tmp
function code(t) tmp = 0.0 if (t <= -0.43) tmp = 0.8333333333333334; elseif (t <= 1.85) tmp = Float64(Float64(Float64(4.0 * Float64(t * Float64(t / Float64(t + 1.0)))) + 1.0) / Float64(2.0 + Float64(4.0 * Float64(t * t)))); else tmp = Float64(Float64(5.0 - Float64(8.0 / t)) / Float64(6.0 - Float64(8.0 / t))); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if (t <= -0.43) tmp = 0.8333333333333334; elseif (t <= 1.85) tmp = ((4.0 * (t * (t / (t + 1.0)))) + 1.0) / (2.0 + (4.0 * (t * t))); else tmp = (5.0 - (8.0 / t)) / (6.0 - (8.0 / t)); end tmp_2 = tmp; end
code[t_] := If[LessEqual[t, -0.43], 0.8333333333333334, If[LessEqual[t, 1.85], N[(N[(N[(4.0 * N[(t * N[(t / N[(t + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / N[(2.0 + N[(4.0 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(5.0 - N[(8.0 / t), $MachinePrecision]), $MachinePrecision] / N[(6.0 - N[(8.0 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.43:\\
\;\;\;\;0.8333333333333334\\
\mathbf{elif}\;t \leq 1.85:\\
\;\;\;\;\frac{4 \cdot \left(t \cdot \frac{t}{t + 1}\right) + 1}{2 + 4 \cdot \left(t \cdot t\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{5 - \frac{8}{t}}{6 - \frac{8}{t}}\\
\end{array}
\end{array}
if t < -0.429999999999999993Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
Taylor expanded in t around -inf 100.0%
mul-1-neg100.0%
unsub-neg100.0%
sub-neg100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
if -0.429999999999999993 < t < 1.8500000000000001Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.5%
if 1.8500000000000001 < t Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 97.4%
associate-*r/97.4%
metadata-eval97.4%
Simplified97.4%
Taylor expanded in t around inf 97.6%
associate-*r/97.6%
metadata-eval97.6%
Simplified97.6%
Final simplification99.2%
(FPCore (t)
:precision binary64
(if (<= t -0.43)
0.8333333333333334
(if (<= t 1.1)
(/ (+ (* 4.0 (* t (/ t (+ t 1.0)))) 1.0) (+ 2.0 (* 4.0 (* t t))))
(/ (- 5.0 (/ (+ 8.0 (/ -12.0 t)) t)) (+ 6.0 (/ (- (/ 12.0 t) 8.0) t))))))
double code(double t) {
double tmp;
if (t <= -0.43) {
tmp = 0.8333333333333334;
} else if (t <= 1.1) {
tmp = ((4.0 * (t * (t / (t + 1.0)))) + 1.0) / (2.0 + (4.0 * (t * t)));
} else {
tmp = (5.0 - ((8.0 + (-12.0 / t)) / t)) / (6.0 + (((12.0 / t) - 8.0) / t));
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-0.43d0)) then
tmp = 0.8333333333333334d0
else if (t <= 1.1d0) then
tmp = ((4.0d0 * (t * (t / (t + 1.0d0)))) + 1.0d0) / (2.0d0 + (4.0d0 * (t * t)))
else
tmp = (5.0d0 - ((8.0d0 + ((-12.0d0) / t)) / t)) / (6.0d0 + (((12.0d0 / t) - 8.0d0) / t))
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if (t <= -0.43) {
tmp = 0.8333333333333334;
} else if (t <= 1.1) {
tmp = ((4.0 * (t * (t / (t + 1.0)))) + 1.0) / (2.0 + (4.0 * (t * t)));
} else {
tmp = (5.0 - ((8.0 + (-12.0 / t)) / t)) / (6.0 + (((12.0 / t) - 8.0) / t));
}
return tmp;
}
def code(t): tmp = 0 if t <= -0.43: tmp = 0.8333333333333334 elif t <= 1.1: tmp = ((4.0 * (t * (t / (t + 1.0)))) + 1.0) / (2.0 + (4.0 * (t * t))) else: tmp = (5.0 - ((8.0 + (-12.0 / t)) / t)) / (6.0 + (((12.0 / t) - 8.0) / t)) return tmp
function code(t) tmp = 0.0 if (t <= -0.43) tmp = 0.8333333333333334; elseif (t <= 1.1) tmp = Float64(Float64(Float64(4.0 * Float64(t * Float64(t / Float64(t + 1.0)))) + 1.0) / Float64(2.0 + Float64(4.0 * Float64(t * t)))); else tmp = Float64(Float64(5.0 - Float64(Float64(8.0 + Float64(-12.0 / t)) / t)) / Float64(6.0 + Float64(Float64(Float64(12.0 / t) - 8.0) / t))); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if (t <= -0.43) tmp = 0.8333333333333334; elseif (t <= 1.1) tmp = ((4.0 * (t * (t / (t + 1.0)))) + 1.0) / (2.0 + (4.0 * (t * t))); else tmp = (5.0 - ((8.0 + (-12.0 / t)) / t)) / (6.0 + (((12.0 / t) - 8.0) / t)); end tmp_2 = tmp; end
code[t_] := If[LessEqual[t, -0.43], 0.8333333333333334, If[LessEqual[t, 1.1], N[(N[(N[(4.0 * N[(t * N[(t / N[(t + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / N[(2.0 + N[(4.0 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(5.0 - N[(N[(8.0 + N[(-12.0 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] / N[(6.0 + N[(N[(N[(12.0 / t), $MachinePrecision] - 8.0), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.43:\\
\;\;\;\;0.8333333333333334\\
\mathbf{elif}\;t \leq 1.1:\\
\;\;\;\;\frac{4 \cdot \left(t \cdot \frac{t}{t + 1}\right) + 1}{2 + 4 \cdot \left(t \cdot t\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{5 - \frac{8 + \frac{-12}{t}}{t}}{6 + \frac{\frac{12}{t} - 8}{t}}\\
\end{array}
\end{array}
if t < -0.429999999999999993Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
Taylor expanded in t around -inf 100.0%
mul-1-neg100.0%
unsub-neg100.0%
sub-neg100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
if -0.429999999999999993 < t < 1.1000000000000001Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.5%
if 1.1000000000000001 < t Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around -inf 98.3%
mul-1-neg98.3%
unsub-neg98.3%
sub-neg98.3%
associate-*r/98.3%
metadata-eval98.3%
distribute-neg-frac98.3%
metadata-eval98.3%
Simplified98.3%
Taylor expanded in t around inf 98.4%
associate--l+98.4%
associate-*r/98.4%
metadata-eval98.4%
unpow298.4%
associate-/r*98.4%
metadata-eval98.4%
associate-*r/98.4%
div-sub98.4%
sub-neg98.4%
associate-*r/98.4%
metadata-eval98.4%
metadata-eval98.4%
distribute-neg-frac98.4%
distribute-neg-in98.4%
+-commutative98.4%
distribute-neg-frac98.4%
sub-neg98.4%
Simplified98.4%
Final simplification99.4%
(FPCore (t) :precision binary64 (if (<= t -0.34) 0.8333333333333334 (if (<= t 2.0) 0.5 (/ (- 5.0 (/ 8.0 t)) (- 6.0 (/ 8.0 t))))))
double code(double t) {
double tmp;
if (t <= -0.34) {
tmp = 0.8333333333333334;
} else if (t <= 2.0) {
tmp = 0.5;
} else {
tmp = (5.0 - (8.0 / t)) / (6.0 - (8.0 / t));
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-0.34d0)) then
tmp = 0.8333333333333334d0
else if (t <= 2.0d0) then
tmp = 0.5d0
else
tmp = (5.0d0 - (8.0d0 / t)) / (6.0d0 - (8.0d0 / t))
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if (t <= -0.34) {
tmp = 0.8333333333333334;
} else if (t <= 2.0) {
tmp = 0.5;
} else {
tmp = (5.0 - (8.0 / t)) / (6.0 - (8.0 / t));
}
return tmp;
}
def code(t): tmp = 0 if t <= -0.34: tmp = 0.8333333333333334 elif t <= 2.0: tmp = 0.5 else: tmp = (5.0 - (8.0 / t)) / (6.0 - (8.0 / t)) return tmp
function code(t) tmp = 0.0 if (t <= -0.34) tmp = 0.8333333333333334; elseif (t <= 2.0) tmp = 0.5; else tmp = Float64(Float64(5.0 - Float64(8.0 / t)) / Float64(6.0 - Float64(8.0 / t))); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if (t <= -0.34) tmp = 0.8333333333333334; elseif (t <= 2.0) tmp = 0.5; else tmp = (5.0 - (8.0 / t)) / (6.0 - (8.0 / t)); end tmp_2 = tmp; end
code[t_] := If[LessEqual[t, -0.34], 0.8333333333333334, If[LessEqual[t, 2.0], 0.5, N[(N[(5.0 - N[(8.0 / t), $MachinePrecision]), $MachinePrecision] / N[(6.0 - N[(8.0 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.34:\\
\;\;\;\;0.8333333333333334\\
\mathbf{elif}\;t \leq 2:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;\frac{5 - \frac{8}{t}}{6 - \frac{8}{t}}\\
\end{array}
\end{array}
if t < -0.340000000000000024Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
Taylor expanded in t around -inf 100.0%
mul-1-neg100.0%
unsub-neg100.0%
sub-neg100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
if -0.340000000000000024 < t < 2Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.5%
neg-mul-199.5%
sub-neg99.5%
Simplified99.5%
Taylor expanded in t around 0 99.5%
if 2 < t Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 97.4%
associate-*r/97.4%
metadata-eval97.4%
Simplified97.4%
Taylor expanded in t around inf 97.6%
associate-*r/97.6%
metadata-eval97.6%
Simplified97.6%
Final simplification99.2%
(FPCore (t) :precision binary64 (if (<= t -0.34) 0.8333333333333334 (if (<= t 1.0) 0.5 0.8333333333333334)))
double code(double t) {
double tmp;
if (t <= -0.34) {
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.34d0)) 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.34) {
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.34: 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.34) 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.34) 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.34], 0.8333333333333334, If[LessEqual[t, 1.0], 0.5, 0.8333333333333334]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.34:\\
\;\;\;\;0.8333333333333334\\
\mathbf{elif}\;t \leq 1:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;0.8333333333333334\\
\end{array}
\end{array}
if t < -0.340000000000000024 or 1 < t Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 98.3%
Taylor expanded in t around -inf 98.3%
mul-1-neg98.3%
unsub-neg98.3%
sub-neg98.3%
distribute-neg-frac98.3%
metadata-eval98.3%
Simplified98.3%
Taylor expanded in t around inf 98.4%
if -0.340000000000000024 < t < 1Initial program 100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around 0 99.5%
Taylor expanded in t around 0 99.5%
neg-mul-199.5%
sub-neg99.5%
Simplified99.5%
Taylor expanded in t around 0 99.5%
Final simplification99.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%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
associate-/l*100.0%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around 0 55.8%
Taylor expanded in t around 0 55.5%
neg-mul-155.5%
sub-neg55.5%
Simplified55.5%
Taylor expanded in t around 0 63.0%
Final simplification63.0%
herbie shell --seed 2024053
(FPCore (t)
:name "Kahan p13 Example 1"
:precision binary64
(/ (+ 1.0 (* (/ (* 2.0 t) (+ 1.0 t)) (/ (* 2.0 t) (+ 1.0 t)))) (+ 2.0 (* (/ (* 2.0 t) (+ 1.0 t)) (/ (* 2.0 t) (+ 1.0 t))))))