
(FPCore (x y z t a) :precision binary64 (- (+ x y) (/ (* (- z t) y) (- a t))))
double code(double x, double y, double z, double t, double a) {
return (x + y) - (((z - t) * y) / (a - 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) - (((z - t) * y) / (a - t))
end function
public static double code(double x, double y, double z, double t, double a) {
return (x + y) - (((z - t) * y) / (a - t));
}
def code(x, y, z, t, a): return (x + y) - (((z - t) * y) / (a - t))
function code(x, y, z, t, a) return Float64(Float64(x + y) - Float64(Float64(Float64(z - t) * y) / Float64(a - t))) end
function tmp = code(x, y, z, t, a) tmp = (x + y) - (((z - t) * y) / (a - t)); end
code[x_, y_, z_, t_, a_] := N[(N[(x + y), $MachinePrecision] - N[(N[(N[(z - t), $MachinePrecision] * y), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 15 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a) :precision binary64 (- (+ x y) (/ (* (- z t) y) (- a t))))
double code(double x, double y, double z, double t, double a) {
return (x + y) - (((z - t) * y) / (a - 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) - (((z - t) * y) / (a - t))
end function
public static double code(double x, double y, double z, double t, double a) {
return (x + y) - (((z - t) * y) / (a - t));
}
def code(x, y, z, t, a): return (x + y) - (((z - t) * y) / (a - t))
function code(x, y, z, t, a) return Float64(Float64(x + y) - Float64(Float64(Float64(z - t) * y) / Float64(a - t))) end
function tmp = code(x, y, z, t, a) tmp = (x + y) - (((z - t) * y) / (a - t)); end
code[x_, y_, z_, t_, a_] := N[(N[(x + y), $MachinePrecision] - N[(N[(N[(z - t), $MachinePrecision] * y), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t}
\end{array}
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ x (* y (+ 1.0 (/ (- z t) (- t a))))))
(t_2 (+ (+ x y) (/ (* y (- z t)) (- t a)))))
(if (<= t_2 -1e-247) t_1 (if (<= t_2 0.0) (+ x (/ (* y (- z a)) t)) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x + (y * (1.0 + ((z - t) / (t - a))));
double t_2 = (x + y) + ((y * (z - t)) / (t - a));
double tmp;
if (t_2 <= -1e-247) {
tmp = t_1;
} else if (t_2 <= 0.0) {
tmp = x + ((y * (z - a)) / t);
} 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) :: t_2
real(8) :: tmp
t_1 = x + (y * (1.0d0 + ((z - t) / (t - a))))
t_2 = (x + y) + ((y * (z - t)) / (t - a))
if (t_2 <= (-1d-247)) then
tmp = t_1
else if (t_2 <= 0.0d0) then
tmp = x + ((y * (z - a)) / t)
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 * (1.0 + ((z - t) / (t - a))));
double t_2 = (x + y) + ((y * (z - t)) / (t - a));
double tmp;
if (t_2 <= -1e-247) {
tmp = t_1;
} else if (t_2 <= 0.0) {
tmp = x + ((y * (z - a)) / t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x + (y * (1.0 + ((z - t) / (t - a)))) t_2 = (x + y) + ((y * (z - t)) / (t - a)) tmp = 0 if t_2 <= -1e-247: tmp = t_1 elif t_2 <= 0.0: tmp = x + ((y * (z - a)) / t) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x + Float64(y * Float64(1.0 + Float64(Float64(z - t) / Float64(t - a))))) t_2 = Float64(Float64(x + y) + Float64(Float64(y * Float64(z - t)) / Float64(t - a))) tmp = 0.0 if (t_2 <= -1e-247) tmp = t_1; elseif (t_2 <= 0.0) tmp = Float64(x + Float64(Float64(y * Float64(z - a)) / t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x + (y * (1.0 + ((z - t) / (t - a)))); t_2 = (x + y) + ((y * (z - t)) / (t - a)); tmp = 0.0; if (t_2 <= -1e-247) tmp = t_1; elseif (t_2 <= 0.0) tmp = x + ((y * (z - a)) / t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x + N[(y * N[(1.0 + N[(N[(z - t), $MachinePrecision] / N[(t - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(x + y), $MachinePrecision] + N[(N[(y * N[(z - t), $MachinePrecision]), $MachinePrecision] / N[(t - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -1e-247], t$95$1, If[LessEqual[t$95$2, 0.0], N[(x + N[(N[(y * N[(z - a), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + y \cdot \left(1 + \frac{z - t}{t - a}\right)\\
t_2 := \left(x + y\right) + \frac{y \cdot \left(z - t\right)}{t - a}\\
\mathbf{if}\;t\_2 \leq -1 \cdot 10^{-247}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 \leq 0:\\
\;\;\;\;x + \frac{y \cdot \left(z - a\right)}{t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (-.f64 (+.f64 x y) (/.f64 (*.f64 (-.f64 z t) y) (-.f64 a t))) < -1e-247 or 0.0 < (-.f64 (+.f64 x y) (/.f64 (*.f64 (-.f64 z t) y) (-.f64 a t))) Initial program 84.7%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval97.0%
Simplified97.0%
if -1e-247 < (-.f64 (+.f64 x y) (/.f64 (*.f64 (-.f64 z t) y) (-.f64 a t))) < 0.0Initial program 4.7%
Taylor expanded in t around inf
Simplified99.7%
Final simplification97.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (* z (/ y a)))))
(if (<= a -1.15e+147)
(+ x y)
(if (<= a -2.62e-11)
t_1
(if (<= a 9.8e-92)
(+ x (* y (/ (- z a) t)))
(if (<= a 1.65e+79) t_1 (+ x y)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - (z * (y / a));
double tmp;
if (a <= -1.15e+147) {
tmp = x + y;
} else if (a <= -2.62e-11) {
tmp = t_1;
} else if (a <= 9.8e-92) {
tmp = x + (y * ((z - a) / t));
} else if (a <= 1.65e+79) {
tmp = t_1;
} else {
tmp = x + y;
}
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 - (z * (y / a))
if (a <= (-1.15d+147)) then
tmp = x + y
else if (a <= (-2.62d-11)) then
tmp = t_1
else if (a <= 9.8d-92) then
tmp = x + (y * ((z - a) / t))
else if (a <= 1.65d+79) then
tmp = t_1
else
tmp = x + y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x - (z * (y / a));
double tmp;
if (a <= -1.15e+147) {
tmp = x + y;
} else if (a <= -2.62e-11) {
tmp = t_1;
} else if (a <= 9.8e-92) {
tmp = x + (y * ((z - a) / t));
} else if (a <= 1.65e+79) {
tmp = t_1;
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - (z * (y / a)) tmp = 0 if a <= -1.15e+147: tmp = x + y elif a <= -2.62e-11: tmp = t_1 elif a <= 9.8e-92: tmp = x + (y * ((z - a) / t)) elif a <= 1.65e+79: tmp = t_1 else: tmp = x + y return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(z * Float64(y / a))) tmp = 0.0 if (a <= -1.15e+147) tmp = Float64(x + y); elseif (a <= -2.62e-11) tmp = t_1; elseif (a <= 9.8e-92) tmp = Float64(x + Float64(y * Float64(Float64(z - a) / t))); elseif (a <= 1.65e+79) tmp = t_1; else tmp = Float64(x + y); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - (z * (y / a)); tmp = 0.0; if (a <= -1.15e+147) tmp = x + y; elseif (a <= -2.62e-11) tmp = t_1; elseif (a <= 9.8e-92) tmp = x + (y * ((z - a) / t)); elseif (a <= 1.65e+79) tmp = t_1; else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(z * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -1.15e+147], N[(x + y), $MachinePrecision], If[LessEqual[a, -2.62e-11], t$95$1, If[LessEqual[a, 9.8e-92], N[(x + N[(y * N[(N[(z - a), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 1.65e+79], t$95$1, N[(x + y), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - z \cdot \frac{y}{a}\\
\mathbf{if}\;a \leq -1.15 \cdot 10^{+147}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq -2.62 \cdot 10^{-11}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 9.8 \cdot 10^{-92}:\\
\;\;\;\;x + y \cdot \frac{z - a}{t}\\
\mathbf{elif}\;a \leq 1.65 \cdot 10^{+79}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -1.15e147 or 1.6500000000000001e79 < a Initial program 77.0%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval95.8%
Simplified95.8%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6484.4%
Simplified84.4%
if -1.15e147 < a < -2.62000000000000006e-11 or 9.8e-92 < a < 1.6500000000000001e79Initial program 87.4%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6489.1%
Simplified89.1%
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6487.9%
Applied egg-rr87.9%
Taylor expanded in x around inf
Simplified78.8%
if -2.62000000000000006e-11 < a < 9.8e-92Initial program 75.2%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval89.9%
Simplified89.9%
Taylor expanded in t around inf
+-lowering-+.f64N/A
associate-/l*N/A
mul-1-negN/A
sub-negN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6487.7%
Simplified87.7%
Final simplification84.3%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (* z (/ y a)))))
(if (<= a -1.15e+147)
(+ x y)
(if (<= a -1.42e-10)
t_1
(if (<= a 2.6e-85)
(+ x (* y (/ z t)))
(if (<= a 2.35e+79) t_1 (+ x y)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - (z * (y / a));
double tmp;
if (a <= -1.15e+147) {
tmp = x + y;
} else if (a <= -1.42e-10) {
tmp = t_1;
} else if (a <= 2.6e-85) {
tmp = x + (y * (z / t));
} else if (a <= 2.35e+79) {
tmp = t_1;
} else {
tmp = x + y;
}
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 - (z * (y / a))
if (a <= (-1.15d+147)) then
tmp = x + y
else if (a <= (-1.42d-10)) then
tmp = t_1
else if (a <= 2.6d-85) then
tmp = x + (y * (z / t))
else if (a <= 2.35d+79) then
tmp = t_1
else
tmp = x + y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = x - (z * (y / a));
double tmp;
if (a <= -1.15e+147) {
tmp = x + y;
} else if (a <= -1.42e-10) {
tmp = t_1;
} else if (a <= 2.6e-85) {
tmp = x + (y * (z / t));
} else if (a <= 2.35e+79) {
tmp = t_1;
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - (z * (y / a)) tmp = 0 if a <= -1.15e+147: tmp = x + y elif a <= -1.42e-10: tmp = t_1 elif a <= 2.6e-85: tmp = x + (y * (z / t)) elif a <= 2.35e+79: tmp = t_1 else: tmp = x + y return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(z * Float64(y / a))) tmp = 0.0 if (a <= -1.15e+147) tmp = Float64(x + y); elseif (a <= -1.42e-10) tmp = t_1; elseif (a <= 2.6e-85) tmp = Float64(x + Float64(y * Float64(z / t))); elseif (a <= 2.35e+79) tmp = t_1; else tmp = Float64(x + y); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - (z * (y / a)); tmp = 0.0; if (a <= -1.15e+147) tmp = x + y; elseif (a <= -1.42e-10) tmp = t_1; elseif (a <= 2.6e-85) tmp = x + (y * (z / t)); elseif (a <= 2.35e+79) tmp = t_1; else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(z * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -1.15e+147], N[(x + y), $MachinePrecision], If[LessEqual[a, -1.42e-10], t$95$1, If[LessEqual[a, 2.6e-85], N[(x + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 2.35e+79], t$95$1, N[(x + y), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - z \cdot \frac{y}{a}\\
\mathbf{if}\;a \leq -1.15 \cdot 10^{+147}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq -1.42 \cdot 10^{-10}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 2.6 \cdot 10^{-85}:\\
\;\;\;\;x + y \cdot \frac{z}{t}\\
\mathbf{elif}\;a \leq 2.35 \cdot 10^{+79}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -1.15e147 or 2.35000000000000011e79 < a Initial program 77.0%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval95.8%
Simplified95.8%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6484.4%
Simplified84.4%
if -1.15e147 < a < -1.42000000000000001e-10 or 2.60000000000000011e-85 < a < 2.35000000000000011e79Initial program 87.3%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6488.9%
Simplified88.9%
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6487.7%
Applied egg-rr87.7%
Taylor expanded in x around inf
Simplified79.8%
if -1.42000000000000001e-10 < a < 2.60000000000000011e-85Initial program 75.5%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval90.0%
Simplified90.0%
Taylor expanded in a around 0
/-lowering-/.f6484.8%
Simplified84.8%
Final simplification83.4%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- x (* y (/ z a)))))
(if (<= a -1.5e+147)
(+ x y)
(if (<= a -3.9e-7)
t_1
(if (<= a 2.6e-85)
(+ x (* y (/ z t)))
(if (<= a 2.6e+79) t_1 (+ x y)))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x - (y * (z / a));
double tmp;
if (a <= -1.5e+147) {
tmp = x + y;
} else if (a <= -3.9e-7) {
tmp = t_1;
} else if (a <= 2.6e-85) {
tmp = x + (y * (z / t));
} else if (a <= 2.6e+79) {
tmp = t_1;
} else {
tmp = x + y;
}
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 * (z / a))
if (a <= (-1.5d+147)) then
tmp = x + y
else if (a <= (-3.9d-7)) then
tmp = t_1
else if (a <= 2.6d-85) then
tmp = x + (y * (z / t))
else if (a <= 2.6d+79) then
tmp = t_1
else
tmp = x + y
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 * (z / a));
double tmp;
if (a <= -1.5e+147) {
tmp = x + y;
} else if (a <= -3.9e-7) {
tmp = t_1;
} else if (a <= 2.6e-85) {
tmp = x + (y * (z / t));
} else if (a <= 2.6e+79) {
tmp = t_1;
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x - (y * (z / a)) tmp = 0 if a <= -1.5e+147: tmp = x + y elif a <= -3.9e-7: tmp = t_1 elif a <= 2.6e-85: tmp = x + (y * (z / t)) elif a <= 2.6e+79: tmp = t_1 else: tmp = x + y return tmp
function code(x, y, z, t, a) t_1 = Float64(x - Float64(y * Float64(z / a))) tmp = 0.0 if (a <= -1.5e+147) tmp = Float64(x + y); elseif (a <= -3.9e-7) tmp = t_1; elseif (a <= 2.6e-85) tmp = Float64(x + Float64(y * Float64(z / t))); elseif (a <= 2.6e+79) tmp = t_1; else tmp = Float64(x + y); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x - (y * (z / a)); tmp = 0.0; if (a <= -1.5e+147) tmp = x + y; elseif (a <= -3.9e-7) tmp = t_1; elseif (a <= 2.6e-85) tmp = x + (y * (z / t)); elseif (a <= 2.6e+79) tmp = t_1; else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x - N[(y * N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -1.5e+147], N[(x + y), $MachinePrecision], If[LessEqual[a, -3.9e-7], t$95$1, If[LessEqual[a, 2.6e-85], N[(x + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 2.6e+79], t$95$1, N[(x + y), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - y \cdot \frac{z}{a}\\
\mathbf{if}\;a \leq -1.5 \cdot 10^{+147}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq -3.9 \cdot 10^{-7}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 2.6 \cdot 10^{-85}:\\
\;\;\;\;x + y \cdot \frac{z}{t}\\
\mathbf{elif}\;a \leq 2.6 \cdot 10^{+79}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -1.49999999999999997e147 or 2.60000000000000015e79 < a Initial program 77.0%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval95.8%
Simplified95.8%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6484.4%
Simplified84.4%
if -1.49999999999999997e147 < a < -3.90000000000000025e-7 or 2.60000000000000011e-85 < a < 2.60000000000000015e79Initial program 86.9%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6488.5%
Simplified88.5%
Taylor expanded in x around inf
Simplified80.2%
if -3.90000000000000025e-7 < a < 2.60000000000000011e-85Initial program 75.9%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval90.2%
Simplified90.2%
Taylor expanded in a around 0
/-lowering-/.f6484.2%
Simplified84.2%
Final simplification83.3%
(FPCore (x y z t a) :precision binary64 (+ x (* y (* z (- (/ -1.0 (- a t)) (/ (+ (/ t (- t a)) -1.0) z))))))
double code(double x, double y, double z, double t, double a) {
return x + (y * (z * ((-1.0 / (a - t)) - (((t / (t - a)) + -1.0) / z))));
}
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 * (((-1.0d0) / (a - t)) - (((t / (t - a)) + (-1.0d0)) / z))))
end function
public static double code(double x, double y, double z, double t, double a) {
return x + (y * (z * ((-1.0 / (a - t)) - (((t / (t - a)) + -1.0) / z))));
}
def code(x, y, z, t, a): return x + (y * (z * ((-1.0 / (a - t)) - (((t / (t - a)) + -1.0) / z))))
function code(x, y, z, t, a) return Float64(x + Float64(y * Float64(z * Float64(Float64(-1.0 / Float64(a - t)) - Float64(Float64(Float64(t / Float64(t - a)) + -1.0) / z))))) end
function tmp = code(x, y, z, t, a) tmp = x + (y * (z * ((-1.0 / (a - t)) - (((t / (t - a)) + -1.0) / z)))); end
code[x_, y_, z_, t_, a_] := N[(x + N[(y * N[(z * N[(N[(-1.0 / N[(a - t), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(t / N[(t - a), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + y \cdot \left(z \cdot \left(\frac{-1}{a - t} - \frac{\frac{t}{t - a} + -1}{z}\right)\right)
\end{array}
Initial program 79.0%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval93.7%
Simplified93.7%
Taylor expanded in z around -inf
associate-*r*N/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified96.2%
Final simplification96.2%
(FPCore (x y z t a) :precision binary64 (if (<= a -6.4e-13) (- (+ x y) (* z (/ y a))) (if (<= a 2.8e-92) (+ x (* y (/ (- z a) t))) (- (+ x y) (/ y (/ a z))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -6.4e-13) {
tmp = (x + y) - (z * (y / a));
} else if (a <= 2.8e-92) {
tmp = x + (y * ((z - a) / t));
} else {
tmp = (x + y) - (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 (a <= (-6.4d-13)) then
tmp = (x + y) - (z * (y / a))
else if (a <= 2.8d-92) then
tmp = x + (y * ((z - a) / t))
else
tmp = (x + y) - (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 (a <= -6.4e-13) {
tmp = (x + y) - (z * (y / a));
} else if (a <= 2.8e-92) {
tmp = x + (y * ((z - a) / t));
} else {
tmp = (x + y) - (y / (a / z));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -6.4e-13: tmp = (x + y) - (z * (y / a)) elif a <= 2.8e-92: tmp = x + (y * ((z - a) / t)) else: tmp = (x + y) - (y / (a / z)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -6.4e-13) tmp = Float64(Float64(x + y) - Float64(z * Float64(y / a))); elseif (a <= 2.8e-92) tmp = Float64(x + Float64(y * Float64(Float64(z - a) / t))); else tmp = Float64(Float64(x + y) - Float64(y / Float64(a / z))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -6.4e-13) tmp = (x + y) - (z * (y / a)); elseif (a <= 2.8e-92) tmp = x + (y * ((z - a) / t)); else tmp = (x + y) - (y / (a / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -6.4e-13], N[(N[(x + y), $MachinePrecision] - N[(z * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 2.8e-92], N[(x + N[(y * N[(N[(z - a), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x + y), $MachinePrecision] - N[(y / N[(a / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -6.4 \cdot 10^{-13}:\\
\;\;\;\;\left(x + y\right) - z \cdot \frac{y}{a}\\
\mathbf{elif}\;a \leq 2.8 \cdot 10^{-92}:\\
\;\;\;\;x + y \cdot \frac{z - a}{t}\\
\mathbf{else}:\\
\;\;\;\;\left(x + y\right) - \frac{y}{\frac{a}{z}}\\
\end{array}
\end{array}
if a < -6.39999999999999999e-13Initial program 81.3%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6492.7%
Simplified92.7%
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6493.0%
Applied egg-rr93.0%
if -6.39999999999999999e-13 < a < 2.8e-92Initial program 75.2%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval89.9%
Simplified89.9%
Taylor expanded in t around inf
+-lowering-+.f64N/A
associate-/l*N/A
mul-1-negN/A
sub-negN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6487.7%
Simplified87.7%
if 2.8e-92 < a Initial program 82.1%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6488.9%
Simplified88.9%
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f6488.9%
Applied egg-rr88.9%
(FPCore (x y z t a) :precision binary64 (if (<= a -1.26e-10) (- (+ x y) (* z (/ y a))) (if (<= a 1.05e-85) (+ x (* y (/ (- z a) t))) (- (+ x y) (* y (/ z a))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.26e-10) {
tmp = (x + y) - (z * (y / a));
} else if (a <= 1.05e-85) {
tmp = x + (y * ((z - a) / t));
} else {
tmp = (x + y) - (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 (a <= (-1.26d-10)) then
tmp = (x + y) - (z * (y / a))
else if (a <= 1.05d-85) then
tmp = x + (y * ((z - a) / t))
else
tmp = (x + y) - (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 (a <= -1.26e-10) {
tmp = (x + y) - (z * (y / a));
} else if (a <= 1.05e-85) {
tmp = x + (y * ((z - a) / t));
} else {
tmp = (x + y) - (y * (z / a));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -1.26e-10: tmp = (x + y) - (z * (y / a)) elif a <= 1.05e-85: tmp = x + (y * ((z - a) / t)) else: tmp = (x + y) - (y * (z / a)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -1.26e-10) tmp = Float64(Float64(x + y) - Float64(z * Float64(y / a))); elseif (a <= 1.05e-85) tmp = Float64(x + Float64(y * Float64(Float64(z - a) / t))); else tmp = Float64(Float64(x + y) - Float64(y * Float64(z / a))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -1.26e-10) tmp = (x + y) - (z * (y / a)); elseif (a <= 1.05e-85) tmp = x + (y * ((z - a) / t)); else tmp = (x + y) - (y * (z / a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -1.26e-10], N[(N[(x + y), $MachinePrecision] - N[(z * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 1.05e-85], N[(x + N[(y * N[(N[(z - a), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x + y), $MachinePrecision] - N[(y * N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.26 \cdot 10^{-10}:\\
\;\;\;\;\left(x + y\right) - z \cdot \frac{y}{a}\\
\mathbf{elif}\;a \leq 1.05 \cdot 10^{-85}:\\
\;\;\;\;x + y \cdot \frac{z - a}{t}\\
\mathbf{else}:\\
\;\;\;\;\left(x + y\right) - y \cdot \frac{z}{a}\\
\end{array}
\end{array}
if a < -1.26000000000000004e-10Initial program 81.3%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6492.7%
Simplified92.7%
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6493.0%
Applied egg-rr93.0%
if -1.26000000000000004e-10 < a < 1.05e-85Initial program 75.2%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval89.9%
Simplified89.9%
Taylor expanded in t around inf
+-lowering-+.f64N/A
associate-/l*N/A
mul-1-negN/A
sub-negN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6487.7%
Simplified87.7%
if 1.05e-85 < a Initial program 82.1%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6488.9%
Simplified88.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- (+ x y) (* y (/ z a)))))
(if (<= a -1.65e-11)
t_1
(if (<= a 7.5e-86) (+ x (* y (/ (- z a) t))) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (x + y) - (y * (z / a));
double tmp;
if (a <= -1.65e-11) {
tmp = t_1;
} else if (a <= 7.5e-86) {
tmp = x + (y * ((z - a) / t));
} 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) - (y * (z / a))
if (a <= (-1.65d-11)) then
tmp = t_1
else if (a <= 7.5d-86) then
tmp = x + (y * ((z - a) / t))
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) - (y * (z / a));
double tmp;
if (a <= -1.65e-11) {
tmp = t_1;
} else if (a <= 7.5e-86) {
tmp = x + (y * ((z - a) / t));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (x + y) - (y * (z / a)) tmp = 0 if a <= -1.65e-11: tmp = t_1 elif a <= 7.5e-86: tmp = x + (y * ((z - a) / t)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(x + y) - Float64(y * Float64(z / a))) tmp = 0.0 if (a <= -1.65e-11) tmp = t_1; elseif (a <= 7.5e-86) tmp = Float64(x + Float64(y * Float64(Float64(z - a) / t))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (x + y) - (y * (z / a)); tmp = 0.0; if (a <= -1.65e-11) tmp = t_1; elseif (a <= 7.5e-86) tmp = x + (y * ((z - a) / t)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(x + y), $MachinePrecision] - N[(y * N[(z / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -1.65e-11], t$95$1, If[LessEqual[a, 7.5e-86], N[(x + N[(y * N[(N[(z - a), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(x + y\right) - y \cdot \frac{z}{a}\\
\mathbf{if}\;a \leq -1.65 \cdot 10^{-11}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 7.5 \cdot 10^{-86}:\\
\;\;\;\;x + y \cdot \frac{z - a}{t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -1.6500000000000001e-11 or 7.50000000000000055e-86 < a Initial program 81.7%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6490.6%
Simplified90.6%
if -1.6500000000000001e-11 < a < 7.50000000000000055e-86Initial program 75.2%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval89.9%
Simplified89.9%
Taylor expanded in t around inf
+-lowering-+.f64N/A
associate-/l*N/A
mul-1-negN/A
sub-negN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6487.7%
Simplified87.7%
(FPCore (x y z t a) :precision binary64 (if (<= a -1.2e-212) (+ x y) (if (<= a 8.5e-259) (/ (* y z) t) (if (<= a 1.2e-95) x (+ x y)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.2e-212) {
tmp = x + y;
} else if (a <= 8.5e-259) {
tmp = (y * z) / t;
} else if (a <= 1.2e-95) {
tmp = x;
} else {
tmp = x + y;
}
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 <= (-1.2d-212)) then
tmp = x + y
else if (a <= 8.5d-259) then
tmp = (y * z) / t
else if (a <= 1.2d-95) then
tmp = x
else
tmp = x + y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.2e-212) {
tmp = x + y;
} else if (a <= 8.5e-259) {
tmp = (y * z) / t;
} else if (a <= 1.2e-95) {
tmp = x;
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -1.2e-212: tmp = x + y elif a <= 8.5e-259: tmp = (y * z) / t elif a <= 1.2e-95: tmp = x else: tmp = x + y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -1.2e-212) tmp = Float64(x + y); elseif (a <= 8.5e-259) tmp = Float64(Float64(y * z) / t); elseif (a <= 1.2e-95) tmp = x; else tmp = Float64(x + y); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -1.2e-212) tmp = x + y; elseif (a <= 8.5e-259) tmp = (y * z) / t; elseif (a <= 1.2e-95) tmp = x; else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -1.2e-212], N[(x + y), $MachinePrecision], If[LessEqual[a, 8.5e-259], N[(N[(y * z), $MachinePrecision] / t), $MachinePrecision], If[LessEqual[a, 1.2e-95], x, N[(x + y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.2 \cdot 10^{-212}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq 8.5 \cdot 10^{-259}:\\
\;\;\;\;\frac{y \cdot z}{t}\\
\mathbf{elif}\;a \leq 1.2 \cdot 10^{-95}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -1.19999999999999995e-212 or 1.2e-95 < a Initial program 82.0%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval94.9%
Simplified94.9%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6467.3%
Simplified67.3%
if -1.19999999999999995e-212 < a < 8.4999999999999994e-259Initial program 76.7%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval89.2%
Simplified89.2%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6490.1%
Simplified90.1%
Taylor expanded in x around 0
/-lowering-/.f64N/A
*-lowering-*.f6460.4%
Simplified60.4%
if 8.4999999999999994e-259 < a < 1.2e-95Initial program 64.6%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval90.8%
Simplified90.8%
Taylor expanded in x around inf
Simplified59.8%
Final simplification65.5%
(FPCore (x y z t a) :precision binary64 (if (<= a -57.0) (+ x y) (if (<= a 2.4e-50) (+ x (* y (/ z t))) (+ x y))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -57.0) {
tmp = x + y;
} else if (a <= 2.4e-50) {
tmp = x + (y * (z / t));
} else {
tmp = x + y;
}
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 <= (-57.0d0)) then
tmp = x + y
else if (a <= 2.4d-50) then
tmp = x + (y * (z / t))
else
tmp = x + y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -57.0) {
tmp = x + y;
} else if (a <= 2.4e-50) {
tmp = x + (y * (z / t));
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -57.0: tmp = x + y elif a <= 2.4e-50: tmp = x + (y * (z / t)) else: tmp = x + y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -57.0) tmp = Float64(x + y); elseif (a <= 2.4e-50) tmp = Float64(x + Float64(y * Float64(z / t))); else tmp = Float64(x + y); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -57.0) tmp = x + y; elseif (a <= 2.4e-50) tmp = x + (y * (z / t)); else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -57.0], N[(x + y), $MachinePrecision], If[LessEqual[a, 2.4e-50], N[(x + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -57:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq 2.4 \cdot 10^{-50}:\\
\;\;\;\;x + y \cdot \frac{z}{t}\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -57 or 2.40000000000000002e-50 < a Initial program 81.0%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval96.7%
Simplified96.7%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6476.1%
Simplified76.1%
if -57 < a < 2.40000000000000002e-50Initial program 76.9%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval90.4%
Simplified90.4%
Taylor expanded in a around 0
/-lowering-/.f6481.1%
Simplified81.1%
Final simplification78.5%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (* y (/ z (- t a))))) (if (<= z -8.5e+149) t_1 (if (<= z 2.7e+220) (+ x y) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = y * (z / (t - a));
double tmp;
if (z <= -8.5e+149) {
tmp = t_1;
} else if (z <= 2.7e+220) {
tmp = x + y;
} 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 (z <= (-8.5d+149)) then
tmp = t_1
else if (z <= 2.7d+220) then
tmp = x + y
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 (z <= -8.5e+149) {
tmp = t_1;
} else if (z <= 2.7e+220) {
tmp = x + y;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = y * (z / (t - a)) tmp = 0 if z <= -8.5e+149: tmp = t_1 elif z <= 2.7e+220: tmp = x + y else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(y * Float64(z / Float64(t - a))) tmp = 0.0 if (z <= -8.5e+149) tmp = t_1; elseif (z <= 2.7e+220) tmp = Float64(x + y); 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 (z <= -8.5e+149) tmp = t_1; elseif (z <= 2.7e+220) tmp = x + y; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(y * N[(z / N[(t - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -8.5e+149], t$95$1, If[LessEqual[z, 2.7e+220], N[(x + y), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \frac{z}{t - a}\\
\mathbf{if}\;z \leq -8.5 \cdot 10^{+149}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 2.7 \cdot 10^{+220}:\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -8.49999999999999956e149 or 2.6999999999999998e220 < z Initial program 76.8%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval93.2%
Simplified93.2%
Taylor expanded in z around -inf
associate-*r*N/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
sub-negN/A
+-lowering-+.f64N/A
Simplified96.1%
Taylor expanded in z around inf
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6463.9%
Simplified63.9%
if -8.49999999999999956e149 < z < 2.6999999999999998e220Initial program 79.6%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval93.8%
Simplified93.8%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6466.5%
Simplified66.5%
Final simplification66.0%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (* z (- 0.0 (/ y a))))) (if (<= z -2e+214) t_1 (if (<= z 1.75e+239) (+ x y) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = z * (0.0 - (y / a));
double tmp;
if (z <= -2e+214) {
tmp = t_1;
} else if (z <= 1.75e+239) {
tmp = x + y;
} 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 * (0.0d0 - (y / a))
if (z <= (-2d+214)) then
tmp = t_1
else if (z <= 1.75d+239) then
tmp = x + y
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 * (0.0 - (y / a));
double tmp;
if (z <= -2e+214) {
tmp = t_1;
} else if (z <= 1.75e+239) {
tmp = x + y;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = z * (0.0 - (y / a)) tmp = 0 if z <= -2e+214: tmp = t_1 elif z <= 1.75e+239: tmp = x + y else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(z * Float64(0.0 - Float64(y / a))) tmp = 0.0 if (z <= -2e+214) tmp = t_1; elseif (z <= 1.75e+239) tmp = Float64(x + y); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = z * (0.0 - (y / a)); tmp = 0.0; if (z <= -2e+214) tmp = t_1; elseif (z <= 1.75e+239) tmp = x + y; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(z * N[(0.0 - N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -2e+214], t$95$1, If[LessEqual[z, 1.75e+239], N[(x + y), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := z \cdot \left(0 - \frac{y}{a}\right)\\
\mathbf{if}\;z \leq -2 \cdot 10^{+214}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 1.75 \cdot 10^{+239}:\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.9999999999999999e214 or 1.7500000000000001e239 < z Initial program 76.6%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6475.5%
Simplified75.5%
Taylor expanded in z around inf
associate-*r/N/A
neg-mul-1N/A
/-lowering-/.f64N/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f6448.3%
Simplified48.3%
sub0-negN/A
distribute-frac-negN/A
associate-*l/N/A
distribute-lft-neg-inN/A
*-lowering-*.f64N/A
neg-lowering-neg.f64N/A
/-lowering-/.f6463.2%
Applied egg-rr63.2%
if -1.9999999999999999e214 < z < 1.7500000000000001e239Initial program 79.4%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval93.2%
Simplified93.2%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6464.1%
Simplified64.1%
Final simplification64.0%
(FPCore (x y z t a) :precision binary64 (if (<= t -2.7e+135) x (if (<= t 9.8e+89) (+ x y) x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -2.7e+135) {
tmp = x;
} else if (t <= 9.8e+89) {
tmp = x + y;
} 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 (t <= (-2.7d+135)) then
tmp = x
else if (t <= 9.8d+89) then
tmp = x + y
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 (t <= -2.7e+135) {
tmp = x;
} else if (t <= 9.8e+89) {
tmp = x + y;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -2.7e+135: tmp = x elif t <= 9.8e+89: tmp = x + y else: tmp = x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -2.7e+135) tmp = x; elseif (t <= 9.8e+89) tmp = Float64(x + y); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (t <= -2.7e+135) tmp = x; elseif (t <= 9.8e+89) tmp = x + y; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -2.7e+135], x, If[LessEqual[t, 9.8e+89], N[(x + y), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -2.7 \cdot 10^{+135}:\\
\;\;\;\;x\\
\mathbf{elif}\;t \leq 9.8 \cdot 10^{+89}:\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if t < -2.69999999999999985e135 or 9.79999999999999992e89 < t Initial program 54.5%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval89.5%
Simplified89.5%
Taylor expanded in x around inf
Simplified73.2%
if -2.69999999999999985e135 < t < 9.79999999999999992e89Initial program 89.2%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval95.4%
Simplified95.4%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6459.9%
Simplified59.9%
Final simplification63.8%
(FPCore (x y z t a) :precision binary64 (if (<= y -6.2e+90) y (if (<= y 1.2e+67) x y)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= -6.2e+90) {
tmp = y;
} else if (y <= 1.2e+67) {
tmp = x;
} else {
tmp = y;
}
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 <= (-6.2d+90)) then
tmp = y
else if (y <= 1.2d+67) then
tmp = x
else
tmp = y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= -6.2e+90) {
tmp = y;
} else if (y <= 1.2e+67) {
tmp = x;
} else {
tmp = y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if y <= -6.2e+90: tmp = y elif y <= 1.2e+67: tmp = x else: tmp = y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (y <= -6.2e+90) tmp = y; elseif (y <= 1.2e+67) tmp = x; else tmp = y; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (y <= -6.2e+90) tmp = y; elseif (y <= 1.2e+67) tmp = x; else tmp = y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[y, -6.2e+90], y, If[LessEqual[y, 1.2e+67], x, y]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -6.2 \cdot 10^{+90}:\\
\;\;\;\;y\\
\mathbf{elif}\;y \leq 1.2 \cdot 10^{+67}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if y < -6.19999999999999977e90 or 1.20000000000000001e67 < y Initial program 60.9%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval89.5%
Simplified89.5%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6442.8%
Simplified42.8%
Taylor expanded in y around inf
Simplified32.4%
if -6.19999999999999977e90 < y < 1.20000000000000001e67Initial program 90.5%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval96.3%
Simplified96.3%
Taylor expanded in x around inf
Simplified66.1%
(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 79.0%
associate--l+N/A
+-lowering-+.f64N/A
associate-*l/N/A
cancel-sign-sub-invN/A
distribute-rgt1-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
+-lowering-+.f64N/A
distribute-frac-neg2N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
unsub-negN/A
remove-double-negN/A
--lowering--.f64N/A
metadata-eval93.7%
Simplified93.7%
Taylor expanded in x around inf
Simplified47.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (- (+ y x) (* (* (- z t) (/ 1.0 (- a t))) y)))
(t_2 (- (+ x y) (/ (* (- z t) y) (- a t)))))
(if (< t_2 -1.3664970889390727e-7)
t_1
(if (< t_2 1.4754293444577233e-239)
(/ (- (* y (- a z)) (* x t)) (- a t))
t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (y + x) - (((z - t) * (1.0 / (a - t))) * y);
double t_2 = (x + y) - (((z - t) * y) / (a - t));
double tmp;
if (t_2 < -1.3664970889390727e-7) {
tmp = t_1;
} else if (t_2 < 1.4754293444577233e-239) {
tmp = ((y * (a - z)) - (x * t)) / (a - t);
} 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) :: t_2
real(8) :: tmp
t_1 = (y + x) - (((z - t) * (1.0d0 / (a - t))) * y)
t_2 = (x + y) - (((z - t) * y) / (a - t))
if (t_2 < (-1.3664970889390727d-7)) then
tmp = t_1
else if (t_2 < 1.4754293444577233d-239) then
tmp = ((y * (a - z)) - (x * t)) / (a - t)
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 + x) - (((z - t) * (1.0 / (a - t))) * y);
double t_2 = (x + y) - (((z - t) * y) / (a - t));
double tmp;
if (t_2 < -1.3664970889390727e-7) {
tmp = t_1;
} else if (t_2 < 1.4754293444577233e-239) {
tmp = ((y * (a - z)) - (x * t)) / (a - t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (y + x) - (((z - t) * (1.0 / (a - t))) * y) t_2 = (x + y) - (((z - t) * y) / (a - t)) tmp = 0 if t_2 < -1.3664970889390727e-7: tmp = t_1 elif t_2 < 1.4754293444577233e-239: tmp = ((y * (a - z)) - (x * t)) / (a - t) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(y + x) - Float64(Float64(Float64(z - t) * Float64(1.0 / Float64(a - t))) * y)) t_2 = Float64(Float64(x + y) - Float64(Float64(Float64(z - t) * y) / Float64(a - t))) tmp = 0.0 if (t_2 < -1.3664970889390727e-7) tmp = t_1; elseif (t_2 < 1.4754293444577233e-239) tmp = Float64(Float64(Float64(y * Float64(a - z)) - Float64(x * t)) / Float64(a - t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (y + x) - (((z - t) * (1.0 / (a - t))) * y); t_2 = (x + y) - (((z - t) * y) / (a - t)); tmp = 0.0; if (t_2 < -1.3664970889390727e-7) tmp = t_1; elseif (t_2 < 1.4754293444577233e-239) tmp = ((y * (a - z)) - (x * t)) / (a - t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(y + x), $MachinePrecision] - N[(N[(N[(z - t), $MachinePrecision] * N[(1.0 / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(x + y), $MachinePrecision] - N[(N[(N[(z - t), $MachinePrecision] * y), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Less[t$95$2, -1.3664970889390727e-7], t$95$1, If[Less[t$95$2, 1.4754293444577233e-239], N[(N[(N[(y * N[(a - z), $MachinePrecision]), $MachinePrecision] - N[(x * t), $MachinePrecision]), $MachinePrecision] / N[(a - t), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(y + x\right) - \left(\left(z - t\right) \cdot \frac{1}{a - t}\right) \cdot y\\
t_2 := \left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t}\\
\mathbf{if}\;t\_2 < -1.3664970889390727 \cdot 10^{-7}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_2 < 1.4754293444577233 \cdot 10^{-239}:\\
\;\;\;\;\frac{y \cdot \left(a - z\right) - x \cdot t}{a - t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
herbie shell --seed 2024158
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTick from plot-0.2.3.4, B"
:precision binary64
:alt
(! :herbie-platform default (if (< (- (+ x y) (/ (* (- z t) y) (- a t))) -13664970889390727/100000000000000000000000) (- (+ y x) (* (* (- z t) (/ 1 (- a t))) y)) (if (< (- (+ x y) (/ (* (- z t) y) (- a t))) 14754293444577233/1000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (/ (- (* y (- a z)) (* x t)) (- a t)) (- (+ y x) (* (* (- z t) (/ 1 (- a t))) y)))))
(- (+ x y) (/ (* (- z t) y) (- a t))))