
(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 12 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 (* (/ y a) (- z t))))
double code(double x, double y, double z, double t, double a) {
return x + ((y / a) * (z - t));
}
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 / a) * (z - t))
end function
public static double code(double x, double y, double z, double t, double a) {
return x + ((y / a) * (z - t));
}
def code(x, y, z, t, a): return x + ((y / a) * (z - t))
function code(x, y, z, t, a) return Float64(x + Float64(Float64(y / a) * Float64(z - t))) end
function tmp = code(x, y, z, t, a) tmp = x + ((y / a) * (z - t)); end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(y / a), $MachinePrecision] * N[(z - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{y}{a} \cdot \left(z - t\right)
\end{array}
Initial program 92.9%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6492.9%
Simplified92.9%
*-commutativeN/A
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6497.2%
Applied egg-rr97.2%
Final simplification97.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (* y (- z t)) a)))
(if (<= t_1 (- INFINITY))
(/ y (/ a (- z t)))
(if (<= t_1 -5e+161)
t_1
(if (<= t_1 1e-188)
(+ x (/ (* y z) a))
(if (<= t_1 1e+103)
(- x (* y (/ t a)))
(if (<= t_1 1e+301) t_1 (* y (/ (- z t) a)))))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (y * (z - t)) / a;
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = y / (a / (z - t));
} else if (t_1 <= -5e+161) {
tmp = t_1;
} else if (t_1 <= 1e-188) {
tmp = x + ((y * z) / a);
} else if (t_1 <= 1e+103) {
tmp = x - (y * (t / a));
} else if (t_1 <= 1e+301) {
tmp = t_1;
} else {
tmp = y * ((z - t) / a);
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (y * (z - t)) / a;
double tmp;
if (t_1 <= -Double.POSITIVE_INFINITY) {
tmp = y / (a / (z - t));
} else if (t_1 <= -5e+161) {
tmp = t_1;
} else if (t_1 <= 1e-188) {
tmp = x + ((y * z) / a);
} else if (t_1 <= 1e+103) {
tmp = x - (y * (t / a));
} else if (t_1 <= 1e+301) {
tmp = t_1;
} else {
tmp = y * ((z - t) / a);
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (y * (z - t)) / a tmp = 0 if t_1 <= -math.inf: tmp = y / (a / (z - t)) elif t_1 <= -5e+161: tmp = t_1 elif t_1 <= 1e-188: tmp = x + ((y * z) / a) elif t_1 <= 1e+103: tmp = x - (y * (t / a)) elif t_1 <= 1e+301: tmp = t_1 else: tmp = y * ((z - t) / a) return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(y * Float64(z - t)) / a) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(y / Float64(a / Float64(z - t))); elseif (t_1 <= -5e+161) tmp = t_1; elseif (t_1 <= 1e-188) tmp = Float64(x + Float64(Float64(y * z) / a)); elseif (t_1 <= 1e+103) tmp = Float64(x - Float64(y * Float64(t / a))); elseif (t_1 <= 1e+301) tmp = t_1; else tmp = Float64(y * Float64(Float64(z - t) / a)); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (y * (z - t)) / a; tmp = 0.0; if (t_1 <= -Inf) tmp = y / (a / (z - t)); elseif (t_1 <= -5e+161) tmp = t_1; elseif (t_1 <= 1e-188) tmp = x + ((y * z) / a); elseif (t_1 <= 1e+103) tmp = x - (y * (t / a)); elseif (t_1 <= 1e+301) tmp = t_1; else tmp = y * ((z - t) / a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(y * N[(z - t), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], N[(y / N[(a / N[(z - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -5e+161], t$95$1, If[LessEqual[t$95$1, 1e-188], N[(x + N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 1e+103], N[(x - N[(y * N[(t / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 1e+301], t$95$1, N[(y * N[(N[(z - t), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{y \cdot \left(z - t\right)}{a}\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;\frac{y}{\frac{a}{z - t}}\\
\mathbf{elif}\;t\_1 \leq -5 \cdot 10^{+161}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_1 \leq 10^{-188}:\\
\;\;\;\;x + \frac{y \cdot z}{a}\\
\mathbf{elif}\;t\_1 \leq 10^{+103}:\\
\;\;\;\;x - y \cdot \frac{t}{a}\\
\mathbf{elif}\;t\_1 \leq 10^{+301}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;y \cdot \frac{z - t}{a}\\
\end{array}
\end{array}
if (/.f64 (*.f64 y (-.f64 z t)) a) < -inf.0Initial program 82.4%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6482.4%
Simplified82.4%
Taylor expanded in x around 0
mul-1-negN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
mul-1-negN/A
sub-negN/A
--lowering--.f6482.4%
Simplified82.4%
*-commutativeN/A
associate-/l*N/A
sub-negN/A
+-commutativeN/A
remove-double-negN/A
distribute-neg-inN/A
sub-negN/A
distribute-lft-neg-inN/A
*-commutativeN/A
associate-/r/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
sub-negN/A
distribute-neg-inN/A
remove-double-negN/A
+-commutativeN/A
sub-negN/A
--lowering--.f6492.7%
Applied egg-rr92.7%
if -inf.0 < (/.f64 (*.f64 y (-.f64 z t)) a) < -4.9999999999999997e161 or 1e103 < (/.f64 (*.f64 y (-.f64 z t)) a) < 1.00000000000000005e301Initial program 99.6%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6499.6%
Simplified99.6%
Taylor expanded in x around 0
mul-1-negN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
mul-1-negN/A
sub-negN/A
--lowering--.f6492.8%
Simplified92.8%
if -4.9999999999999997e161 < (/.f64 (*.f64 y (-.f64 z t)) a) < 9.9999999999999995e-189Initial program 99.9%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6499.9%
Simplified99.9%
Taylor expanded in t around 0
sub-negN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6490.6%
Simplified90.6%
if 9.9999999999999995e-189 < (/.f64 (*.f64 y (-.f64 z t)) a) < 1e103Initial program 99.8%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6499.8%
Simplified99.8%
*-commutativeN/A
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6499.8%
Applied egg-rr99.8%
Taylor expanded in t around inf
Simplified88.4%
*-commutativeN/A
frac-2negN/A
frac-2negN/A
*-lft-identityN/A
associate-*l/N/A
associate-*r*N/A
div-invN/A
*-lowering-*.f64N/A
/-lowering-/.f6488.4%
Applied egg-rr88.4%
if 1.00000000000000005e301 < (/.f64 (*.f64 y (-.f64 z t)) a) Initial program 77.1%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6477.1%
Simplified77.1%
Taylor expanded in x around 0
mul-1-negN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
mul-1-negN/A
sub-negN/A
--lowering--.f6477.1%
Simplified77.1%
associate-/l*N/A
*-commutativeN/A
clear-numN/A
sub-negN/A
+-commutativeN/A
remove-double-negN/A
distribute-neg-inN/A
sub-negN/A
distribute-neg-frac2N/A
clear-numN/A
distribute-neg-frac2N/A
remove-double-divN/A
*-lowering-*.f64N/A
Applied egg-rr97.5%
Final simplification92.0%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (* y (- z t)) a)) (t_2 (/ (- z t) (/ a y))))
(if (<= t_1 -1e+161)
t_2
(if (<= t_1 1e-188)
(+ x (/ (* y z) a))
(if (<= t_1 1e+103) (- x (* y (/ t a))) t_2)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (y * (z - t)) / a;
double t_2 = (z - t) / (a / y);
double tmp;
if (t_1 <= -1e+161) {
tmp = t_2;
} else if (t_1 <= 1e-188) {
tmp = x + ((y * z) / a);
} else if (t_1 <= 1e+103) {
tmp = x - (y * (t / a));
} else {
tmp = t_2;
}
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) :: t_2
real(8) :: tmp
t_1 = (y * (z - t)) / a
t_2 = (z - t) / (a / y)
if (t_1 <= (-1d+161)) then
tmp = t_2
else if (t_1 <= 1d-188) then
tmp = x + ((y * z) / a)
else if (t_1 <= 1d+103) then
tmp = x - (y * (t / a))
else
tmp = t_2
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (y * (z - t)) / a;
double t_2 = (z - t) / (a / y);
double tmp;
if (t_1 <= -1e+161) {
tmp = t_2;
} else if (t_1 <= 1e-188) {
tmp = x + ((y * z) / a);
} else if (t_1 <= 1e+103) {
tmp = x - (y * (t / a));
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (y * (z - t)) / a t_2 = (z - t) / (a / y) tmp = 0 if t_1 <= -1e+161: tmp = t_2 elif t_1 <= 1e-188: tmp = x + ((y * z) / a) elif t_1 <= 1e+103: tmp = x - (y * (t / a)) else: tmp = t_2 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(y * Float64(z - t)) / a) t_2 = Float64(Float64(z - t) / Float64(a / y)) tmp = 0.0 if (t_1 <= -1e+161) tmp = t_2; elseif (t_1 <= 1e-188) tmp = Float64(x + Float64(Float64(y * z) / a)); elseif (t_1 <= 1e+103) tmp = Float64(x - Float64(y * Float64(t / a))); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (y * (z - t)) / a; t_2 = (z - t) / (a / y); tmp = 0.0; if (t_1 <= -1e+161) tmp = t_2; elseif (t_1 <= 1e-188) tmp = x + ((y * z) / a); elseif (t_1 <= 1e+103) tmp = x - (y * (t / a)); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(y * N[(z - t), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]}, Block[{t$95$2 = N[(N[(z - t), $MachinePrecision] / N[(a / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -1e+161], t$95$2, If[LessEqual[t$95$1, 1e-188], N[(x + N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 1e+103], N[(x - N[(y * N[(t / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$2]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{y \cdot \left(z - t\right)}{a}\\
t_2 := \frac{z - t}{\frac{a}{y}}\\
\mathbf{if}\;t\_1 \leq -1 \cdot 10^{+161}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 10^{-188}:\\
\;\;\;\;x + \frac{y \cdot z}{a}\\
\mathbf{elif}\;t\_1 \leq 10^{+103}:\\
\;\;\;\;x - y \cdot \frac{t}{a}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if (/.f64 (*.f64 y (-.f64 z t)) a) < -1e161 or 1e103 < (/.f64 (*.f64 y (-.f64 z t)) a) Initial program 86.0%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6486.0%
Simplified86.0%
Taylor expanded in x around 0
mul-1-negN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
mul-1-negN/A
sub-negN/A
--lowering--.f6484.0%
Simplified84.0%
*-commutativeN/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
sub-negN/A
+-commutativeN/A
remove-double-negN/A
distribute-neg-inN/A
sub-negN/A
/-lowering-/.f64N/A
sub-negN/A
distribute-neg-inN/A
remove-double-negN/A
+-commutativeN/A
sub-negN/A
--lowering--.f64N/A
/-lowering-/.f6492.1%
Applied egg-rr92.1%
if -1e161 < (/.f64 (*.f64 y (-.f64 z t)) a) < 9.9999999999999995e-189Initial program 99.9%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6499.9%
Simplified99.9%
Taylor expanded in t around 0
sub-negN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6490.5%
Simplified90.5%
if 9.9999999999999995e-189 < (/.f64 (*.f64 y (-.f64 z t)) a) < 1e103Initial program 99.8%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6499.8%
Simplified99.8%
*-commutativeN/A
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6499.8%
Applied egg-rr99.8%
Taylor expanded in t around inf
Simplified88.4%
*-commutativeN/A
frac-2negN/A
frac-2negN/A
*-lft-identityN/A
associate-*l/N/A
associate-*r*N/A
div-invN/A
*-lowering-*.f64N/A
/-lowering-/.f6488.4%
Applied egg-rr88.4%
Final simplification91.0%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (* y (- z t)))) (if (<= t_1 -5e+304) (/ y (/ a (- z t))) (+ x (/ t_1 a)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y * (z - t);
double tmp;
if (t_1 <= -5e+304) {
tmp = y / (a / (z - t));
} else {
tmp = x + (t_1 / 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) :: t_1
real(8) :: tmp
t_1 = y * (z - t)
if (t_1 <= (-5d+304)) then
tmp = y / (a / (z - t))
else
tmp = x + (t_1 / a)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = y * (z - t);
double tmp;
if (t_1 <= -5e+304) {
tmp = y / (a / (z - t));
} else {
tmp = x + (t_1 / a);
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y * (z - t) tmp = 0 if t_1 <= -5e+304: tmp = y / (a / (z - t)) else: tmp = x + (t_1 / a) return tmp
function code(x, y, z, t, a) t_1 = Float64(y * Float64(z - t)) tmp = 0.0 if (t_1 <= -5e+304) tmp = Float64(y / Float64(a / Float64(z - t))); else tmp = Float64(x + Float64(t_1 / a)); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y * (z - t); tmp = 0.0; if (t_1 <= -5e+304) tmp = y / (a / (z - t)); else tmp = x + (t_1 / a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y * N[(z - t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -5e+304], N[(y / N[(a / N[(z - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(t$95$1 / a), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \left(z - t\right)\\
\mathbf{if}\;t\_1 \leq -5 \cdot 10^{+304}:\\
\;\;\;\;\frac{y}{\frac{a}{z - t}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{t\_1}{a}\\
\end{array}
\end{array}
if (*.f64 y (-.f64 z t)) < -4.9999999999999997e304Initial program 64.0%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6464.0%
Simplified64.0%
Taylor expanded in x around 0
mul-1-negN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
mul-1-negN/A
sub-negN/A
--lowering--.f6464.0%
Simplified64.0%
*-commutativeN/A
associate-/l*N/A
sub-negN/A
+-commutativeN/A
remove-double-negN/A
distribute-neg-inN/A
sub-negN/A
distribute-lft-neg-inN/A
*-commutativeN/A
associate-/r/N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
sub-negN/A
distribute-neg-inN/A
remove-double-negN/A
+-commutativeN/A
sub-negN/A
--lowering--.f6487.8%
Applied egg-rr87.8%
if -4.9999999999999997e304 < (*.f64 y (-.f64 z t)) Initial program 98.1%
(FPCore (x y z t a) :precision binary64 (if (<= t -1.6e+34) (- x (/ y (/ a t))) (if (<= t 16600000000.0) (+ x (/ (* y z) a)) (- x (* (/ y a) t)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -1.6e+34) {
tmp = x - (y / (a / t));
} else if (t <= 16600000000.0) {
tmp = x + ((y * z) / a);
} else {
tmp = x - ((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 (t <= (-1.6d+34)) then
tmp = x - (y / (a / t))
else if (t <= 16600000000.0d0) then
tmp = x + ((y * z) / a)
else
tmp = x - ((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 (t <= -1.6e+34) {
tmp = x - (y / (a / t));
} else if (t <= 16600000000.0) {
tmp = x + ((y * z) / a);
} else {
tmp = x - ((y / a) * t);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -1.6e+34: tmp = x - (y / (a / t)) elif t <= 16600000000.0: tmp = x + ((y * z) / a) else: tmp = x - ((y / a) * t) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -1.6e+34) tmp = Float64(x - Float64(y / Float64(a / t))); elseif (t <= 16600000000.0) tmp = Float64(x + Float64(Float64(y * z) / a)); else tmp = Float64(x - Float64(Float64(y / a) * t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -1.6e+34) tmp = x - (y / (a / t)); elseif (t <= 16600000000.0) tmp = x + ((y * z) / a); else tmp = x - ((y / a) * t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -1.6e+34], N[(x - N[(y / N[(a / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 16600000000.0], N[(x + N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], N[(x - N[(N[(y / a), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.6 \cdot 10^{+34}:\\
\;\;\;\;x - \frac{y}{\frac{a}{t}}\\
\mathbf{elif}\;t \leq 16600000000:\\
\;\;\;\;x + \frac{y \cdot z}{a}\\
\mathbf{else}:\\
\;\;\;\;x - \frac{y}{a} \cdot t\\
\end{array}
\end{array}
if t < -1.5999999999999999e34Initial program 90.3%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6490.3%
Simplified90.3%
--lowering--.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
--lowering--.f6497.9%
Applied egg-rr97.9%
Taylor expanded in t around inf
Simplified85.6%
if -1.5999999999999999e34 < t < 1.66e10Initial program 94.4%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6494.4%
Simplified94.4%
Taylor expanded in t around 0
sub-negN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6486.9%
Simplified86.9%
if 1.66e10 < t Initial program 92.0%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6492.0%
Simplified92.0%
*-commutativeN/A
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6498.9%
Applied egg-rr98.9%
Taylor expanded in t around inf
Simplified84.4%
(FPCore (x y z t a) :precision binary64 (if (<= t -1.3e+33) (- x (* y (/ t a))) (if (<= t 102000000000.0) (+ x (/ (* y z) a)) (- x (* (/ y a) t)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -1.3e+33) {
tmp = x - (y * (t / a));
} else if (t <= 102000000000.0) {
tmp = x + ((y * z) / a);
} else {
tmp = x - ((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 (t <= (-1.3d+33)) then
tmp = x - (y * (t / a))
else if (t <= 102000000000.0d0) then
tmp = x + ((y * z) / a)
else
tmp = x - ((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 (t <= -1.3e+33) {
tmp = x - (y * (t / a));
} else if (t <= 102000000000.0) {
tmp = x + ((y * z) / a);
} else {
tmp = x - ((y / a) * t);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -1.3e+33: tmp = x - (y * (t / a)) elif t <= 102000000000.0: tmp = x + ((y * z) / a) else: tmp = x - ((y / a) * t) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -1.3e+33) tmp = Float64(x - Float64(y * Float64(t / a))); elseif (t <= 102000000000.0) tmp = Float64(x + Float64(Float64(y * z) / a)); else tmp = Float64(x - Float64(Float64(y / a) * t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -1.3e+33) tmp = x - (y * (t / a)); elseif (t <= 102000000000.0) tmp = x + ((y * z) / a); else tmp = x - ((y / a) * t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -1.3e+33], N[(x - N[(y * N[(t / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 102000000000.0], N[(x + N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], N[(x - N[(N[(y / a), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.3 \cdot 10^{+33}:\\
\;\;\;\;x - y \cdot \frac{t}{a}\\
\mathbf{elif}\;t \leq 102000000000:\\
\;\;\;\;x + \frac{y \cdot z}{a}\\
\mathbf{else}:\\
\;\;\;\;x - \frac{y}{a} \cdot t\\
\end{array}
\end{array}
if t < -1.2999999999999999e33Initial program 90.3%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6490.3%
Simplified90.3%
*-commutativeN/A
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6499.0%
Applied egg-rr99.0%
Taylor expanded in t around inf
Simplified84.8%
*-commutativeN/A
frac-2negN/A
frac-2negN/A
*-lft-identityN/A
associate-*l/N/A
associate-*r*N/A
div-invN/A
*-lowering-*.f64N/A
/-lowering-/.f6485.6%
Applied egg-rr85.6%
if -1.2999999999999999e33 < t < 1.02e11Initial program 94.4%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6494.4%
Simplified94.4%
Taylor expanded in t around 0
sub-negN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6486.9%
Simplified86.9%
if 1.02e11 < t Initial program 92.0%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6492.0%
Simplified92.0%
*-commutativeN/A
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6498.9%
Applied egg-rr98.9%
Taylor expanded in t around inf
Simplified84.4%
Final simplification85.9%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (- x (* y (/ t a))))) (if (<= t -7.5e+34) t_1 (if (<= t 210000000.0) (+ x (/ (* y z) a)) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y * (t / a));
double tmp;
if (t <= -7.5e+34) {
tmp = t_1;
} else if (t <= 210000000.0) {
tmp = x + ((y * z) / 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 = x - (y * (t / a))
if (t <= (-7.5d+34)) then
tmp = t_1
else if (t <= 210000000.0d0) then
tmp = x + ((y * z) / 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 = x - (y * (t / a));
double tmp;
if (t <= -7.5e+34) {
tmp = t_1;
} else if (t <= 210000000.0) {
tmp = x + ((y * z) / a);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - (y * (t / a)) tmp = 0 if t <= -7.5e+34: tmp = t_1 elif t <= 210000000.0: tmp = x + ((y * z) / a) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(y * Float64(t / a))) tmp = 0.0 if (t <= -7.5e+34) tmp = t_1; elseif (t <= 210000000.0) tmp = Float64(x + Float64(Float64(y * z) / a)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - (y * (t / a)); tmp = 0.0; if (t <= -7.5e+34) tmp = t_1; elseif (t <= 210000000.0) tmp = x + ((y * z) / a); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(y * N[(t / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -7.5e+34], t$95$1, If[LessEqual[t, 210000000.0], N[(x + N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - y \cdot \frac{t}{a}\\
\mathbf{if}\;t \leq -7.5 \cdot 10^{+34}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 210000000:\\
\;\;\;\;x + \frac{y \cdot z}{a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -7.49999999999999976e34 or 2.1e8 < t Initial program 91.3%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6491.3%
Simplified91.3%
*-commutativeN/A
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6498.9%
Applied egg-rr98.9%
Taylor expanded in t around inf
Simplified84.5%
*-commutativeN/A
frac-2negN/A
frac-2negN/A
*-lft-identityN/A
associate-*l/N/A
associate-*r*N/A
div-invN/A
*-lowering-*.f64N/A
/-lowering-/.f6484.6%
Applied egg-rr84.6%
if -7.49999999999999976e34 < t < 2.1e8Initial program 94.4%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6494.4%
Simplified94.4%
Taylor expanded in t around 0
sub-negN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6486.9%
Simplified86.9%
Final simplification85.8%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (* y (/ (- z t) a)))) (if (<= y -5.6e+20) t_1 (if (<= y 9e-14) (+ x (/ (* y z) a)) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y * ((z - t) / a);
double tmp;
if (y <= -5.6e+20) {
tmp = t_1;
} else if (y <= 9e-14) {
tmp = x + ((y * z) / 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 = y * ((z - t) / a)
if (y <= (-5.6d+20)) then
tmp = t_1
else if (y <= 9d-14) then
tmp = x + ((y * z) / 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 = y * ((z - t) / a);
double tmp;
if (y <= -5.6e+20) {
tmp = t_1;
} else if (y <= 9e-14) {
tmp = x + ((y * z) / a);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y * ((z - t) / a) tmp = 0 if y <= -5.6e+20: tmp = t_1 elif y <= 9e-14: tmp = x + ((y * z) / a) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(y * Float64(Float64(z - t) / a)) tmp = 0.0 if (y <= -5.6e+20) tmp = t_1; elseif (y <= 9e-14) tmp = Float64(x + Float64(Float64(y * z) / a)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = y * ((z - t) / a); tmp = 0.0; if (y <= -5.6e+20) tmp = t_1; elseif (y <= 9e-14) tmp = x + ((y * z) / a); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y * N[(N[(z - t), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -5.6e+20], t$95$1, If[LessEqual[y, 9e-14], N[(x + N[(N[(y * z), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \frac{z - t}{a}\\
\mathbf{if}\;y \leq -5.6 \cdot 10^{+20}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 9 \cdot 10^{-14}:\\
\;\;\;\;x + \frac{y \cdot z}{a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -5.6e20 or 8.9999999999999995e-14 < y Initial program 86.9%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6486.9%
Simplified86.9%
Taylor expanded in x around 0
mul-1-negN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
mul-1-negN/A
sub-negN/A
--lowering--.f6470.1%
Simplified70.1%
associate-/l*N/A
*-commutativeN/A
clear-numN/A
sub-negN/A
+-commutativeN/A
remove-double-negN/A
distribute-neg-inN/A
sub-negN/A
distribute-neg-frac2N/A
clear-numN/A
distribute-neg-frac2N/A
remove-double-divN/A
*-lowering-*.f64N/A
Applied egg-rr79.4%
if -5.6e20 < y < 8.9999999999999995e-14Initial program 99.9%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6499.9%
Simplified99.9%
Taylor expanded in t around 0
sub-negN/A
mul-1-negN/A
remove-double-negN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6483.8%
Simplified83.8%
Final simplification81.4%
(FPCore (x y z t a) :precision binary64 (if (<= x -7.4e+58) x (if (<= x 3.6e+93) (* y (/ (- z t) a)) x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (x <= -7.4e+58) {
tmp = x;
} else if (x <= 3.6e+93) {
tmp = y * ((z - t) / a);
} else {
tmp = x;
}
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 <= (-7.4d+58)) then
tmp = x
else if (x <= 3.6d+93) then
tmp = y * ((z - t) / a)
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (x <= -7.4e+58) {
tmp = x;
} else if (x <= 3.6e+93) {
tmp = y * ((z - t) / a);
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if x <= -7.4e+58: tmp = x elif x <= 3.6e+93: tmp = y * ((z - t) / a) else: tmp = x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (x <= -7.4e+58) tmp = x; elseif (x <= 3.6e+93) tmp = Float64(y * Float64(Float64(z - t) / a)); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (x <= -7.4e+58) tmp = x; elseif (x <= 3.6e+93) tmp = y * ((z - t) / a); else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[x, -7.4e+58], x, If[LessEqual[x, 3.6e+93], N[(y * N[(N[(z - t), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -7.4 \cdot 10^{+58}:\\
\;\;\;\;x\\
\mathbf{elif}\;x \leq 3.6 \cdot 10^{+93}:\\
\;\;\;\;y \cdot \frac{z - t}{a}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if x < -7.4000000000000004e58 or 3.5999999999999999e93 < x Initial program 93.1%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6493.1%
Simplified93.1%
Taylor expanded in x around inf
Simplified63.8%
if -7.4000000000000004e58 < x < 3.5999999999999999e93Initial program 92.8%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6492.8%
Simplified92.8%
Taylor expanded in x around 0
mul-1-negN/A
distribute-neg-fracN/A
/-lowering-/.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
mul-1-negN/A
sub-negN/A
--lowering--.f6473.1%
Simplified73.1%
associate-/l*N/A
*-commutativeN/A
clear-numN/A
sub-negN/A
+-commutativeN/A
remove-double-negN/A
distribute-neg-inN/A
sub-negN/A
distribute-neg-frac2N/A
clear-numN/A
distribute-neg-frac2N/A
remove-double-divN/A
*-lowering-*.f64N/A
Applied egg-rr74.2%
Final simplification70.4%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (/ z (/ a y)))) (if (<= z -2.3e+30) t_1 (if (<= z 5.9e+63) x t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = z / (a / y);
double tmp;
if (z <= -2.3e+30) {
tmp = t_1;
} else if (z <= 5.9e+63) {
tmp = x;
} 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 / (a / y)
if (z <= (-2.3d+30)) then
tmp = t_1
else if (z <= 5.9d+63) then
tmp = x
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 / (a / y);
double tmp;
if (z <= -2.3e+30) {
tmp = t_1;
} else if (z <= 5.9e+63) {
tmp = x;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = z / (a / y) tmp = 0 if z <= -2.3e+30: tmp = t_1 elif z <= 5.9e+63: tmp = x else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(z / Float64(a / y)) tmp = 0.0 if (z <= -2.3e+30) tmp = t_1; elseif (z <= 5.9e+63) tmp = x; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = z / (a / y); tmp = 0.0; if (z <= -2.3e+30) tmp = t_1; elseif (z <= 5.9e+63) tmp = x; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(z / N[(a / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -2.3e+30], t$95$1, If[LessEqual[z, 5.9e+63], x, t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z}{\frac{a}{y}}\\
\mathbf{if}\;z \leq -2.3 \cdot 10^{+30}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 5.9 \cdot 10^{+63}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -2.3e30 or 5.90000000000000029e63 < z Initial program 89.0%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6489.0%
Simplified89.0%
Taylor expanded in z around inf
/-lowering-/.f64N/A
*-lowering-*.f6458.7%
Simplified58.7%
*-commutativeN/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f6463.3%
Applied egg-rr63.3%
if -2.3e30 < z < 5.90000000000000029e63Initial program 96.4%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6496.4%
Simplified96.4%
Taylor expanded in x around inf
Simplified49.3%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (* (/ y a) z))) (if (<= z -9e+36) t_1 (if (<= z 1.72e+62) x t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (y / a) * z;
double tmp;
if (z <= -9e+36) {
tmp = t_1;
} else if (z <= 1.72e+62) {
tmp = x;
} 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 = (y / a) * z
if (z <= (-9d+36)) then
tmp = t_1
else if (z <= 1.72d+62) then
tmp = x
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 = (y / a) * z;
double tmp;
if (z <= -9e+36) {
tmp = t_1;
} else if (z <= 1.72e+62) {
tmp = x;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (y / a) * z tmp = 0 if z <= -9e+36: tmp = t_1 elif z <= 1.72e+62: tmp = x else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(y / a) * z) tmp = 0.0 if (z <= -9e+36) tmp = t_1; elseif (z <= 1.72e+62) tmp = x; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (y / a) * z; tmp = 0.0; if (z <= -9e+36) tmp = t_1; elseif (z <= 1.72e+62) tmp = x; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(y / a), $MachinePrecision] * z), $MachinePrecision]}, If[LessEqual[z, -9e+36], t$95$1, If[LessEqual[z, 1.72e+62], x, t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{y}{a} \cdot z\\
\mathbf{if}\;z \leq -9 \cdot 10^{+36}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 1.72 \cdot 10^{+62}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -8.99999999999999994e36 or 1.7200000000000001e62 < z Initial program 89.0%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6489.0%
Simplified89.0%
Taylor expanded in z around inf
/-lowering-/.f64N/A
*-lowering-*.f6458.7%
Simplified58.7%
*-commutativeN/A
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6463.3%
Applied egg-rr63.3%
if -8.99999999999999994e36 < z < 1.7200000000000001e62Initial program 96.4%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6496.4%
Simplified96.4%
Taylor expanded in x around inf
Simplified49.3%
(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 92.9%
*-commutativeN/A
associate-/l*N/A
cancel-sign-subN/A
sub0-negN/A
associate-+l-N/A
neg-sub0N/A
+-commutativeN/A
sub-negN/A
*-commutativeN/A
--lowering--.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6492.9%
Simplified92.9%
Taylor expanded in x around inf
Simplified36.9%
(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 2024138
(FPCore (x y z t a)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, E"
:precision binary64
:alt
(! :herbie-platform default (if (< y -430450648655599/4000000000000000000000000) (+ x (/ 1 (/ (/ a (- z t)) y))) (if (< y 2894426862792089/10000000000000000000000000000000000000000000000000000000000000000) (+ x (/ (* y (- z t)) a)) (+ x (/ y (/ a (- z t)))))))
(+ x (/ (* y (- z t)) a)))