
(FPCore (x y z t a) :precision binary64 (+ x (/ (* y (- z t)) a)))
double code(double x, double y, double z, double t, double a) {
return x + ((y * (z - t)) / a);
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = x + ((y * (z - t)) / a)
end function
public static double code(double x, double y, double z, double t, double a) {
return x + ((y * (z - t)) / a);
}
def code(x, y, z, t, a): return x + ((y * (z - t)) / a)
function code(x, y, z, t, a) return Float64(x + Float64(Float64(y * Float64(z - t)) / a)) end
function tmp = code(x, y, z, t, a) tmp = x + ((y * (z - t)) / a); end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(y * N[(z - t), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{y \cdot \left(z - t\right)}{a}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 21 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a) :precision binary64 (+ x (/ (* y (- z t)) a)))
double code(double x, double y, double z, double t, double a) {
return x + ((y * (z - t)) / a);
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = x + ((y * (z - t)) / a)
end function
public static double code(double x, double y, double z, double t, double a) {
return x + ((y * (z - t)) / a);
}
def code(x, y, z, t, a): return x + ((y * (z - t)) / a)
function code(x, y, z, t, a) return Float64(x + Float64(Float64(y * Float64(z - t)) / a)) end
function tmp = code(x, y, z, t, a) tmp = x + ((y * (z - t)) / a); end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(y * N[(z - t), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{y \cdot \left(z - t\right)}{a}
\end{array}
(FPCore (x y z t a) :precision binary64 (+ x (* (- z t) (/ y a))))
double code(double x, double y, double z, double t, double a) {
return x + ((z - t) * (y / a));
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = x + ((z - t) * (y / a))
end function
public static double code(double x, double y, double z, double t, double a) {
return x + ((z - t) * (y / a));
}
def code(x, y, z, t, a): return x + ((z - t) * (y / a))
function code(x, y, z, t, a) return Float64(x + Float64(Float64(z - t) * Float64(y / a))) end
function tmp = code(x, y, z, t, a) tmp = x + ((z - t) * (y / a)); end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(z - t), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(z - t\right) \cdot \frac{y}{a}
\end{array}
Initial program 95.1%
associate-/l*93.9%
Simplified93.9%
Taylor expanded in y around 0 95.1%
associate-*l/97.1%
*-commutative97.1%
Simplified97.1%
Final simplification97.1%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (/ (* (- z t) y) a))) (if (<= t_1 1e+138) (+ x (* z (/ y a))) t_1)))
double code(double x, double y, double z, double t, double a) {
double t_1 = ((z - t) * y) / a;
double tmp;
if (t_1 <= 1e+138) {
tmp = x + (z * (y / a));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = ((z - t) * y) / a
if (t_1 <= 1d+138) then
tmp = x + (z * (y / a))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = ((z - t) * y) / a;
double tmp;
if (t_1 <= 1e+138) {
tmp = x + (z * (y / a));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = ((z - t) * y) / a tmp = 0 if t_1 <= 1e+138: tmp = x + (z * (y / a)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(Float64(z - t) * y) / a) tmp = 0.0 if (t_1 <= 1e+138) tmp = Float64(x + Float64(z * Float64(y / a))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = ((z - t) * y) / a; tmp = 0.0; if (t_1 <= 1e+138) tmp = x + (z * (y / a)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(N[(z - t), $MachinePrecision] * y), $MachinePrecision] / a), $MachinePrecision]}, If[LessEqual[t$95$1, 1e+138], N[(x + N[(z * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\left(z - t\right) \cdot y}{a}\\
\mathbf{if}\;t\_1 \leq 10^{+138}:\\
\;\;\;\;x + z \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (/.f64 (*.f64 y (-.f64 z t)) a) < 1e138Initial program 96.1%
associate-/l*94.6%
Simplified94.6%
Taylor expanded in t around 0 78.1%
+-commutative78.1%
associate-/l*76.9%
Simplified76.9%
clear-num28.7%
un-div-inv28.7%
Applied egg-rr77.4%
associate-/r/80.0%
Applied egg-rr80.0%
if 1e138 < (/.f64 (*.f64 y (-.f64 z t)) a) Initial program 91.5%
associate-/l*91.5%
Simplified91.5%
Taylor expanded in x around inf 73.2%
associate-/l*78.3%
Simplified78.3%
Taylor expanded in x around 0 89.5%
Final simplification82.1%
(FPCore (x y z t a) :precision binary64 (if (<= x -1.6e+109) x (* y (/ (- z t) a))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (x <= -1.6e+109) {
tmp = x;
} else {
tmp = y * ((z - t) / a);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (x <= (-1.6d+109)) then
tmp = x
else
tmp = y * ((z - t) / a)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (x <= -1.6e+109) {
tmp = x;
} else {
tmp = y * ((z - t) / a);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if x <= -1.6e+109: tmp = x else: tmp = y * ((z - t) / a) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (x <= -1.6e+109) tmp = x; else tmp = Float64(y * Float64(Float64(z - t) / a)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (x <= -1.6e+109) tmp = x; else tmp = y * ((z - t) / a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[x, -1.6e+109], x, N[(y * N[(N[(z - t), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.6 \cdot 10^{+109}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;y \cdot \frac{z - t}{a}\\
\end{array}
\end{array}
if x < -1.6000000000000001e109Initial program 97.2%
associate-/l*91.6%
Simplified91.6%
Taylor expanded in x around inf 69.8%
if -1.6000000000000001e109 < x Initial program 94.7%
associate-/l*94.3%
Simplified94.3%
Taylor expanded in x around inf 82.7%
associate-/l*81.5%
Simplified81.5%
Taylor expanded in x around 0 59.8%
associate-/l*60.6%
Simplified60.6%
Final simplification61.9%
(FPCore (x y z t a) :precision binary64 (if (<= a -6.5e+104) x (* (- z t) (/ y a))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -6.5e+104) {
tmp = x;
} else {
tmp = (z - t) * (y / a);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (a <= (-6.5d+104)) then
tmp = x
else
tmp = (z - t) * (y / a)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -6.5e+104) {
tmp = x;
} else {
tmp = (z - t) * (y / a);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -6.5e+104: tmp = x else: tmp = (z - t) * (y / a) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -6.5e+104) tmp = x; else tmp = Float64(Float64(z - t) * Float64(y / a)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -6.5e+104) tmp = x; else tmp = (z - t) * (y / a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -6.5e+104], x, N[(N[(z - t), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -6.5 \cdot 10^{+104}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;\left(z - t\right) \cdot \frac{y}{a}\\
\end{array}
\end{array}
if a < -6.5000000000000005e104Initial program 83.1%
associate-/l*99.9%
Simplified99.9%
Taylor expanded in x around inf 64.8%
if -6.5000000000000005e104 < a Initial program 97.6%
associate-/l*92.7%
Simplified92.7%
Taylor expanded in x around inf 85.4%
associate-/l*81.9%
Simplified81.9%
Taylor expanded in x around 0 62.2%
associate-/l*59.6%
Simplified59.6%
Taylor expanded in y around 0 62.2%
associate-*l/98.1%
*-commutative98.1%
Simplified62.7%
Final simplification63.0%
(FPCore (x y z t a) :precision binary64 (if (<= t -8.5e+56) (* (- z t) (/ y a)) (+ x (/ (* z y) a))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -8.5e+56) {
tmp = (z - t) * (y / a);
} else {
tmp = x + ((z * y) / a);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (t <= (-8.5d+56)) then
tmp = (z - t) * (y / a)
else
tmp = x + ((z * y) / a)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -8.5e+56) {
tmp = (z - t) * (y / a);
} else {
tmp = x + ((z * y) / a);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -8.5e+56: tmp = (z - t) * (y / a) else: tmp = x + ((z * y) / a) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -8.5e+56) tmp = Float64(Float64(z - t) * Float64(y / a)); else tmp = Float64(x + Float64(Float64(z * y) / a)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -8.5e+56) tmp = (z - t) * (y / a); else tmp = x + ((z * y) / a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -8.5e+56], N[(N[(z - t), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z * y), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -8.5 \cdot 10^{+56}:\\
\;\;\;\;\left(z - t\right) \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{z \cdot y}{a}\\
\end{array}
\end{array}
if t < -8.4999999999999998e56Initial program 94.6%
associate-/l*89.2%
Simplified89.2%
Taylor expanded in x around inf 78.8%
associate-/l*70.8%
Simplified70.8%
Taylor expanded in x around 0 69.7%
associate-/l*67.7%
Simplified67.7%
Taylor expanded in y around 0 69.7%
associate-*l/100.0%
*-commutative100.0%
Simplified75.1%
if -8.4999999999999998e56 < t Initial program 95.2%
associate-/l*94.7%
Simplified94.7%
Taylor expanded in z around inf 76.5%
Final simplification76.3%
(FPCore (x y z t a) :precision binary64 (if (<= t -4.2e+65) (* (- z t) (/ y a)) (+ x (* y (/ z a)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -4.2e+65) {
tmp = (z - t) * (y / a);
} else {
tmp = x + (y * (z / a));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (t <= (-4.2d+65)) then
tmp = (z - t) * (y / a)
else
tmp = x + (y * (z / a))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -4.2e+65) {
tmp = (z - t) * (y / a);
} else {
tmp = x + (y * (z / a));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -4.2e+65: tmp = (z - t) * (y / a) else: tmp = x + (y * (z / a)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -4.2e+65) tmp = Float64(Float64(z - t) * Float64(y / a)); else tmp = Float64(x + Float64(y * Float64(z / a))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -4.2e+65) tmp = (z - t) * (y / a); else tmp = x + (y * (z / a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -4.2e+65], N[(N[(z - t), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision], N[(x + N[(y * N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -4.2 \cdot 10^{+65}:\\
\;\;\;\;\left(z - t\right) \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;x + y \cdot \frac{z}{a}\\
\end{array}
\end{array}
if t < -4.19999999999999983e65Initial program 94.6%
associate-/l*89.2%
Simplified89.2%
Taylor expanded in x around inf 78.8%
associate-/l*70.8%
Simplified70.8%
Taylor expanded in x around 0 69.7%
associate-/l*67.7%
Simplified67.7%
Taylor expanded in y around 0 69.7%
associate-*l/100.0%
*-commutative100.0%
Simplified75.1%
if -4.19999999999999983e65 < t Initial program 95.2%
associate-/l*94.7%
Simplified94.7%
Taylor expanded in t around 0 76.5%
+-commutative76.5%
associate-/l*76.2%
Simplified76.2%
Final simplification76.1%
(FPCore (x y z t a) :precision binary64 (if (<= t -1.52e+66) (* (- z t) (/ y a)) (+ x (* z (/ y a)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -1.52e+66) {
tmp = (z - t) * (y / a);
} else {
tmp = x + (z * (y / a));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (t <= (-1.52d+66)) then
tmp = (z - t) * (y / a)
else
tmp = x + (z * (y / a))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -1.52e+66) {
tmp = (z - t) * (y / a);
} else {
tmp = x + (z * (y / a));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -1.52e+66: tmp = (z - t) * (y / a) else: tmp = x + (z * (y / a)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -1.52e+66) tmp = Float64(Float64(z - t) * Float64(y / a)); else tmp = Float64(x + Float64(z * Float64(y / a))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -1.52e+66) tmp = (z - t) * (y / a); else tmp = x + (z * (y / a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -1.52e+66], N[(N[(z - t), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision], N[(x + N[(z * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.52 \cdot 10^{+66}:\\
\;\;\;\;\left(z - t\right) \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;x + z \cdot \frac{y}{a}\\
\end{array}
\end{array}
if t < -1.52000000000000004e66Initial program 94.6%
associate-/l*89.2%
Simplified89.2%
Taylor expanded in x around inf 78.8%
associate-/l*70.8%
Simplified70.8%
Taylor expanded in x around 0 69.7%
associate-/l*67.7%
Simplified67.7%
Taylor expanded in y around 0 69.7%
associate-*l/100.0%
*-commutative100.0%
Simplified75.1%
if -1.52000000000000004e66 < t Initial program 95.2%
associate-/l*94.7%
Simplified94.7%
Taylor expanded in t around 0 76.5%
+-commutative76.5%
associate-/l*76.2%
Simplified76.2%
clear-num35.4%
un-div-inv35.5%
Applied egg-rr76.2%
associate-/r/77.3%
Applied egg-rr77.3%
Final simplification77.0%
(FPCore (x y z t a) :precision binary64 (if (<= t -3e+55) (* (- z t) (/ y a)) (+ x (/ y (/ a z)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -3e+55) {
tmp = (z - t) * (y / a);
} else {
tmp = x + (y / (a / z));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (t <= (-3d+55)) then
tmp = (z - t) * (y / a)
else
tmp = x + (y / (a / z))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -3e+55) {
tmp = (z - t) * (y / a);
} else {
tmp = x + (y / (a / z));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -3e+55: tmp = (z - t) * (y / a) else: tmp = x + (y / (a / z)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -3e+55) tmp = Float64(Float64(z - t) * Float64(y / a)); else tmp = Float64(x + Float64(y / Float64(a / z))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -3e+55) tmp = (z - t) * (y / a); else tmp = x + (y / (a / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -3e+55], N[(N[(z - t), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision], N[(x + N[(y / N[(a / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -3 \cdot 10^{+55}:\\
\;\;\;\;\left(z - t\right) \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{\frac{a}{z}}\\
\end{array}
\end{array}
if t < -3.00000000000000017e55Initial program 94.6%
associate-/l*89.2%
Simplified89.2%
Taylor expanded in x around inf 78.8%
associate-/l*70.8%
Simplified70.8%
Taylor expanded in x around 0 69.7%
associate-/l*67.7%
Simplified67.7%
Taylor expanded in y around 0 69.7%
associate-*l/100.0%
*-commutative100.0%
Simplified75.1%
if -3.00000000000000017e55 < t Initial program 95.2%
associate-/l*94.7%
Simplified94.7%
Taylor expanded in t around 0 76.5%
+-commutative76.5%
associate-/l*76.2%
Simplified76.2%
clear-num35.4%
un-div-inv35.5%
Applied egg-rr76.2%
Final simplification76.1%
(FPCore (x y z t a) :precision binary64 (if (<= t -2.25e+55) (* (- z t) (/ y a)) (+ x (/ (* z y) a))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -2.25e+55) {
tmp = (z - t) * (y / a);
} else {
tmp = x + ((z * y) / a);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (t <= (-2.25d+55)) then
tmp = (z - t) * (y / a)
else
tmp = x + ((z * y) / a)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -2.25e+55) {
tmp = (z - t) * (y / a);
} else {
tmp = x + ((z * y) / a);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -2.25e+55: tmp = (z - t) * (y / a) else: tmp = x + ((z * y) / a) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -2.25e+55) tmp = Float64(Float64(z - t) * Float64(y / a)); else tmp = Float64(x + Float64(Float64(z * y) / a)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -2.25e+55) tmp = (z - t) * (y / a); else tmp = x + ((z * y) / a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -2.25e+55], N[(N[(z - t), $MachinePrecision] * N[(y / a), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z * y), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -2.25 \cdot 10^{+55}:\\
\;\;\;\;\left(z - t\right) \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{z \cdot y}{a}\\
\end{array}
\end{array}
if t < -2.24999999999999999e55Initial program 94.6%
associate-/l*89.2%
Simplified89.2%
Taylor expanded in x around inf 78.8%
associate-/l*70.8%
Simplified70.8%
Taylor expanded in x around 0 69.7%
associate-/l*67.7%
Simplified67.7%
Taylor expanded in y around 0 69.7%
associate-*l/100.0%
*-commutative100.0%
Simplified75.1%
if -2.24999999999999999e55 < t Initial program 95.2%
associate-/l*94.7%
Simplified94.7%
Taylor expanded in t around 0 76.5%
+-commutative76.5%
associate-/l*76.2%
Simplified76.2%
*-commutative35.5%
associate-*l/34.5%
Applied egg-rr76.5%
Final simplification76.3%
(FPCore (x y z t a) :precision binary64 (if (<= t 2.05e+82) (+ x (* z (/ y a))) (- x (* t (/ y a)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 2.05e+82) {
tmp = x + (z * (y / a));
} else {
tmp = x - (t * (y / a));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (t <= 2.05d+82) then
tmp = x + (z * (y / a))
else
tmp = x - (t * (y / a))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 2.05e+82) {
tmp = x + (z * (y / a));
} else {
tmp = x - (t * (y / a));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= 2.05e+82: tmp = x + (z * (y / a)) else: tmp = x - (t * (y / a)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= 2.05e+82) tmp = Float64(x + Float64(z * Float64(y / a))); else tmp = Float64(x - Float64(t * Float64(y / a))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= 2.05e+82) tmp = x + (z * (y / a)); else tmp = x - (t * (y / a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, 2.05e+82], N[(x + N[(z * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 2.05 \cdot 10^{+82}:\\
\;\;\;\;x + z \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;x - t \cdot \frac{y}{a}\\
\end{array}
\end{array}
if t < 2.04999999999999998e82Initial program 96.1%
associate-/l*94.7%
Simplified94.7%
Taylor expanded in t around 0 79.0%
+-commutative79.0%
associate-/l*78.5%
Simplified78.5%
clear-num37.3%
un-div-inv36.9%
Applied egg-rr78.5%
associate-/r/80.7%
Applied egg-rr80.7%
if 2.04999999999999998e82 < t Initial program 91.4%
associate-/l*91.4%
Simplified91.4%
Taylor expanded in y around 0 91.4%
associate-*l/96.3%
*-commutative96.3%
Simplified96.3%
Taylor expanded in z around 0 84.8%
associate-*l/78.5%
*-commutative78.5%
neg-mul-178.5%
sub-neg78.5%
*-commutative78.5%
associate-*l/84.8%
associate-*r/83.6%
Simplified83.6%
Final simplification81.3%
(FPCore (x y z t a) :precision binary64 (if (<= t 1.45e+80) (+ x (* z (/ y a))) (- x (* y (/ t a)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 1.45e+80) {
tmp = x + (z * (y / a));
} else {
tmp = x - (y * (t / a));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (t <= 1.45d+80) then
tmp = x + (z * (y / a))
else
tmp = x - (y * (t / a))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 1.45e+80) {
tmp = x + (z * (y / a));
} else {
tmp = x - (y * (t / a));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= 1.45e+80: tmp = x + (z * (y / a)) else: tmp = x - (y * (t / a)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= 1.45e+80) tmp = Float64(x + Float64(z * Float64(y / a))); else tmp = Float64(x - Float64(y * Float64(t / a))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= 1.45e+80) tmp = x + (z * (y / a)); else tmp = x - (y * (t / a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, 1.45e+80], N[(x + N[(z * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(y * N[(t / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 1.45 \cdot 10^{+80}:\\
\;\;\;\;x + z \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;x - y \cdot \frac{t}{a}\\
\end{array}
\end{array}
if t < 1.44999999999999993e80Initial program 96.1%
associate-/l*94.7%
Simplified94.7%
Taylor expanded in t around 0 79.0%
+-commutative79.0%
associate-/l*78.5%
Simplified78.5%
clear-num37.3%
un-div-inv36.9%
Applied egg-rr78.5%
associate-/r/80.7%
Applied egg-rr80.7%
if 1.44999999999999993e80 < t Initial program 91.4%
associate-/l*91.4%
Simplified91.4%
Taylor expanded in z around 0 84.8%
mul-1-neg84.8%
unsub-neg84.8%
*-commutative84.8%
associate-/l*78.5%
Simplified78.5%
Final simplification80.2%
(FPCore (x y z t a) :precision binary64 (if (<= t 1.65e+79) (+ x (* z (/ y a))) (- x (/ (* t y) a))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 1.65e+79) {
tmp = x + (z * (y / a));
} else {
tmp = x - ((t * y) / a);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (t <= 1.65d+79) then
tmp = x + (z * (y / a))
else
tmp = x - ((t * y) / a)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= 1.65e+79) {
tmp = x + (z * (y / a));
} else {
tmp = x - ((t * y) / a);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= 1.65e+79: tmp = x + (z * (y / a)) else: tmp = x - ((t * y) / a) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= 1.65e+79) tmp = Float64(x + Float64(z * Float64(y / a))); else tmp = Float64(x - Float64(Float64(t * y) / a)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= 1.65e+79) tmp = x + (z * (y / a)); else tmp = x - ((t * y) / a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, 1.65e+79], N[(x + N[(z * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(N[(t * y), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 1.65 \cdot 10^{+79}:\\
\;\;\;\;x + z \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;x - \frac{t \cdot y}{a}\\
\end{array}
\end{array}
if t < 1.6500000000000001e79Initial program 96.1%
associate-/l*94.7%
Simplified94.7%
Taylor expanded in t around 0 79.0%
+-commutative79.0%
associate-/l*78.5%
Simplified78.5%
clear-num37.3%
un-div-inv36.9%
Applied egg-rr78.5%
associate-/r/80.7%
Applied egg-rr80.7%
if 1.6500000000000001e79 < t Initial program 91.4%
associate-/l*91.4%
Simplified91.4%
Taylor expanded in z around 0 84.8%
+-commutative84.8%
associate-*r/84.8%
mul-1-neg84.8%
distribute-lft-neg-out84.8%
*-commutative84.8%
Simplified84.8%
Final simplification81.6%
(FPCore (x y z t a) :precision binary64 (if (<= y 7.8e+15) x (* y (/ (- t) a))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= 7.8e+15) {
tmp = x;
} else {
tmp = y * (-t / a);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (y <= 7.8d+15) then
tmp = x
else
tmp = y * (-t / a)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= 7.8e+15) {
tmp = x;
} else {
tmp = y * (-t / a);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if y <= 7.8e+15: tmp = x else: tmp = y * (-t / a) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (y <= 7.8e+15) tmp = x; else tmp = Float64(y * Float64(Float64(-t) / a)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (y <= 7.8e+15) tmp = x; else tmp = y * (-t / a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[y, 7.8e+15], x, N[(y * N[((-t) / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 7.8 \cdot 10^{+15}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;y \cdot \frac{-t}{a}\\
\end{array}
\end{array}
if y < 7.8e15Initial program 96.4%
associate-/l*92.1%
Simplified92.1%
Taylor expanded in x around inf 50.2%
if 7.8e15 < y Initial program 90.9%
associate-/l*99.8%
Simplified99.8%
Taylor expanded in x around inf 72.5%
associate-/l*77.3%
Simplified77.3%
Taylor expanded in t around inf 58.2%
associate-*r/58.2%
neg-mul-158.2%
distribute-rgt-neg-in58.2%
associate-*l/58.2%
Simplified58.2%
Final simplification52.2%
(FPCore (x y z t a) :precision binary64 (if (<= y 510000000.0) x (/ y (/ a (- t)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= 510000000.0) {
tmp = x;
} else {
tmp = y / (a / -t);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (y <= 510000000.0d0) then
tmp = x
else
tmp = y / (a / -t)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= 510000000.0) {
tmp = x;
} else {
tmp = y / (a / -t);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if y <= 510000000.0: tmp = x else: tmp = y / (a / -t) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (y <= 510000000.0) tmp = x; else tmp = Float64(y / Float64(a / Float64(-t))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (y <= 510000000.0) tmp = x; else tmp = y / (a / -t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[y, 510000000.0], x, N[(y / N[(a / (-t)), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 510000000:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{\frac{a}{-t}}\\
\end{array}
\end{array}
if y < 5.1e8Initial program 96.4%
associate-/l*92.1%
Simplified92.1%
Taylor expanded in x around inf 50.2%
if 5.1e8 < y Initial program 90.9%
associate-/l*99.8%
Simplified99.8%
Taylor expanded in x around inf 72.5%
associate-/l*77.3%
Simplified77.3%
Taylor expanded in t around inf 58.2%
mul-1-neg58.2%
*-commutative58.2%
associate-*l/59.8%
associate-/r/58.3%
distribute-neg-frac258.3%
Simplified58.3%
Final simplification52.2%
(FPCore (x y z t a) :precision binary64 (if (<= y 2900000000.0) x (/ (* y (- t)) a)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= 2900000000.0) {
tmp = x;
} else {
tmp = (y * -t) / a;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (y <= 2900000000.0d0) then
tmp = x
else
tmp = (y * -t) / a
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= 2900000000.0) {
tmp = x;
} else {
tmp = (y * -t) / a;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if y <= 2900000000.0: tmp = x else: tmp = (y * -t) / a return tmp
function code(x, y, z, t, a) tmp = 0.0 if (y <= 2900000000.0) tmp = x; else tmp = Float64(Float64(y * Float64(-t)) / a); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (y <= 2900000000.0) tmp = x; else tmp = (y * -t) / a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[y, 2900000000.0], x, N[(N[(y * (-t)), $MachinePrecision] / a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 2900000000:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;\frac{y \cdot \left(-t\right)}{a}\\
\end{array}
\end{array}
if y < 2.9e9Initial program 96.4%
associate-/l*92.1%
Simplified92.1%
Taylor expanded in x around inf 50.2%
if 2.9e9 < y Initial program 90.9%
associate-/l*99.8%
Simplified99.8%
Taylor expanded in x around inf 72.5%
associate-/l*77.3%
Simplified77.3%
Taylor expanded in x around 0 78.7%
associate-/l*86.2%
Simplified86.2%
Taylor expanded in z around 0 58.2%
associate-*r/58.2%
mul-1-neg58.2%
distribute-lft-neg-out58.2%
*-commutative58.2%
Simplified58.2%
Final simplification52.2%
(FPCore (x y z t a) :precision binary64 (if (<= y 3.5e+111) x (* y (/ z a))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= 3.5e+111) {
tmp = x;
} else {
tmp = y * (z / a);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (y <= 3.5d+111) then
tmp = x
else
tmp = y * (z / a)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= 3.5e+111) {
tmp = x;
} else {
tmp = y * (z / a);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if y <= 3.5e+111: tmp = x else: tmp = y * (z / a) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (y <= 3.5e+111) tmp = x; else tmp = Float64(y * Float64(z / a)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (y <= 3.5e+111) tmp = x; else tmp = y * (z / a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[y, 3.5e+111], x, N[(y * N[(z / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 3.5 \cdot 10^{+111}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;y \cdot \frac{z}{a}\\
\end{array}
\end{array}
if y < 3.5000000000000002e111Initial program 96.3%
associate-/l*92.9%
Simplified92.9%
Taylor expanded in x around inf 48.1%
if 3.5000000000000002e111 < y Initial program 88.2%
associate-/l*99.8%
Simplified99.8%
Taylor expanded in x around inf 74.0%
associate-/l*78.8%
Simplified78.8%
Taylor expanded in z around inf 41.8%
associate-*r/48.8%
Simplified48.8%
Final simplification48.2%
(FPCore (x y z t a) :precision binary64 (if (<= y 2.2e+111) x (/ y (/ a z))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= 2.2e+111) {
tmp = x;
} else {
tmp = y / (a / z);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (y <= 2.2d+111) then
tmp = x
else
tmp = y / (a / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= 2.2e+111) {
tmp = x;
} else {
tmp = y / (a / z);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if y <= 2.2e+111: tmp = x else: tmp = y / (a / z) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (y <= 2.2e+111) tmp = x; else tmp = Float64(y / Float64(a / z)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (y <= 2.2e+111) tmp = x; else tmp = y / (a / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[y, 2.2e+111], x, N[(y / N[(a / z), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 2.2 \cdot 10^{+111}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{\frac{a}{z}}\\
\end{array}
\end{array}
if y < 2.19999999999999999e111Initial program 96.3%
associate-/l*92.9%
Simplified92.9%
Taylor expanded in x around inf 48.1%
if 2.19999999999999999e111 < y Initial program 88.2%
associate-/l*99.8%
Simplified99.8%
Taylor expanded in x around inf 74.0%
associate-/l*78.8%
Simplified78.8%
Taylor expanded in z around inf 41.8%
associate-*r/48.8%
Simplified48.8%
clear-num48.7%
un-div-inv48.7%
Applied egg-rr48.7%
Final simplification48.2%
(FPCore (x y z t a) :precision binary64 (if (<= y 1.85e+111) x (/ (* z y) a)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= 1.85e+111) {
tmp = x;
} else {
tmp = (z * y) / a;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (y <= 1.85d+111) then
tmp = x
else
tmp = (z * y) / a
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= 1.85e+111) {
tmp = x;
} else {
tmp = (z * y) / a;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if y <= 1.85e+111: tmp = x else: tmp = (z * y) / a return tmp
function code(x, y, z, t, a) tmp = 0.0 if (y <= 1.85e+111) tmp = x; else tmp = Float64(Float64(z * y) / a); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (y <= 1.85e+111) tmp = x; else tmp = (z * y) / a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[y, 1.85e+111], x, N[(N[(z * y), $MachinePrecision] / a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 1.85 \cdot 10^{+111}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;\frac{z \cdot y}{a}\\
\end{array}
\end{array}
if y < 1.8500000000000001e111Initial program 96.3%
associate-/l*92.9%
Simplified92.9%
Taylor expanded in x around inf 48.1%
if 1.8500000000000001e111 < y Initial program 88.2%
associate-/l*99.8%
Simplified99.8%
Taylor expanded in x around inf 74.0%
associate-/l*78.8%
Simplified78.8%
Taylor expanded in z around inf 41.8%
Final simplification47.1%
(FPCore (x y z t a) :precision binary64 (if (<= y 2.4e+111) x (/ (* z y) a)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= 2.4e+111) {
tmp = x;
} else {
tmp = (z * y) / a;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (y <= 2.4d+111) then
tmp = x
else
tmp = (z * y) / a
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= 2.4e+111) {
tmp = x;
} else {
tmp = (z * y) / a;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if y <= 2.4e+111: tmp = x else: tmp = (z * y) / a return tmp
function code(x, y, z, t, a) tmp = 0.0 if (y <= 2.4e+111) tmp = x; else tmp = Float64(Float64(z * y) / a); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (y <= 2.4e+111) tmp = x; else tmp = (z * y) / a; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[y, 2.4e+111], x, N[(N[(z * y), $MachinePrecision] / a), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 2.4 \cdot 10^{+111}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;\frac{z \cdot y}{a}\\
\end{array}
\end{array}
if y < 2.40000000000000006e111Initial program 96.3%
associate-/l*92.9%
Simplified92.9%
Taylor expanded in x around inf 48.1%
if 2.40000000000000006e111 < y Initial program 88.2%
associate-/l*99.8%
Simplified99.8%
Taylor expanded in x around inf 74.0%
associate-/l*78.8%
Simplified78.8%
Taylor expanded in z around inf 41.8%
associate-*r/48.8%
Simplified48.8%
*-commutative48.8%
associate-*l/41.8%
Applied egg-rr41.8%
Final simplification47.1%
(FPCore (x y z t a) :precision binary64 (+ x (* y (/ (- z t) a))))
double code(double x, double y, double z, double t, double a) {
return x + (y * ((z - t) / a));
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = x + (y * ((z - t) / a))
end function
public static double code(double x, double y, double z, double t, double a) {
return x + (y * ((z - t) / a));
}
def code(x, y, z, t, a): return x + (y * ((z - t) / a))
function code(x, y, z, t, a) return Float64(x + Float64(y * Float64(Float64(z - t) / a))) end
function tmp = code(x, y, z, t, a) tmp = x + (y * ((z - t) / a)); end
code[x_, y_, z_, t_, a_] := N[(x + N[(y * N[(N[(z - t), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + y \cdot \frac{z - t}{a}
\end{array}
Initial program 95.1%
associate-/l*93.9%
Simplified93.9%
Final simplification93.9%
(FPCore (x y z t a) :precision binary64 x)
double code(double x, double y, double z, double t, double a) {
return x;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = x
end function
public static double code(double x, double y, double z, double t, double a) {
return x;
}
def code(x, y, z, t, a): return x
function code(x, y, z, t, a) return x end
function tmp = code(x, y, z, t, a) tmp = x; end
code[x_, y_, z_, t_, a_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 95.1%
associate-/l*93.9%
Simplified93.9%
Taylor expanded in x around inf 41.6%
Final simplification41.6%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ a (- z t))))
(if (< y -1.0761266216389975e-10)
(+ x (/ 1.0 (/ t_1 y)))
(if (< y 2.894426862792089e-49)
(+ x (/ (* y (- z t)) a))
(+ x (/ y t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = a / (z - t);
double tmp;
if (y < -1.0761266216389975e-10) {
tmp = x + (1.0 / (t_1 / y));
} else if (y < 2.894426862792089e-49) {
tmp = x + ((y * (z - t)) / a);
} else {
tmp = x + (y / t_1);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = a / (z - t)
if (y < (-1.0761266216389975d-10)) then
tmp = x + (1.0d0 / (t_1 / y))
else if (y < 2.894426862792089d-49) then
tmp = x + ((y * (z - t)) / a)
else
tmp = x + (y / t_1)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = a / (z - t);
double tmp;
if (y < -1.0761266216389975e-10) {
tmp = x + (1.0 / (t_1 / y));
} else if (y < 2.894426862792089e-49) {
tmp = x + ((y * (z - t)) / a);
} else {
tmp = x + (y / t_1);
}
return tmp;
}
def code(x, y, z, t, a): t_1 = a / (z - t) tmp = 0 if y < -1.0761266216389975e-10: tmp = x + (1.0 / (t_1 / y)) elif y < 2.894426862792089e-49: tmp = x + ((y * (z - t)) / a) else: tmp = x + (y / t_1) return tmp
function code(x, y, z, t, a) t_1 = Float64(a / Float64(z - t)) tmp = 0.0 if (y < -1.0761266216389975e-10) tmp = Float64(x + Float64(1.0 / Float64(t_1 / y))); elseif (y < 2.894426862792089e-49) tmp = Float64(x + Float64(Float64(y * Float64(z - t)) / a)); else tmp = Float64(x + Float64(y / t_1)); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = a / (z - t); tmp = 0.0; if (y < -1.0761266216389975e-10) tmp = x + (1.0 / (t_1 / y)); elseif (y < 2.894426862792089e-49) tmp = x + ((y * (z - t)) / a); else tmp = x + (y / t_1); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(a / N[(z - t), $MachinePrecision]), $MachinePrecision]}, If[Less[y, -1.0761266216389975e-10], N[(x + N[(1.0 / N[(t$95$1 / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Less[y, 2.894426862792089e-49], N[(x + N[(N[(y * N[(z - t), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], N[(x + N[(y / t$95$1), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{a}{z - t}\\
\mathbf{if}\;y < -1.0761266216389975 \cdot 10^{-10}:\\
\;\;\;\;x + \frac{1}{\frac{t\_1}{y}}\\
\mathbf{elif}\;y < 2.894426862792089 \cdot 10^{-49}:\\
\;\;\;\;x + \frac{y \cdot \left(z - t\right)}{a}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{t\_1}\\
\end{array}
\end{array}
herbie shell --seed 2024066
(FPCore (x y z t a)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, E"
:precision binary64
:alt
(if (< y -1.0761266216389975e-10) (+ x (/ 1.0 (/ (/ a (- z t)) y))) (if (< y 2.894426862792089e-49) (+ x (/ (* y (- z t)) a)) (+ x (/ y (/ a (- z t))))))
(+ x (/ (* y (- z t)) a)))