
(FPCore (x y z t a) :precision binary64 (+ x (/ (* (- y z) t) (- a z))))
double code(double x, double y, double z, double t, double a) {
return x + (((y - z) * t) / (a - z));
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = x + (((y - z) * t) / (a - z))
end function
public static double code(double x, double y, double z, double t, double a) {
return x + (((y - z) * t) / (a - z));
}
def code(x, y, z, t, a): return x + (((y - z) * t) / (a - z))
function code(x, y, z, t, a) return Float64(x + Float64(Float64(Float64(y - z) * t) / Float64(a - z))) end
function tmp = code(x, y, z, t, a) tmp = x + (((y - z) * t) / (a - z)); end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(N[(y - z), $MachinePrecision] * t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{\left(y - z\right) \cdot t}{a - z}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 12 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a) :precision binary64 (+ x (/ (* (- y z) t) (- a z))))
double code(double x, double y, double z, double t, double a) {
return x + (((y - z) * t) / (a - z));
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = x + (((y - z) * t) / (a - z))
end function
public static double code(double x, double y, double z, double t, double a) {
return x + (((y - z) * t) / (a - z));
}
def code(x, y, z, t, a): return x + (((y - z) * t) / (a - z))
function code(x, y, z, t, a) return Float64(x + Float64(Float64(Float64(y - z) * t) / Float64(a - z))) end
function tmp = code(x, y, z, t, a) tmp = x + (((y - z) * t) / (a - z)); end
code[x_, y_, z_, t_, a_] := N[(x + N[(N[(N[(y - z), $MachinePrecision] * t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{\left(y - z\right) \cdot t}{a - z}
\end{array}
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (* (- y z) t) (- a z))))
(if (or (<= t_1 (- INFINITY)) (not (<= t_1 2e+157)))
(+ x (/ (- y z) (/ (- a z) t)))
(+ t_1 x))))
double code(double x, double y, double z, double t, double a) {
double t_1 = ((y - z) * t) / (a - z);
double tmp;
if ((t_1 <= -((double) INFINITY)) || !(t_1 <= 2e+157)) {
tmp = x + ((y - z) / ((a - z) / t));
} else {
tmp = t_1 + x;
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a) {
double t_1 = ((y - z) * t) / (a - z);
double tmp;
if ((t_1 <= -Double.POSITIVE_INFINITY) || !(t_1 <= 2e+157)) {
tmp = x + ((y - z) / ((a - z) / t));
} else {
tmp = t_1 + x;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = ((y - z) * t) / (a - z) tmp = 0 if (t_1 <= -math.inf) or not (t_1 <= 2e+157): tmp = x + ((y - z) / ((a - z) / t)) else: tmp = t_1 + x return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(Float64(y - z) * t) / Float64(a - z)) tmp = 0.0 if ((t_1 <= Float64(-Inf)) || !(t_1 <= 2e+157)) tmp = Float64(x + Float64(Float64(y - z) / Float64(Float64(a - z) / t))); else tmp = Float64(t_1 + x); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = ((y - z) * t) / (a - z); tmp = 0.0; if ((t_1 <= -Inf) || ~((t_1 <= 2e+157))) tmp = x + ((y - z) / ((a - z) / t)); else tmp = t_1 + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(N[(y - z), $MachinePrecision] * t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[t$95$1, (-Infinity)], N[Not[LessEqual[t$95$1, 2e+157]], $MachinePrecision]], N[(x + N[(N[(y - z), $MachinePrecision] / N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 + x), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\left(y - z\right) \cdot t}{a - z}\\
\mathbf{if}\;t\_1 \leq -\infty \lor \neg \left(t\_1 \leq 2 \cdot 10^{+157}\right):\\
\;\;\;\;x + \frac{y - z}{\frac{a - z}{t}}\\
\mathbf{else}:\\
\;\;\;\;t\_1 + x\\
\end{array}
\end{array}
if (/.f64 (*.f64 (-.f64 y z) t) (-.f64 a z)) < -inf.0 or 1.99999999999999997e157 < (/.f64 (*.f64 (-.f64 y z) t) (-.f64 a z)) Initial program 41.0%
associate-/l*99.7%
Simplified99.7%
clear-num99.6%
un-div-inv99.8%
Applied egg-rr99.8%
if -inf.0 < (/.f64 (*.f64 (-.f64 y z) t) (-.f64 a z)) < 1.99999999999999997e157Initial program 99.8%
Final simplification99.8%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (* (- y z) t) (- a z))))
(if (or (<= t_1 (- INFINITY)) (not (<= t_1 2e+203)))
(+ x (* (- y z) (/ t (- a z))))
(+ t_1 x))))
double code(double x, double y, double z, double t, double a) {
double t_1 = ((y - z) * t) / (a - z);
double tmp;
if ((t_1 <= -((double) INFINITY)) || !(t_1 <= 2e+203)) {
tmp = x + ((y - z) * (t / (a - z)));
} else {
tmp = t_1 + x;
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a) {
double t_1 = ((y - z) * t) / (a - z);
double tmp;
if ((t_1 <= -Double.POSITIVE_INFINITY) || !(t_1 <= 2e+203)) {
tmp = x + ((y - z) * (t / (a - z)));
} else {
tmp = t_1 + x;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = ((y - z) * t) / (a - z) tmp = 0 if (t_1 <= -math.inf) or not (t_1 <= 2e+203): tmp = x + ((y - z) * (t / (a - z))) else: tmp = t_1 + x return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(Float64(y - z) * t) / Float64(a - z)) tmp = 0.0 if ((t_1 <= Float64(-Inf)) || !(t_1 <= 2e+203)) tmp = Float64(x + Float64(Float64(y - z) * Float64(t / Float64(a - z)))); else tmp = Float64(t_1 + x); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = ((y - z) * t) / (a - z); tmp = 0.0; if ((t_1 <= -Inf) || ~((t_1 <= 2e+203))) tmp = x + ((y - z) * (t / (a - z))); else tmp = t_1 + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(N[(y - z), $MachinePrecision] * t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[t$95$1, (-Infinity)], N[Not[LessEqual[t$95$1, 2e+203]], $MachinePrecision]], N[(x + N[(N[(y - z), $MachinePrecision] * N[(t / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 + x), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\left(y - z\right) \cdot t}{a - z}\\
\mathbf{if}\;t\_1 \leq -\infty \lor \neg \left(t\_1 \leq 2 \cdot 10^{+203}\right):\\
\;\;\;\;x + \left(y - z\right) \cdot \frac{t}{a - z}\\
\mathbf{else}:\\
\;\;\;\;t\_1 + x\\
\end{array}
\end{array}
if (/.f64 (*.f64 (-.f64 y z) t) (-.f64 a z)) < -inf.0 or 2e203 < (/.f64 (*.f64 (-.f64 y z) t) (-.f64 a z)) Initial program 35.4%
associate-/l*99.7%
Simplified99.7%
if -inf.0 < (/.f64 (*.f64 (-.f64 y z) t) (-.f64 a z)) < 2e203Initial program 99.8%
Final simplification99.8%
(FPCore (x y z t a)
:precision binary64
(if (<= z -1.05e+238)
(+ t x)
(if (or (<= z -2.8e-97) (not (<= z 7000.0)))
(+ x (* z (/ t (- z a))))
(+ x (* t (/ y a))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -1.05e+238) {
tmp = t + x;
} else if ((z <= -2.8e-97) || !(z <= 7000.0)) {
tmp = x + (z * (t / (z - a)));
} else {
tmp = x + (t * (y / a));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (z <= (-1.05d+238)) then
tmp = t + x
else if ((z <= (-2.8d-97)) .or. (.not. (z <= 7000.0d0))) then
tmp = x + (z * (t / (z - a)))
else
tmp = x + (t * (y / a))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -1.05e+238) {
tmp = t + x;
} else if ((z <= -2.8e-97) || !(z <= 7000.0)) {
tmp = x + (z * (t / (z - a)));
} else {
tmp = x + (t * (y / a));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -1.05e+238: tmp = t + x elif (z <= -2.8e-97) or not (z <= 7000.0): tmp = x + (z * (t / (z - a))) else: tmp = x + (t * (y / a)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -1.05e+238) tmp = Float64(t + x); elseif ((z <= -2.8e-97) || !(z <= 7000.0)) tmp = Float64(x + Float64(z * Float64(t / Float64(z - a)))); else tmp = Float64(x + Float64(t * Float64(y / a))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -1.05e+238) tmp = t + x; elseif ((z <= -2.8e-97) || ~((z <= 7000.0))) tmp = x + (z * (t / (z - a))); else tmp = x + (t * (y / a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -1.05e+238], N[(t + x), $MachinePrecision], If[Or[LessEqual[z, -2.8e-97], N[Not[LessEqual[z, 7000.0]], $MachinePrecision]], N[(x + N[(z * N[(t / N[(z - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.05 \cdot 10^{+238}:\\
\;\;\;\;t + x\\
\mathbf{elif}\;z \leq -2.8 \cdot 10^{-97} \lor \neg \left(z \leq 7000\right):\\
\;\;\;\;x + z \cdot \frac{t}{z - a}\\
\mathbf{else}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\end{array}
\end{array}
if z < -1.05000000000000004e238Initial program 84.6%
associate-/l*68.2%
Simplified68.2%
Taylor expanded in z around inf 93.0%
if -1.05000000000000004e238 < z < -2.8000000000000002e-97 or 7e3 < z Initial program 78.3%
associate-/l*97.3%
Simplified97.3%
Taylor expanded in y around 0 64.7%
associate-*r/64.7%
mul-1-neg64.7%
distribute-rgt-neg-out64.7%
associate-*l/80.4%
*-commutative80.4%
distribute-lft-neg-out80.4%
distribute-rgt-neg-in80.4%
distribute-frac-neg280.4%
neg-sub080.4%
sub-neg80.4%
+-commutative80.4%
associate--r+80.4%
neg-sub080.4%
remove-double-neg80.4%
Simplified80.4%
if -2.8000000000000002e-97 < z < 7e3Initial program 93.4%
associate-/l*96.6%
Simplified96.6%
Taylor expanded in z around 0 76.3%
+-commutative76.3%
associate-/l*81.2%
Simplified81.2%
Final simplification81.6%
(FPCore (x y z t a)
:precision binary64
(if (<= z -9.8e+135)
(+ t x)
(if (<= z -69000000000000.0)
(- x (* t (/ y z)))
(if (<= z 9.2e+42) (+ x (* t (/ y a))) (+ t x)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -9.8e+135) {
tmp = t + x;
} else if (z <= -69000000000000.0) {
tmp = x - (t * (y / z));
} else if (z <= 9.2e+42) {
tmp = x + (t * (y / a));
} else {
tmp = t + 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 (z <= (-9.8d+135)) then
tmp = t + x
else if (z <= (-69000000000000.0d0)) then
tmp = x - (t * (y / z))
else if (z <= 9.2d+42) then
tmp = x + (t * (y / a))
else
tmp = t + x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -9.8e+135) {
tmp = t + x;
} else if (z <= -69000000000000.0) {
tmp = x - (t * (y / z));
} else if (z <= 9.2e+42) {
tmp = x + (t * (y / a));
} else {
tmp = t + x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -9.8e+135: tmp = t + x elif z <= -69000000000000.0: tmp = x - (t * (y / z)) elif z <= 9.2e+42: tmp = x + (t * (y / a)) else: tmp = t + x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -9.8e+135) tmp = Float64(t + x); elseif (z <= -69000000000000.0) tmp = Float64(x - Float64(t * Float64(y / z))); elseif (z <= 9.2e+42) tmp = Float64(x + Float64(t * Float64(y / a))); else tmp = Float64(t + x); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -9.8e+135) tmp = t + x; elseif (z <= -69000000000000.0) tmp = x - (t * (y / z)); elseif (z <= 9.2e+42) tmp = x + (t * (y / a)); else tmp = t + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -9.8e+135], N[(t + x), $MachinePrecision], If[LessEqual[z, -69000000000000.0], N[(x - N[(t * N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 9.2e+42], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t + x), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -9.8 \cdot 10^{+135}:\\
\;\;\;\;t + x\\
\mathbf{elif}\;z \leq -69000000000000:\\
\;\;\;\;x - t \cdot \frac{y}{z}\\
\mathbf{elif}\;z \leq 9.2 \cdot 10^{+42}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\mathbf{else}:\\
\;\;\;\;t + x\\
\end{array}
\end{array}
if z < -9.8000000000000002e135 or 9.2e42 < z Initial program 71.5%
associate-/l*90.7%
Simplified90.7%
Taylor expanded in z around inf 80.5%
if -9.8000000000000002e135 < z < -6.9e13Initial program 89.8%
Taylor expanded in y around inf 75.5%
Taylor expanded in a around 0 71.8%
mul-1-neg71.8%
unsub-neg71.8%
associate-/l*76.9%
Simplified76.9%
if -6.9e13 < z < 9.2e42Initial program 92.1%
associate-/l*97.5%
Simplified97.5%
Taylor expanded in z around 0 70.8%
+-commutative70.8%
associate-/l*75.8%
Simplified75.8%
Final simplification77.8%
(FPCore (x y z t a) :precision binary64 (if (or (<= z -500.0) (not (<= z 1.02e+43))) (+ t (- x (* t (/ y z)))) (+ x (/ y (/ (- a z) t)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((z <= -500.0) || !(z <= 1.02e+43)) {
tmp = t + (x - (t * (y / z)));
} else {
tmp = x + (y / ((a - z) / t));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((z <= (-500.0d0)) .or. (.not. (z <= 1.02d+43))) then
tmp = t + (x - (t * (y / z)))
else
tmp = x + (y / ((a - z) / t))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((z <= -500.0) || !(z <= 1.02e+43)) {
tmp = t + (x - (t * (y / z)));
} else {
tmp = x + (y / ((a - z) / t));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (z <= -500.0) or not (z <= 1.02e+43): tmp = t + (x - (t * (y / z))) else: tmp = x + (y / ((a - z) / t)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((z <= -500.0) || !(z <= 1.02e+43)) tmp = Float64(t + Float64(x - Float64(t * Float64(y / z)))); else tmp = Float64(x + Float64(y / Float64(Float64(a - z) / t))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((z <= -500.0) || ~((z <= 1.02e+43))) tmp = t + (x - (t * (y / z))); else tmp = x + (y / ((a - z) / t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[z, -500.0], N[Not[LessEqual[z, 1.02e+43]], $MachinePrecision]], N[(t + N[(x - N[(t * N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y / N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -500 \lor \neg \left(z \leq 1.02 \cdot 10^{+43}\right):\\
\;\;\;\;t + \left(x - t \cdot \frac{y}{z}\right)\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{\frac{a - z}{t}}\\
\end{array}
\end{array}
if z < -500 or 1.02e43 < z Initial program 76.1%
associate-/l*92.8%
Simplified92.8%
clear-num92.6%
un-div-inv92.7%
Applied egg-rr92.7%
Taylor expanded in a around 0 87.6%
neg-mul-187.6%
+-commutative87.6%
sub-neg87.6%
Simplified87.6%
Taylor expanded in a around 0 83.2%
cancel-sign-sub-inv83.2%
mul-1-neg83.2%
associate-*r/87.4%
sub-neg87.4%
metadata-eval87.4%
*-lft-identity87.4%
Simplified87.4%
if -500 < z < 1.02e43Initial program 93.3%
associate-/l*97.4%
Simplified97.4%
clear-num97.4%
un-div-inv97.8%
Applied egg-rr97.8%
Taylor expanded in y around inf 89.1%
Final simplification88.2%
(FPCore (x y z t a) :precision binary64 (if (or (<= y -1.55e-7) (not (<= y 4.4e+74))) (+ x (/ y (/ (- a z) t))) (+ x (* z (/ t (- z a))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((y <= -1.55e-7) || !(y <= 4.4e+74)) {
tmp = x + (y / ((a - z) / t));
} else {
tmp = x + (z * (t / (z - a)));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((y <= (-1.55d-7)) .or. (.not. (y <= 4.4d+74))) then
tmp = x + (y / ((a - z) / t))
else
tmp = x + (z * (t / (z - a)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((y <= -1.55e-7) || !(y <= 4.4e+74)) {
tmp = x + (y / ((a - z) / t));
} else {
tmp = x + (z * (t / (z - a)));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (y <= -1.55e-7) or not (y <= 4.4e+74): tmp = x + (y / ((a - z) / t)) else: tmp = x + (z * (t / (z - a))) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((y <= -1.55e-7) || !(y <= 4.4e+74)) tmp = Float64(x + Float64(y / Float64(Float64(a - z) / t))); else tmp = Float64(x + Float64(z * Float64(t / Float64(z - a)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((y <= -1.55e-7) || ~((y <= 4.4e+74))) tmp = x + (y / ((a - z) / t)); else tmp = x + (z * (t / (z - a))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[y, -1.55e-7], N[Not[LessEqual[y, 4.4e+74]], $MachinePrecision]], N[(x + N[(y / N[(N[(a - z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(z * N[(t / N[(z - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.55 \cdot 10^{-7} \lor \neg \left(y \leq 4.4 \cdot 10^{+74}\right):\\
\;\;\;\;x + \frac{y}{\frac{a - z}{t}}\\
\mathbf{else}:\\
\;\;\;\;x + z \cdot \frac{t}{z - a}\\
\end{array}
\end{array}
if y < -1.55e-7 or 4.4000000000000002e74 < y Initial program 83.2%
associate-/l*95.2%
Simplified95.2%
clear-num95.1%
un-div-inv95.2%
Applied egg-rr95.2%
Taylor expanded in y around inf 87.3%
if -1.55e-7 < y < 4.4000000000000002e74Initial program 84.6%
associate-/l*94.7%
Simplified94.7%
Taylor expanded in y around 0 74.4%
associate-*r/74.4%
mul-1-neg74.4%
distribute-rgt-neg-out74.4%
associate-*l/87.9%
*-commutative87.9%
distribute-lft-neg-out87.9%
distribute-rgt-neg-in87.9%
distribute-frac-neg287.9%
neg-sub087.9%
sub-neg87.9%
+-commutative87.9%
associate--r+87.9%
neg-sub087.9%
remove-double-neg87.9%
Simplified87.9%
Final simplification87.6%
(FPCore (x y z t a) :precision binary64 (if (or (<= z -16000000000.0) (not (<= z 1.2e+43))) (+ t x) (+ x (* t (/ y a)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((z <= -16000000000.0) || !(z <= 1.2e+43)) {
tmp = t + x;
} else {
tmp = x + (t * (y / a));
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((z <= (-16000000000.0d0)) .or. (.not. (z <= 1.2d+43))) then
tmp = t + x
else
tmp = x + (t * (y / a))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((z <= -16000000000.0) || !(z <= 1.2e+43)) {
tmp = t + x;
} else {
tmp = x + (t * (y / a));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (z <= -16000000000.0) or not (z <= 1.2e+43): tmp = t + x else: tmp = x + (t * (y / a)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((z <= -16000000000.0) || !(z <= 1.2e+43)) tmp = Float64(t + x); else tmp = Float64(x + Float64(t * Float64(y / a))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((z <= -16000000000.0) || ~((z <= 1.2e+43))) tmp = t + x; else tmp = x + (t * (y / a)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[z, -16000000000.0], N[Not[LessEqual[z, 1.2e+43]], $MachinePrecision]], N[(t + x), $MachinePrecision], N[(x + N[(t * N[(y / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -16000000000 \lor \neg \left(z \leq 1.2 \cdot 10^{+43}\right):\\
\;\;\;\;t + x\\
\mathbf{else}:\\
\;\;\;\;x + t \cdot \frac{y}{a}\\
\end{array}
\end{array}
if z < -1.6e10 or 1.20000000000000012e43 < z Initial program 76.7%
associate-/l*92.6%
Simplified92.6%
Taylor expanded in z around inf 76.0%
if -1.6e10 < z < 1.20000000000000012e43Initial program 92.1%
associate-/l*97.5%
Simplified97.5%
Taylor expanded in z around 0 70.9%
+-commutative70.9%
associate-/l*76.0%
Simplified76.0%
Final simplification76.0%
(FPCore (x y z t a) :precision binary64 (if (or (<= z -54.0) (not (<= z 1.08e+43))) (+ t x) (+ x (/ (* y t) a))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((z <= -54.0) || !(z <= 1.08e+43)) {
tmp = t + x;
} else {
tmp = x + ((y * t) / a);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((z <= (-54.0d0)) .or. (.not. (z <= 1.08d+43))) then
tmp = t + x
else
tmp = x + ((y * t) / a)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((z <= -54.0) || !(z <= 1.08e+43)) {
tmp = t + x;
} else {
tmp = x + ((y * t) / a);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (z <= -54.0) or not (z <= 1.08e+43): tmp = t + x else: tmp = x + ((y * t) / a) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((z <= -54.0) || !(z <= 1.08e+43)) tmp = Float64(t + x); else tmp = Float64(x + Float64(Float64(y * t) / a)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((z <= -54.0) || ~((z <= 1.08e+43))) tmp = t + x; else tmp = x + ((y * t) / a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[z, -54.0], N[Not[LessEqual[z, 1.08e+43]], $MachinePrecision]], N[(t + x), $MachinePrecision], N[(x + N[(N[(y * t), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -54 \lor \neg \left(z \leq 1.08 \cdot 10^{+43}\right):\\
\;\;\;\;t + x\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y \cdot t}{a}\\
\end{array}
\end{array}
if z < -54 or 1.08e43 < z Initial program 76.1%
associate-/l*92.8%
Simplified92.8%
Taylor expanded in z around inf 75.4%
if -54 < z < 1.08e43Initial program 93.3%
associate-/l*97.4%
Simplified97.4%
Taylor expanded in z around 0 72.3%
Final simplification74.0%
(FPCore (x y z t a) :precision binary64 (if (or (<= x -3.8e-183) (not (<= x 4e-32))) (+ t x) (* t (- 1.0 (/ y z)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((x <= -3.8e-183) || !(x <= 4e-32)) {
tmp = t + x;
} else {
tmp = t * (1.0 - (y / 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 ((x <= (-3.8d-183)) .or. (.not. (x <= 4d-32))) then
tmp = t + x
else
tmp = t * (1.0d0 - (y / z))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((x <= -3.8e-183) || !(x <= 4e-32)) {
tmp = t + x;
} else {
tmp = t * (1.0 - (y / z));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (x <= -3.8e-183) or not (x <= 4e-32): tmp = t + x else: tmp = t * (1.0 - (y / z)) return tmp
function code(x, y, z, t, a) tmp = 0.0 if ((x <= -3.8e-183) || !(x <= 4e-32)) tmp = Float64(t + x); else tmp = Float64(t * Float64(1.0 - Float64(y / z))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((x <= -3.8e-183) || ~((x <= 4e-32))) tmp = t + x; else tmp = t * (1.0 - (y / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[Or[LessEqual[x, -3.8e-183], N[Not[LessEqual[x, 4e-32]], $MachinePrecision]], N[(t + x), $MachinePrecision], N[(t * N[(1.0 - N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -3.8 \cdot 10^{-183} \lor \neg \left(x \leq 4 \cdot 10^{-32}\right):\\
\;\;\;\;t + x\\
\mathbf{else}:\\
\;\;\;\;t \cdot \left(1 - \frac{y}{z}\right)\\
\end{array}
\end{array}
if x < -3.7999999999999996e-183 or 4.00000000000000022e-32 < x Initial program 83.1%
associate-/l*97.4%
Simplified97.4%
Taylor expanded in z around inf 74.0%
if -3.7999999999999996e-183 < x < 4.00000000000000022e-32Initial program 85.6%
associate-/l*90.1%
Simplified90.1%
clear-num89.7%
un-div-inv90.4%
Applied egg-rr90.4%
Taylor expanded in a around 0 88.0%
neg-mul-188.0%
+-commutative88.0%
sub-neg88.0%
Simplified88.0%
Taylor expanded in a around 0 66.5%
cancel-sign-sub-inv66.5%
mul-1-neg66.5%
associate-*r/68.8%
sub-neg68.8%
metadata-eval68.8%
*-lft-identity68.8%
Simplified68.8%
Taylor expanded in x around 0 55.1%
*-rgt-identity55.1%
associate-*r/57.3%
distribute-lft-out--57.3%
Simplified57.3%
Final simplification68.4%
(FPCore (x y z t a) :precision binary64 (if (<= z -1.1e+223) (+ t (- x (* t (/ y z)))) (+ x (* (- y z) (/ t (- a z))))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (z <= -1.1e+223) {
tmp = t + (x - (t * (y / z)));
} else {
tmp = x + ((y - z) * (t / (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 (z <= (-1.1d+223)) then
tmp = t + (x - (t * (y / z)))
else
tmp = x + ((y - z) * (t / (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 (z <= -1.1e+223) {
tmp = t + (x - (t * (y / z)));
} else {
tmp = x + ((y - z) * (t / (a - z)));
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if z <= -1.1e+223: tmp = t + (x - (t * (y / z))) else: tmp = x + ((y - z) * (t / (a - z))) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (z <= -1.1e+223) tmp = Float64(t + Float64(x - Float64(t * Float64(y / z)))); else tmp = Float64(x + Float64(Float64(y - z) * Float64(t / Float64(a - z)))); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (z <= -1.1e+223) tmp = t + (x - (t * (y / z))); else tmp = x + ((y - z) * (t / (a - z))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[z, -1.1e+223], N[(t + N[(x - N[(t * N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(y - z), $MachinePrecision] * N[(t / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.1 \cdot 10^{+223}:\\
\;\;\;\;t + \left(x - t \cdot \frac{y}{z}\right)\\
\mathbf{else}:\\
\;\;\;\;x + \left(y - z\right) \cdot \frac{t}{a - z}\\
\end{array}
\end{array}
if z < -1.1e223Initial program 75.0%
associate-/l*73.7%
Simplified73.7%
clear-num72.8%
un-div-inv72.8%
Applied egg-rr72.8%
Taylor expanded in a around 0 59.5%
neg-mul-159.5%
+-commutative59.5%
sub-neg59.5%
Simplified59.5%
Taylor expanded in a around 0 91.8%
cancel-sign-sub-inv91.8%
mul-1-neg91.8%
associate-*r/99.9%
sub-neg99.9%
metadata-eval99.9%
*-lft-identity99.9%
Simplified99.9%
if -1.1e223 < z Initial program 84.9%
associate-/l*97.0%
Simplified97.0%
Final simplification97.3%
(FPCore (x y z t a) :precision binary64 (if (<= a -1.62e+156) x (+ t x)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (a <= -1.62e+156) {
tmp = x;
} else {
tmp = t + 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 (a <= (-1.62d+156)) then
tmp = x
else
tmp = t + x
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.62e+156) {
tmp = x;
} else {
tmp = t + x;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if a <= -1.62e+156: tmp = x else: tmp = t + x return tmp
function code(x, y, z, t, a) tmp = 0.0 if (a <= -1.62e+156) tmp = x; else tmp = Float64(t + x); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (a <= -1.62e+156) tmp = x; else tmp = t + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[a, -1.62e+156], x, N[(t + x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq -1.62 \cdot 10^{+156}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;t + x\\
\end{array}
\end{array}
if a < -1.62000000000000006e156Initial program 79.5%
associate-/l*94.2%
Simplified94.2%
Taylor expanded in x around inf 68.9%
if -1.62000000000000006e156 < a Initial program 84.6%
associate-/l*95.0%
Simplified95.0%
Taylor expanded in z around inf 63.3%
Final simplification64.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 84.0%
associate-/l*94.9%
Simplified94.9%
Taylor expanded in x around inf 49.9%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (+ x (* (/ (- y z) (- a z)) t))))
(if (< t -1.0682974490174067e-39)
t_1
(if (< t 3.9110949887586375e-141) (+ x (/ (* (- y z) t) (- a z))) t_1))))
double code(double x, double y, double z, double t, double a) {
double t_1 = x + (((y - z) / (a - z)) * t);
double tmp;
if (t < -1.0682974490174067e-39) {
tmp = t_1;
} else if (t < 3.9110949887586375e-141) {
tmp = x + (((y - z) * t) / (a - z));
} 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) / (a - z)) * t)
if (t < (-1.0682974490174067d-39)) then
tmp = t_1
else if (t < 3.9110949887586375d-141) then
tmp = x + (((y - z) * t) / (a - z))
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) / (a - z)) * t);
double tmp;
if (t < -1.0682974490174067e-39) {
tmp = t_1;
} else if (t < 3.9110949887586375e-141) {
tmp = x + (((y - z) * t) / (a - z));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = x + (((y - z) / (a - z)) * t) tmp = 0 if t < -1.0682974490174067e-39: tmp = t_1 elif t < 3.9110949887586375e-141: tmp = x + (((y - z) * t) / (a - z)) else: tmp = t_1 return tmp
function code(x, y, z, t, a) t_1 = Float64(x + Float64(Float64(Float64(y - z) / Float64(a - z)) * t)) tmp = 0.0 if (t < -1.0682974490174067e-39) tmp = t_1; elseif (t < 3.9110949887586375e-141) tmp = Float64(x + Float64(Float64(Float64(y - z) * t) / Float64(a - z))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = x + (((y - z) / (a - z)) * t); tmp = 0.0; if (t < -1.0682974490174067e-39) tmp = t_1; elseif (t < 3.9110949887586375e-141) tmp = x + (((y - z) * t) / (a - z)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(x + N[(N[(N[(y - z), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]}, If[Less[t, -1.0682974490174067e-39], t$95$1, If[Less[t, 3.9110949887586375e-141], N[(x + N[(N[(N[(y - z), $MachinePrecision] * t), $MachinePrecision] / N[(a - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{y - z}{a - z} \cdot t\\
\mathbf{if}\;t < -1.0682974490174067 \cdot 10^{-39}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t < 3.9110949887586375 \cdot 10^{-141}:\\
\;\;\;\;x + \frac{\left(y - z\right) \cdot t}{a - z}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
herbie shell --seed 2024139
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTick from plot-0.2.3.4, A"
:precision binary64
:alt
(! :herbie-platform default (if (< t -10682974490174067/10000000000000000000000000000000000000000000000000000000) (+ x (* (/ (- y z) (- a z)) t)) (if (< t 312887599100691/80000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (+ x (/ (* (- y z) t) (- a z))) (+ x (* (/ (- y z) (- a z)) t)))))
(+ x (/ (* (- y z) t) (- a z))))