
(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 7 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 99.6%
associate-/l*99.6%
associate-/l*99.6%
swap-sqr99.6%
metadata-eval99.6%
associate-/l*99.6%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (t)
:precision binary64
(if (<= t -0.35)
(+
0.8333333333333334
(/
(-
(/ (+ 0.037037037037037035 (/ 0.04938271604938271 t)) t)
0.2222222222222222)
t))
(if (<= t 0.49)
(/ (+ (* 4.0 (* t t)) 1.0) 2.0)
(/ (- 5.0 (/ (+ (/ -12.0 t) 8.0) t)) (+ (/ (+ -8.0 (/ 12.0 t)) t) 6.0)))))
double code(double t) {
double tmp;
if (t <= -0.35) {
tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t);
} else if (t <= 0.49) {
tmp = ((4.0 * (t * t)) + 1.0) / 2.0;
} else {
tmp = (5.0 - (((-12.0 / t) + 8.0) / t)) / (((-8.0 + (12.0 / t)) / t) + 6.0);
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-0.35d0)) then
tmp = 0.8333333333333334d0 + ((((0.037037037037037035d0 + (0.04938271604938271d0 / t)) / t) - 0.2222222222222222d0) / t)
else if (t <= 0.49d0) then
tmp = ((4.0d0 * (t * t)) + 1.0d0) / 2.0d0
else
tmp = (5.0d0 - ((((-12.0d0) / t) + 8.0d0) / t)) / ((((-8.0d0) + (12.0d0 / t)) / t) + 6.0d0)
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if (t <= -0.35) {
tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t);
} else if (t <= 0.49) {
tmp = ((4.0 * (t * t)) + 1.0) / 2.0;
} else {
tmp = (5.0 - (((-12.0 / t) + 8.0) / t)) / (((-8.0 + (12.0 / t)) / t) + 6.0);
}
return tmp;
}
def code(t): tmp = 0 if t <= -0.35: tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t) elif t <= 0.49: tmp = ((4.0 * (t * t)) + 1.0) / 2.0 else: tmp = (5.0 - (((-12.0 / t) + 8.0) / t)) / (((-8.0 + (12.0 / t)) / t) + 6.0) return tmp
function code(t) tmp = 0.0 if (t <= -0.35) tmp = Float64(0.8333333333333334 + Float64(Float64(Float64(Float64(0.037037037037037035 + Float64(0.04938271604938271 / t)) / t) - 0.2222222222222222) / t)); elseif (t <= 0.49) tmp = Float64(Float64(Float64(4.0 * Float64(t * t)) + 1.0) / 2.0); else tmp = Float64(Float64(5.0 - Float64(Float64(Float64(-12.0 / t) + 8.0) / t)) / Float64(Float64(Float64(-8.0 + Float64(12.0 / t)) / t) + 6.0)); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if (t <= -0.35) tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t); elseif (t <= 0.49) tmp = ((4.0 * (t * t)) + 1.0) / 2.0; else tmp = (5.0 - (((-12.0 / t) + 8.0) / t)) / (((-8.0 + (12.0 / t)) / t) + 6.0); end tmp_2 = tmp; end
code[t_] := If[LessEqual[t, -0.35], N[(0.8333333333333334 + N[(N[(N[(N[(0.037037037037037035 + N[(0.04938271604938271 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] - 0.2222222222222222), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.49], N[(N[(N[(4.0 * N[(t * t), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / 2.0), $MachinePrecision], N[(N[(5.0 - N[(N[(N[(-12.0 / t), $MachinePrecision] + 8.0), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(-8.0 + N[(12.0 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + 6.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.35:\\
\;\;\;\;0.8333333333333334 + \frac{\frac{0.037037037037037035 + \frac{0.04938271604938271}{t}}{t} - 0.2222222222222222}{t}\\
\mathbf{elif}\;t \leq 0.49:\\
\;\;\;\;\frac{4 \cdot \left(t \cdot t\right) + 1}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{5 - \frac{\frac{-12}{t} + 8}{t}}{\frac{-8 + \frac{12}{t}}{t} + 6}\\
\end{array}
\end{array}
if t < -0.34999999999999998Initial 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.1%
mul-1-neg97.1%
unsub-neg97.1%
mul-1-neg97.1%
unsub-neg97.1%
sub-neg97.1%
associate-*r/97.1%
metadata-eval97.1%
distribute-neg-frac97.1%
metadata-eval97.1%
Simplified97.1%
Taylor expanded in t around -inf 97.5%
mul-1-neg97.5%
unsub-neg97.5%
mul-1-neg97.5%
unsub-neg97.5%
associate-*r/97.5%
metadata-eval97.5%
Simplified97.5%
if -0.34999999999999998 < t < 0.48999999999999999Initial 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.6%
Taylor expanded in t around 0 99.6%
Taylor expanded in t around 0 99.6%
if 0.48999999999999999 < t Initial program 98.6%
associate-/l*98.6%
associate-/l*98.6%
swap-sqr98.6%
metadata-eval98.6%
associate-/l*98.6%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
sub-neg100.0%
+-commutative100.0%
associate-+r+100.0%
+-commutative100.0%
sub-neg100.0%
associate--r-100.0%
associate-*r/100.0%
metadata-eval100.0%
unpow2100.0%
associate-/r*100.0%
metadata-eval100.0%
associate-*r/100.0%
div-sub100.0%
unsub-neg100.0%
mul-1-neg100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
sub-neg100.0%
+-commutative100.0%
associate-+r+100.0%
+-commutative100.0%
sub-neg100.0%
associate--r-100.0%
associate-*r/100.0%
metadata-eval100.0%
unpow2100.0%
associate-/r*100.0%
metadata-eval100.0%
associate-*r/100.0%
div-sub100.0%
sub-neg100.0%
associate-*r/100.0%
metadata-eval100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Simplified100.0%
Final simplification99.1%
(FPCore (t)
:precision binary64
(if (or (<= t -0.195) (not (<= t 0.17)))
(-
0.8333333333333334
(/ (+ 0.2222222222222222 (/ -0.037037037037037035 t)) t))
(/ (+ (* 4.0 t) 1.0) 2.0)))
double code(double t) {
double tmp;
if ((t <= -0.195) || !(t <= 0.17)) {
tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t);
} else {
tmp = ((4.0 * t) + 1.0) / 2.0;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if ((t <= (-0.195d0)) .or. (.not. (t <= 0.17d0))) then
tmp = 0.8333333333333334d0 - ((0.2222222222222222d0 + ((-0.037037037037037035d0) / t)) / t)
else
tmp = ((4.0d0 * t) + 1.0d0) / 2.0d0
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if ((t <= -0.195) || !(t <= 0.17)) {
tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t);
} else {
tmp = ((4.0 * t) + 1.0) / 2.0;
}
return tmp;
}
def code(t): tmp = 0 if (t <= -0.195) or not (t <= 0.17): tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t) else: tmp = ((4.0 * t) + 1.0) / 2.0 return tmp
function code(t) tmp = 0.0 if ((t <= -0.195) || !(t <= 0.17)) tmp = Float64(0.8333333333333334 - Float64(Float64(0.2222222222222222 + Float64(-0.037037037037037035 / t)) / t)); else tmp = Float64(Float64(Float64(4.0 * t) + 1.0) / 2.0); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if ((t <= -0.195) || ~((t <= 0.17))) tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t); else tmp = ((4.0 * t) + 1.0) / 2.0; end tmp_2 = tmp; end
code[t_] := If[Or[LessEqual[t, -0.195], N[Not[LessEqual[t, 0.17]], $MachinePrecision]], N[(0.8333333333333334 - N[(N[(0.2222222222222222 + N[(-0.037037037037037035 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], N[(N[(N[(4.0 * t), $MachinePrecision] + 1.0), $MachinePrecision] / 2.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.195 \lor \neg \left(t \leq 0.17\right):\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222 + \frac{-0.037037037037037035}{t}}{t}\\
\mathbf{else}:\\
\;\;\;\;\frac{4 \cdot t + 1}{2}\\
\end{array}
\end{array}
if t < -0.19500000000000001 or 0.170000000000000012 < t Initial program 99.3%
associate-/l*99.3%
associate-/l*99.3%
swap-sqr99.3%
metadata-eval99.3%
associate-/l*99.3%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 98.3%
sub-neg98.3%
+-commutative98.3%
associate-+r+98.3%
+-commutative98.3%
sub-neg98.3%
associate--r-98.3%
associate-*r/98.3%
metadata-eval98.3%
unpow298.3%
associate-/r*98.3%
metadata-eval98.3%
associate-*r/98.3%
div-sub98.3%
unsub-neg98.3%
mul-1-neg98.3%
+-commutative98.3%
Simplified98.3%
Taylor expanded in t around -inf 98.5%
mul-1-neg98.5%
unsub-neg98.5%
sub-neg98.5%
associate-*r/98.5%
metadata-eval98.5%
distribute-neg-frac98.5%
metadata-eval98.5%
Simplified98.5%
if -0.19500000000000001 < t < 0.170000000000000012Initial 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.6%
Taylor expanded in t around 0 99.6%
Taylor expanded in t around inf 98.2%
*-commutative98.2%
Simplified98.2%
Final simplification98.4%
(FPCore (t)
:precision binary64
(if (or (<= t -0.62) (not (<= t 0.215)))
(-
0.8333333333333334
(/ (+ 0.2222222222222222 (/ -0.037037037037037035 t)) t))
(/ (+ (* 4.0 (* t t)) 1.0) 2.0)))
double code(double t) {
double tmp;
if ((t <= -0.62) || !(t <= 0.215)) {
tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t);
} else {
tmp = ((4.0 * (t * t)) + 1.0) / 2.0;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if ((t <= (-0.62d0)) .or. (.not. (t <= 0.215d0))) then
tmp = 0.8333333333333334d0 - ((0.2222222222222222d0 + ((-0.037037037037037035d0) / t)) / t)
else
tmp = ((4.0d0 * (t * t)) + 1.0d0) / 2.0d0
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if ((t <= -0.62) || !(t <= 0.215)) {
tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t);
} else {
tmp = ((4.0 * (t * t)) + 1.0) / 2.0;
}
return tmp;
}
def code(t): tmp = 0 if (t <= -0.62) or not (t <= 0.215): tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t) else: tmp = ((4.0 * (t * t)) + 1.0) / 2.0 return tmp
function code(t) tmp = 0.0 if ((t <= -0.62) || !(t <= 0.215)) tmp = Float64(0.8333333333333334 - Float64(Float64(0.2222222222222222 + Float64(-0.037037037037037035 / t)) / t)); else tmp = Float64(Float64(Float64(4.0 * Float64(t * t)) + 1.0) / 2.0); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if ((t <= -0.62) || ~((t <= 0.215))) tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t); else tmp = ((4.0 * (t * t)) + 1.0) / 2.0; end tmp_2 = tmp; end
code[t_] := If[Or[LessEqual[t, -0.62], N[Not[LessEqual[t, 0.215]], $MachinePrecision]], N[(0.8333333333333334 - N[(N[(0.2222222222222222 + N[(-0.037037037037037035 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], N[(N[(N[(4.0 * N[(t * t), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / 2.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.62 \lor \neg \left(t \leq 0.215\right):\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222 + \frac{-0.037037037037037035}{t}}{t}\\
\mathbf{else}:\\
\;\;\;\;\frac{4 \cdot \left(t \cdot t\right) + 1}{2}\\
\end{array}
\end{array}
if t < -0.619999999999999996 or 0.214999999999999997 < t Initial program 99.3%
associate-/l*99.3%
associate-/l*99.3%
swap-sqr99.3%
metadata-eval99.3%
associate-/l*99.3%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 98.3%
sub-neg98.3%
+-commutative98.3%
associate-+r+98.3%
+-commutative98.3%
sub-neg98.3%
associate--r-98.3%
associate-*r/98.3%
metadata-eval98.3%
unpow298.3%
associate-/r*98.3%
metadata-eval98.3%
associate-*r/98.3%
div-sub98.3%
unsub-neg98.3%
mul-1-neg98.3%
+-commutative98.3%
Simplified98.3%
Taylor expanded in t around -inf 98.5%
mul-1-neg98.5%
unsub-neg98.5%
sub-neg98.5%
associate-*r/98.5%
metadata-eval98.5%
distribute-neg-frac98.5%
metadata-eval98.5%
Simplified98.5%
if -0.619999999999999996 < t < 0.214999999999999997Initial 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.6%
Taylor expanded in t around 0 99.6%
Taylor expanded in t around 0 99.6%
Final simplification99.0%
(FPCore (t)
:precision binary64
(if (<= t -0.35)
(+
0.8333333333333334
(/
(-
(/ (+ 0.037037037037037035 (/ 0.04938271604938271 t)) t)
0.2222222222222222)
t))
(if (<= t 0.215)
(/ (+ (* 4.0 (* t t)) 1.0) 2.0)
(-
0.8333333333333334
(/ (+ 0.2222222222222222 (/ -0.037037037037037035 t)) t)))))
double code(double t) {
double tmp;
if (t <= -0.35) {
tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t);
} else if (t <= 0.215) {
tmp = ((4.0 * (t * t)) + 1.0) / 2.0;
} else {
tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t);
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-0.35d0)) then
tmp = 0.8333333333333334d0 + ((((0.037037037037037035d0 + (0.04938271604938271d0 / t)) / t) - 0.2222222222222222d0) / t)
else if (t <= 0.215d0) then
tmp = ((4.0d0 * (t * t)) + 1.0d0) / 2.0d0
else
tmp = 0.8333333333333334d0 - ((0.2222222222222222d0 + ((-0.037037037037037035d0) / t)) / t)
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if (t <= -0.35) {
tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t);
} else if (t <= 0.215) {
tmp = ((4.0 * (t * t)) + 1.0) / 2.0;
} else {
tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t);
}
return tmp;
}
def code(t): tmp = 0 if t <= -0.35: tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t) elif t <= 0.215: tmp = ((4.0 * (t * t)) + 1.0) / 2.0 else: tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t) return tmp
function code(t) tmp = 0.0 if (t <= -0.35) tmp = Float64(0.8333333333333334 + Float64(Float64(Float64(Float64(0.037037037037037035 + Float64(0.04938271604938271 / t)) / t) - 0.2222222222222222) / t)); elseif (t <= 0.215) tmp = Float64(Float64(Float64(4.0 * Float64(t * t)) + 1.0) / 2.0); else tmp = Float64(0.8333333333333334 - Float64(Float64(0.2222222222222222 + Float64(-0.037037037037037035 / t)) / t)); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if (t <= -0.35) tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t); elseif (t <= 0.215) tmp = ((4.0 * (t * t)) + 1.0) / 2.0; else tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t); end tmp_2 = tmp; end
code[t_] := If[LessEqual[t, -0.35], N[(0.8333333333333334 + N[(N[(N[(N[(0.037037037037037035 + N[(0.04938271604938271 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] - 0.2222222222222222), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.215], N[(N[(N[(4.0 * N[(t * t), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / 2.0), $MachinePrecision], N[(0.8333333333333334 - N[(N[(0.2222222222222222 + N[(-0.037037037037037035 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.35:\\
\;\;\;\;0.8333333333333334 + \frac{\frac{0.037037037037037035 + \frac{0.04938271604938271}{t}}{t} - 0.2222222222222222}{t}\\
\mathbf{elif}\;t \leq 0.215:\\
\;\;\;\;\frac{4 \cdot \left(t \cdot t\right) + 1}{2}\\
\mathbf{else}:\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222 + \frac{-0.037037037037037035}{t}}{t}\\
\end{array}
\end{array}
if t < -0.34999999999999998Initial 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.1%
mul-1-neg97.1%
unsub-neg97.1%
mul-1-neg97.1%
unsub-neg97.1%
sub-neg97.1%
associate-*r/97.1%
metadata-eval97.1%
distribute-neg-frac97.1%
metadata-eval97.1%
Simplified97.1%
Taylor expanded in t around -inf 97.5%
mul-1-neg97.5%
unsub-neg97.5%
mul-1-neg97.5%
unsub-neg97.5%
associate-*r/97.5%
metadata-eval97.5%
Simplified97.5%
if -0.34999999999999998 < t < 0.214999999999999997Initial 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.6%
Taylor expanded in t around 0 99.6%
Taylor expanded in t around 0 99.6%
if 0.214999999999999997 < t Initial program 98.6%
associate-/l*98.6%
associate-/l*98.6%
swap-sqr98.6%
metadata-eval98.6%
associate-/l*98.6%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
sub-neg100.0%
+-commutative100.0%
associate-+r+100.0%
+-commutative100.0%
sub-neg100.0%
associate--r-100.0%
associate-*r/100.0%
metadata-eval100.0%
unpow2100.0%
associate-/r*100.0%
metadata-eval100.0%
associate-*r/100.0%
div-sub100.0%
unsub-neg100.0%
mul-1-neg100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in t around -inf 100.0%
mul-1-neg100.0%
unsub-neg100.0%
sub-neg100.0%
associate-*r/100.0%
metadata-eval100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Simplified100.0%
Final simplification99.1%
(FPCore (t) :precision binary64 (if (or (<= t -0.18) (not (<= t 0.38))) (- 0.8333333333333334 (/ 0.2222222222222222 t)) (/ (+ (* 4.0 t) 1.0) 2.0)))
double code(double t) {
double tmp;
if ((t <= -0.18) || !(t <= 0.38)) {
tmp = 0.8333333333333334 - (0.2222222222222222 / t);
} else {
tmp = ((4.0 * t) + 1.0) / 2.0;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if ((t <= (-0.18d0)) .or. (.not. (t <= 0.38d0))) then
tmp = 0.8333333333333334d0 - (0.2222222222222222d0 / t)
else
tmp = ((4.0d0 * t) + 1.0d0) / 2.0d0
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if ((t <= -0.18) || !(t <= 0.38)) {
tmp = 0.8333333333333334 - (0.2222222222222222 / t);
} else {
tmp = ((4.0 * t) + 1.0) / 2.0;
}
return tmp;
}
def code(t): tmp = 0 if (t <= -0.18) or not (t <= 0.38): tmp = 0.8333333333333334 - (0.2222222222222222 / t) else: tmp = ((4.0 * t) + 1.0) / 2.0 return tmp
function code(t) tmp = 0.0 if ((t <= -0.18) || !(t <= 0.38)) tmp = Float64(0.8333333333333334 - Float64(0.2222222222222222 / t)); else tmp = Float64(Float64(Float64(4.0 * t) + 1.0) / 2.0); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if ((t <= -0.18) || ~((t <= 0.38))) tmp = 0.8333333333333334 - (0.2222222222222222 / t); else tmp = ((4.0 * t) + 1.0) / 2.0; end tmp_2 = tmp; end
code[t_] := If[Or[LessEqual[t, -0.18], N[Not[LessEqual[t, 0.38]], $MachinePrecision]], N[(0.8333333333333334 - N[(0.2222222222222222 / t), $MachinePrecision]), $MachinePrecision], N[(N[(N[(4.0 * t), $MachinePrecision] + 1.0), $MachinePrecision] / 2.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.18 \lor \neg \left(t \leq 0.38\right):\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222}{t}\\
\mathbf{else}:\\
\;\;\;\;\frac{4 \cdot t + 1}{2}\\
\end{array}
\end{array}
if t < -0.17999999999999999 or 0.38 < t Initial program 99.3%
associate-/l*99.3%
associate-/l*99.3%
swap-sqr99.3%
metadata-eval99.3%
associate-/l*99.3%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 97.6%
associate-*r/97.6%
metadata-eval97.6%
Simplified97.6%
Taylor expanded in t around inf 98.0%
associate-*r/98.0%
metadata-eval98.0%
Simplified98.0%
if -0.17999999999999999 < 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 0 99.6%
Taylor expanded in t around 0 99.6%
Taylor expanded in t around inf 98.2%
*-commutative98.2%
Simplified98.2%
Final simplification98.1%
(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 99.6%
associate-/l*99.6%
associate-/l*99.6%
swap-sqr99.6%
metadata-eval99.6%
associate-/l*99.6%
associate-/l*100.0%
swap-sqr100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 61.1%
Taylor expanded in t around inf 62.2%
Final simplification62.2%
herbie shell --seed 2024066
(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))))))