
(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 16 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 (* z (+ (/ y (- t a)) (/ (* y (+ (/ t (- a t)) 1.0)) z))))))
(if (<= z -3.4e-113)
t_1
(if (<= z 2.1e-165)
(+ x (* y (+ (* (/ -1.0 (- a t)) (- z t)) 1.0)))
t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x + (z * ((y / (t - a)) + ((y * ((t / (a - t)) + 1.0)) / z)));
double tmp;
if (z <= -3.4e-113) {
tmp = t_1;
} else if (z <= 2.1e-165) {
tmp = x + (y * (((-1.0 / (a - t)) * (z - t)) + 1.0));
} 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 + (z * ((y / (t - a)) + ((y * ((t / (a - t)) + 1.0d0)) / z)))
if (z <= (-3.4d-113)) then
tmp = t_1
else if (z <= 2.1d-165) then
tmp = x + (y * ((((-1.0d0) / (a - t)) * (z - t)) + 1.0d0))
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 + (z * ((y / (t - a)) + ((y * ((t / (a - t)) + 1.0)) / z)));
double tmp;
if (z <= -3.4e-113) {
tmp = t_1;
} else if (z <= 2.1e-165) {
tmp = x + (y * (((-1.0 / (a - t)) * (z - t)) + 1.0));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x + (z * ((y / (t - a)) + ((y * ((t / (a - t)) + 1.0)) / z))) tmp = 0 if z <= -3.4e-113: tmp = t_1 elif z <= 2.1e-165: tmp = x + (y * (((-1.0 / (a - t)) * (z - t)) + 1.0)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x + Float64(z * Float64(Float64(y / Float64(t - a)) + Float64(Float64(y * Float64(Float64(t / Float64(a - t)) + 1.0)) / z)))) tmp = 0.0 if (z <= -3.4e-113) tmp = t_1; elseif (z <= 2.1e-165) tmp = Float64(x + Float64(y * Float64(Float64(Float64(-1.0 / Float64(a - t)) * Float64(z - t)) + 1.0))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x + (z * ((y / (t - a)) + ((y * ((t / (a - t)) + 1.0)) / z))); tmp = 0.0; if (z <= -3.4e-113) tmp = t_1; elseif (z <= 2.1e-165) tmp = x + (y * (((-1.0 / (a - t)) * (z - t)) + 1.0)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x + N[(z * N[(N[(y / N[(t - a), $MachinePrecision]), $MachinePrecision] + N[(N[(y * N[(N[(t / N[(a - t), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -3.4e-113], t$95$1, If[LessEqual[z, 2.1e-165], N[(x + N[(y * N[(N[(N[(-1.0 / N[(a - t), $MachinePrecision]), $MachinePrecision] * N[(z - t), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + z \cdot \left(\frac{y}{t - a} + \frac{y \cdot \left(\frac{t}{a - t} + 1\right)}{z}\right)\\
\mathbf{if}\;z \leq -3.4 \cdot 10^{-113}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 2.1 \cdot 10^{-165}:\\
\;\;\;\;x + y \cdot \left(\frac{-1}{a - t} \cdot \left(z - t\right) + 1\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -3.4000000000000002e-113 or 2.09999999999999995e-165 < z Initial program 74.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-eval85.9%
Simplified85.9%
Taylor expanded in z around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f6496.1%
Simplified96.1%
if -3.4000000000000002e-113 < z < 2.09999999999999995e-165Initial 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.0%
Simplified93.0%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f6493.0%
Applied egg-rr93.0%
Final simplification95.2%
(FPCore (x y z t a)
:precision binary64
(if (<= a -7e-10)
(+ x (* y (+ (/ (- z t) (- t a)) 1.0)))
(if (<= a 2.3e+57)
(+ x (* y (/ z (- t a))))
(+ x (* y (+ (* (/ -1.0 (- a t)) (- z t)) 1.0))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -7e-10) {
tmp = x + (y * (((z - t) / (t - a)) + 1.0));
} else if (a <= 2.3e+57) {
tmp = x + (y * (z / (t - a)));
} else {
tmp = x + (y * (((-1.0 / (a - t)) * (z - t)) + 1.0));
}
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 <= (-7d-10)) then
tmp = x + (y * (((z - t) / (t - a)) + 1.0d0))
else if (a <= 2.3d+57) then
tmp = x + (y * (z / (t - a)))
else
tmp = x + (y * ((((-1.0d0) / (a - t)) * (z - t)) + 1.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -7e-10) {
tmp = x + (y * (((z - t) / (t - a)) + 1.0));
} else if (a <= 2.3e+57) {
tmp = x + (y * (z / (t - a)));
} else {
tmp = x + (y * (((-1.0 / (a - t)) * (z - t)) + 1.0));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -7e-10: tmp = x + (y * (((z - t) / (t - a)) + 1.0)) elif a <= 2.3e+57: tmp = x + (y * (z / (t - a))) else: tmp = x + (y * (((-1.0 / (a - t)) * (z - t)) + 1.0)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -7e-10) tmp = Float64(x + Float64(y * Float64(Float64(Float64(z - t) / Float64(t - a)) + 1.0))); elseif (a <= 2.3e+57) tmp = Float64(x + Float64(y * Float64(z / Float64(t - a)))); else tmp = Float64(x + Float64(y * Float64(Float64(Float64(-1.0 / Float64(a - t)) * Float64(z - t)) + 1.0))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -7e-10) tmp = x + (y * (((z - t) / (t - a)) + 1.0)); elseif (a <= 2.3e+57) tmp = x + (y * (z / (t - a))); else tmp = x + (y * (((-1.0 / (a - t)) * (z - t)) + 1.0)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -7e-10], N[(x + N[(y * N[(N[(N[(z - t), $MachinePrecision] / N[(t - a), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 2.3e+57], N[(x + N[(y * N[(z / N[(t - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y * N[(N[(N[(-1.0 / N[(a - t), $MachinePrecision]), $MachinePrecision] * N[(z - t), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -7 \cdot 10^{-10}:\\
\;\;\;\;x + y \cdot \left(\frac{z - t}{t - a} + 1\right)\\
\mathbf{elif}\;a \leq 2.3 \cdot 10^{+57}:\\
\;\;\;\;x + y \cdot \frac{z}{t - a}\\
\mathbf{else}:\\
\;\;\;\;x + y \cdot \left(\frac{-1}{a - t} \cdot \left(z - t\right) + 1\right)\\
\end{array}
\end{array}
if a < -6.99999999999999961e-10Initial program 79.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.9%
Simplified89.9%
if -6.99999999999999961e-10 < a < 2.2999999999999999e57Initial program 73.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-eval83.9%
Simplified83.9%
Taylor expanded in z around inf
/-lowering-/.f64N/A
--lowering--.f6494.1%
Simplified94.1%
if 2.2999999999999999e57 < a Initial program 77.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-eval95.2%
Simplified95.2%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f6495.3%
Applied egg-rr95.3%
Final simplification93.6%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (+ x (* y (+ (/ (- z t) (- t a)) 1.0))))) (if (<= a -5.3e-11) t_1 (if (<= a 7e+56) (+ x (* y (/ z (- t a)))) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x + (y * (((z - t) / (t - a)) + 1.0));
double tmp;
if (a <= -5.3e-11) {
tmp = t_1;
} else if (a <= 7e+56) {
tmp = x + (y * (z / (t - 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 * (((z - t) / (t - a)) + 1.0d0))
if (a <= (-5.3d-11)) then
tmp = t_1
else if (a <= 7d+56) then
tmp = x + (y * (z / (t - 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 * (((z - t) / (t - a)) + 1.0));
double tmp;
if (a <= -5.3e-11) {
tmp = t_1;
} else if (a <= 7e+56) {
tmp = x + (y * (z / (t - a)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x + (y * (((z - t) / (t - a)) + 1.0)) tmp = 0 if a <= -5.3e-11: tmp = t_1 elif a <= 7e+56: tmp = x + (y * (z / (t - a))) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x + Float64(y * Float64(Float64(Float64(z - t) / Float64(t - a)) + 1.0))) tmp = 0.0 if (a <= -5.3e-11) tmp = t_1; elseif (a <= 7e+56) tmp = Float64(x + Float64(y * Float64(z / Float64(t - a)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x + (y * (((z - t) / (t - a)) + 1.0)); tmp = 0.0; if (a <= -5.3e-11) tmp = t_1; elseif (a <= 7e+56) tmp = x + (y * (z / (t - 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[(N[(N[(z - t), $MachinePrecision] / N[(t - a), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -5.3e-11], t$95$1, If[LessEqual[a, 7e+56], N[(x + N[(y * N[(z / N[(t - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + y \cdot \left(\frac{z - t}{t - a} + 1\right)\\
\mathbf{if}\;a \leq -5.3 \cdot 10^{-11}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 7 \cdot 10^{+56}:\\
\;\;\;\;x + y \cdot \frac{z}{t - a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -5.2999999999999998e-11 or 6.99999999999999999e56 < a Initial program 78.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-eval92.9%
Simplified92.9%
if -5.2999999999999998e-11 < a < 6.99999999999999999e56Initial program 73.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-eval83.9%
Simplified83.9%
Taylor expanded in z around inf
/-lowering-/.f64N/A
--lowering--.f6494.1%
Simplified94.1%
(FPCore (x y z t a) :precision binary64 (if (<= a -5.7e+85) (+ x (+ y (* t (/ y (- a t))))) (if (<= a 9.5e+57) (+ x (* y (/ z (- t a)))) (- (+ x y) (/ y (/ a z))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -5.7e+85) {
tmp = x + (y + (t * (y / (a - t))));
} else if (a <= 9.5e+57) {
tmp = x + (y * (z / (t - a)));
} 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 <= (-5.7d+85)) then
tmp = x + (y + (t * (y / (a - t))))
else if (a <= 9.5d+57) then
tmp = x + (y * (z / (t - a)))
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 <= -5.7e+85) {
tmp = x + (y + (t * (y / (a - t))));
} else if (a <= 9.5e+57) {
tmp = x + (y * (z / (t - a)));
} else {
tmp = (x + y) - (y / (a / z));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -5.7e+85: tmp = x + (y + (t * (y / (a - t)))) elif a <= 9.5e+57: tmp = x + (y * (z / (t - a))) else: tmp = (x + y) - (y / (a / z)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -5.7e+85) tmp = Float64(x + Float64(y + Float64(t * Float64(y / Float64(a - t))))); elseif (a <= 9.5e+57) tmp = Float64(x + Float64(y * Float64(z / Float64(t - a)))); 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 <= -5.7e+85) tmp = x + (y + (t * (y / (a - t)))); elseif (a <= 9.5e+57) tmp = x + (y * (z / (t - a))); else tmp = (x + y) - (y / (a / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -5.7e+85], N[(x + N[(y + N[(t * N[(y / N[(a - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 9.5e+57], N[(x + N[(y * N[(z / N[(t - a), $MachinePrecision]), $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 -5.7 \cdot 10^{+85}:\\
\;\;\;\;x + \left(y + t \cdot \frac{y}{a - t}\right)\\
\mathbf{elif}\;a \leq 9.5 \cdot 10^{+57}:\\
\;\;\;\;x + y \cdot \frac{z}{t - a}\\
\mathbf{else}:\\
\;\;\;\;\left(x + y\right) - \frac{y}{\frac{a}{z}}\\
\end{array}
\end{array}
if a < -5.7000000000000002e85Initial program 81.8%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
--lowering--.f6481.8%
Applied egg-rr81.8%
Taylor expanded in z around 0
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
mul-1-negN/A
remove-double-negN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6490.6%
Simplified90.6%
if -5.7000000000000002e85 < a < 9.4999999999999997e57Initial program 74.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-eval83.8%
Simplified83.8%
Taylor expanded in z around inf
/-lowering-/.f64N/A
--lowering--.f6491.2%
Simplified91.2%
if 9.4999999999999997e57 < a Initial program 77.4%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6491.4%
Simplified91.4%
--lowering--.f64N/A
+-lowering-+.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f6491.5%
Applied egg-rr91.5%
(FPCore (x y z t a) :precision binary64 (if (<= a -1.35e+79) (+ x (* y (+ (/ t (- a t)) 1.0))) (if (<= a 6.5e+58) (+ x (* y (/ z (- t a)))) (- (+ x y) (/ y (/ a z))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.35e+79) {
tmp = x + (y * ((t / (a - t)) + 1.0));
} else if (a <= 6.5e+58) {
tmp = x + (y * (z / (t - a)));
} 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 <= (-1.35d+79)) then
tmp = x + (y * ((t / (a - t)) + 1.0d0))
else if (a <= 6.5d+58) then
tmp = x + (y * (z / (t - a)))
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 <= -1.35e+79) {
tmp = x + (y * ((t / (a - t)) + 1.0));
} else if (a <= 6.5e+58) {
tmp = x + (y * (z / (t - a)));
} else {
tmp = (x + y) - (y / (a / z));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -1.35e+79: tmp = x + (y * ((t / (a - t)) + 1.0)) elif a <= 6.5e+58: tmp = x + (y * (z / (t - a))) else: tmp = (x + y) - (y / (a / z)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -1.35e+79) tmp = Float64(x + Float64(y * Float64(Float64(t / Float64(a - t)) + 1.0))); elseif (a <= 6.5e+58) tmp = Float64(x + Float64(y * Float64(z / Float64(t - a)))); 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 <= -1.35e+79) tmp = x + (y * ((t / (a - t)) + 1.0)); elseif (a <= 6.5e+58) tmp = x + (y * (z / (t - a))); else tmp = (x + y) - (y / (a / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -1.35e+79], N[(x + N[(y * N[(N[(t / N[(a - t), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 6.5e+58], N[(x + N[(y * N[(z / N[(t - a), $MachinePrecision]), $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 -1.35 \cdot 10^{+79}:\\
\;\;\;\;x + y \cdot \left(\frac{t}{a - t} + 1\right)\\
\mathbf{elif}\;a \leq 6.5 \cdot 10^{+58}:\\
\;\;\;\;x + y \cdot \frac{z}{t - a}\\
\mathbf{else}:\\
\;\;\;\;\left(x + y\right) - \frac{y}{\frac{a}{z}}\\
\end{array}
\end{array}
if a < -1.35e79Initial program 81.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.7%
Simplified93.7%
Taylor expanded in z around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f6490.5%
Simplified90.5%
if -1.35e79 < a < 6.49999999999999998e58Initial program 74.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-eval83.8%
Simplified83.8%
Taylor expanded in z around inf
/-lowering-/.f64N/A
--lowering--.f6491.2%
Simplified91.2%
if 6.49999999999999998e58 < a Initial program 77.4%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6491.4%
Simplified91.4%
--lowering--.f64N/A
+-lowering-+.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f6491.5%
Applied egg-rr91.5%
Final simplification91.2%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (- (+ x y) (/ y (/ a z))))) (if (<= a -2.3e+20) t_1 (if (<= a 4.6e+57) (+ x (* y (/ z (- t a)))) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (x + y) - (y / (a / z));
double tmp;
if (a <= -2.3e+20) {
tmp = t_1;
} else if (a <= 4.6e+57) {
tmp = x + (y * (z / (t - 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) - (y / (a / z))
if (a <= (-2.3d+20)) then
tmp = t_1
else if (a <= 4.6d+57) then
tmp = x + (y * (z / (t - 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) - (y / (a / z));
double tmp;
if (a <= -2.3e+20) {
tmp = t_1;
} else if (a <= 4.6e+57) {
tmp = x + (y * (z / (t - a)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (x + y) - (y / (a / z)) tmp = 0 if a <= -2.3e+20: tmp = t_1 elif a <= 4.6e+57: tmp = x + (y * (z / (t - a))) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(x + y) - Float64(y / Float64(a / z))) tmp = 0.0 if (a <= -2.3e+20) tmp = t_1; elseif (a <= 4.6e+57) tmp = Float64(x + Float64(y * Float64(z / Float64(t - a)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (x + y) - (y / (a / z)); tmp = 0.0; if (a <= -2.3e+20) tmp = t_1; elseif (a <= 4.6e+57) tmp = x + (y * (z / (t - a))); 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[(a / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -2.3e+20], t$95$1, If[LessEqual[a, 4.6e+57], N[(x + N[(y * N[(z / N[(t - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(x + y\right) - \frac{y}{\frac{a}{z}}\\
\mathbf{if}\;a \leq -2.3 \cdot 10^{+20}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 4.6 \cdot 10^{+57}:\\
\;\;\;\;x + y \cdot \frac{z}{t - a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -2.3e20 or 4.5999999999999998e57 < a Initial program 78.4%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6488.4%
Simplified88.4%
--lowering--.f64N/A
+-lowering-+.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f6488.4%
Applied egg-rr88.4%
if -2.3e20 < a < 4.5999999999999998e57Initial program 74.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-eval83.8%
Simplified83.8%
Taylor expanded in z around inf
/-lowering-/.f64N/A
--lowering--.f6493.1%
Simplified93.1%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (- (+ x y) (* y (/ z a))))) (if (<= a -6.6e+20) t_1 (if (<= a 9e+57) (+ x (* y (/ z (- t a)))) 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 <= -6.6e+20) {
tmp = t_1;
} else if (a <= 9e+57) {
tmp = x + (y * (z / (t - 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) - (y * (z / a))
if (a <= (-6.6d+20)) then
tmp = t_1
else if (a <= 9d+57) then
tmp = x + (y * (z / (t - 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) - (y * (z / a));
double tmp;
if (a <= -6.6e+20) {
tmp = t_1;
} else if (a <= 9e+57) {
tmp = x + (y * (z / (t - a)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (x + y) - (y * (z / a)) tmp = 0 if a <= -6.6e+20: tmp = t_1 elif a <= 9e+57: tmp = x + (y * (z / (t - a))) 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 <= -6.6e+20) tmp = t_1; elseif (a <= 9e+57) tmp = Float64(x + Float64(y * Float64(z / Float64(t - a)))); 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 <= -6.6e+20) tmp = t_1; elseif (a <= 9e+57) tmp = x + (y * (z / (t - a))); 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, -6.6e+20], t$95$1, If[LessEqual[a, 9e+57], N[(x + N[(y * N[(z / N[(t - a), $MachinePrecision]), $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 -6.6 \cdot 10^{+20}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 9 \cdot 10^{+57}:\\
\;\;\;\;x + y \cdot \frac{z}{t - a}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -6.6e20 or 8.99999999999999991e57 < a Initial program 78.4%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6488.4%
Simplified88.4%
if -6.6e20 < a < 8.99999999999999991e57Initial program 74.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-eval83.8%
Simplified83.8%
Taylor expanded in z around inf
/-lowering-/.f64N/A
--lowering--.f6493.1%
Simplified93.1%
(FPCore (x y z t a) :precision binary64 (if (<= a -2.8e+76) (+ x y) (if (<= a 6e+58) (+ x (* y (/ z (- t a)))) (+ x y))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -2.8e+76) {
tmp = x + y;
} else if (a <= 6e+58) {
tmp = x + (y * (z / (t - a)));
} 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 <= (-2.8d+76)) then
tmp = x + y
else if (a <= 6d+58) then
tmp = x + (y * (z / (t - a)))
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 <= -2.8e+76) {
tmp = x + y;
} else if (a <= 6e+58) {
tmp = x + (y * (z / (t - a)));
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -2.8e+76: tmp = x + y elif a <= 6e+58: tmp = x + (y * (z / (t - a))) else: tmp = x + y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -2.8e+76) tmp = Float64(x + y); elseif (a <= 6e+58) tmp = Float64(x + Float64(y * Float64(z / Float64(t - a)))); else tmp = Float64(x + y); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -2.8e+76) tmp = x + y; elseif (a <= 6e+58) tmp = x + (y * (z / (t - a))); else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -2.8e+76], N[(x + y), $MachinePrecision], If[LessEqual[a, 6e+58], N[(x + N[(y * N[(z / N[(t - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -2.8 \cdot 10^{+76}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq 6 \cdot 10^{+58}:\\
\;\;\;\;x + y \cdot \frac{z}{t - a}\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -2.7999999999999999e76 or 6.0000000000000005e58 < a Initial program 78.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-eval94.7%
Simplified94.7%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6483.7%
Simplified83.7%
if -2.7999999999999999e76 < a < 6.0000000000000005e58Initial program 74.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-eval83.8%
Simplified83.8%
Taylor expanded in z around inf
/-lowering-/.f64N/A
--lowering--.f6491.2%
Simplified91.2%
Final simplification88.4%
(FPCore (x y z t a) :precision binary64 (if (<= a -1.3e+20) (+ x y) (if (<= a 3.8e+57) (+ x (* z (/ y t))) (+ x y))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.3e+20) {
tmp = x + y;
} else if (a <= 3.8e+57) {
tmp = x + (z * (y / 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 <= (-1.3d+20)) then
tmp = x + y
else if (a <= 3.8d+57) then
tmp = x + (z * (y / 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 <= -1.3e+20) {
tmp = x + y;
} else if (a <= 3.8e+57) {
tmp = x + (z * (y / t));
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -1.3e+20: tmp = x + y elif a <= 3.8e+57: tmp = x + (z * (y / t)) else: tmp = x + y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -1.3e+20) tmp = Float64(x + y); elseif (a <= 3.8e+57) tmp = Float64(x + Float64(z * Float64(y / t))); else tmp = Float64(x + y); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -1.3e+20) tmp = x + y; elseif (a <= 3.8e+57) tmp = x + (z * (y / t)); else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -1.3e+20], N[(x + y), $MachinePrecision], If[LessEqual[a, 3.8e+57], N[(x + N[(z * N[(y / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.3 \cdot 10^{+20}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq 3.8 \cdot 10^{+57}:\\
\;\;\;\;x + z \cdot \frac{y}{t}\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -1.3e20 or 3.7999999999999999e57 < a Initial program 78.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.5%
Simplified93.5%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6480.5%
Simplified80.5%
if -1.3e20 < a < 3.7999999999999999e57Initial program 74.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-eval83.8%
Simplified83.8%
Taylor expanded in z around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
--lowering--.f6494.8%
Simplified94.8%
Taylor expanded in a around 0
/-lowering-/.f6484.9%
Simplified84.9%
Final simplification83.1%
(FPCore (x y z t a) :precision binary64 (if (<= a -2.8e+19) (+ x y) (if (<= a 2.9e+57) (+ x (* y (/ z t))) (+ x y))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -2.8e+19) {
tmp = x + y;
} else if (a <= 2.9e+57) {
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 <= (-2.8d+19)) then
tmp = x + y
else if (a <= 2.9d+57) 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 <= -2.8e+19) {
tmp = x + y;
} else if (a <= 2.9e+57) {
tmp = x + (y * (z / t));
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -2.8e+19: tmp = x + y elif a <= 2.9e+57: tmp = x + (y * (z / t)) else: tmp = x + y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -2.8e+19) tmp = Float64(x + y); elseif (a <= 2.9e+57) 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 <= -2.8e+19) tmp = x + y; elseif (a <= 2.9e+57) tmp = x + (y * (z / t)); else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -2.8e+19], N[(x + y), $MachinePrecision], If[LessEqual[a, 2.9e+57], N[(x + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -2.8 \cdot 10^{+19}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq 2.9 \cdot 10^{+57}:\\
\;\;\;\;x + y \cdot \frac{z}{t}\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -2.8e19 or 2.9000000000000002e57 < a Initial program 78.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.5%
Simplified93.5%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6480.5%
Simplified80.5%
if -2.8e19 < a < 2.9000000000000002e57Initial program 74.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-eval83.8%
Simplified83.8%
Taylor expanded in a around 0
/-lowering-/.f6483.7%
Simplified83.7%
Final simplification82.4%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (* z (/ y (- t a))))) (if (<= z -2.05e+18) t_1 (if (<= z 2.1e+78) (+ x y) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = z * (y / (t - a));
double tmp;
if (z <= -2.05e+18) {
tmp = t_1;
} else if (z <= 2.1e+78) {
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 * (y / (t - a))
if (z <= (-2.05d+18)) then
tmp = t_1
else if (z <= 2.1d+78) 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 * (y / (t - a));
double tmp;
if (z <= -2.05e+18) {
tmp = t_1;
} else if (z <= 2.1e+78) {
tmp = x + y;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = z * (y / (t - a)) tmp = 0 if z <= -2.05e+18: tmp = t_1 elif z <= 2.1e+78: tmp = x + y else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(z * Float64(y / Float64(t - a))) tmp = 0.0 if (z <= -2.05e+18) tmp = t_1; elseif (z <= 2.1e+78) tmp = Float64(x + y); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = z * (y / (t - a)); tmp = 0.0; if (z <= -2.05e+18) tmp = t_1; elseif (z <= 2.1e+78) 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[(y / N[(t - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -2.05e+18], t$95$1, If[LessEqual[z, 2.1e+78], N[(x + y), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := z \cdot \frac{y}{t - a}\\
\mathbf{if}\;z \leq -2.05 \cdot 10^{+18}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 2.1 \cdot 10^{+78}:\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -2.05e18 or 2.1000000000000001e78 < z Initial program 79.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-eval87.6%
Simplified87.6%
Taylor expanded in z around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6448.1%
Simplified48.1%
associate-*l/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6458.1%
Applied egg-rr58.1%
if -2.05e18 < z < 2.1000000000000001e78Initial program 72.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-eval88.0%
Simplified88.0%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6476.3%
Simplified76.3%
Final simplification67.8%
(FPCore (x y z t a) :precision binary64 (if (<= a -6.2e-225) (+ x y) (if (<= a 1.05e-36) (/ z (/ t y)) (+ x y))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -6.2e-225) {
tmp = x + y;
} else if (a <= 1.05e-36) {
tmp = z / (t / y);
} 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 <= (-6.2d-225)) then
tmp = x + y
else if (a <= 1.05d-36) then
tmp = z / (t / y)
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 <= -6.2e-225) {
tmp = x + y;
} else if (a <= 1.05e-36) {
tmp = z / (t / y);
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -6.2e-225: tmp = x + y elif a <= 1.05e-36: tmp = z / (t / y) else: tmp = x + y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -6.2e-225) tmp = Float64(x + y); elseif (a <= 1.05e-36) tmp = Float64(z / Float64(t / y)); else tmp = Float64(x + y); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -6.2e-225) tmp = x + y; elseif (a <= 1.05e-36) tmp = z / (t / y); else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -6.2e-225], N[(x + y), $MachinePrecision], If[LessEqual[a, 1.05e-36], N[(z / N[(t / y), $MachinePrecision]), $MachinePrecision], N[(x + y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -6.2 \cdot 10^{-225}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq 1.05 \cdot 10^{-36}:\\
\;\;\;\;\frac{z}{\frac{t}{y}}\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -6.19999999999999993e-225 or 1.04999999999999995e-36 < a Initial program 77.1%
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-eval91.1%
Simplified91.1%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6469.9%
Simplified69.9%
if -6.19999999999999993e-225 < a < 1.04999999999999995e-36Initial program 72.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-eval80.5%
Simplified80.5%
Taylor expanded in z around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6456.1%
Simplified56.1%
*-commutativeN/A
associate-*r/N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
--lowering--.f6462.2%
Applied egg-rr62.2%
Taylor expanded in t around inf
/-lowering-/.f6455.5%
Simplified55.5%
Final simplification65.5%
(FPCore (x y z t a) :precision binary64 (if (<= a -8.2e-225) (+ x y) (if (<= a 1.05e-36) (* y (/ z t)) (+ x y))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -8.2e-225) {
tmp = x + y;
} else if (a <= 1.05e-36) {
tmp = 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 <= (-8.2d-225)) then
tmp = x + y
else if (a <= 1.05d-36) then
tmp = 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 <= -8.2e-225) {
tmp = x + y;
} else if (a <= 1.05e-36) {
tmp = y * (z / t);
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -8.2e-225: tmp = x + y elif a <= 1.05e-36: tmp = y * (z / t) else: tmp = x + y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -8.2e-225) tmp = Float64(x + y); elseif (a <= 1.05e-36) tmp = 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 <= -8.2e-225) tmp = x + y; elseif (a <= 1.05e-36) tmp = y * (z / t); else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -8.2e-225], N[(x + y), $MachinePrecision], If[LessEqual[a, 1.05e-36], N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision], N[(x + y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -8.2 \cdot 10^{-225}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq 1.05 \cdot 10^{-36}:\\
\;\;\;\;y \cdot \frac{z}{t}\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -8.20000000000000044e-225 or 1.04999999999999995e-36 < a Initial program 77.1%
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-eval91.1%
Simplified91.1%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6469.9%
Simplified69.9%
if -8.20000000000000044e-225 < a < 1.04999999999999995e-36Initial program 72.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-eval80.5%
Simplified80.5%
Taylor expanded in z around inf
/-lowering-/.f64N/A
*-lowering-*.f64N/A
--lowering--.f6456.1%
Simplified56.1%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6459.9%
Applied egg-rr59.9%
Taylor expanded in t around inf
/-lowering-/.f6453.2%
Simplified53.2%
Final simplification64.7%
(FPCore (x y z t a) :precision binary64 (if (<= a -6.1e-38) (+ x y) (if (<= a 2.2e+57) x (+ x y))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -6.1e-38) {
tmp = x + y;
} else if (a <= 2.2e+57) {
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 <= (-6.1d-38)) then
tmp = x + y
else if (a <= 2.2d+57) 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 <= -6.1e-38) {
tmp = x + y;
} else if (a <= 2.2e+57) {
tmp = x;
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -6.1e-38: tmp = x + y elif a <= 2.2e+57: tmp = x else: tmp = x + y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -6.1e-38) tmp = Float64(x + y); elseif (a <= 2.2e+57) 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 <= -6.1e-38) tmp = x + y; elseif (a <= 2.2e+57) tmp = x; else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -6.1e-38], N[(x + y), $MachinePrecision], If[LessEqual[a, 2.2e+57], x, N[(x + y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -6.1 \cdot 10^{-38}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq 2.2 \cdot 10^{+57}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -6.09999999999999972e-38 or 2.2000000000000001e57 < a Initial program 78.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-eval92.3%
Simplified92.3%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6476.9%
Simplified76.9%
if -6.09999999999999972e-38 < a < 2.2000000000000001e57Initial program 73.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-eval84.2%
Simplified84.2%
Taylor expanded in x around inf
Simplified49.4%
Final simplification61.7%
(FPCore (x y z t a) :precision binary64 (if (<= x -2.8e-61) x (if (<= x 1.35e-210) y x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (x <= -2.8e-61) {
tmp = x;
} else if (x <= 1.35e-210) {
tmp = 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 (x <= (-2.8d-61)) then
tmp = x
else if (x <= 1.35d-210) then
tmp = 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 (x <= -2.8e-61) {
tmp = x;
} else if (x <= 1.35e-210) {
tmp = y;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if x <= -2.8e-61: tmp = x elif x <= 1.35e-210: tmp = y else: tmp = x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (x <= -2.8e-61) tmp = x; elseif (x <= 1.35e-210) tmp = y; else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (x <= -2.8e-61) tmp = x; elseif (x <= 1.35e-210) tmp = y; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[x, -2.8e-61], x, If[LessEqual[x, 1.35e-210], y, x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -2.8 \cdot 10^{-61}:\\
\;\;\;\;x\\
\mathbf{elif}\;x \leq 1.35 \cdot 10^{-210}:\\
\;\;\;\;y\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if x < -2.8000000000000001e-61 or 1.34999999999999996e-210 < x Initial program 78.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-eval92.3%
Simplified92.3%
Taylor expanded in x around inf
Simplified59.8%
if -2.8000000000000001e-61 < x < 1.34999999999999996e-210Initial program 67.3%
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-eval74.6%
Simplified74.6%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6434.8%
Simplified34.8%
Taylor expanded in y around inf
Simplified30.0%
(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 75.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-eval87.8%
Simplified87.8%
Taylor expanded in x around inf
Simplified47.1%
(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 2024138
(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))))