
(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 12 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 (+ x (* y (+ (/ (- z t) (- t a)) 1.0))))
double code(double x, double y, double z, double t, double a) {
return x + (y * (((z - t) / (t - a)) + 1.0));
}
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) / (t - a)) + 1.0d0))
end function
public static double code(double x, double y, double z, double t, double a) {
return x + (y * (((z - t) / (t - a)) + 1.0));
}
def code(x, y, z, t, a): return x + (y * (((z - t) / (t - a)) + 1.0))
function code(x, y, z, t, a) return Float64(x + Float64(y * Float64(Float64(Float64(z - t) / Float64(t - a)) + 1.0))) end
function tmp = code(x, y, z, t, a) tmp = x + (y * (((z - t) / (t - a)) + 1.0)); end
code[x_, y_, z_, t_, a_] := N[(x + N[(y * N[(N[(N[(z - t), $MachinePrecision] / N[(t - a), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + y \cdot \left(\frac{z - t}{t - a} + 1\right)
\end{array}
Initial program 82.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-eval92.4%
Simplified92.4%
(FPCore (x y z t a) :precision binary64 (if (<= a -1.6e-59) (- (+ x y) (/ y (/ a z))) (if (<= a 2.8e-34) (+ x (* (/ y t) (- z a))) (- (+ x y) (* y (/ z a))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.6e-59) {
tmp = (x + y) - (y / (a / z));
} else if (a <= 2.8e-34) {
tmp = x + ((y / t) * (z - a));
} 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.6d-59)) then
tmp = (x + y) - (y / (a / z))
else if (a <= 2.8d-34) then
tmp = x + ((y / t) * (z - a))
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.6e-59) {
tmp = (x + y) - (y / (a / z));
} else if (a <= 2.8e-34) {
tmp = x + ((y / t) * (z - a));
} else {
tmp = (x + y) - (y * (z / a));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -1.6e-59: tmp = (x + y) - (y / (a / z)) elif a <= 2.8e-34: tmp = x + ((y / t) * (z - a)) else: tmp = (x + y) - (y * (z / a)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -1.6e-59) tmp = Float64(Float64(x + y) - Float64(y / Float64(a / z))); elseif (a <= 2.8e-34) tmp = Float64(x + Float64(Float64(y / t) * Float64(z - a))); 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.6e-59) tmp = (x + y) - (y / (a / z)); elseif (a <= 2.8e-34) tmp = x + ((y / t) * (z - a)); else tmp = (x + y) - (y * (z / a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -1.6e-59], N[(N[(x + y), $MachinePrecision] - N[(y / N[(a / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 2.8e-34], N[(x + N[(N[(y / t), $MachinePrecision] * N[(z - a), $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.6 \cdot 10^{-59}:\\
\;\;\;\;\left(x + y\right) - \frac{y}{\frac{a}{z}}\\
\mathbf{elif}\;a \leq 2.8 \cdot 10^{-34}:\\
\;\;\;\;x + \frac{y}{t} \cdot \left(z - a\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x + y\right) - y \cdot \frac{z}{a}\\
\end{array}
\end{array}
if a < -1.6e-59Initial program 86.6%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6488.5%
Simplified88.5%
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f6488.5%
Applied egg-rr88.5%
if -1.6e-59 < a < 2.79999999999999997e-34Initial program 79.3%
Taylor expanded in t around inf
Simplified84.1%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
/-lowering-/.f6484.6%
Applied egg-rr84.6%
if 2.79999999999999997e-34 < a Initial program 84.4%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6490.1%
Simplified90.1%
Final simplification87.2%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (- (+ x y) (* y (/ z a))))) (if (<= a -3.1e-61) t_1 (if (<= a 2.2e-33) (+ x (* (/ y t) (- z 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 <= -3.1e-61) {
tmp = t_1;
} else if (a <= 2.2e-33) {
tmp = x + ((y / t) * (z - a));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = (x + y) - (y * (z / a))
if (a <= (-3.1d-61)) then
tmp = t_1
else if (a <= 2.2d-33) then
tmp = x + ((y / t) * (z - a))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (x + y) - (y * (z / a));
double tmp;
if (a <= -3.1e-61) {
tmp = t_1;
} else if (a <= 2.2e-33) {
tmp = x + ((y / t) * (z - 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 <= -3.1e-61: tmp = t_1 elif a <= 2.2e-33: tmp = x + ((y / t) * (z - 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 <= -3.1e-61) tmp = t_1; elseif (a <= 2.2e-33) tmp = Float64(x + Float64(Float64(y / t) * Float64(z - 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 <= -3.1e-61) tmp = t_1; elseif (a <= 2.2e-33) tmp = x + ((y / t) * (z - 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, -3.1e-61], t$95$1, If[LessEqual[a, 2.2e-33], N[(x + N[(N[(y / t), $MachinePrecision] * N[(z - a), $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 -3.1 \cdot 10^{-61}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 2.2 \cdot 10^{-33}:\\
\;\;\;\;x + \frac{y}{t} \cdot \left(z - a\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -3.09999999999999995e-61 or 2.20000000000000005e-33 < a Initial program 85.6%
Taylor expanded in t around 0
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6489.2%
Simplified89.2%
if -3.09999999999999995e-61 < a < 2.20000000000000005e-33Initial program 79.3%
Taylor expanded in t around inf
Simplified84.1%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
/-lowering-/.f6484.6%
Applied egg-rr84.6%
Final simplification87.2%
(FPCore (x y z t a) :precision binary64 (if (<= a -3.1e+29) (+ x y) (if (<= a 4.3e-17) (+ x (* (/ y t) (- z a))) (+ x y))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -3.1e+29) {
tmp = x + y;
} else if (a <= 4.3e-17) {
tmp = x + ((y / t) * (z - 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 <= (-3.1d+29)) then
tmp = x + y
else if (a <= 4.3d-17) then
tmp = x + ((y / t) * (z - 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 <= -3.1e+29) {
tmp = x + y;
} else if (a <= 4.3e-17) {
tmp = x + ((y / t) * (z - a));
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -3.1e+29: tmp = x + y elif a <= 4.3e-17: tmp = x + ((y / t) * (z - a)) else: tmp = x + y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -3.1e+29) tmp = Float64(x + y); elseif (a <= 4.3e-17) tmp = Float64(x + Float64(Float64(y / t) * Float64(z - a))); else tmp = Float64(x + y); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -3.1e+29) tmp = x + y; elseif (a <= 4.3e-17) tmp = x + ((y / t) * (z - a)); else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -3.1e+29], N[(x + y), $MachinePrecision], If[LessEqual[a, 4.3e-17], N[(x + N[(N[(y / t), $MachinePrecision] * N[(z - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -3.1 \cdot 10^{+29}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq 4.3 \cdot 10^{-17}:\\
\;\;\;\;x + \frac{y}{t} \cdot \left(z - a\right)\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -3.0999999999999999e29 or 4.30000000000000023e-17 < a Initial program 84.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.0%
Simplified96.0%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6482.8%
Simplified82.8%
if -3.0999999999999999e29 < a < 4.30000000000000023e-17Initial program 81.3%
Taylor expanded in t around inf
Simplified79.5%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
/-lowering-/.f6480.0%
Applied egg-rr80.0%
Final simplification81.3%
(FPCore (x y z t a) :precision binary64 (if (<= a -2.7e+29) (+ x y) (if (<= a 4.5e-16) (+ x (* y (/ (- z a) t))) (+ x y))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -2.7e+29) {
tmp = x + y;
} else if (a <= 4.5e-16) {
tmp = x + (y * ((z - a) / 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.7d+29)) then
tmp = x + y
else if (a <= 4.5d-16) then
tmp = x + (y * ((z - a) / 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.7e+29) {
tmp = x + y;
} else if (a <= 4.5e-16) {
tmp = x + (y * ((z - a) / t));
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -2.7e+29: tmp = x + y elif a <= 4.5e-16: tmp = x + (y * ((z - a) / t)) else: tmp = x + y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -2.7e+29) tmp = Float64(x + y); elseif (a <= 4.5e-16) tmp = Float64(x + Float64(y * Float64(Float64(z - a) / 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.7e+29) tmp = x + y; elseif (a <= 4.5e-16) tmp = x + (y * ((z - a) / t)); else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -2.7e+29], N[(x + y), $MachinePrecision], If[LessEqual[a, 4.5e-16], N[(x + N[(y * N[(N[(z - a), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -2.7 \cdot 10^{+29}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq 4.5 \cdot 10^{-16}:\\
\;\;\;\;x + y \cdot \frac{z - a}{t}\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -2.7e29 or 4.5000000000000002e-16 < a Initial program 84.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.0%
Simplified96.0%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6482.8%
Simplified82.8%
if -2.7e29 < a < 4.5000000000000002e-16Initial program 81.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-eval89.2%
Simplified89.2%
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--.f6478.2%
Simplified78.2%
Final simplification80.3%
(FPCore (x y z t a) :precision binary64 (if (<= a -3.1e+29) (+ x y) (if (<= a 1.18e-70) (+ x (* z (/ y t))) (+ x y))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -3.1e+29) {
tmp = x + y;
} else if (a <= 1.18e-70) {
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 <= (-3.1d+29)) then
tmp = x + y
else if (a <= 1.18d-70) 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 <= -3.1e+29) {
tmp = x + y;
} else if (a <= 1.18e-70) {
tmp = x + (z * (y / t));
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -3.1e+29: tmp = x + y elif a <= 1.18e-70: tmp = x + (z * (y / t)) else: tmp = x + y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -3.1e+29) tmp = Float64(x + y); elseif (a <= 1.18e-70) 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 <= -3.1e+29) tmp = x + y; elseif (a <= 1.18e-70) tmp = x + (z * (y / t)); else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -3.1e+29], N[(x + y), $MachinePrecision], If[LessEqual[a, 1.18e-70], N[(x + N[(z * N[(y / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -3.1 \cdot 10^{+29}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq 1.18 \cdot 10^{-70}:\\
\;\;\;\;x + z \cdot \frac{y}{t}\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -3.0999999999999999e29 or 1.18e-70 < a Initial program 84.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-eval94.7%
Simplified94.7%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6479.7%
Simplified79.7%
if -3.0999999999999999e29 < a < 1.18e-70Initial program 81.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-eval90.0%
Simplified90.0%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6478.4%
Simplified78.4%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f6479.2%
Applied egg-rr79.2%
Final simplification79.5%
(FPCore (x y z t a) :precision binary64 (if (<= a -3.2e+29) (+ x y) (if (<= a 4.5e-71) (+ x (* y (/ z t))) (+ x y))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -3.2e+29) {
tmp = x + y;
} else if (a <= 4.5e-71) {
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 <= (-3.2d+29)) then
tmp = x + y
else if (a <= 4.5d-71) 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 <= -3.2e+29) {
tmp = x + y;
} else if (a <= 4.5e-71) {
tmp = x + (y * (z / t));
} else {
tmp = x + y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -3.2e+29: tmp = x + y elif a <= 4.5e-71: tmp = x + (y * (z / t)) else: tmp = x + y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -3.2e+29) tmp = Float64(x + y); elseif (a <= 4.5e-71) 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 <= -3.2e+29) tmp = x + y; elseif (a <= 4.5e-71) tmp = x + (y * (z / t)); else tmp = x + y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -3.2e+29], N[(x + y), $MachinePrecision], If[LessEqual[a, 4.5e-71], N[(x + N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -3.2 \cdot 10^{+29}:\\
\;\;\;\;x + y\\
\mathbf{elif}\;a \leq 4.5 \cdot 10^{-71}:\\
\;\;\;\;x + y \cdot \frac{z}{t}\\
\mathbf{else}:\\
\;\;\;\;x + y\\
\end{array}
\end{array}
if a < -3.19999999999999987e29 or 4.5000000000000002e-71 < a Initial program 84.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-eval94.7%
Simplified94.7%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6479.7%
Simplified79.7%
if -3.19999999999999987e29 < a < 4.5000000000000002e-71Initial program 81.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-eval90.0%
Simplified90.0%
Taylor expanded in a around 0
/-lowering-/.f6477.0%
Simplified77.0%
Final simplification78.4%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (* y (/ z (- t a))))) (if (<= z -5.5e+118) t_1 (if (<= z 2.3e+179) (+ 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 <= -5.5e+118) {
tmp = t_1;
} else if (z <= 2.3e+179) {
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 <= (-5.5d+118)) then
tmp = t_1
else if (z <= 2.3d+179) 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 <= -5.5e+118) {
tmp = t_1;
} else if (z <= 2.3e+179) {
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 <= -5.5e+118: tmp = t_1 elif z <= 2.3e+179: 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 <= -5.5e+118) tmp = t_1; elseif (z <= 2.3e+179) 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 <= -5.5e+118) tmp = t_1; elseif (z <= 2.3e+179) 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, -5.5e+118], t$95$1, If[LessEqual[z, 2.3e+179], 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 -5.5 \cdot 10^{+118}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 2.3 \cdot 10^{+179}:\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -5.5000000000000003e118 or 2.29999999999999994e179 < z Initial program 89.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-eval91.7%
Simplified91.7%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
--lowering--.f6491.6%
Applied egg-rr91.6%
Taylor expanded in z around inf
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f6470.4%
Simplified70.4%
if -5.5000000000000003e118 < z < 2.29999999999999994e179Initial program 80.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-eval92.6%
Simplified92.6%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6472.1%
Simplified72.1%
Final simplification71.7%
(FPCore (x y z t a) :precision binary64 (let* ((t_1 (/ (* y z) t))) (if (<= z -4.8e+119) t_1 (if (<= z 6e+178) (+ x y) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (y * z) / t;
double tmp;
if (z <= -4.8e+119) {
tmp = t_1;
} else if (z <= 6e+178) {
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
if (z <= (-4.8d+119)) then
tmp = t_1
else if (z <= 6d+178) 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;
double tmp;
if (z <= -4.8e+119) {
tmp = t_1;
} else if (z <= 6e+178) {
tmp = x + y;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (y * z) / t tmp = 0 if z <= -4.8e+119: tmp = t_1 elif z <= 6e+178: tmp = x + y else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(y * z) / t) tmp = 0.0 if (z <= -4.8e+119) tmp = t_1; elseif (z <= 6e+178) 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; tmp = 0.0; if (z <= -4.8e+119) tmp = t_1; elseif (z <= 6e+178) tmp = x + y; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(y * z), $MachinePrecision] / t), $MachinePrecision]}, If[LessEqual[z, -4.8e+119], t$95$1, If[LessEqual[z, 6e+178], N[(x + y), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{y \cdot z}{t}\\
\mathbf{if}\;z \leq -4.8 \cdot 10^{+119}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 6 \cdot 10^{+178}:\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -4.8e119 or 6.00000000000000031e178 < z Initial program 89.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-eval91.7%
Simplified91.7%
Taylor expanded in a around 0
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6469.3%
Simplified69.3%
Taylor expanded in x around 0
/-lowering-/.f64N/A
*-lowering-*.f6454.9%
Simplified54.9%
if -4.8e119 < z < 6.00000000000000031e178Initial program 80.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-eval92.6%
Simplified92.6%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6472.1%
Simplified72.1%
Final simplification68.2%
(FPCore (x y z t a) :precision binary64 (if (<= t -2.1e+204) x (if (<= t 2.05e+120) (+ x y) x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (t <= -2.1e+204) {
tmp = x;
} else if (t <= 2.05e+120) {
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.1d+204)) then
tmp = x
else if (t <= 2.05d+120) 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.1e+204) {
tmp = x;
} else if (t <= 2.05e+120) {
tmp = x + y;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if t <= -2.1e+204: tmp = x elif t <= 2.05e+120: tmp = x + y else: tmp = x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (t <= -2.1e+204) tmp = x; elseif (t <= 2.05e+120) 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.1e+204) tmp = x; elseif (t <= 2.05e+120) tmp = x + y; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[t, -2.1e+204], x, If[LessEqual[t, 2.05e+120], N[(x + y), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -2.1 \cdot 10^{+204}:\\
\;\;\;\;x\\
\mathbf{elif}\;t \leq 2.05 \cdot 10^{+120}:\\
\;\;\;\;x + y\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if t < -2.1e204 or 2.05e120 < t Initial program 61.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-eval88.3%
Simplified88.3%
Taylor expanded in x around inf
Simplified74.3%
if -2.1e204 < t < 2.05e120Initial program 89.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.8%
Simplified93.8%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6463.1%
Simplified63.1%
Final simplification65.9%
(FPCore (x y z t a) :precision binary64 (if (<= y -6e+165) y (if (<= y 8.8e+119) x y)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (y <= -6e+165) {
tmp = y;
} else if (y <= 8.8e+119) {
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 <= (-6d+165)) then
tmp = y
else if (y <= 8.8d+119) 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 <= -6e+165) {
tmp = y;
} else if (y <= 8.8e+119) {
tmp = x;
} else {
tmp = y;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if y <= -6e+165: tmp = y elif y <= 8.8e+119: tmp = x else: tmp = y return tmp
function code(x, y, z, t, a) tmp = 0.0 if (y <= -6e+165) tmp = y; elseif (y <= 8.8e+119) tmp = x; else tmp = y; end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (y <= -6e+165) tmp = y; elseif (y <= 8.8e+119) tmp = x; else tmp = y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[y, -6e+165], y, If[LessEqual[y, 8.8e+119], x, y]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -6 \cdot 10^{+165}:\\
\;\;\;\;y\\
\mathbf{elif}\;y \leq 8.8 \cdot 10^{+119}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;y\\
\end{array}
\end{array}
if y < -5.99999999999999981e165 or 8.8000000000000005e119 < y Initial program 67.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-eval87.6%
Simplified87.6%
Taylor expanded in a around inf
+-commutativeN/A
+-lowering-+.f6448.8%
Simplified48.8%
Taylor expanded in y around inf
Simplified47.2%
if -5.99999999999999981e165 < y < 8.8000000000000005e119Initial program 87.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-eval93.8%
Simplified93.8%
Taylor expanded in x around inf
Simplified61.7%
(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 82.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-eval92.4%
Simplified92.4%
Taylor expanded in x around inf
Simplified50.4%
(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 2024152
(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))))