
(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 9 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))))
(/
(- 1.0 (* 4.0 (* (/ t (+ -1.0 (+ t 2.0))) (/ t (- -1.0 t)))))
(+ 2.0 (* 4.0 (* t_1 t_1))))))
double code(double t) {
double t_1 = t / (t + 1.0);
return (1.0 - (4.0 * ((t / (-1.0 + (t + 2.0))) * (t / (-1.0 - t))))) / (2.0 + (4.0 * (t_1 * t_1)));
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: t_1
t_1 = t / (t + 1.0d0)
code = (1.0d0 - (4.0d0 * ((t / ((-1.0d0) + (t + 2.0d0))) * (t / ((-1.0d0) - t))))) / (2.0d0 + (4.0d0 * (t_1 * t_1)))
end function
public static double code(double t) {
double t_1 = t / (t + 1.0);
return (1.0 - (4.0 * ((t / (-1.0 + (t + 2.0))) * (t / (-1.0 - t))))) / (2.0 + (4.0 * (t_1 * t_1)));
}
def code(t): t_1 = t / (t + 1.0) return (1.0 - (4.0 * ((t / (-1.0 + (t + 2.0))) * (t / (-1.0 - t))))) / (2.0 + (4.0 * (t_1 * t_1)))
function code(t) t_1 = Float64(t / Float64(t + 1.0)) return Float64(Float64(1.0 - Float64(4.0 * Float64(Float64(t / Float64(-1.0 + Float64(t + 2.0))) * Float64(t / Float64(-1.0 - t))))) / Float64(2.0 + Float64(4.0 * Float64(t_1 * t_1)))) end
function tmp = code(t) t_1 = t / (t + 1.0); tmp = (1.0 - (4.0 * ((t / (-1.0 + (t + 2.0))) * (t / (-1.0 - t))))) / (2.0 + (4.0 * (t_1 * t_1))); end
code[t_] := Block[{t$95$1 = N[(t / N[(t + 1.0), $MachinePrecision]), $MachinePrecision]}, N[(N[(1.0 - N[(4.0 * N[(N[(t / N[(-1.0 + N[(t + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t / N[(-1.0 - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[(4.0 * N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{t + 1}\\
\frac{1 - 4 \cdot \left(\frac{t}{-1 + \left(t + 2\right)} \cdot \frac{t}{-1 - t}\right)}{2 + 4 \cdot \left(t\_1 \cdot t\_1\right)}
\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%
expm1-log1p-u73.5%
expm1-undefine73.5%
+-commutative73.5%
Applied egg-rr73.5%
sub-neg73.5%
metadata-eval73.5%
+-commutative73.5%
log1p-undefine73.5%
rem-exp-log100.0%
+-commutative100.0%
associate-+r+100.0%
metadata-eval100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (t) :precision binary64 (let* ((t_1 (/ t (+ t 1.0))) (t_2 (* 4.0 (* t_1 t_1)))) (/ (+ t_2 1.0) (+ 2.0 t_2))))
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) / (2.0 + t_2);
}
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) / (2.0d0 + t_2)
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) / (2.0 + t_2);
}
def code(t): t_1 = t / (t + 1.0) t_2 = 4.0 * (t_1 * t_1) return (t_2 + 1.0) / (2.0 + t_2)
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(2.0 + t_2)) end
function tmp = code(t) t_1 = t / (t + 1.0); t_2 = 4.0 * (t_1 * t_1); tmp = (t_2 + 1.0) / (2.0 + t_2); 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[(2.0 + t$95$2), $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}{2 + t\_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 (or (<= t -0.52) (not (<= t 0.58)))
(+
0.8333333333333334
(/
(-
(/ (+ 0.037037037037037035 (/ 0.04938271604938271 t)) t)
0.2222222222222222)
t))
(+ (* (* t t) (+ (* t -2.0) 1.0)) 0.5)))
double code(double t) {
double tmp;
if ((t <= -0.52) || !(t <= 0.58)) {
tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t);
} else {
tmp = ((t * t) * ((t * -2.0) + 1.0)) + 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.58d0))) then
tmp = 0.8333333333333334d0 + ((((0.037037037037037035d0 + (0.04938271604938271d0 / t)) / t) - 0.2222222222222222d0) / t)
else
tmp = ((t * t) * ((t * (-2.0d0)) + 1.0d0)) + 0.5d0
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if ((t <= -0.52) || !(t <= 0.58)) {
tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t);
} else {
tmp = ((t * t) * ((t * -2.0) + 1.0)) + 0.5;
}
return tmp;
}
def code(t): tmp = 0 if (t <= -0.52) or not (t <= 0.58): tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t) else: tmp = ((t * t) * ((t * -2.0) + 1.0)) + 0.5 return tmp
function code(t) tmp = 0.0 if ((t <= -0.52) || !(t <= 0.58)) tmp = Float64(0.8333333333333334 + Float64(Float64(Float64(Float64(0.037037037037037035 + Float64(0.04938271604938271 / t)) / t) - 0.2222222222222222) / t)); else tmp = Float64(Float64(Float64(t * t) * Float64(Float64(t * -2.0) + 1.0)) + 0.5); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if ((t <= -0.52) || ~((t <= 0.58))) tmp = 0.8333333333333334 + ((((0.037037037037037035 + (0.04938271604938271 / t)) / t) - 0.2222222222222222) / t); else tmp = ((t * t) * ((t * -2.0) + 1.0)) + 0.5; end tmp_2 = tmp; end
code[t_] := If[Or[LessEqual[t, -0.52], N[Not[LessEqual[t, 0.58]], $MachinePrecision]], N[(0.8333333333333334 + N[(N[(N[(N[(0.037037037037037035 + N[(0.04938271604938271 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] - 0.2222222222222222), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], N[(N[(N[(t * t), $MachinePrecision] * N[(N[(t * -2.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + 0.5), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.52 \lor \neg \left(t \leq 0.58\right):\\
\;\;\;\;0.8333333333333334 + \frac{\frac{0.037037037037037035 + \frac{0.04938271604938271}{t}}{t} - 0.2222222222222222}{t}\\
\mathbf{else}:\\
\;\;\;\;\left(t \cdot t\right) \cdot \left(t \cdot -2 + 1\right) + 0.5\\
\end{array}
\end{array}
if t < -0.52000000000000002 or 0.57999999999999996 < 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 99.1%
mul-1-neg99.1%
unsub-neg99.1%
mul-1-neg99.1%
unsub-neg99.1%
associate-*r/99.1%
metadata-eval99.1%
Simplified99.1%
if -0.52000000000000002 < t < 0.57999999999999996Initial 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 100.0%
+-commutative100.0%
*-commutative100.0%
Simplified100.0%
unpow2100.0%
Applied egg-rr100.0%
Final simplification99.5%
(FPCore (t)
:precision binary64
(if (or (<= t -0.6) (not (<= t 0.44)))
(-
0.8333333333333334
(/ (+ 0.2222222222222222 (/ -0.037037037037037035 t)) t))
(+ (* (* t t) (+ (* t -2.0) 1.0)) 0.5)))
double code(double t) {
double tmp;
if ((t <= -0.6) || !(t <= 0.44)) {
tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t);
} else {
tmp = ((t * t) * ((t * -2.0) + 1.0)) + 0.5;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if ((t <= (-0.6d0)) .or. (.not. (t <= 0.44d0))) then
tmp = 0.8333333333333334d0 - ((0.2222222222222222d0 + ((-0.037037037037037035d0) / t)) / t)
else
tmp = ((t * t) * ((t * (-2.0d0)) + 1.0d0)) + 0.5d0
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if ((t <= -0.6) || !(t <= 0.44)) {
tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t);
} else {
tmp = ((t * t) * ((t * -2.0) + 1.0)) + 0.5;
}
return tmp;
}
def code(t): tmp = 0 if (t <= -0.6) or not (t <= 0.44): tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t) else: tmp = ((t * t) * ((t * -2.0) + 1.0)) + 0.5 return tmp
function code(t) tmp = 0.0 if ((t <= -0.6) || !(t <= 0.44)) tmp = Float64(0.8333333333333334 - Float64(Float64(0.2222222222222222 + Float64(-0.037037037037037035 / t)) / t)); else tmp = Float64(Float64(Float64(t * t) * Float64(Float64(t * -2.0) + 1.0)) + 0.5); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if ((t <= -0.6) || ~((t <= 0.44))) tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t); else tmp = ((t * t) * ((t * -2.0) + 1.0)) + 0.5; end tmp_2 = tmp; end
code[t_] := If[Or[LessEqual[t, -0.6], N[Not[LessEqual[t, 0.44]], $MachinePrecision]], N[(0.8333333333333334 - N[(N[(0.2222222222222222 + N[(-0.037037037037037035 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], N[(N[(N[(t * t), $MachinePrecision] * N[(N[(t * -2.0), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + 0.5), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.6 \lor \neg \left(t \leq 0.44\right):\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222 + \frac{-0.037037037037037035}{t}}{t}\\
\mathbf{else}:\\
\;\;\;\;\left(t \cdot t\right) \cdot \left(t \cdot -2 + 1\right) + 0.5\\
\end{array}
\end{array}
if t < -0.599999999999999978 or 0.440000000000000002 < 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 99.5%
mul-1-neg99.5%
unsub-neg99.5%
sub-neg99.5%
associate-*r/99.5%
metadata-eval99.5%
distribute-neg-frac99.5%
metadata-eval99.5%
Simplified99.5%
if -0.599999999999999978 < t < 0.440000000000000002Initial 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.4%
+-commutative99.4%
*-commutative99.4%
Simplified99.4%
unpow299.4%
Applied egg-rr99.4%
Final simplification99.4%
(FPCore (t)
:precision binary64
(if (or (<= t -0.81) (not (<= t 0.235)))
(-
0.8333333333333334
(/ (+ 0.2222222222222222 (/ -0.037037037037037035 t)) t))
(+ (* t t) 0.5)))
double code(double t) {
double tmp;
if ((t <= -0.81) || !(t <= 0.235)) {
tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / 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.81d0)) .or. (.not. (t <= 0.235d0))) then
tmp = 0.8333333333333334d0 - ((0.2222222222222222d0 + ((-0.037037037037037035d0) / t)) / t)
else
tmp = (t * t) + 0.5d0
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if ((t <= -0.81) || !(t <= 0.235)) {
tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t);
} else {
tmp = (t * t) + 0.5;
}
return tmp;
}
def code(t): tmp = 0 if (t <= -0.81) or not (t <= 0.235): tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t) else: tmp = (t * t) + 0.5 return tmp
function code(t) tmp = 0.0 if ((t <= -0.81) || !(t <= 0.235)) tmp = Float64(0.8333333333333334 - Float64(Float64(0.2222222222222222 + Float64(-0.037037037037037035 / t)) / t)); else tmp = Float64(Float64(t * t) + 0.5); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if ((t <= -0.81) || ~((t <= 0.235))) tmp = 0.8333333333333334 - ((0.2222222222222222 + (-0.037037037037037035 / t)) / t); else tmp = (t * t) + 0.5; end tmp_2 = tmp; end
code[t_] := If[Or[LessEqual[t, -0.81], N[Not[LessEqual[t, 0.235]], $MachinePrecision]], N[(0.8333333333333334 - N[(N[(0.2222222222222222 + N[(-0.037037037037037035 / t), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], N[(N[(t * t), $MachinePrecision] + 0.5), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -0.81 \lor \neg \left(t \leq 0.235\right):\\
\;\;\;\;0.8333333333333334 - \frac{0.2222222222222222 + \frac{-0.037037037037037035}{t}}{t}\\
\mathbf{else}:\\
\;\;\;\;t \cdot t + 0.5\\
\end{array}
\end{array}
if t < -0.81000000000000005 or 0.23499999999999999 < 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 99.5%
mul-1-neg99.5%
unsub-neg99.5%
sub-neg99.5%
associate-*r/99.5%
metadata-eval99.5%
distribute-neg-frac99.5%
metadata-eval99.5%
Simplified99.5%
if -0.81000000000000005 < t < 0.23499999999999999Initial 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.4%
+-commutative99.4%
*-commutative99.4%
Simplified99.4%
unpow299.4%
Applied egg-rr99.4%
Taylor expanded in t around 0 99.2%
Final simplification99.3%
(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*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 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-/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.4%
+-commutative99.4%
*-commutative99.4%
Simplified99.4%
unpow299.4%
Applied egg-rr99.4%
Taylor expanded in t around 0 99.2%
Final simplification99.2%
(FPCore (t) :precision binary64 (if (or (<= t -0.49) (not (<= t 0.68))) (- 0.8333333333333334 (/ 0.2222222222222222 t)) 0.5))
double code(double t) {
double tmp;
if ((t <= -0.49) || !(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.49d0)) .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.49) || !(t <= 0.68)) {
tmp = 0.8333333333333334 - (0.2222222222222222 / t);
} else {
tmp = 0.5;
}
return tmp;
}
def code(t): tmp = 0 if (t <= -0.49) 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.49) || !(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.49) || ~((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.49], 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.49 \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.48999999999999999 or 0.680000000000000049 < 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%
associate-*r/98.6%
metadata-eval98.6%
Simplified98.6%
if -0.48999999999999999 < t < 0.680000000000000049Initial 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.3%
Final simplification98.9%
(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%
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 96.9%
if -0.330000000000000016 < 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.3%
(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 58.3%
herbie shell --seed 2024132
(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))))))