
(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 4.0) (+ 1.0 t)) (+ 1.0 t))))) (/ (+ 1.0 t_1) (+ t_1 2.0))))
double code(double t) {
double t_1 = t * (((t * 4.0) / (1.0 + t)) / (1.0 + t));
return (1.0 + t_1) / (t_1 + 2.0);
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
t_1 = t * (((t * 4.0d0) / (1.0d0 + t)) / (1.0d0 + t))
code = (1.0d0 + t_1) / (t_1 + 2.0d0)
end function
public static double code(double t) {
double t_1 = t * (((t * 4.0) / (1.0 + t)) / (1.0 + t));
return (1.0 + t_1) / (t_1 + 2.0);
}
def code(t): t_1 = t * (((t * 4.0) / (1.0 + t)) / (1.0 + t)) return (1.0 + t_1) / (t_1 + 2.0)
function code(t) t_1 = Float64(t * Float64(Float64(Float64(t * 4.0) / Float64(1.0 + t)) / Float64(1.0 + t))) return Float64(Float64(1.0 + t_1) / Float64(t_1 + 2.0)) end
function tmp = code(t) t_1 = t * (((t * 4.0) / (1.0 + t)) / (1.0 + t)); tmp = (1.0 + t_1) / (t_1 + 2.0); end
code[t_] := Block[{t$95$1 = N[(t * N[(N[(N[(t * 4.0), $MachinePrecision] / N[(1.0 + t), $MachinePrecision]), $MachinePrecision] / N[(1.0 + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, N[(N[(1.0 + t$95$1), $MachinePrecision] / N[(t$95$1 + 2.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \frac{\frac{t \cdot 4}{1 + t}}{1 + t}\\
\frac{1 + t_1}{t_1 + 2}
\end{array}
\end{array}
Initial program 100.0%
associate-/l*100.0%
associate-*r/100.0%
associate-/r/100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*r/100.0%
*-commutative100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*l*100.0%
metadata-eval100.0%
associate-/l*100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (t)
:precision binary64
(if (<= t -0.8)
(- 0.8333333333333334 (/ 0.2222222222222222 t))
(if (<= t 0.235)
(+ (* t t) 0.5)
(+
0.8333333333333334
(+ (/ 0.037037037037037035 (* t t)) (/ -0.2222222222222222 t))))))
double code(double t) {
double tmp;
if (t <= -0.8) {
tmp = 0.8333333333333334 - (0.2222222222222222 / t);
} else if (t <= 0.235) {
tmp = (t * t) + 0.5;
} else {
tmp = 0.8333333333333334 + ((0.037037037037037035 / (t * t)) + (-0.2222222222222222 / t));
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-0.8d0)) then
tmp = 0.8333333333333334d0 - (0.2222222222222222d0 / t)
else if (t <= 0.235d0) then
tmp = (t * t) + 0.5d0
else
tmp = 0.8333333333333334d0 + ((0.037037037037037035d0 / (t * t)) + ((-0.2222222222222222d0) / t))
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if (t <= -0.8) {
tmp = 0.8333333333333334 - (0.2222222222222222 / t);
} else if (t <= 0.235) {
tmp = (t * t) + 0.5;
} else {
tmp = 0.8333333333333334 + ((0.037037037037037035 / (t * t)) + (-0.2222222222222222 / t));
}
return tmp;
}
def code(t): tmp = 0 if t <= -0.8: tmp = 0.8333333333333334 - (0.2222222222222222 / t) elif t <= 0.235: tmp = (t * t) + 0.5 else: tmp = 0.8333333333333334 + ((0.037037037037037035 / (t * t)) + (-0.2222222222222222 / t)) return tmp
function code(t) tmp = 0.0 if (t <= -0.8) tmp = Float64(0.8333333333333334 - Float64(0.2222222222222222 / t)); elseif (t <= 0.235) tmp = Float64(Float64(t * t) + 0.5); else tmp = Float64(0.8333333333333334 + Float64(Float64(0.037037037037037035 / Float64(t * t)) + Float64(-0.2222222222222222 / t))); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if (t <= -0.8) tmp = 0.8333333333333334 - (0.2222222222222222 / t); elseif (t <= 0.235) tmp = (t * t) + 0.5; else tmp = 0.8333333333333334 + ((0.037037037037037035 / (t * t)) + (-0.2222222222222222 / t)); end tmp_2 = tmp; end
code[t_] := If[LessEqual[t, -0.8], N[(0.8333333333333334 - N[(0.2222222222222222 / t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.235], N[(N[(t * t), $MachinePrecision] + 0.5), $MachinePrecision], N[(0.8333333333333334 + N[(N[(0.037037037037037035 / N[(t * t), $MachinePrecision]), $MachinePrecision] + N[(-0.2222222222222222 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.8:\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222}{t}\\
\mathbf{elif}\;t \leq 0.235:\\
\;\;\;\;t \cdot t + 0.5\\
\mathbf{else}:\\
\;\;\;\;0.8333333333333334 + \left(\frac{0.037037037037037035}{t \cdot t} + \frac{-0.2222222222222222}{t}\right)\\
\end{array}
\end{array}
if t < -0.80000000000000004Initial program 100.0%
associate-/l*99.9%
associate-*r/99.9%
associate-/r/100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*r/100.0%
*-commutative100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*l*100.0%
metadata-eval100.0%
associate-/l*99.9%
Simplified100.0%
Taylor expanded in t around inf 99.9%
associate-*r/99.9%
metadata-eval99.9%
Simplified99.9%
if -0.80000000000000004 < t < 0.23499999999999999Initial program 100.0%
associate-/l*100.0%
associate-*r/100.0%
associate-/r/100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*r/100.0%
*-commutative100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*l*100.0%
metadata-eval100.0%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in t around 0 99.1%
+-commutative99.1%
unpow299.1%
Simplified99.1%
if 0.23499999999999999 < t Initial program 100.0%
associate-/l*99.9%
associate-*r/99.9%
associate-/r/100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*r/100.0%
*-commutative100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*l*100.0%
metadata-eval100.0%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in t around inf 98.8%
sub-neg98.8%
+-commutative98.8%
associate-+l+98.8%
associate-*r/98.8%
metadata-eval98.8%
unpow298.8%
associate-*r/98.8%
metadata-eval98.8%
distribute-neg-frac98.8%
metadata-eval98.8%
Simplified98.8%
Final simplification99.2%
(FPCore (t)
:precision binary64
(if (<= t -0.55)
(/ (- -5.0 (/ -8.0 t)) (- -6.0 (/ -8.0 t)))
(if (<= t 0.235)
(+ (* t t) 0.5)
(+
0.8333333333333334
(+ (/ 0.037037037037037035 (* t t)) (/ -0.2222222222222222 t))))))
double code(double t) {
double tmp;
if (t <= -0.55) {
tmp = (-5.0 - (-8.0 / t)) / (-6.0 - (-8.0 / t));
} else if (t <= 0.235) {
tmp = (t * t) + 0.5;
} else {
tmp = 0.8333333333333334 + ((0.037037037037037035 / (t * t)) + (-0.2222222222222222 / t));
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-0.55d0)) then
tmp = ((-5.0d0) - ((-8.0d0) / t)) / ((-6.0d0) - ((-8.0d0) / t))
else if (t <= 0.235d0) then
tmp = (t * t) + 0.5d0
else
tmp = 0.8333333333333334d0 + ((0.037037037037037035d0 / (t * t)) + ((-0.2222222222222222d0) / t))
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if (t <= -0.55) {
tmp = (-5.0 - (-8.0 / t)) / (-6.0 - (-8.0 / t));
} else if (t <= 0.235) {
tmp = (t * t) + 0.5;
} else {
tmp = 0.8333333333333334 + ((0.037037037037037035 / (t * t)) + (-0.2222222222222222 / t));
}
return tmp;
}
def code(t): tmp = 0 if t <= -0.55: tmp = (-5.0 - (-8.0 / t)) / (-6.0 - (-8.0 / t)) elif t <= 0.235: tmp = (t * t) + 0.5 else: tmp = 0.8333333333333334 + ((0.037037037037037035 / (t * t)) + (-0.2222222222222222 / t)) return tmp
function code(t) tmp = 0.0 if (t <= -0.55) tmp = Float64(Float64(-5.0 - Float64(-8.0 / t)) / Float64(-6.0 - Float64(-8.0 / t))); elseif (t <= 0.235) tmp = Float64(Float64(t * t) + 0.5); else tmp = Float64(0.8333333333333334 + Float64(Float64(0.037037037037037035 / Float64(t * t)) + Float64(-0.2222222222222222 / t))); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if (t <= -0.55) tmp = (-5.0 - (-8.0 / t)) / (-6.0 - (-8.0 / t)); elseif (t <= 0.235) tmp = (t * t) + 0.5; else tmp = 0.8333333333333334 + ((0.037037037037037035 / (t * t)) + (-0.2222222222222222 / t)); end tmp_2 = tmp; end
code[t_] := If[LessEqual[t, -0.55], N[(N[(-5.0 - N[(-8.0 / t), $MachinePrecision]), $MachinePrecision] / N[(-6.0 - N[(-8.0 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.235], N[(N[(t * t), $MachinePrecision] + 0.5), $MachinePrecision], N[(0.8333333333333334 + N[(N[(0.037037037037037035 / N[(t * t), $MachinePrecision]), $MachinePrecision] + N[(-0.2222222222222222 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.55:\\
\;\;\;\;\frac{-5 - \frac{-8}{t}}{-6 - \frac{-8}{t}}\\
\mathbf{elif}\;t \leq 0.235:\\
\;\;\;\;t \cdot t + 0.5\\
\mathbf{else}:\\
\;\;\;\;0.8333333333333334 + \left(\frac{0.037037037037037035}{t \cdot t} + \frac{-0.2222222222222222}{t}\right)\\
\end{array}
\end{array}
if t < -0.55000000000000004Initial program 100.0%
associate-/l*99.9%
associate-*r/99.9%
associate-/r/100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*r/100.0%
*-commutative100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*l*100.0%
metadata-eval100.0%
associate-/l*99.9%
Simplified100.0%
Taylor expanded in t around inf 100.0%
associate-*r/100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in t around inf 100.0%
associate-*r/100.0%
metadata-eval100.0%
Simplified100.0%
add-log-exp99.9%
*-un-lft-identity99.9%
log-prod99.9%
metadata-eval99.9%
add-log-exp100.0%
associate-+r-100.0%
sub-neg100.0%
metadata-eval100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
associate-+r-100.0%
sub-neg100.0%
metadata-eval100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Applied egg-rr100.0%
+-lft-identity100.0%
*-lft-identity100.0%
metadata-eval100.0%
times-frac100.0%
neg-mul-1100.0%
neg-mul-1100.0%
distribute-neg-in100.0%
unsub-neg100.0%
metadata-eval100.0%
distribute-neg-in100.0%
unsub-neg100.0%
metadata-eval100.0%
Simplified100.0%
if -0.55000000000000004 < t < 0.23499999999999999Initial program 100.0%
associate-/l*100.0%
associate-*r/100.0%
associate-/r/100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*r/100.0%
*-commutative100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*l*100.0%
metadata-eval100.0%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in t around 0 99.1%
+-commutative99.1%
unpow299.1%
Simplified99.1%
if 0.23499999999999999 < t Initial program 100.0%
associate-/l*99.9%
associate-*r/99.9%
associate-/r/100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*r/100.0%
*-commutative100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*l*100.0%
metadata-eval100.0%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in t around inf 98.8%
sub-neg98.8%
+-commutative98.8%
associate-+l+98.8%
associate-*r/98.8%
metadata-eval98.8%
unpow298.8%
associate-*r/98.8%
metadata-eval98.8%
distribute-neg-frac98.8%
metadata-eval98.8%
Simplified98.8%
Final simplification99.2%
(FPCore (t) :precision binary64 (if (or (<= t -0.8) (not (<= t 0.58))) (- 0.8333333333333334 (/ 0.2222222222222222 t)) (+ (* t t) 0.5)))
double code(double t) {
double tmp;
if ((t <= -0.8) || !(t <= 0.58)) {
tmp = 0.8333333333333334 - (0.2222222222222222 / t);
} else {
tmp = (t * t) + 0.5;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if ((t <= (-0.8d0)) .or. (.not. (t <= 0.58d0))) then
tmp = 0.8333333333333334d0 - (0.2222222222222222d0 / t)
else
tmp = (t * t) + 0.5d0
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if ((t <= -0.8) || !(t <= 0.58)) {
tmp = 0.8333333333333334 - (0.2222222222222222 / t);
} else {
tmp = (t * t) + 0.5;
}
return tmp;
}
def code(t): tmp = 0 if (t <= -0.8) or not (t <= 0.58): tmp = 0.8333333333333334 - (0.2222222222222222 / t) else: tmp = (t * t) + 0.5 return tmp
function code(t) tmp = 0.0 if ((t <= -0.8) || !(t <= 0.58)) tmp = Float64(0.8333333333333334 - Float64(0.2222222222222222 / t)); else tmp = Float64(Float64(t * t) + 0.5); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if ((t <= -0.8) || ~((t <= 0.58))) tmp = 0.8333333333333334 - (0.2222222222222222 / t); else tmp = (t * t) + 0.5; end tmp_2 = tmp; end
code[t_] := If[Or[LessEqual[t, -0.8], N[Not[LessEqual[t, 0.58]], $MachinePrecision]], N[(0.8333333333333334 - N[(0.2222222222222222 / t), $MachinePrecision]), $MachinePrecision], N[(N[(t * t), $MachinePrecision] + 0.5), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.8 \lor \neg \left(t \leq 0.58\right):\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222}{t}\\
\mathbf{else}:\\
\;\;\;\;t \cdot t + 0.5\\
\end{array}
\end{array}
if t < -0.80000000000000004 or 0.57999999999999996 < t Initial program 100.0%
associate-/l*99.9%
associate-*r/99.9%
associate-/r/100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*r/100.0%
*-commutative100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*l*100.0%
metadata-eval100.0%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in t around inf 99.2%
associate-*r/99.2%
metadata-eval99.2%
Simplified99.2%
if -0.80000000000000004 < t < 0.57999999999999996Initial program 100.0%
associate-/l*100.0%
associate-*r/100.0%
associate-/r/100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*r/100.0%
*-commutative100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*l*100.0%
metadata-eval100.0%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in t around 0 99.1%
+-commutative99.1%
unpow299.1%
Simplified99.1%
Final simplification99.1%
(FPCore (t) :precision binary64 (if (<= t -0.9) 0.8333333333333334 (if (<= t 0.58) (+ (* t t) 0.5) 0.8333333333333334)))
double code(double t) {
double tmp;
if (t <= -0.9) {
tmp = 0.8333333333333334;
} else if (t <= 0.58) {
tmp = (t * t) + 0.5;
} else {
tmp = 0.8333333333333334;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-0.9d0)) then
tmp = 0.8333333333333334d0
else if (t <= 0.58d0) then
tmp = (t * t) + 0.5d0
else
tmp = 0.8333333333333334d0
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if (t <= -0.9) {
tmp = 0.8333333333333334;
} else if (t <= 0.58) {
tmp = (t * t) + 0.5;
} else {
tmp = 0.8333333333333334;
}
return tmp;
}
def code(t): tmp = 0 if t <= -0.9: tmp = 0.8333333333333334 elif t <= 0.58: tmp = (t * t) + 0.5 else: tmp = 0.8333333333333334 return tmp
function code(t) tmp = 0.0 if (t <= -0.9) tmp = 0.8333333333333334; elseif (t <= 0.58) tmp = Float64(Float64(t * t) + 0.5); else tmp = 0.8333333333333334; end return tmp end
function tmp_2 = code(t) tmp = 0.0; if (t <= -0.9) tmp = 0.8333333333333334; elseif (t <= 0.58) tmp = (t * t) + 0.5; else tmp = 0.8333333333333334; end tmp_2 = tmp; end
code[t_] := If[LessEqual[t, -0.9], 0.8333333333333334, If[LessEqual[t, 0.58], N[(N[(t * t), $MachinePrecision] + 0.5), $MachinePrecision], 0.8333333333333334]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.9:\\
\;\;\;\;0.8333333333333334\\
\mathbf{elif}\;t \leq 0.58:\\
\;\;\;\;t \cdot t + 0.5\\
\mathbf{else}:\\
\;\;\;\;0.8333333333333334\\
\end{array}
\end{array}
if t < -0.900000000000000022 or 0.57999999999999996 < t Initial program 100.0%
associate-/l*99.9%
associate-*r/99.9%
associate-/r/100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*r/100.0%
*-commutative100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*l*100.0%
metadata-eval100.0%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in t around inf 98.7%
if -0.900000000000000022 < t < 0.57999999999999996Initial program 100.0%
associate-/l*100.0%
associate-*r/100.0%
associate-/r/100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*r/100.0%
*-commutative100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*l*100.0%
metadata-eval100.0%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in t around 0 99.1%
+-commutative99.1%
unpow299.1%
Simplified99.1%
Final simplification98.9%
(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*99.9%
associate-*r/99.9%
associate-/r/100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*r/100.0%
*-commutative100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*l*100.0%
metadata-eval100.0%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in t around inf 98.7%
if -0.340000000000000024 < t < 1Initial program 100.0%
associate-/l*100.0%
associate-*r/100.0%
associate-/r/100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*r/100.0%
*-commutative100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*l*100.0%
metadata-eval100.0%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in t around 0 98.8%
Final simplification98.8%
(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-*r/100.0%
associate-/r/100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*r/100.0%
*-commutative100.0%
associate-*l/100.0%
*-commutative100.0%
associate-*l*100.0%
metadata-eval100.0%
associate-/l*100.0%
Simplified100.0%
Taylor expanded in t around 0 55.5%
Final simplification55.5%
herbie shell --seed 2023221
(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))))))