
(FPCore (x y z t) :precision binary64 (+ (- x (/ y (* z 3.0))) (/ t (* (* z 3.0) y))))
double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x - (y / (z * 3.0d0))) + (t / ((z * 3.0d0) * y))
end function
public static double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y));
}
def code(x, y, z, t): return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y))
function code(x, y, z, t) return Float64(Float64(x - Float64(y / Float64(z * 3.0))) + Float64(t / Float64(Float64(z * 3.0) * y))) end
function tmp = code(x, y, z, t) tmp = (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y)); end
code[x_, y_, z_, t_] := N[(N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t / N[(N[(z * 3.0), $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x - \frac{y}{z \cdot 3}\right) + \frac{t}{\left(z \cdot 3\right) \cdot y}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 18 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t) :precision binary64 (+ (- x (/ y (* z 3.0))) (/ t (* (* z 3.0) y))))
double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x - (y / (z * 3.0d0))) + (t / ((z * 3.0d0) * y))
end function
public static double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y));
}
def code(x, y, z, t): return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y))
function code(x, y, z, t) return Float64(Float64(x - Float64(y / Float64(z * 3.0))) + Float64(t / Float64(Float64(z * 3.0) * y))) end
function tmp = code(x, y, z, t) tmp = (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y)); end
code[x_, y_, z_, t_] := N[(N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t / N[(N[(z * 3.0), $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x - \frac{y}{z \cdot 3}\right) + \frac{t}{\left(z \cdot 3\right) \cdot y}
\end{array}
(FPCore (x y z t) :precision binary64 (+ (+ (/ t (* z (* y 3.0))) x) (/ y (* z -3.0))))
double code(double x, double y, double z, double t) {
return ((t / (z * (y * 3.0))) + x) + (y / (z * -3.0));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = ((t / (z * (y * 3.0d0))) + x) + (y / (z * (-3.0d0)))
end function
public static double code(double x, double y, double z, double t) {
return ((t / (z * (y * 3.0))) + x) + (y / (z * -3.0));
}
def code(x, y, z, t): return ((t / (z * (y * 3.0))) + x) + (y / (z * -3.0))
function code(x, y, z, t) return Float64(Float64(Float64(t / Float64(z * Float64(y * 3.0))) + x) + Float64(y / Float64(z * -3.0))) end
function tmp = code(x, y, z, t) tmp = ((t / (z * (y * 3.0))) + x) + (y / (z * -3.0)); end
code[x_, y_, z_, t_] := N[(N[(N[(t / N[(z * N[(y * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision] + N[(y / N[(z * -3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\frac{t}{z \cdot \left(y \cdot 3\right)} + x\right) + \frac{y}{z \cdot -3}
\end{array}
Initial program 97.1%
+-commutative97.1%
associate-+r-97.1%
sub-neg97.1%
associate-*l*97.4%
*-commutative97.4%
distribute-frac-neg297.4%
distribute-rgt-neg-in97.4%
metadata-eval97.4%
Simplified97.4%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (* 0.3333333333333333 (/ (/ t z) y)))
(t_2 (+ x (/ (* y -0.3333333333333333) z))))
(if (<= y -5.6e-99)
(- x (* 0.3333333333333333 (/ y z)))
(if (<= y -1.06e-113)
t_1
(if (<= y -2.3e-156)
t_2
(if (<= y 3.2e-5)
t_1
(if (or (<= y 6.2e+27) (not (<= y 4.4e+49)))
t_2
(/ t (* y (* z 3.0))))))))))
double code(double x, double y, double z, double t) {
double t_1 = 0.3333333333333333 * ((t / z) / y);
double t_2 = x + ((y * -0.3333333333333333) / z);
double tmp;
if (y <= -5.6e-99) {
tmp = x - (0.3333333333333333 * (y / z));
} else if (y <= -1.06e-113) {
tmp = t_1;
} else if (y <= -2.3e-156) {
tmp = t_2;
} else if (y <= 3.2e-5) {
tmp = t_1;
} else if ((y <= 6.2e+27) || !(y <= 4.4e+49)) {
tmp = t_2;
} else {
tmp = t / (y * (z * 3.0));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = 0.3333333333333333d0 * ((t / z) / y)
t_2 = x + ((y * (-0.3333333333333333d0)) / z)
if (y <= (-5.6d-99)) then
tmp = x - (0.3333333333333333d0 * (y / z))
else if (y <= (-1.06d-113)) then
tmp = t_1
else if (y <= (-2.3d-156)) then
tmp = t_2
else if (y <= 3.2d-5) then
tmp = t_1
else if ((y <= 6.2d+27) .or. (.not. (y <= 4.4d+49))) then
tmp = t_2
else
tmp = t / (y * (z * 3.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = 0.3333333333333333 * ((t / z) / y);
double t_2 = x + ((y * -0.3333333333333333) / z);
double tmp;
if (y <= -5.6e-99) {
tmp = x - (0.3333333333333333 * (y / z));
} else if (y <= -1.06e-113) {
tmp = t_1;
} else if (y <= -2.3e-156) {
tmp = t_2;
} else if (y <= 3.2e-5) {
tmp = t_1;
} else if ((y <= 6.2e+27) || !(y <= 4.4e+49)) {
tmp = t_2;
} else {
tmp = t / (y * (z * 3.0));
}
return tmp;
}
def code(x, y, z, t): t_1 = 0.3333333333333333 * ((t / z) / y) t_2 = x + ((y * -0.3333333333333333) / z) tmp = 0 if y <= -5.6e-99: tmp = x - (0.3333333333333333 * (y / z)) elif y <= -1.06e-113: tmp = t_1 elif y <= -2.3e-156: tmp = t_2 elif y <= 3.2e-5: tmp = t_1 elif (y <= 6.2e+27) or not (y <= 4.4e+49): tmp = t_2 else: tmp = t / (y * (z * 3.0)) return tmp
function code(x, y, z, t) t_1 = Float64(0.3333333333333333 * Float64(Float64(t / z) / y)) t_2 = Float64(x + Float64(Float64(y * -0.3333333333333333) / z)) tmp = 0.0 if (y <= -5.6e-99) tmp = Float64(x - Float64(0.3333333333333333 * Float64(y / z))); elseif (y <= -1.06e-113) tmp = t_1; elseif (y <= -2.3e-156) tmp = t_2; elseif (y <= 3.2e-5) tmp = t_1; elseif ((y <= 6.2e+27) || !(y <= 4.4e+49)) tmp = t_2; else tmp = Float64(t / Float64(y * Float64(z * 3.0))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = 0.3333333333333333 * ((t / z) / y); t_2 = x + ((y * -0.3333333333333333) / z); tmp = 0.0; if (y <= -5.6e-99) tmp = x - (0.3333333333333333 * (y / z)); elseif (y <= -1.06e-113) tmp = t_1; elseif (y <= -2.3e-156) tmp = t_2; elseif (y <= 3.2e-5) tmp = t_1; elseif ((y <= 6.2e+27) || ~((y <= 4.4e+49))) tmp = t_2; else tmp = t / (y * (z * 3.0)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(0.3333333333333333 * N[(N[(t / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -5.6e-99], N[(x - N[(0.3333333333333333 * N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -1.06e-113], t$95$1, If[LessEqual[y, -2.3e-156], t$95$2, If[LessEqual[y, 3.2e-5], t$95$1, If[Or[LessEqual[y, 6.2e+27], N[Not[LessEqual[y, 4.4e+49]], $MachinePrecision]], t$95$2, N[(t / N[(y * N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 0.3333333333333333 \cdot \frac{\frac{t}{z}}{y}\\
t_2 := x + \frac{y \cdot -0.3333333333333333}{z}\\
\mathbf{if}\;y \leq -5.6 \cdot 10^{-99}:\\
\;\;\;\;x - 0.3333333333333333 \cdot \frac{y}{z}\\
\mathbf{elif}\;y \leq -1.06 \cdot 10^{-113}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -2.3 \cdot 10^{-156}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq 3.2 \cdot 10^{-5}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 6.2 \cdot 10^{+27} \lor \neg \left(y \leq 4.4 \cdot 10^{+49}\right):\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;\frac{t}{y \cdot \left(z \cdot 3\right)}\\
\end{array}
\end{array}
if y < -5.6000000000000001e-99Initial program 97.6%
Taylor expanded in t around 0 85.2%
if -5.6000000000000001e-99 < y < -1.05999999999999995e-113 or -2.3e-156 < y < 3.19999999999999986e-5Initial program 95.6%
+-commutative95.6%
associate-+r-95.6%
sub-neg95.6%
associate-*l*96.6%
*-commutative96.6%
distribute-frac-neg296.6%
distribute-rgt-neg-in96.6%
metadata-eval96.6%
Simplified96.6%
*-un-lft-identity96.6%
times-frac96.6%
Applied egg-rr96.6%
*-commutative96.6%
associate-*l/96.6%
associate-*r/96.6%
*-rgt-identity96.6%
Simplified96.6%
Taylor expanded in t around inf 71.0%
*-commutative71.0%
associate-/r*71.2%
Simplified71.2%
if -1.05999999999999995e-113 < y < -2.3e-156 or 3.19999999999999986e-5 < y < 6.19999999999999992e27 or 4.4000000000000001e49 < y Initial program 99.9%
sub-neg99.9%
associate-+l+99.9%
+-commutative99.9%
remove-double-neg99.9%
distribute-frac-neg99.9%
distribute-neg-in99.9%
remove-double-neg99.9%
sub-neg99.9%
neg-mul-199.9%
times-frac99.8%
distribute-frac-neg99.8%
neg-mul-199.8%
*-commutative99.8%
associate-/l*99.8%
*-commutative99.8%
Simplified99.7%
Taylor expanded in z around 0 99.8%
associate-*r/99.8%
distribute-lft-out--99.8%
clear-num99.7%
*-un-lft-identity99.7%
distribute-lft-out--99.7%
times-frac99.9%
metadata-eval99.9%
Applied egg-rr99.9%
Taylor expanded in t around 0 94.1%
associate-*r/94.1%
*-commutative94.1%
Simplified94.1%
if 6.19999999999999992e27 < y < 4.4000000000000001e49Initial program 85.7%
+-commutative85.7%
associate-+r-85.7%
sub-neg85.7%
associate-*l*85.5%
*-commutative85.5%
distribute-frac-neg285.5%
distribute-rgt-neg-in85.5%
metadata-eval85.5%
Simplified85.5%
*-un-lft-identity85.5%
times-frac85.3%
Applied egg-rr85.3%
*-commutative85.3%
associate-*l/85.5%
associate-*r/85.5%
*-rgt-identity85.5%
Simplified85.5%
Taylor expanded in t around inf 79.0%
*-commutative79.0%
associate-/r*79.0%
Simplified79.0%
associate-/l/79.0%
associate-*r/78.9%
associate-*l/79.0%
associate-/l/79.0%
*-commutative79.0%
clear-num79.0%
un-div-inv79.0%
div-inv79.2%
clear-num79.0%
div-inv79.2%
metadata-eval79.2%
Applied egg-rr79.2%
Final simplification81.8%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (* 0.3333333333333333 (/ (/ t z) y)))
(t_2 (+ x (/ (* y -0.3333333333333333) z))))
(if (<= y -7e-99)
(- x (* 0.3333333333333333 (/ y z)))
(if (<= y -3.3e-114)
t_1
(if (<= y -1.85e-156)
t_2
(if (<= y 3.1e-5)
t_1
(if (or (<= y 6.2e+27) (not (<= y 4.4e+49)))
t_2
(* 0.3333333333333333 (/ t (* z y))))))))))
double code(double x, double y, double z, double t) {
double t_1 = 0.3333333333333333 * ((t / z) / y);
double t_2 = x + ((y * -0.3333333333333333) / z);
double tmp;
if (y <= -7e-99) {
tmp = x - (0.3333333333333333 * (y / z));
} else if (y <= -3.3e-114) {
tmp = t_1;
} else if (y <= -1.85e-156) {
tmp = t_2;
} else if (y <= 3.1e-5) {
tmp = t_1;
} else if ((y <= 6.2e+27) || !(y <= 4.4e+49)) {
tmp = t_2;
} else {
tmp = 0.3333333333333333 * (t / (z * y));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = 0.3333333333333333d0 * ((t / z) / y)
t_2 = x + ((y * (-0.3333333333333333d0)) / z)
if (y <= (-7d-99)) then
tmp = x - (0.3333333333333333d0 * (y / z))
else if (y <= (-3.3d-114)) then
tmp = t_1
else if (y <= (-1.85d-156)) then
tmp = t_2
else if (y <= 3.1d-5) then
tmp = t_1
else if ((y <= 6.2d+27) .or. (.not. (y <= 4.4d+49))) then
tmp = t_2
else
tmp = 0.3333333333333333d0 * (t / (z * y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = 0.3333333333333333 * ((t / z) / y);
double t_2 = x + ((y * -0.3333333333333333) / z);
double tmp;
if (y <= -7e-99) {
tmp = x - (0.3333333333333333 * (y / z));
} else if (y <= -3.3e-114) {
tmp = t_1;
} else if (y <= -1.85e-156) {
tmp = t_2;
} else if (y <= 3.1e-5) {
tmp = t_1;
} else if ((y <= 6.2e+27) || !(y <= 4.4e+49)) {
tmp = t_2;
} else {
tmp = 0.3333333333333333 * (t / (z * y));
}
return tmp;
}
def code(x, y, z, t): t_1 = 0.3333333333333333 * ((t / z) / y) t_2 = x + ((y * -0.3333333333333333) / z) tmp = 0 if y <= -7e-99: tmp = x - (0.3333333333333333 * (y / z)) elif y <= -3.3e-114: tmp = t_1 elif y <= -1.85e-156: tmp = t_2 elif y <= 3.1e-5: tmp = t_1 elif (y <= 6.2e+27) or not (y <= 4.4e+49): tmp = t_2 else: tmp = 0.3333333333333333 * (t / (z * y)) return tmp
function code(x, y, z, t) t_1 = Float64(0.3333333333333333 * Float64(Float64(t / z) / y)) t_2 = Float64(x + Float64(Float64(y * -0.3333333333333333) / z)) tmp = 0.0 if (y <= -7e-99) tmp = Float64(x - Float64(0.3333333333333333 * Float64(y / z))); elseif (y <= -3.3e-114) tmp = t_1; elseif (y <= -1.85e-156) tmp = t_2; elseif (y <= 3.1e-5) tmp = t_1; elseif ((y <= 6.2e+27) || !(y <= 4.4e+49)) tmp = t_2; else tmp = Float64(0.3333333333333333 * Float64(t / Float64(z * y))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = 0.3333333333333333 * ((t / z) / y); t_2 = x + ((y * -0.3333333333333333) / z); tmp = 0.0; if (y <= -7e-99) tmp = x - (0.3333333333333333 * (y / z)); elseif (y <= -3.3e-114) tmp = t_1; elseif (y <= -1.85e-156) tmp = t_2; elseif (y <= 3.1e-5) tmp = t_1; elseif ((y <= 6.2e+27) || ~((y <= 4.4e+49))) tmp = t_2; else tmp = 0.3333333333333333 * (t / (z * y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(0.3333333333333333 * N[(N[(t / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -7e-99], N[(x - N[(0.3333333333333333 * N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -3.3e-114], t$95$1, If[LessEqual[y, -1.85e-156], t$95$2, If[LessEqual[y, 3.1e-5], t$95$1, If[Or[LessEqual[y, 6.2e+27], N[Not[LessEqual[y, 4.4e+49]], $MachinePrecision]], t$95$2, N[(0.3333333333333333 * N[(t / N[(z * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 0.3333333333333333 \cdot \frac{\frac{t}{z}}{y}\\
t_2 := x + \frac{y \cdot -0.3333333333333333}{z}\\
\mathbf{if}\;y \leq -7 \cdot 10^{-99}:\\
\;\;\;\;x - 0.3333333333333333 \cdot \frac{y}{z}\\
\mathbf{elif}\;y \leq -3.3 \cdot 10^{-114}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -1.85 \cdot 10^{-156}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq 3.1 \cdot 10^{-5}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 6.2 \cdot 10^{+27} \lor \neg \left(y \leq 4.4 \cdot 10^{+49}\right):\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{t}{z \cdot y}\\
\end{array}
\end{array}
if y < -6.9999999999999997e-99Initial program 97.6%
Taylor expanded in t around 0 85.2%
if -6.9999999999999997e-99 < y < -3.30000000000000035e-114 or -1.85e-156 < y < 3.10000000000000014e-5Initial program 95.6%
+-commutative95.6%
associate-+r-95.6%
sub-neg95.6%
associate-*l*96.6%
*-commutative96.6%
distribute-frac-neg296.6%
distribute-rgt-neg-in96.6%
metadata-eval96.6%
Simplified96.6%
*-un-lft-identity96.6%
times-frac96.6%
Applied egg-rr96.6%
*-commutative96.6%
associate-*l/96.6%
associate-*r/96.6%
*-rgt-identity96.6%
Simplified96.6%
Taylor expanded in t around inf 71.0%
*-commutative71.0%
associate-/r*71.2%
Simplified71.2%
if -3.30000000000000035e-114 < y < -1.85e-156 or 3.10000000000000014e-5 < y < 6.19999999999999992e27 or 4.4000000000000001e49 < y Initial program 99.9%
sub-neg99.9%
associate-+l+99.9%
+-commutative99.9%
remove-double-neg99.9%
distribute-frac-neg99.9%
distribute-neg-in99.9%
remove-double-neg99.9%
sub-neg99.9%
neg-mul-199.9%
times-frac99.8%
distribute-frac-neg99.8%
neg-mul-199.8%
*-commutative99.8%
associate-/l*99.8%
*-commutative99.8%
Simplified99.7%
Taylor expanded in z around 0 99.8%
associate-*r/99.8%
distribute-lft-out--99.8%
clear-num99.7%
*-un-lft-identity99.7%
distribute-lft-out--99.7%
times-frac99.9%
metadata-eval99.9%
Applied egg-rr99.9%
Taylor expanded in t around 0 94.1%
associate-*r/94.1%
*-commutative94.1%
Simplified94.1%
if 6.19999999999999992e27 < y < 4.4000000000000001e49Initial program 85.7%
+-commutative85.7%
associate-+r-85.7%
sub-neg85.7%
associate-*l*85.5%
*-commutative85.5%
distribute-frac-neg285.5%
distribute-rgt-neg-in85.5%
metadata-eval85.5%
Simplified85.5%
*-un-lft-identity85.5%
times-frac85.3%
Applied egg-rr85.3%
*-commutative85.3%
associate-*l/85.5%
associate-*r/85.5%
*-rgt-identity85.5%
Simplified85.5%
Taylor expanded in t around inf 79.0%
Final simplification81.8%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (* 0.3333333333333333 (/ (/ t z) y)))
(t_2 (+ x (/ (* y -0.3333333333333333) z))))
(if (<= y -2.4e-99)
t_2
(if (<= y -1.65e-114)
t_1
(if (<= y -4.2e-158)
t_2
(if (<= y 3.1e-5)
t_1
(if (or (<= y 6.2e+27) (not (<= y 4.4e+49)))
t_2
(* 0.3333333333333333 (/ t (* z y))))))))))
double code(double x, double y, double z, double t) {
double t_1 = 0.3333333333333333 * ((t / z) / y);
double t_2 = x + ((y * -0.3333333333333333) / z);
double tmp;
if (y <= -2.4e-99) {
tmp = t_2;
} else if (y <= -1.65e-114) {
tmp = t_1;
} else if (y <= -4.2e-158) {
tmp = t_2;
} else if (y <= 3.1e-5) {
tmp = t_1;
} else if ((y <= 6.2e+27) || !(y <= 4.4e+49)) {
tmp = t_2;
} else {
tmp = 0.3333333333333333 * (t / (z * y));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = 0.3333333333333333d0 * ((t / z) / y)
t_2 = x + ((y * (-0.3333333333333333d0)) / z)
if (y <= (-2.4d-99)) then
tmp = t_2
else if (y <= (-1.65d-114)) then
tmp = t_1
else if (y <= (-4.2d-158)) then
tmp = t_2
else if (y <= 3.1d-5) then
tmp = t_1
else if ((y <= 6.2d+27) .or. (.not. (y <= 4.4d+49))) then
tmp = t_2
else
tmp = 0.3333333333333333d0 * (t / (z * y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = 0.3333333333333333 * ((t / z) / y);
double t_2 = x + ((y * -0.3333333333333333) / z);
double tmp;
if (y <= -2.4e-99) {
tmp = t_2;
} else if (y <= -1.65e-114) {
tmp = t_1;
} else if (y <= -4.2e-158) {
tmp = t_2;
} else if (y <= 3.1e-5) {
tmp = t_1;
} else if ((y <= 6.2e+27) || !(y <= 4.4e+49)) {
tmp = t_2;
} else {
tmp = 0.3333333333333333 * (t / (z * y));
}
return tmp;
}
def code(x, y, z, t): t_1 = 0.3333333333333333 * ((t / z) / y) t_2 = x + ((y * -0.3333333333333333) / z) tmp = 0 if y <= -2.4e-99: tmp = t_2 elif y <= -1.65e-114: tmp = t_1 elif y <= -4.2e-158: tmp = t_2 elif y <= 3.1e-5: tmp = t_1 elif (y <= 6.2e+27) or not (y <= 4.4e+49): tmp = t_2 else: tmp = 0.3333333333333333 * (t / (z * y)) return tmp
function code(x, y, z, t) t_1 = Float64(0.3333333333333333 * Float64(Float64(t / z) / y)) t_2 = Float64(x + Float64(Float64(y * -0.3333333333333333) / z)) tmp = 0.0 if (y <= -2.4e-99) tmp = t_2; elseif (y <= -1.65e-114) tmp = t_1; elseif (y <= -4.2e-158) tmp = t_2; elseif (y <= 3.1e-5) tmp = t_1; elseif ((y <= 6.2e+27) || !(y <= 4.4e+49)) tmp = t_2; else tmp = Float64(0.3333333333333333 * Float64(t / Float64(z * y))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = 0.3333333333333333 * ((t / z) / y); t_2 = x + ((y * -0.3333333333333333) / z); tmp = 0.0; if (y <= -2.4e-99) tmp = t_2; elseif (y <= -1.65e-114) tmp = t_1; elseif (y <= -4.2e-158) tmp = t_2; elseif (y <= 3.1e-5) tmp = t_1; elseif ((y <= 6.2e+27) || ~((y <= 4.4e+49))) tmp = t_2; else tmp = 0.3333333333333333 * (t / (z * y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(0.3333333333333333 * N[(N[(t / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.4e-99], t$95$2, If[LessEqual[y, -1.65e-114], t$95$1, If[LessEqual[y, -4.2e-158], t$95$2, If[LessEqual[y, 3.1e-5], t$95$1, If[Or[LessEqual[y, 6.2e+27], N[Not[LessEqual[y, 4.4e+49]], $MachinePrecision]], t$95$2, N[(0.3333333333333333 * N[(t / N[(z * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 0.3333333333333333 \cdot \frac{\frac{t}{z}}{y}\\
t_2 := x + \frac{y \cdot -0.3333333333333333}{z}\\
\mathbf{if}\;y \leq -2.4 \cdot 10^{-99}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -1.65 \cdot 10^{-114}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -4.2 \cdot 10^{-158}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq 3.1 \cdot 10^{-5}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 6.2 \cdot 10^{+27} \lor \neg \left(y \leq 4.4 \cdot 10^{+49}\right):\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{t}{z \cdot y}\\
\end{array}
\end{array}
if y < -2.4e-99 or -1.65000000000000017e-114 < y < -4.19999999999999983e-158 or 3.10000000000000014e-5 < y < 6.19999999999999992e27 or 4.4000000000000001e49 < y Initial program 98.5%
sub-neg98.5%
associate-+l+98.5%
+-commutative98.5%
remove-double-neg98.5%
distribute-frac-neg98.5%
distribute-neg-in98.5%
remove-double-neg98.5%
sub-neg98.5%
neg-mul-198.5%
times-frac99.1%
distribute-frac-neg99.1%
neg-mul-199.1%
*-commutative99.1%
associate-/l*99.1%
*-commutative99.1%
Simplified99.0%
Taylor expanded in z around 0 99.1%
associate-*r/99.1%
distribute-lft-out--99.1%
clear-num99.1%
*-un-lft-identity99.1%
distribute-lft-out--99.1%
times-frac99.2%
metadata-eval99.2%
Applied egg-rr99.2%
Taylor expanded in t around 0 88.9%
associate-*r/88.9%
*-commutative88.9%
Simplified88.9%
if -2.4e-99 < y < -1.65000000000000017e-114 or -4.19999999999999983e-158 < y < 3.10000000000000014e-5Initial program 95.6%
+-commutative95.6%
associate-+r-95.6%
sub-neg95.6%
associate-*l*96.6%
*-commutative96.6%
distribute-frac-neg296.6%
distribute-rgt-neg-in96.6%
metadata-eval96.6%
Simplified96.6%
*-un-lft-identity96.6%
times-frac96.6%
Applied egg-rr96.6%
*-commutative96.6%
associate-*l/96.6%
associate-*r/96.6%
*-rgt-identity96.6%
Simplified96.6%
Taylor expanded in t around inf 71.0%
*-commutative71.0%
associate-/r*71.2%
Simplified71.2%
if 6.19999999999999992e27 < y < 4.4000000000000001e49Initial program 85.7%
+-commutative85.7%
associate-+r-85.7%
sub-neg85.7%
associate-*l*85.5%
*-commutative85.5%
distribute-frac-neg285.5%
distribute-rgt-neg-in85.5%
metadata-eval85.5%
Simplified85.5%
*-un-lft-identity85.5%
times-frac85.3%
Applied egg-rr85.3%
*-commutative85.3%
associate-*l/85.5%
associate-*r/85.5%
*-rgt-identity85.5%
Simplified85.5%
Taylor expanded in t around inf 79.0%
Final simplification81.8%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (* 0.3333333333333333 (/ (/ t z) y))))
(if (<= y -7.2e+80)
(/ (/ y -3.0) z)
(if (<= y -9.2e-99)
x
(if (<= y -1.1e-114)
t_1
(if (<= y -2.2e-157) x (if (<= y 5e+53) t_1 (/ y (* z -3.0)))))))))
double code(double x, double y, double z, double t) {
double t_1 = 0.3333333333333333 * ((t / z) / y);
double tmp;
if (y <= -7.2e+80) {
tmp = (y / -3.0) / z;
} else if (y <= -9.2e-99) {
tmp = x;
} else if (y <= -1.1e-114) {
tmp = t_1;
} else if (y <= -2.2e-157) {
tmp = x;
} else if (y <= 5e+53) {
tmp = t_1;
} else {
tmp = y / (z * -3.0);
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = 0.3333333333333333d0 * ((t / z) / y)
if (y <= (-7.2d+80)) then
tmp = (y / (-3.0d0)) / z
else if (y <= (-9.2d-99)) then
tmp = x
else if (y <= (-1.1d-114)) then
tmp = t_1
else if (y <= (-2.2d-157)) then
tmp = x
else if (y <= 5d+53) then
tmp = t_1
else
tmp = y / (z * (-3.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = 0.3333333333333333 * ((t / z) / y);
double tmp;
if (y <= -7.2e+80) {
tmp = (y / -3.0) / z;
} else if (y <= -9.2e-99) {
tmp = x;
} else if (y <= -1.1e-114) {
tmp = t_1;
} else if (y <= -2.2e-157) {
tmp = x;
} else if (y <= 5e+53) {
tmp = t_1;
} else {
tmp = y / (z * -3.0);
}
return tmp;
}
def code(x, y, z, t): t_1 = 0.3333333333333333 * ((t / z) / y) tmp = 0 if y <= -7.2e+80: tmp = (y / -3.0) / z elif y <= -9.2e-99: tmp = x elif y <= -1.1e-114: tmp = t_1 elif y <= -2.2e-157: tmp = x elif y <= 5e+53: tmp = t_1 else: tmp = y / (z * -3.0) return tmp
function code(x, y, z, t) t_1 = Float64(0.3333333333333333 * Float64(Float64(t / z) / y)) tmp = 0.0 if (y <= -7.2e+80) tmp = Float64(Float64(y / -3.0) / z); elseif (y <= -9.2e-99) tmp = x; elseif (y <= -1.1e-114) tmp = t_1; elseif (y <= -2.2e-157) tmp = x; elseif (y <= 5e+53) tmp = t_1; else tmp = Float64(y / Float64(z * -3.0)); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = 0.3333333333333333 * ((t / z) / y); tmp = 0.0; if (y <= -7.2e+80) tmp = (y / -3.0) / z; elseif (y <= -9.2e-99) tmp = x; elseif (y <= -1.1e-114) tmp = t_1; elseif (y <= -2.2e-157) tmp = x; elseif (y <= 5e+53) tmp = t_1; else tmp = y / (z * -3.0); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(0.3333333333333333 * N[(N[(t / z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -7.2e+80], N[(N[(y / -3.0), $MachinePrecision] / z), $MachinePrecision], If[LessEqual[y, -9.2e-99], x, If[LessEqual[y, -1.1e-114], t$95$1, If[LessEqual[y, -2.2e-157], x, If[LessEqual[y, 5e+53], t$95$1, N[(y / N[(z * -3.0), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 0.3333333333333333 \cdot \frac{\frac{t}{z}}{y}\\
\mathbf{if}\;y \leq -7.2 \cdot 10^{+80}:\\
\;\;\;\;\frac{\frac{y}{-3}}{z}\\
\mathbf{elif}\;y \leq -9.2 \cdot 10^{-99}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq -1.1 \cdot 10^{-114}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -2.2 \cdot 10^{-157}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 5 \cdot 10^{+53}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{z \cdot -3}\\
\end{array}
\end{array}
if y < -7.1999999999999999e80Initial program 97.7%
+-commutative97.7%
associate-+r-97.7%
sub-neg97.7%
associate-*l*97.7%
*-commutative97.7%
distribute-frac-neg297.7%
distribute-rgt-neg-in97.7%
metadata-eval97.7%
Simplified97.7%
*-un-lft-identity97.7%
times-frac97.8%
Applied egg-rr97.8%
*-commutative97.8%
associate-*l/97.6%
associate-*r/97.6%
*-rgt-identity97.6%
Simplified97.6%
Taylor expanded in y around inf 67.3%
*-commutative67.3%
associate-*l/67.3%
associate-*r/67.4%
Simplified67.4%
clear-num67.4%
un-div-inv67.4%
div-inv67.3%
metadata-eval67.3%
Applied egg-rr67.3%
*-commutative67.3%
associate-/r*67.4%
Simplified67.4%
if -7.1999999999999999e80 < y < -9.1999999999999994e-99 or -1.10000000000000006e-114 < y < -2.2000000000000001e-157Initial program 98.0%
+-commutative98.0%
associate-+r-98.0%
sub-neg98.0%
associate-*l*98.0%
*-commutative98.0%
distribute-frac-neg298.0%
distribute-rgt-neg-in98.0%
metadata-eval98.0%
Simplified98.0%
Taylor expanded in z around inf 52.3%
if -9.1999999999999994e-99 < y < -1.10000000000000006e-114 or -2.2000000000000001e-157 < y < 5.0000000000000004e53Initial program 95.2%
+-commutative95.2%
associate-+r-95.2%
sub-neg95.2%
associate-*l*96.1%
*-commutative96.1%
distribute-frac-neg296.1%
distribute-rgt-neg-in96.1%
metadata-eval96.1%
Simplified96.1%
*-un-lft-identity96.1%
times-frac96.1%
Applied egg-rr96.1%
*-commutative96.1%
associate-*l/96.1%
associate-*r/96.1%
*-rgt-identity96.1%
Simplified96.1%
Taylor expanded in t around inf 67.9%
*-commutative67.9%
associate-/r*68.1%
Simplified68.1%
if 5.0000000000000004e53 < y Initial program 99.8%
+-commutative99.8%
associate-+r-99.8%
sub-neg99.8%
associate-*l*99.8%
*-commutative99.8%
distribute-frac-neg299.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
*-un-lft-identity99.8%
times-frac99.7%
Applied egg-rr99.7%
*-commutative99.7%
associate-*l/99.7%
associate-*r/99.8%
*-rgt-identity99.8%
Simplified99.8%
Taylor expanded in y around inf 72.4%
*-commutative72.4%
associate-*l/72.4%
associate-*r/72.3%
Simplified72.3%
clear-num72.3%
un-div-inv72.4%
div-inv72.4%
metadata-eval72.4%
Applied egg-rr72.4%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (* 0.3333333333333333 (/ t (* z y)))))
(if (<= y -1.85e+82)
(/ (/ y -3.0) z)
(if (<= y -7.8e-99)
x
(if (<= y -1.38e-114)
t_1
(if (<= y -3.2e-156) x (if (<= y 1.7e+54) t_1 (/ y (* z -3.0)))))))))
double code(double x, double y, double z, double t) {
double t_1 = 0.3333333333333333 * (t / (z * y));
double tmp;
if (y <= -1.85e+82) {
tmp = (y / -3.0) / z;
} else if (y <= -7.8e-99) {
tmp = x;
} else if (y <= -1.38e-114) {
tmp = t_1;
} else if (y <= -3.2e-156) {
tmp = x;
} else if (y <= 1.7e+54) {
tmp = t_1;
} else {
tmp = y / (z * -3.0);
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = 0.3333333333333333d0 * (t / (z * y))
if (y <= (-1.85d+82)) then
tmp = (y / (-3.0d0)) / z
else if (y <= (-7.8d-99)) then
tmp = x
else if (y <= (-1.38d-114)) then
tmp = t_1
else if (y <= (-3.2d-156)) then
tmp = x
else if (y <= 1.7d+54) then
tmp = t_1
else
tmp = y / (z * (-3.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = 0.3333333333333333 * (t / (z * y));
double tmp;
if (y <= -1.85e+82) {
tmp = (y / -3.0) / z;
} else if (y <= -7.8e-99) {
tmp = x;
} else if (y <= -1.38e-114) {
tmp = t_1;
} else if (y <= -3.2e-156) {
tmp = x;
} else if (y <= 1.7e+54) {
tmp = t_1;
} else {
tmp = y / (z * -3.0);
}
return tmp;
}
def code(x, y, z, t): t_1 = 0.3333333333333333 * (t / (z * y)) tmp = 0 if y <= -1.85e+82: tmp = (y / -3.0) / z elif y <= -7.8e-99: tmp = x elif y <= -1.38e-114: tmp = t_1 elif y <= -3.2e-156: tmp = x elif y <= 1.7e+54: tmp = t_1 else: tmp = y / (z * -3.0) return tmp
function code(x, y, z, t) t_1 = Float64(0.3333333333333333 * Float64(t / Float64(z * y))) tmp = 0.0 if (y <= -1.85e+82) tmp = Float64(Float64(y / -3.0) / z); elseif (y <= -7.8e-99) tmp = x; elseif (y <= -1.38e-114) tmp = t_1; elseif (y <= -3.2e-156) tmp = x; elseif (y <= 1.7e+54) tmp = t_1; else tmp = Float64(y / Float64(z * -3.0)); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = 0.3333333333333333 * (t / (z * y)); tmp = 0.0; if (y <= -1.85e+82) tmp = (y / -3.0) / z; elseif (y <= -7.8e-99) tmp = x; elseif (y <= -1.38e-114) tmp = t_1; elseif (y <= -3.2e-156) tmp = x; elseif (y <= 1.7e+54) tmp = t_1; else tmp = y / (z * -3.0); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(0.3333333333333333 * N[(t / N[(z * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.85e+82], N[(N[(y / -3.0), $MachinePrecision] / z), $MachinePrecision], If[LessEqual[y, -7.8e-99], x, If[LessEqual[y, -1.38e-114], t$95$1, If[LessEqual[y, -3.2e-156], x, If[LessEqual[y, 1.7e+54], t$95$1, N[(y / N[(z * -3.0), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 0.3333333333333333 \cdot \frac{t}{z \cdot y}\\
\mathbf{if}\;y \leq -1.85 \cdot 10^{+82}:\\
\;\;\;\;\frac{\frac{y}{-3}}{z}\\
\mathbf{elif}\;y \leq -7.8 \cdot 10^{-99}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq -1.38 \cdot 10^{-114}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -3.2 \cdot 10^{-156}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 1.7 \cdot 10^{+54}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{z \cdot -3}\\
\end{array}
\end{array}
if y < -1.8500000000000001e82Initial program 97.7%
+-commutative97.7%
associate-+r-97.7%
sub-neg97.7%
associate-*l*97.7%
*-commutative97.7%
distribute-frac-neg297.7%
distribute-rgt-neg-in97.7%
metadata-eval97.7%
Simplified97.7%
*-un-lft-identity97.7%
times-frac97.8%
Applied egg-rr97.8%
*-commutative97.8%
associate-*l/97.6%
associate-*r/97.6%
*-rgt-identity97.6%
Simplified97.6%
Taylor expanded in y around inf 67.3%
*-commutative67.3%
associate-*l/67.3%
associate-*r/67.4%
Simplified67.4%
clear-num67.4%
un-div-inv67.4%
div-inv67.3%
metadata-eval67.3%
Applied egg-rr67.3%
*-commutative67.3%
associate-/r*67.4%
Simplified67.4%
if -1.8500000000000001e82 < y < -7.79999999999999975e-99 or -1.38e-114 < y < -3.19999999999999982e-156Initial program 98.0%
+-commutative98.0%
associate-+r-98.0%
sub-neg98.0%
associate-*l*98.0%
*-commutative98.0%
distribute-frac-neg298.0%
distribute-rgt-neg-in98.0%
metadata-eval98.0%
Simplified98.0%
Taylor expanded in z around inf 52.3%
if -7.79999999999999975e-99 < y < -1.38e-114 or -3.19999999999999982e-156 < y < 1.7e54Initial program 95.2%
+-commutative95.2%
associate-+r-95.2%
sub-neg95.2%
associate-*l*96.1%
*-commutative96.1%
distribute-frac-neg296.1%
distribute-rgt-neg-in96.1%
metadata-eval96.1%
Simplified96.1%
*-un-lft-identity96.1%
times-frac96.1%
Applied egg-rr96.1%
*-commutative96.1%
associate-*l/96.1%
associate-*r/96.1%
*-rgt-identity96.1%
Simplified96.1%
Taylor expanded in t around inf 67.9%
if 1.7e54 < y Initial program 99.8%
+-commutative99.8%
associate-+r-99.8%
sub-neg99.8%
associate-*l*99.8%
*-commutative99.8%
distribute-frac-neg299.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
*-un-lft-identity99.8%
times-frac99.7%
Applied egg-rr99.7%
*-commutative99.7%
associate-*l/99.7%
associate-*r/99.8%
*-rgt-identity99.8%
Simplified99.8%
Taylor expanded in y around inf 72.4%
*-commutative72.4%
associate-*l/72.4%
associate-*r/72.3%
Simplified72.3%
clear-num72.3%
un-div-inv72.4%
div-inv72.4%
metadata-eval72.4%
Applied egg-rr72.4%
Final simplification65.5%
(FPCore (x y z t)
:precision binary64
(if (<= y -1.65e+24)
(- x (* 0.3333333333333333 (/ y z)))
(if (<= y 4.4e+49)
(+ x (* 0.3333333333333333 (* (/ 1.0 y) (/ t z))))
(+ x (/ (* y -0.3333333333333333) z)))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -1.65e+24) {
tmp = x - (0.3333333333333333 * (y / z));
} else if (y <= 4.4e+49) {
tmp = x + (0.3333333333333333 * ((1.0 / y) * (t / z)));
} else {
tmp = x + ((y * -0.3333333333333333) / z);
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (y <= (-1.65d+24)) then
tmp = x - (0.3333333333333333d0 * (y / z))
else if (y <= 4.4d+49) then
tmp = x + (0.3333333333333333d0 * ((1.0d0 / y) * (t / z)))
else
tmp = x + ((y * (-0.3333333333333333d0)) / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -1.65e+24) {
tmp = x - (0.3333333333333333 * (y / z));
} else if (y <= 4.4e+49) {
tmp = x + (0.3333333333333333 * ((1.0 / y) * (t / z)));
} else {
tmp = x + ((y * -0.3333333333333333) / z);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -1.65e+24: tmp = x - (0.3333333333333333 * (y / z)) elif y <= 4.4e+49: tmp = x + (0.3333333333333333 * ((1.0 / y) * (t / z))) else: tmp = x + ((y * -0.3333333333333333) / z) return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -1.65e+24) tmp = Float64(x - Float64(0.3333333333333333 * Float64(y / z))); elseif (y <= 4.4e+49) tmp = Float64(x + Float64(0.3333333333333333 * Float64(Float64(1.0 / y) * Float64(t / z)))); else tmp = Float64(x + Float64(Float64(y * -0.3333333333333333) / z)); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -1.65e+24) tmp = x - (0.3333333333333333 * (y / z)); elseif (y <= 4.4e+49) tmp = x + (0.3333333333333333 * ((1.0 / y) * (t / z))); else tmp = x + ((y * -0.3333333333333333) / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -1.65e+24], N[(x - N[(0.3333333333333333 * N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.4e+49], N[(x + N[(0.3333333333333333 * N[(N[(1.0 / y), $MachinePrecision] * N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.65 \cdot 10^{+24}:\\
\;\;\;\;x - 0.3333333333333333 \cdot \frac{y}{z}\\
\mathbf{elif}\;y \leq 4.4 \cdot 10^{+49}:\\
\;\;\;\;x + 0.3333333333333333 \cdot \left(\frac{1}{y} \cdot \frac{t}{z}\right)\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y \cdot -0.3333333333333333}{z}\\
\end{array}
\end{array}
if y < -1.6499999999999999e24Initial program 98.2%
Taylor expanded in t around 0 91.9%
if -1.6499999999999999e24 < y < 4.4000000000000001e49Initial program 95.7%
sub-neg95.7%
associate-+l+95.7%
+-commutative95.7%
remove-double-neg95.7%
distribute-frac-neg95.7%
distribute-neg-in95.7%
remove-double-neg95.7%
sub-neg95.7%
neg-mul-195.7%
times-frac94.0%
distribute-frac-neg94.0%
neg-mul-194.0%
*-commutative94.0%
associate-/l*93.9%
*-commutative93.9%
Simplified94.5%
Taylor expanded in t around inf 90.9%
*-un-lft-identity90.9%
times-frac91.0%
Applied egg-rr91.0%
if 4.4000000000000001e49 < y Initial program 99.8%
sub-neg99.8%
associate-+l+99.8%
+-commutative99.8%
remove-double-neg99.8%
distribute-frac-neg99.8%
distribute-neg-in99.8%
remove-double-neg99.8%
sub-neg99.8%
neg-mul-199.8%
times-frac99.8%
distribute-frac-neg99.8%
neg-mul-199.8%
*-commutative99.8%
associate-/l*99.7%
*-commutative99.7%
Simplified99.7%
Taylor expanded in z around 0 99.7%
associate-*r/99.7%
distribute-lft-out--99.7%
clear-num99.7%
*-un-lft-identity99.7%
distribute-lft-out--99.7%
times-frac99.8%
metadata-eval99.8%
Applied egg-rr99.8%
Taylor expanded in t around 0 99.7%
associate-*r/99.7%
*-commutative99.7%
Simplified99.7%
(FPCore (x y z t)
:precision binary64
(if (<= y -5.8e+24)
(- x (* 0.3333333333333333 (/ y z)))
(if (<= y 4.6e+50)
(+ x (* 0.3333333333333333 (/ t (* z y))))
(+ x (/ (* y -0.3333333333333333) z)))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -5.8e+24) {
tmp = x - (0.3333333333333333 * (y / z));
} else if (y <= 4.6e+50) {
tmp = x + (0.3333333333333333 * (t / (z * y)));
} else {
tmp = x + ((y * -0.3333333333333333) / z);
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (y <= (-5.8d+24)) then
tmp = x - (0.3333333333333333d0 * (y / z))
else if (y <= 4.6d+50) then
tmp = x + (0.3333333333333333d0 * (t / (z * y)))
else
tmp = x + ((y * (-0.3333333333333333d0)) / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -5.8e+24) {
tmp = x - (0.3333333333333333 * (y / z));
} else if (y <= 4.6e+50) {
tmp = x + (0.3333333333333333 * (t / (z * y)));
} else {
tmp = x + ((y * -0.3333333333333333) / z);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -5.8e+24: tmp = x - (0.3333333333333333 * (y / z)) elif y <= 4.6e+50: tmp = x + (0.3333333333333333 * (t / (z * y))) else: tmp = x + ((y * -0.3333333333333333) / z) return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -5.8e+24) tmp = Float64(x - Float64(0.3333333333333333 * Float64(y / z))); elseif (y <= 4.6e+50) tmp = Float64(x + Float64(0.3333333333333333 * Float64(t / Float64(z * y)))); else tmp = Float64(x + Float64(Float64(y * -0.3333333333333333) / z)); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -5.8e+24) tmp = x - (0.3333333333333333 * (y / z)); elseif (y <= 4.6e+50) tmp = x + (0.3333333333333333 * (t / (z * y))); else tmp = x + ((y * -0.3333333333333333) / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -5.8e+24], N[(x - N[(0.3333333333333333 * N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.6e+50], N[(x + N[(0.3333333333333333 * N[(t / N[(z * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5.8 \cdot 10^{+24}:\\
\;\;\;\;x - 0.3333333333333333 \cdot \frac{y}{z}\\
\mathbf{elif}\;y \leq 4.6 \cdot 10^{+50}:\\
\;\;\;\;x + 0.3333333333333333 \cdot \frac{t}{z \cdot y}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y \cdot -0.3333333333333333}{z}\\
\end{array}
\end{array}
if y < -5.79999999999999958e24Initial program 98.2%
Taylor expanded in t around 0 91.9%
if -5.79999999999999958e24 < y < 4.59999999999999994e50Initial program 95.7%
sub-neg95.7%
associate-+l+95.7%
+-commutative95.7%
remove-double-neg95.7%
distribute-frac-neg95.7%
distribute-neg-in95.7%
remove-double-neg95.7%
sub-neg95.7%
neg-mul-195.7%
times-frac94.0%
distribute-frac-neg94.0%
neg-mul-194.0%
*-commutative94.0%
associate-/l*93.9%
*-commutative93.9%
Simplified94.5%
Taylor expanded in t around inf 90.9%
if 4.59999999999999994e50 < y Initial program 99.8%
sub-neg99.8%
associate-+l+99.8%
+-commutative99.8%
remove-double-neg99.8%
distribute-frac-neg99.8%
distribute-neg-in99.8%
remove-double-neg99.8%
sub-neg99.8%
neg-mul-199.8%
times-frac99.8%
distribute-frac-neg99.8%
neg-mul-199.8%
*-commutative99.8%
associate-/l*99.7%
*-commutative99.7%
Simplified99.7%
Taylor expanded in z around 0 99.7%
associate-*r/99.7%
distribute-lft-out--99.7%
clear-num99.7%
*-un-lft-identity99.7%
distribute-lft-out--99.7%
times-frac99.8%
metadata-eval99.8%
Applied egg-rr99.8%
Taylor expanded in t around 0 99.7%
associate-*r/99.7%
*-commutative99.7%
Simplified99.7%
Final simplification92.8%
(FPCore (x y z t)
:precision binary64
(if (<= x -1.25e+112)
(- x (* y (/ 0.3333333333333333 z)))
(if (<= x 1.38e+181)
(* 0.3333333333333333 (/ (- (/ t y) y) z))
(+ x (/ (* y -0.3333333333333333) z)))))
double code(double x, double y, double z, double t) {
double tmp;
if (x <= -1.25e+112) {
tmp = x - (y * (0.3333333333333333 / z));
} else if (x <= 1.38e+181) {
tmp = 0.3333333333333333 * (((t / y) - y) / z);
} else {
tmp = x + ((y * -0.3333333333333333) / z);
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (x <= (-1.25d+112)) then
tmp = x - (y * (0.3333333333333333d0 / z))
else if (x <= 1.38d+181) then
tmp = 0.3333333333333333d0 * (((t / y) - y) / z)
else
tmp = x + ((y * (-0.3333333333333333d0)) / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (x <= -1.25e+112) {
tmp = x - (y * (0.3333333333333333 / z));
} else if (x <= 1.38e+181) {
tmp = 0.3333333333333333 * (((t / y) - y) / z);
} else {
tmp = x + ((y * -0.3333333333333333) / z);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if x <= -1.25e+112: tmp = x - (y * (0.3333333333333333 / z)) elif x <= 1.38e+181: tmp = 0.3333333333333333 * (((t / y) - y) / z) else: tmp = x + ((y * -0.3333333333333333) / z) return tmp
function code(x, y, z, t) tmp = 0.0 if (x <= -1.25e+112) tmp = Float64(x - Float64(y * Float64(0.3333333333333333 / z))); elseif (x <= 1.38e+181) tmp = Float64(0.3333333333333333 * Float64(Float64(Float64(t / y) - y) / z)); else tmp = Float64(x + Float64(Float64(y * -0.3333333333333333) / z)); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (x <= -1.25e+112) tmp = x - (y * (0.3333333333333333 / z)); elseif (x <= 1.38e+181) tmp = 0.3333333333333333 * (((t / y) - y) / z); else tmp = x + ((y * -0.3333333333333333) / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[x, -1.25e+112], N[(x - N[(y * N[(0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.38e+181], N[(0.3333333333333333 * N[(N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.25 \cdot 10^{+112}:\\
\;\;\;\;x - y \cdot \frac{0.3333333333333333}{z}\\
\mathbf{elif}\;x \leq 1.38 \cdot 10^{+181}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{\frac{t}{y} - y}{z}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y \cdot -0.3333333333333333}{z}\\
\end{array}
\end{array}
if x < -1.25e112Initial program 99.9%
sub-neg99.9%
associate-+l+99.9%
+-commutative99.9%
remove-double-neg99.9%
distribute-frac-neg99.9%
distribute-neg-in99.9%
remove-double-neg99.9%
sub-neg99.9%
neg-mul-199.9%
times-frac97.8%
distribute-frac-neg97.8%
neg-mul-197.8%
*-commutative97.8%
associate-/l*97.8%
*-commutative97.8%
Simplified97.8%
Taylor expanded in t around 0 95.2%
neg-mul-195.2%
Simplified95.2%
if -1.25e112 < x < 1.3799999999999999e181Initial program 96.2%
Taylor expanded in z around 0 77.3%
distribute-lft-out--77.3%
associate-*r/77.3%
*-commutative77.3%
Applied egg-rr77.3%
if 1.3799999999999999e181 < x Initial program 100.0%
sub-neg100.0%
associate-+l+100.0%
+-commutative100.0%
remove-double-neg100.0%
distribute-frac-neg100.0%
distribute-neg-in100.0%
remove-double-neg100.0%
sub-neg100.0%
neg-mul-1100.0%
times-frac94.8%
distribute-frac-neg94.8%
neg-mul-194.8%
*-commutative94.8%
associate-/l*94.7%
*-commutative94.7%
Simplified94.7%
Taylor expanded in z around 0 94.8%
associate-*r/94.8%
distribute-lft-out--94.8%
clear-num94.8%
*-un-lft-identity94.8%
distribute-lft-out--94.8%
times-frac94.8%
metadata-eval94.8%
Applied egg-rr94.8%
Taylor expanded in t around 0 82.4%
associate-*r/82.4%
*-commutative82.4%
Simplified82.4%
Final simplification80.6%
(FPCore (x y z t)
:precision binary64
(if (<= x -1.25e+112)
(- x (* y (/ 0.3333333333333333 z)))
(if (<= x 1.38e+181)
(* (- (/ t y) y) (/ 0.3333333333333333 z))
(+ x (/ (* y -0.3333333333333333) z)))))
double code(double x, double y, double z, double t) {
double tmp;
if (x <= -1.25e+112) {
tmp = x - (y * (0.3333333333333333 / z));
} else if (x <= 1.38e+181) {
tmp = ((t / y) - y) * (0.3333333333333333 / z);
} else {
tmp = x + ((y * -0.3333333333333333) / z);
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (x <= (-1.25d+112)) then
tmp = x - (y * (0.3333333333333333d0 / z))
else if (x <= 1.38d+181) then
tmp = ((t / y) - y) * (0.3333333333333333d0 / z)
else
tmp = x + ((y * (-0.3333333333333333d0)) / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (x <= -1.25e+112) {
tmp = x - (y * (0.3333333333333333 / z));
} else if (x <= 1.38e+181) {
tmp = ((t / y) - y) * (0.3333333333333333 / z);
} else {
tmp = x + ((y * -0.3333333333333333) / z);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if x <= -1.25e+112: tmp = x - (y * (0.3333333333333333 / z)) elif x <= 1.38e+181: tmp = ((t / y) - y) * (0.3333333333333333 / z) else: tmp = x + ((y * -0.3333333333333333) / z) return tmp
function code(x, y, z, t) tmp = 0.0 if (x <= -1.25e+112) tmp = Float64(x - Float64(y * Float64(0.3333333333333333 / z))); elseif (x <= 1.38e+181) tmp = Float64(Float64(Float64(t / y) - y) * Float64(0.3333333333333333 / z)); else tmp = Float64(x + Float64(Float64(y * -0.3333333333333333) / z)); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (x <= -1.25e+112) tmp = x - (y * (0.3333333333333333 / z)); elseif (x <= 1.38e+181) tmp = ((t / y) - y) * (0.3333333333333333 / z); else tmp = x + ((y * -0.3333333333333333) / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[x, -1.25e+112], N[(x - N[(y * N[(0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.38e+181], N[(N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision] * N[(0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.25 \cdot 10^{+112}:\\
\;\;\;\;x - y \cdot \frac{0.3333333333333333}{z}\\
\mathbf{elif}\;x \leq 1.38 \cdot 10^{+181}:\\
\;\;\;\;\left(\frac{t}{y} - y\right) \cdot \frac{0.3333333333333333}{z}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y \cdot -0.3333333333333333}{z}\\
\end{array}
\end{array}
if x < -1.25e112Initial program 99.9%
sub-neg99.9%
associate-+l+99.9%
+-commutative99.9%
remove-double-neg99.9%
distribute-frac-neg99.9%
distribute-neg-in99.9%
remove-double-neg99.9%
sub-neg99.9%
neg-mul-199.9%
times-frac97.8%
distribute-frac-neg97.8%
neg-mul-197.8%
*-commutative97.8%
associate-/l*97.8%
*-commutative97.8%
Simplified97.8%
Taylor expanded in t around 0 95.2%
neg-mul-195.2%
Simplified95.2%
if -1.25e112 < x < 1.3799999999999999e181Initial program 96.2%
+-commutative96.2%
associate-+r-96.2%
sub-neg96.2%
associate-*l*96.6%
*-commutative96.6%
distribute-frac-neg296.6%
distribute-rgt-neg-in96.6%
metadata-eval96.6%
Simplified96.6%
*-un-lft-identity96.6%
times-frac96.7%
Applied egg-rr96.7%
*-commutative96.7%
associate-*l/96.6%
associate-*r/96.6%
*-rgt-identity96.6%
Simplified96.6%
Taylor expanded in x around 0 76.6%
metadata-eval76.6%
distribute-lft-neg-in76.6%
+-commutative76.6%
associate-/r*76.8%
fma-define76.8%
fma-neg76.8%
distribute-lft-out--76.8%
div-sub77.3%
associate-*r/77.3%
associate-*l/77.2%
Simplified77.2%
if 1.3799999999999999e181 < x Initial program 100.0%
sub-neg100.0%
associate-+l+100.0%
+-commutative100.0%
remove-double-neg100.0%
distribute-frac-neg100.0%
distribute-neg-in100.0%
remove-double-neg100.0%
sub-neg100.0%
neg-mul-1100.0%
times-frac94.8%
distribute-frac-neg94.8%
neg-mul-194.8%
*-commutative94.8%
associate-/l*94.7%
*-commutative94.7%
Simplified94.7%
Taylor expanded in z around 0 94.8%
associate-*r/94.8%
distribute-lft-out--94.8%
clear-num94.8%
*-un-lft-identity94.8%
distribute-lft-out--94.8%
times-frac94.8%
metadata-eval94.8%
Applied egg-rr94.8%
Taylor expanded in t around 0 82.4%
associate-*r/82.4%
*-commutative82.4%
Simplified82.4%
Final simplification80.5%
(FPCore (x y z t) :precision binary64 (if (or (<= y -5.5e+80) (not (<= y 1.52e+81))) (/ (/ y -3.0) z) x))
double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -5.5e+80) || !(y <= 1.52e+81)) {
tmp = (y / -3.0) / z;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if ((y <= (-5.5d+80)) .or. (.not. (y <= 1.52d+81))) then
tmp = (y / (-3.0d0)) / z
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -5.5e+80) || !(y <= 1.52e+81)) {
tmp = (y / -3.0) / z;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (y <= -5.5e+80) or not (y <= 1.52e+81): tmp = (y / -3.0) / z else: tmp = x return tmp
function code(x, y, z, t) tmp = 0.0 if ((y <= -5.5e+80) || !(y <= 1.52e+81)) tmp = Float64(Float64(y / -3.0) / z); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((y <= -5.5e+80) || ~((y <= 1.52e+81))) tmp = (y / -3.0) / z; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[y, -5.5e+80], N[Not[LessEqual[y, 1.52e+81]], $MachinePrecision]], N[(N[(y / -3.0), $MachinePrecision] / z), $MachinePrecision], x]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5.5 \cdot 10^{+80} \lor \neg \left(y \leq 1.52 \cdot 10^{+81}\right):\\
\;\;\;\;\frac{\frac{y}{-3}}{z}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -5.49999999999999967e80 or 1.51999999999999999e81 < y Initial program 98.6%
+-commutative98.6%
associate-+r-98.6%
sub-neg98.6%
associate-*l*98.6%
*-commutative98.6%
distribute-frac-neg298.6%
distribute-rgt-neg-in98.6%
metadata-eval98.6%
Simplified98.6%
*-un-lft-identity98.6%
times-frac98.6%
Applied egg-rr98.6%
*-commutative98.6%
associate-*l/98.6%
associate-*r/98.6%
*-rgt-identity98.6%
Simplified98.6%
Taylor expanded in y around inf 72.9%
*-commutative72.9%
associate-*l/72.9%
associate-*r/72.9%
Simplified72.9%
clear-num72.9%
un-div-inv72.9%
div-inv72.9%
metadata-eval72.9%
Applied egg-rr72.9%
*-commutative72.9%
associate-/r*73.0%
Simplified73.0%
if -5.49999999999999967e80 < y < 1.51999999999999999e81Initial program 96.3%
+-commutative96.3%
associate-+r-96.3%
sub-neg96.3%
associate-*l*96.8%
*-commutative96.8%
distribute-frac-neg296.8%
distribute-rgt-neg-in96.8%
metadata-eval96.8%
Simplified96.8%
Taylor expanded in z around inf 35.7%
Final simplification48.2%
(FPCore (x y z t) :precision binary64 (if (or (<= y -5.1e+82) (not (<= y 3.7e+80))) (* (/ y z) -0.3333333333333333) x))
double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -5.1e+82) || !(y <= 3.7e+80)) {
tmp = (y / z) * -0.3333333333333333;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if ((y <= (-5.1d+82)) .or. (.not. (y <= 3.7d+80))) then
tmp = (y / z) * (-0.3333333333333333d0)
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -5.1e+82) || !(y <= 3.7e+80)) {
tmp = (y / z) * -0.3333333333333333;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (y <= -5.1e+82) or not (y <= 3.7e+80): tmp = (y / z) * -0.3333333333333333 else: tmp = x return tmp
function code(x, y, z, t) tmp = 0.0 if ((y <= -5.1e+82) || !(y <= 3.7e+80)) tmp = Float64(Float64(y / z) * -0.3333333333333333); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((y <= -5.1e+82) || ~((y <= 3.7e+80))) tmp = (y / z) * -0.3333333333333333; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[y, -5.1e+82], N[Not[LessEqual[y, 3.7e+80]], $MachinePrecision]], N[(N[(y / z), $MachinePrecision] * -0.3333333333333333), $MachinePrecision], x]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5.1 \cdot 10^{+82} \lor \neg \left(y \leq 3.7 \cdot 10^{+80}\right):\\
\;\;\;\;\frac{y}{z} \cdot -0.3333333333333333\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -5.1000000000000003e82 or 3.69999999999999996e80 < y Initial program 98.6%
+-commutative98.6%
associate-+r-98.6%
sub-neg98.6%
associate-*l*98.6%
*-commutative98.6%
distribute-frac-neg298.6%
distribute-rgt-neg-in98.6%
metadata-eval98.6%
Simplified98.6%
*-un-lft-identity98.6%
times-frac98.6%
Applied egg-rr98.6%
*-commutative98.6%
associate-*l/98.6%
associate-*r/98.6%
*-rgt-identity98.6%
Simplified98.6%
Taylor expanded in y around inf 72.9%
if -5.1000000000000003e82 < y < 3.69999999999999996e80Initial program 96.3%
+-commutative96.3%
associate-+r-96.3%
sub-neg96.3%
associate-*l*96.8%
*-commutative96.8%
distribute-frac-neg296.8%
distribute-rgt-neg-in96.8%
metadata-eval96.8%
Simplified96.8%
Taylor expanded in z around inf 35.7%
Final simplification48.2%
(FPCore (x y z t) :precision binary64 (if (<= y -8e+86) (* y (/ -0.3333333333333333 z)) (if (<= y 1.32e+82) x (/ y (* z -3.0)))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -8e+86) {
tmp = y * (-0.3333333333333333 / z);
} else if (y <= 1.32e+82) {
tmp = x;
} else {
tmp = y / (z * -3.0);
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (y <= (-8d+86)) then
tmp = y * ((-0.3333333333333333d0) / z)
else if (y <= 1.32d+82) then
tmp = x
else
tmp = y / (z * (-3.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -8e+86) {
tmp = y * (-0.3333333333333333 / z);
} else if (y <= 1.32e+82) {
tmp = x;
} else {
tmp = y / (z * -3.0);
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -8e+86: tmp = y * (-0.3333333333333333 / z) elif y <= 1.32e+82: tmp = x else: tmp = y / (z * -3.0) return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -8e+86) tmp = Float64(y * Float64(-0.3333333333333333 / z)); elseif (y <= 1.32e+82) tmp = x; else tmp = Float64(y / Float64(z * -3.0)); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -8e+86) tmp = y * (-0.3333333333333333 / z); elseif (y <= 1.32e+82) tmp = x; else tmp = y / (z * -3.0); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -8e+86], N[(y * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.32e+82], x, N[(y / N[(z * -3.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -8 \cdot 10^{+86}:\\
\;\;\;\;y \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 1.32 \cdot 10^{+82}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{z \cdot -3}\\
\end{array}
\end{array}
if y < -8.0000000000000001e86Initial program 97.7%
+-commutative97.7%
associate-+r-97.7%
sub-neg97.7%
associate-*l*97.7%
*-commutative97.7%
distribute-frac-neg297.7%
distribute-rgt-neg-in97.7%
metadata-eval97.7%
Simplified97.7%
*-un-lft-identity97.7%
times-frac97.8%
Applied egg-rr97.8%
*-commutative97.8%
associate-*l/97.6%
associate-*r/97.6%
*-rgt-identity97.6%
Simplified97.6%
Taylor expanded in y around inf 67.3%
*-commutative67.3%
associate-*l/67.3%
associate-*r/67.4%
Simplified67.4%
if -8.0000000000000001e86 < y < 1.32e82Initial program 96.3%
+-commutative96.3%
associate-+r-96.3%
sub-neg96.3%
associate-*l*96.8%
*-commutative96.8%
distribute-frac-neg296.8%
distribute-rgt-neg-in96.8%
metadata-eval96.8%
Simplified96.8%
Taylor expanded in z around inf 35.7%
if 1.32e82 < y Initial program 99.8%
+-commutative99.8%
associate-+r-99.8%
sub-neg99.8%
associate-*l*99.8%
*-commutative99.8%
distribute-frac-neg299.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
*-un-lft-identity99.8%
times-frac99.6%
Applied egg-rr99.6%
*-commutative99.6%
associate-*l/99.7%
associate-*r/99.8%
*-rgt-identity99.8%
Simplified99.8%
Taylor expanded in y around inf 79.6%
*-commutative79.6%
associate-*l/79.6%
associate-*r/79.5%
Simplified79.5%
clear-num79.5%
un-div-inv79.7%
div-inv79.7%
metadata-eval79.7%
Applied egg-rr79.7%
(FPCore (x y z t) :precision binary64 (if (<= y -1.85e+80) (* y (/ -0.3333333333333333 z)) (if (<= y 4.5e+83) x (* (/ y z) -0.3333333333333333))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -1.85e+80) {
tmp = y * (-0.3333333333333333 / z);
} else if (y <= 4.5e+83) {
tmp = x;
} else {
tmp = (y / z) * -0.3333333333333333;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (y <= (-1.85d+80)) then
tmp = y * ((-0.3333333333333333d0) / z)
else if (y <= 4.5d+83) then
tmp = x
else
tmp = (y / z) * (-0.3333333333333333d0)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -1.85e+80) {
tmp = y * (-0.3333333333333333 / z);
} else if (y <= 4.5e+83) {
tmp = x;
} else {
tmp = (y / z) * -0.3333333333333333;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -1.85e+80: tmp = y * (-0.3333333333333333 / z) elif y <= 4.5e+83: tmp = x else: tmp = (y / z) * -0.3333333333333333 return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -1.85e+80) tmp = Float64(y * Float64(-0.3333333333333333 / z)); elseif (y <= 4.5e+83) tmp = x; else tmp = Float64(Float64(y / z) * -0.3333333333333333); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -1.85e+80) tmp = y * (-0.3333333333333333 / z); elseif (y <= 4.5e+83) tmp = x; else tmp = (y / z) * -0.3333333333333333; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -1.85e+80], N[(y * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.5e+83], x, N[(N[(y / z), $MachinePrecision] * -0.3333333333333333), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.85 \cdot 10^{+80}:\\
\;\;\;\;y \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 4.5 \cdot 10^{+83}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{z} \cdot -0.3333333333333333\\
\end{array}
\end{array}
if y < -1.84999999999999998e80Initial program 97.7%
+-commutative97.7%
associate-+r-97.7%
sub-neg97.7%
associate-*l*97.7%
*-commutative97.7%
distribute-frac-neg297.7%
distribute-rgt-neg-in97.7%
metadata-eval97.7%
Simplified97.7%
*-un-lft-identity97.7%
times-frac97.8%
Applied egg-rr97.8%
*-commutative97.8%
associate-*l/97.6%
associate-*r/97.6%
*-rgt-identity97.6%
Simplified97.6%
Taylor expanded in y around inf 67.3%
*-commutative67.3%
associate-*l/67.3%
associate-*r/67.4%
Simplified67.4%
if -1.84999999999999998e80 < y < 4.4999999999999999e83Initial program 96.3%
+-commutative96.3%
associate-+r-96.3%
sub-neg96.3%
associate-*l*96.8%
*-commutative96.8%
distribute-frac-neg296.8%
distribute-rgt-neg-in96.8%
metadata-eval96.8%
Simplified96.8%
Taylor expanded in z around inf 35.7%
if 4.4999999999999999e83 < y Initial program 99.8%
+-commutative99.8%
associate-+r-99.8%
sub-neg99.8%
associate-*l*99.8%
*-commutative99.8%
distribute-frac-neg299.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
*-un-lft-identity99.8%
times-frac99.6%
Applied egg-rr99.6%
*-commutative99.6%
associate-*l/99.7%
associate-*r/99.8%
*-rgt-identity99.8%
Simplified99.8%
Taylor expanded in y around inf 79.6%
Final simplification48.2%
(FPCore (x y z t) :precision binary64 (+ x (/ 1.0 (* 3.0 (/ z (- (/ t y) y))))))
double code(double x, double y, double z, double t) {
return x + (1.0 / (3.0 * (z / ((t / y) - y))));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = x + (1.0d0 / (3.0d0 * (z / ((t / y) - y))))
end function
public static double code(double x, double y, double z, double t) {
return x + (1.0 / (3.0 * (z / ((t / y) - y))));
}
def code(x, y, z, t): return x + (1.0 / (3.0 * (z / ((t / y) - y))))
function code(x, y, z, t) return Float64(x + Float64(1.0 / Float64(3.0 * Float64(z / Float64(Float64(t / y) - y))))) end
function tmp = code(x, y, z, t) tmp = x + (1.0 / (3.0 * (z / ((t / y) - y)))); end
code[x_, y_, z_, t_] := N[(x + N[(1.0 / N[(3.0 * N[(z / N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{1}{3 \cdot \frac{z}{\frac{t}{y} - y}}
\end{array}
Initial program 97.1%
sub-neg97.1%
associate-+l+97.1%
+-commutative97.1%
remove-double-neg97.1%
distribute-frac-neg97.1%
distribute-neg-in97.1%
remove-double-neg97.1%
sub-neg97.1%
neg-mul-197.1%
times-frac96.5%
distribute-frac-neg96.5%
neg-mul-196.5%
*-commutative96.5%
associate-/l*96.4%
*-commutative96.4%
Simplified96.7%
Taylor expanded in z around 0 96.8%
associate-*r/96.8%
distribute-lft-out--96.8%
clear-num96.8%
*-un-lft-identity96.8%
distribute-lft-out--96.8%
times-frac96.9%
metadata-eval96.9%
Applied egg-rr96.9%
(FPCore (x y z t) :precision binary64 (+ x (/ (- (/ t y) y) (* z 3.0))))
double code(double x, double y, double z, double t) {
return x + (((t / y) - y) / (z * 3.0));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = x + (((t / y) - y) / (z * 3.0d0))
end function
public static double code(double x, double y, double z, double t) {
return x + (((t / y) - y) / (z * 3.0));
}
def code(x, y, z, t): return x + (((t / y) - y) / (z * 3.0))
function code(x, y, z, t) return Float64(x + Float64(Float64(Float64(t / y) - y) / Float64(z * 3.0))) end
function tmp = code(x, y, z, t) tmp = x + (((t / y) - y) / (z * 3.0)); end
code[x_, y_, z_, t_] := N[(x + N[(N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision] / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{\frac{t}{y} - y}{z \cdot 3}
\end{array}
Initial program 97.1%
sub-neg97.1%
associate-+l+97.1%
+-commutative97.1%
remove-double-neg97.1%
distribute-frac-neg97.1%
distribute-neg-in97.1%
remove-double-neg97.1%
sub-neg97.1%
neg-mul-197.1%
times-frac96.5%
distribute-frac-neg96.5%
neg-mul-196.5%
*-commutative96.5%
associate-/l*96.4%
*-commutative96.4%
Simplified96.7%
*-commutative96.7%
clear-num96.7%
div-inv96.8%
metadata-eval96.8%
un-div-inv96.9%
Applied egg-rr96.9%
(FPCore (x y z t) :precision binary64 (+ x (* 0.3333333333333333 (/ (- (/ t y) y) z))))
double code(double x, double y, double z, double t) {
return x + (0.3333333333333333 * (((t / y) - y) / z));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = x + (0.3333333333333333d0 * (((t / y) - y) / z))
end function
public static double code(double x, double y, double z, double t) {
return x + (0.3333333333333333 * (((t / y) - y) / z));
}
def code(x, y, z, t): return x + (0.3333333333333333 * (((t / y) - y) / z))
function code(x, y, z, t) return Float64(x + Float64(0.3333333333333333 * Float64(Float64(Float64(t / y) - y) / z))) end
function tmp = code(x, y, z, t) tmp = x + (0.3333333333333333 * (((t / y) - y) / z)); end
code[x_, y_, z_, t_] := N[(x + N[(0.3333333333333333 * N[(N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + 0.3333333333333333 \cdot \frac{\frac{t}{y} - y}{z}
\end{array}
Initial program 97.1%
sub-neg97.1%
associate-+l+97.1%
+-commutative97.1%
remove-double-neg97.1%
distribute-frac-neg97.1%
distribute-neg-in97.1%
remove-double-neg97.1%
sub-neg97.1%
neg-mul-197.1%
times-frac96.5%
distribute-frac-neg96.5%
neg-mul-196.5%
*-commutative96.5%
associate-/l*96.4%
*-commutative96.4%
Simplified96.7%
Taylor expanded in z around 0 96.8%
(FPCore (x y z t) :precision binary64 x)
double code(double x, double y, double z, double t) {
return x;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = x
end function
public static double code(double x, double y, double z, double t) {
return x;
}
def code(x, y, z, t): return x
function code(x, y, z, t) return x end
function tmp = code(x, y, z, t) tmp = x; end
code[x_, y_, z_, t_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 97.1%
+-commutative97.1%
associate-+r-97.1%
sub-neg97.1%
associate-*l*97.4%
*-commutative97.4%
distribute-frac-neg297.4%
distribute-rgt-neg-in97.4%
metadata-eval97.4%
Simplified97.4%
Taylor expanded in z around inf 32.4%
(FPCore (x y z t) :precision binary64 (+ (- x (/ y (* z 3.0))) (/ (/ t (* z 3.0)) y)))
double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + ((t / (z * 3.0)) / y);
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x - (y / (z * 3.0d0))) + ((t / (z * 3.0d0)) / y)
end function
public static double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + ((t / (z * 3.0)) / y);
}
def code(x, y, z, t): return (x - (y / (z * 3.0))) + ((t / (z * 3.0)) / y)
function code(x, y, z, t) return Float64(Float64(x - Float64(y / Float64(z * 3.0))) + Float64(Float64(t / Float64(z * 3.0)) / y)) end
function tmp = code(x, y, z, t) tmp = (x - (y / (z * 3.0))) + ((t / (z * 3.0)) / y); end
code[x_, y_, z_, t_] := N[(N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(t / N[(z * 3.0), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x - \frac{y}{z \cdot 3}\right) + \frac{\frac{t}{z \cdot 3}}{y}
\end{array}
herbie shell --seed 2024085
(FPCore (x y z t)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, H"
:precision binary64
:alt
(+ (- x (/ y (* z 3.0))) (/ (/ t (* z 3.0)) y))
(+ (- x (/ y (* z 3.0))) (/ t (* (* z 3.0) y))))