
(FPCore (x y z t) :precision binary64 (+ (- x (/ y (* z 3.0))) (/ t (* (* z 3.0) y))))
double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x - (y / (z * 3.0d0))) + (t / ((z * 3.0d0) * y))
end function
public static double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y));
}
def code(x, y, z, t): return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y))
function code(x, y, z, t) return Float64(Float64(x - Float64(y / Float64(z * 3.0))) + Float64(t / Float64(Float64(z * 3.0) * y))) end
function tmp = code(x, y, z, t) tmp = (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y)); end
code[x_, y_, z_, t_] := N[(N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t / N[(N[(z * 3.0), $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x - \frac{y}{z \cdot 3}\right) + \frac{t}{\left(z \cdot 3\right) \cdot y}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 20 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t) :precision binary64 (+ (- x (/ y (* z 3.0))) (/ t (* (* z 3.0) y))))
double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x - (y / (z * 3.0d0))) + (t / ((z * 3.0d0) * y))
end function
public static double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y));
}
def code(x, y, z, t): return (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y))
function code(x, y, z, t) return Float64(Float64(x - Float64(y / Float64(z * 3.0))) + Float64(t / Float64(Float64(z * 3.0) * y))) end
function tmp = code(x, y, z, t) tmp = (x - (y / (z * 3.0))) + (t / ((z * 3.0) * y)); end
code[x_, y_, z_, t_] := N[(N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t / N[(N[(z * 3.0), $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x - \frac{y}{z \cdot 3}\right) + \frac{t}{\left(z \cdot 3\right) \cdot y}
\end{array}
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (/ t y) y)))
(if (<= y -2.5e-99)
(+ x (/ (/ t_1 z) 3.0))
(if (<= y 2.4e-222)
(+ (/ (/ t z) (* y 3.0)) x)
(+ x (/ 1.0 (/ (* z 3.0) t_1)))))))
double code(double x, double y, double z, double t) {
double t_1 = (t / y) - y;
double tmp;
if (y <= -2.5e-99) {
tmp = x + ((t_1 / z) / 3.0);
} else if (y <= 2.4e-222) {
tmp = ((t / z) / (y * 3.0)) + x;
} else {
tmp = x + (1.0 / ((z * 3.0) / t_1));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (t / y) - y
if (y <= (-2.5d-99)) then
tmp = x + ((t_1 / z) / 3.0d0)
else if (y <= 2.4d-222) then
tmp = ((t / z) / (y * 3.0d0)) + x
else
tmp = x + (1.0d0 / ((z * 3.0d0) / t_1))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (t / y) - y;
double tmp;
if (y <= -2.5e-99) {
tmp = x + ((t_1 / z) / 3.0);
} else if (y <= 2.4e-222) {
tmp = ((t / z) / (y * 3.0)) + x;
} else {
tmp = x + (1.0 / ((z * 3.0) / t_1));
}
return tmp;
}
def code(x, y, z, t): t_1 = (t / y) - y tmp = 0 if y <= -2.5e-99: tmp = x + ((t_1 / z) / 3.0) elif y <= 2.4e-222: tmp = ((t / z) / (y * 3.0)) + x else: tmp = x + (1.0 / ((z * 3.0) / t_1)) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(t / y) - y) tmp = 0.0 if (y <= -2.5e-99) tmp = Float64(x + Float64(Float64(t_1 / z) / 3.0)); elseif (y <= 2.4e-222) tmp = Float64(Float64(Float64(t / z) / Float64(y * 3.0)) + x); else tmp = Float64(x + Float64(1.0 / Float64(Float64(z * 3.0) / t_1))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (t / y) - y; tmp = 0.0; if (y <= -2.5e-99) tmp = x + ((t_1 / z) / 3.0); elseif (y <= 2.4e-222) tmp = ((t / z) / (y * 3.0)) + x; else tmp = x + (1.0 / ((z * 3.0) / t_1)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision]}, If[LessEqual[y, -2.5e-99], N[(x + N[(N[(t$95$1 / z), $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.4e-222], N[(N[(N[(t / z), $MachinePrecision] / N[(y * 3.0), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision], N[(x + N[(1.0 / N[(N[(z * 3.0), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{y} - y\\
\mathbf{if}\;y \leq -2.5 \cdot 10^{-99}:\\
\;\;\;\;x + \frac{\frac{t\_1}{z}}{3}\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{-222}:\\
\;\;\;\;\frac{\frac{t}{z}}{y \cdot 3} + x\\
\mathbf{else}:\\
\;\;\;\;x + \frac{1}{\frac{z \cdot 3}{t\_1}}\\
\end{array}
\end{array}
if y < -2.49999999999999985e-99Initial program 95.8%
sub-neg95.8%
associate-+l+95.8%
+-commutative95.8%
remove-double-neg95.8%
distribute-frac-neg95.8%
distribute-neg-in95.8%
remove-double-neg95.8%
sub-neg95.8%
neg-mul-195.8%
times-frac98.7%
distribute-frac-neg98.7%
neg-mul-198.7%
*-commutative98.7%
associate-/l*98.7%
*-commutative98.7%
Simplified98.6%
Taylor expanded in z around 0 98.7%
*-commutative98.7%
metadata-eval98.7%
div-inv98.8%
Applied egg-rr98.8%
if -2.49999999999999985e-99 < y < 2.39999999999999993e-222Initial program 80.9%
sub-neg80.9%
associate-+l+80.9%
+-commutative80.9%
remove-double-neg80.9%
distribute-frac-neg80.9%
distribute-neg-in80.9%
remove-double-neg80.9%
sub-neg80.9%
neg-mul-180.9%
times-frac85.7%
distribute-frac-neg85.7%
neg-mul-185.7%
*-commutative85.7%
associate-/l*85.7%
*-commutative85.7%
Simplified85.9%
Taylor expanded in y around 0 80.9%
+-commutative80.9%
Simplified80.9%
*-commutative80.9%
associate-/l/99.7%
metadata-eval99.7%
div-inv99.7%
associate-/r*99.8%
Applied egg-rr99.8%
if 2.39999999999999993e-222 < y Initial program 95.8%
sub-neg95.8%
associate-+l+95.8%
+-commutative95.8%
remove-double-neg95.8%
distribute-frac-neg95.8%
distribute-neg-in95.8%
remove-double-neg95.8%
sub-neg95.8%
neg-mul-195.8%
times-frac98.8%
distribute-frac-neg98.8%
neg-mul-198.8%
*-commutative98.8%
associate-/l*98.7%
*-commutative98.7%
Simplified98.8%
Taylor expanded in z around 0 98.7%
metadata-eval98.7%
times-frac98.8%
*-un-lft-identity98.8%
*-commutative98.8%
clear-num98.9%
Applied egg-rr98.9%
Final simplification99.1%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ x (* 0.3333333333333333 (/ t (* z y))))))
(if (<= (* z 3.0) -1e+139)
t_1
(if (<= (* z 3.0) 2e-10)
(* (/ (- (/ t y) y) z) 0.3333333333333333)
(if (<= (* z 3.0) 5e+122) t_1 (+ x (/ (* y -0.3333333333333333) z)))))))
double code(double x, double y, double z, double t) {
double t_1 = x + (0.3333333333333333 * (t / (z * y)));
double tmp;
if ((z * 3.0) <= -1e+139) {
tmp = t_1;
} else if ((z * 3.0) <= 2e-10) {
tmp = (((t / y) - y) / z) * 0.3333333333333333;
} else if ((z * 3.0) <= 5e+122) {
tmp = t_1;
} else {
tmp = x + ((y * -0.3333333333333333) / z);
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = x + (0.3333333333333333d0 * (t / (z * y)))
if ((z * 3.0d0) <= (-1d+139)) then
tmp = t_1
else if ((z * 3.0d0) <= 2d-10) then
tmp = (((t / y) - y) / z) * 0.3333333333333333d0
else if ((z * 3.0d0) <= 5d+122) then
tmp = t_1
else
tmp = x + ((y * (-0.3333333333333333d0)) / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = x + (0.3333333333333333 * (t / (z * y)));
double tmp;
if ((z * 3.0) <= -1e+139) {
tmp = t_1;
} else if ((z * 3.0) <= 2e-10) {
tmp = (((t / y) - y) / z) * 0.3333333333333333;
} else if ((z * 3.0) <= 5e+122) {
tmp = t_1;
} else {
tmp = x + ((y * -0.3333333333333333) / z);
}
return tmp;
}
def code(x, y, z, t): t_1 = x + (0.3333333333333333 * (t / (z * y))) tmp = 0 if (z * 3.0) <= -1e+139: tmp = t_1 elif (z * 3.0) <= 2e-10: tmp = (((t / y) - y) / z) * 0.3333333333333333 elif (z * 3.0) <= 5e+122: tmp = t_1 else: tmp = x + ((y * -0.3333333333333333) / z) return tmp
function code(x, y, z, t) t_1 = Float64(x + Float64(0.3333333333333333 * Float64(t / Float64(z * y)))) tmp = 0.0 if (Float64(z * 3.0) <= -1e+139) tmp = t_1; elseif (Float64(z * 3.0) <= 2e-10) tmp = Float64(Float64(Float64(Float64(t / y) - y) / z) * 0.3333333333333333); elseif (Float64(z * 3.0) <= 5e+122) tmp = t_1; else tmp = Float64(x + Float64(Float64(y * -0.3333333333333333) / z)); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = x + (0.3333333333333333 * (t / (z * y))); tmp = 0.0; if ((z * 3.0) <= -1e+139) tmp = t_1; elseif ((z * 3.0) <= 2e-10) tmp = (((t / y) - y) / z) * 0.3333333333333333; elseif ((z * 3.0) <= 5e+122) tmp = t_1; else tmp = x + ((y * -0.3333333333333333) / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x + N[(0.3333333333333333 * N[(t / N[(z * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(z * 3.0), $MachinePrecision], -1e+139], t$95$1, If[LessEqual[N[(z * 3.0), $MachinePrecision], 2e-10], N[(N[(N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision] / z), $MachinePrecision] * 0.3333333333333333), $MachinePrecision], If[LessEqual[N[(z * 3.0), $MachinePrecision], 5e+122], t$95$1, N[(x + N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + 0.3333333333333333 \cdot \frac{t}{z \cdot y}\\
\mathbf{if}\;z \cdot 3 \leq -1 \cdot 10^{+139}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \cdot 3 \leq 2 \cdot 10^{-10}:\\
\;\;\;\;\frac{\frac{t}{y} - y}{z} \cdot 0.3333333333333333\\
\mathbf{elif}\;z \cdot 3 \leq 5 \cdot 10^{+122}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y \cdot -0.3333333333333333}{z}\\
\end{array}
\end{array}
if (*.f64 z 3) < -1.00000000000000003e139 or 2.00000000000000007e-10 < (*.f64 z 3) < 4.99999999999999989e122Initial program 98.3%
sub-neg98.3%
associate-+l+98.3%
+-commutative98.3%
remove-double-neg98.3%
distribute-frac-neg98.3%
distribute-neg-in98.3%
remove-double-neg98.3%
sub-neg98.3%
neg-mul-198.3%
times-frac89.0%
distribute-frac-neg89.0%
neg-mul-189.0%
*-commutative89.0%
associate-/l*89.0%
*-commutative89.0%
Simplified89.0%
Taylor expanded in y around 0 84.0%
+-commutative84.0%
Simplified84.0%
if -1.00000000000000003e139 < (*.f64 z 3) < 2.00000000000000007e-10Initial program 89.3%
sub-neg89.3%
associate-+l+89.3%
+-commutative89.3%
remove-double-neg89.3%
distribute-frac-neg89.3%
distribute-neg-in89.3%
remove-double-neg89.3%
sub-neg89.3%
neg-mul-189.3%
times-frac99.2%
distribute-frac-neg99.2%
neg-mul-199.2%
*-commutative99.2%
associate-/l*99.1%
*-commutative99.1%
Simplified99.1%
Taylor expanded in z around 0 99.1%
Taylor expanded in x around 0 91.1%
if 4.99999999999999989e122 < (*.f64 z 3) Initial program 96.6%
sub-neg96.6%
associate-+l+96.6%
+-commutative96.6%
remove-double-neg96.6%
distribute-frac-neg96.6%
distribute-neg-in96.6%
remove-double-neg96.6%
sub-neg96.6%
neg-mul-196.6%
times-frac90.2%
distribute-frac-neg90.2%
neg-mul-190.2%
*-commutative90.2%
associate-/l*90.3%
*-commutative90.3%
Simplified90.3%
Taylor expanded in t around 0 87.2%
+-commutative87.2%
Simplified87.2%
associate-*r/87.2%
Applied egg-rr87.2%
Final simplification89.0%
(FPCore (x y z t) :precision binary64 (if (or (<= (* z 3.0) -5e+83) (not (<= (* z 3.0) 5e+122))) (+ x (/ (* y -0.3333333333333333) z)) (* (/ (- (/ t y) y) z) 0.3333333333333333)))
double code(double x, double y, double z, double t) {
double tmp;
if (((z * 3.0) <= -5e+83) || !((z * 3.0) <= 5e+122)) {
tmp = x + ((y * -0.3333333333333333) / z);
} else {
tmp = (((t / y) - y) / z) * 0.3333333333333333;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (((z * 3.0d0) <= (-5d+83)) .or. (.not. ((z * 3.0d0) <= 5d+122))) then
tmp = x + ((y * (-0.3333333333333333d0)) / z)
else
tmp = (((t / y) - y) / z) * 0.3333333333333333d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (((z * 3.0) <= -5e+83) || !((z * 3.0) <= 5e+122)) {
tmp = x + ((y * -0.3333333333333333) / z);
} else {
tmp = (((t / y) - y) / z) * 0.3333333333333333;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if ((z * 3.0) <= -5e+83) or not ((z * 3.0) <= 5e+122): tmp = x + ((y * -0.3333333333333333) / z) else: tmp = (((t / y) - y) / z) * 0.3333333333333333 return tmp
function code(x, y, z, t) tmp = 0.0 if ((Float64(z * 3.0) <= -5e+83) || !(Float64(z * 3.0) <= 5e+122)) tmp = Float64(x + Float64(Float64(y * -0.3333333333333333) / z)); else tmp = Float64(Float64(Float64(Float64(t / y) - y) / z) * 0.3333333333333333); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (((z * 3.0) <= -5e+83) || ~(((z * 3.0) <= 5e+122))) tmp = x + ((y * -0.3333333333333333) / z); else tmp = (((t / y) - y) / z) * 0.3333333333333333; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[N[(z * 3.0), $MachinePrecision], -5e+83], N[Not[LessEqual[N[(z * 3.0), $MachinePrecision], 5e+122]], $MachinePrecision]], N[(x + N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision] / z), $MachinePrecision] * 0.3333333333333333), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \cdot 3 \leq -5 \cdot 10^{+83} \lor \neg \left(z \cdot 3 \leq 5 \cdot 10^{+122}\right):\\
\;\;\;\;x + \frac{y \cdot -0.3333333333333333}{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{t}{y} - y}{z} \cdot 0.3333333333333333\\
\end{array}
\end{array}
if (*.f64 z 3) < -5.00000000000000029e83 or 4.99999999999999989e122 < (*.f64 z 3) Initial program 97.4%
sub-neg97.4%
associate-+l+97.4%
+-commutative97.4%
remove-double-neg97.4%
distribute-frac-neg97.4%
distribute-neg-in97.4%
remove-double-neg97.4%
sub-neg97.4%
neg-mul-197.4%
times-frac89.8%
distribute-frac-neg89.8%
neg-mul-189.8%
*-commutative89.8%
associate-/l*89.8%
*-commutative89.8%
Simplified89.8%
Taylor expanded in t around 0 79.4%
+-commutative79.4%
Simplified79.4%
associate-*r/79.5%
Applied egg-rr79.5%
if -5.00000000000000029e83 < (*.f64 z 3) < 4.99999999999999989e122Initial program 90.2%
sub-neg90.2%
associate-+l+90.2%
+-commutative90.2%
remove-double-neg90.2%
distribute-frac-neg90.2%
distribute-neg-in90.2%
remove-double-neg90.2%
sub-neg90.2%
neg-mul-190.2%
times-frac98.2%
distribute-frac-neg98.2%
neg-mul-198.2%
*-commutative98.2%
associate-/l*98.2%
*-commutative98.2%
Simplified98.2%
Taylor expanded in z around 0 98.2%
Taylor expanded in x around 0 86.8%
Final simplification84.6%
(FPCore (x y z t) :precision binary64 (if (or (<= y -2.5e-99) (not (<= y 1.22e-220))) (+ x (* (/ (- (/ t y) y) z) 0.3333333333333333)) (+ (/ (/ t z) (* y 3.0)) x)))
double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -2.5e-99) || !(y <= 1.22e-220)) {
tmp = x + ((((t / y) - y) / z) * 0.3333333333333333);
} else {
tmp = ((t / z) / (y * 3.0)) + x;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if ((y <= (-2.5d-99)) .or. (.not. (y <= 1.22d-220))) then
tmp = x + ((((t / y) - y) / z) * 0.3333333333333333d0)
else
tmp = ((t / z) / (y * 3.0d0)) + x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -2.5e-99) || !(y <= 1.22e-220)) {
tmp = x + ((((t / y) - y) / z) * 0.3333333333333333);
} else {
tmp = ((t / z) / (y * 3.0)) + x;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (y <= -2.5e-99) or not (y <= 1.22e-220): tmp = x + ((((t / y) - y) / z) * 0.3333333333333333) else: tmp = ((t / z) / (y * 3.0)) + x return tmp
function code(x, y, z, t) tmp = 0.0 if ((y <= -2.5e-99) || !(y <= 1.22e-220)) tmp = Float64(x + Float64(Float64(Float64(Float64(t / y) - y) / z) * 0.3333333333333333)); else tmp = Float64(Float64(Float64(t / z) / Float64(y * 3.0)) + x); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((y <= -2.5e-99) || ~((y <= 1.22e-220))) tmp = x + ((((t / y) - y) / z) * 0.3333333333333333); else tmp = ((t / z) / (y * 3.0)) + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[y, -2.5e-99], N[Not[LessEqual[y, 1.22e-220]], $MachinePrecision]], N[(x + N[(N[(N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision] / z), $MachinePrecision] * 0.3333333333333333), $MachinePrecision]), $MachinePrecision], N[(N[(N[(t / z), $MachinePrecision] / N[(y * 3.0), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.5 \cdot 10^{-99} \lor \neg \left(y \leq 1.22 \cdot 10^{-220}\right):\\
\;\;\;\;x + \frac{\frac{t}{y} - y}{z} \cdot 0.3333333333333333\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{t}{z}}{y \cdot 3} + x\\
\end{array}
\end{array}
if y < -2.49999999999999985e-99 or 1.22000000000000002e-220 < y Initial program 96.3%
sub-neg96.3%
associate-+l+96.3%
+-commutative96.3%
remove-double-neg96.3%
distribute-frac-neg96.3%
distribute-neg-in96.3%
remove-double-neg96.3%
sub-neg96.3%
neg-mul-196.3%
times-frac98.8%
distribute-frac-neg98.8%
neg-mul-198.8%
*-commutative98.8%
associate-/l*98.7%
*-commutative98.7%
Simplified98.7%
Taylor expanded in z around 0 98.7%
if -2.49999999999999985e-99 < y < 1.22000000000000002e-220Initial program 79.6%
sub-neg79.6%
associate-+l+79.6%
+-commutative79.6%
remove-double-neg79.6%
distribute-frac-neg79.6%
distribute-neg-in79.6%
remove-double-neg79.6%
sub-neg79.6%
neg-mul-179.6%
times-frac85.9%
distribute-frac-neg85.9%
neg-mul-185.9%
*-commutative85.9%
associate-/l*85.9%
*-commutative85.9%
Simplified86.1%
Taylor expanded in y around 0 79.7%
+-commutative79.7%
Simplified79.7%
*-commutative79.7%
associate-/l/99.7%
metadata-eval99.7%
div-inv99.7%
associate-/r*99.8%
Applied egg-rr99.8%
Final simplification99.0%
(FPCore (x y z t) :precision binary64 (if (or (<= y -3.9e-98) (not (<= y 1.8e-219))) (+ x (/ (* (- (/ t y) y) 0.3333333333333333) z)) (+ (/ (/ t z) (* y 3.0)) x)))
double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -3.9e-98) || !(y <= 1.8e-219)) {
tmp = x + ((((t / y) - y) * 0.3333333333333333) / z);
} else {
tmp = ((t / z) / (y * 3.0)) + x;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if ((y <= (-3.9d-98)) .or. (.not. (y <= 1.8d-219))) then
tmp = x + ((((t / y) - y) * 0.3333333333333333d0) / z)
else
tmp = ((t / z) / (y * 3.0d0)) + x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -3.9e-98) || !(y <= 1.8e-219)) {
tmp = x + ((((t / y) - y) * 0.3333333333333333) / z);
} else {
tmp = ((t / z) / (y * 3.0)) + x;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (y <= -3.9e-98) or not (y <= 1.8e-219): tmp = x + ((((t / y) - y) * 0.3333333333333333) / z) else: tmp = ((t / z) / (y * 3.0)) + x return tmp
function code(x, y, z, t) tmp = 0.0 if ((y <= -3.9e-98) || !(y <= 1.8e-219)) tmp = Float64(x + Float64(Float64(Float64(Float64(t / y) - y) * 0.3333333333333333) / z)); else tmp = Float64(Float64(Float64(t / z) / Float64(y * 3.0)) + x); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((y <= -3.9e-98) || ~((y <= 1.8e-219))) tmp = x + ((((t / y) - y) * 0.3333333333333333) / z); else tmp = ((t / z) / (y * 3.0)) + x; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[y, -3.9e-98], N[Not[LessEqual[y, 1.8e-219]], $MachinePrecision]], N[(x + N[(N[(N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision] * 0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], N[(N[(N[(t / z), $MachinePrecision] / N[(y * 3.0), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -3.9 \cdot 10^{-98} \lor \neg \left(y \leq 1.8 \cdot 10^{-219}\right):\\
\;\;\;\;x + \frac{\left(\frac{t}{y} - y\right) \cdot 0.3333333333333333}{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{t}{z}}{y \cdot 3} + x\\
\end{array}
\end{array}
if y < -3.89999999999999971e-98 or 1.79999999999999987e-219 < y Initial program 96.3%
sub-neg96.3%
associate-+l+96.3%
+-commutative96.3%
remove-double-neg96.3%
distribute-frac-neg96.3%
distribute-neg-in96.3%
remove-double-neg96.3%
sub-neg96.3%
neg-mul-196.3%
times-frac98.8%
distribute-frac-neg98.8%
neg-mul-198.8%
*-commutative98.8%
associate-/l*98.7%
*-commutative98.7%
Simplified98.7%
associate-*l/98.8%
Applied egg-rr98.8%
if -3.89999999999999971e-98 < y < 1.79999999999999987e-219Initial program 79.6%
sub-neg79.6%
associate-+l+79.6%
+-commutative79.6%
remove-double-neg79.6%
distribute-frac-neg79.6%
distribute-neg-in79.6%
remove-double-neg79.6%
sub-neg79.6%
neg-mul-179.6%
times-frac85.9%
distribute-frac-neg85.9%
neg-mul-185.9%
*-commutative85.9%
associate-/l*85.9%
*-commutative85.9%
Simplified86.1%
Taylor expanded in y around 0 79.7%
+-commutative79.7%
Simplified79.7%
*-commutative79.7%
associate-/l/99.7%
metadata-eval99.7%
div-inv99.7%
associate-/r*99.8%
Applied egg-rr99.8%
Final simplification99.0%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (/ t y) y)))
(if (<= y -4.6e-98)
(+ x (* (/ t_1 z) 0.3333333333333333))
(if (<= y 4.8e-220)
(+ (/ (/ t z) (* y 3.0)) x)
(+ x (* t_1 (/ 0.3333333333333333 z)))))))
double code(double x, double y, double z, double t) {
double t_1 = (t / y) - y;
double tmp;
if (y <= -4.6e-98) {
tmp = x + ((t_1 / z) * 0.3333333333333333);
} else if (y <= 4.8e-220) {
tmp = ((t / z) / (y * 3.0)) + x;
} else {
tmp = x + (t_1 * (0.3333333333333333 / z));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (t / y) - y
if (y <= (-4.6d-98)) then
tmp = x + ((t_1 / z) * 0.3333333333333333d0)
else if (y <= 4.8d-220) then
tmp = ((t / z) / (y * 3.0d0)) + x
else
tmp = x + (t_1 * (0.3333333333333333d0 / z))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (t / y) - y;
double tmp;
if (y <= -4.6e-98) {
tmp = x + ((t_1 / z) * 0.3333333333333333);
} else if (y <= 4.8e-220) {
tmp = ((t / z) / (y * 3.0)) + x;
} else {
tmp = x + (t_1 * (0.3333333333333333 / z));
}
return tmp;
}
def code(x, y, z, t): t_1 = (t / y) - y tmp = 0 if y <= -4.6e-98: tmp = x + ((t_1 / z) * 0.3333333333333333) elif y <= 4.8e-220: tmp = ((t / z) / (y * 3.0)) + x else: tmp = x + (t_1 * (0.3333333333333333 / z)) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(t / y) - y) tmp = 0.0 if (y <= -4.6e-98) tmp = Float64(x + Float64(Float64(t_1 / z) * 0.3333333333333333)); elseif (y <= 4.8e-220) tmp = Float64(Float64(Float64(t / z) / Float64(y * 3.0)) + x); else tmp = Float64(x + Float64(t_1 * Float64(0.3333333333333333 / z))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (t / y) - y; tmp = 0.0; if (y <= -4.6e-98) tmp = x + ((t_1 / z) * 0.3333333333333333); elseif (y <= 4.8e-220) tmp = ((t / z) / (y * 3.0)) + x; else tmp = x + (t_1 * (0.3333333333333333 / z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision]}, If[LessEqual[y, -4.6e-98], N[(x + N[(N[(t$95$1 / z), $MachinePrecision] * 0.3333333333333333), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 4.8e-220], N[(N[(N[(t / z), $MachinePrecision] / N[(y * 3.0), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision], N[(x + N[(t$95$1 * N[(0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{y} - y\\
\mathbf{if}\;y \leq -4.6 \cdot 10^{-98}:\\
\;\;\;\;x + \frac{t\_1}{z} \cdot 0.3333333333333333\\
\mathbf{elif}\;y \leq 4.8 \cdot 10^{-220}:\\
\;\;\;\;\frac{\frac{t}{z}}{y \cdot 3} + x\\
\mathbf{else}:\\
\;\;\;\;x + t\_1 \cdot \frac{0.3333333333333333}{z}\\
\end{array}
\end{array}
if y < -4.60000000000000001e-98Initial program 95.8%
sub-neg95.8%
associate-+l+95.8%
+-commutative95.8%
remove-double-neg95.8%
distribute-frac-neg95.8%
distribute-neg-in95.8%
remove-double-neg95.8%
sub-neg95.8%
neg-mul-195.8%
times-frac98.7%
distribute-frac-neg98.7%
neg-mul-198.7%
*-commutative98.7%
associate-/l*98.7%
*-commutative98.7%
Simplified98.6%
Taylor expanded in z around 0 98.7%
if -4.60000000000000001e-98 < y < 4.8000000000000003e-220Initial program 79.6%
sub-neg79.6%
associate-+l+79.6%
+-commutative79.6%
remove-double-neg79.6%
distribute-frac-neg79.6%
distribute-neg-in79.6%
remove-double-neg79.6%
sub-neg79.6%
neg-mul-179.6%
times-frac85.9%
distribute-frac-neg85.9%
neg-mul-185.9%
*-commutative85.9%
associate-/l*85.9%
*-commutative85.9%
Simplified86.1%
Taylor expanded in y around 0 79.7%
+-commutative79.7%
Simplified79.7%
*-commutative79.7%
associate-/l/99.7%
metadata-eval99.7%
div-inv99.7%
associate-/r*99.8%
Applied egg-rr99.8%
if 4.8000000000000003e-220 < y Initial program 96.7%
sub-neg96.7%
associate-+l+96.7%
+-commutative96.7%
remove-double-neg96.7%
distribute-frac-neg96.7%
distribute-neg-in96.7%
remove-double-neg96.7%
sub-neg96.7%
neg-mul-196.7%
times-frac98.8%
distribute-frac-neg98.8%
neg-mul-198.8%
*-commutative98.8%
associate-/l*98.7%
*-commutative98.7%
Simplified98.8%
Final simplification99.0%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (/ t y) y)))
(if (<= y -1.9e-98)
(+ x (/ t_1 (* z 3.0)))
(if (<= y 5.8e-219)
(+ (/ (/ t z) (* y 3.0)) x)
(+ x (/ (* t_1 0.3333333333333333) z))))))
double code(double x, double y, double z, double t) {
double t_1 = (t / y) - y;
double tmp;
if (y <= -1.9e-98) {
tmp = x + (t_1 / (z * 3.0));
} else if (y <= 5.8e-219) {
tmp = ((t / z) / (y * 3.0)) + x;
} else {
tmp = x + ((t_1 * 0.3333333333333333) / z);
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (t / y) - y
if (y <= (-1.9d-98)) then
tmp = x + (t_1 / (z * 3.0d0))
else if (y <= 5.8d-219) then
tmp = ((t / z) / (y * 3.0d0)) + x
else
tmp = x + ((t_1 * 0.3333333333333333d0) / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (t / y) - y;
double tmp;
if (y <= -1.9e-98) {
tmp = x + (t_1 / (z * 3.0));
} else if (y <= 5.8e-219) {
tmp = ((t / z) / (y * 3.0)) + x;
} else {
tmp = x + ((t_1 * 0.3333333333333333) / z);
}
return tmp;
}
def code(x, y, z, t): t_1 = (t / y) - y tmp = 0 if y <= -1.9e-98: tmp = x + (t_1 / (z * 3.0)) elif y <= 5.8e-219: tmp = ((t / z) / (y * 3.0)) + x else: tmp = x + ((t_1 * 0.3333333333333333) / z) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(t / y) - y) tmp = 0.0 if (y <= -1.9e-98) tmp = Float64(x + Float64(t_1 / Float64(z * 3.0))); elseif (y <= 5.8e-219) tmp = Float64(Float64(Float64(t / z) / Float64(y * 3.0)) + x); else tmp = Float64(x + Float64(Float64(t_1 * 0.3333333333333333) / z)); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (t / y) - y; tmp = 0.0; if (y <= -1.9e-98) tmp = x + (t_1 / (z * 3.0)); elseif (y <= 5.8e-219) tmp = ((t / z) / (y * 3.0)) + x; else tmp = x + ((t_1 * 0.3333333333333333) / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision]}, If[LessEqual[y, -1.9e-98], N[(x + N[(t$95$1 / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 5.8e-219], N[(N[(N[(t / z), $MachinePrecision] / N[(y * 3.0), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision], N[(x + N[(N[(t$95$1 * 0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{y} - y\\
\mathbf{if}\;y \leq -1.9 \cdot 10^{-98}:\\
\;\;\;\;x + \frac{t\_1}{z \cdot 3}\\
\mathbf{elif}\;y \leq 5.8 \cdot 10^{-219}:\\
\;\;\;\;\frac{\frac{t}{z}}{y \cdot 3} + x\\
\mathbf{else}:\\
\;\;\;\;x + \frac{t\_1 \cdot 0.3333333333333333}{z}\\
\end{array}
\end{array}
if y < -1.9000000000000002e-98Initial program 95.8%
sub-neg95.8%
associate-+l+95.8%
+-commutative95.8%
remove-double-neg95.8%
distribute-frac-neg95.8%
distribute-neg-in95.8%
remove-double-neg95.8%
sub-neg95.8%
neg-mul-195.8%
times-frac98.7%
distribute-frac-neg98.7%
neg-mul-198.7%
*-commutative98.7%
associate-/l*98.7%
*-commutative98.7%
Simplified98.6%
*-commutative98.6%
clear-num98.6%
div-inv98.7%
metadata-eval98.7%
un-div-inv98.7%
Applied egg-rr98.7%
if -1.9000000000000002e-98 < y < 5.79999999999999968e-219Initial program 79.6%
sub-neg79.6%
associate-+l+79.6%
+-commutative79.6%
remove-double-neg79.6%
distribute-frac-neg79.6%
distribute-neg-in79.6%
remove-double-neg79.6%
sub-neg79.6%
neg-mul-179.6%
times-frac85.9%
distribute-frac-neg85.9%
neg-mul-185.9%
*-commutative85.9%
associate-/l*85.9%
*-commutative85.9%
Simplified86.1%
Taylor expanded in y around 0 79.7%
+-commutative79.7%
Simplified79.7%
*-commutative79.7%
associate-/l/99.7%
metadata-eval99.7%
div-inv99.7%
associate-/r*99.8%
Applied egg-rr99.8%
if 5.79999999999999968e-219 < y Initial program 96.7%
sub-neg96.7%
associate-+l+96.7%
+-commutative96.7%
remove-double-neg96.7%
distribute-frac-neg96.7%
distribute-neg-in96.7%
remove-double-neg96.7%
sub-neg96.7%
neg-mul-196.7%
times-frac98.8%
distribute-frac-neg98.8%
neg-mul-198.8%
*-commutative98.8%
associate-/l*98.7%
*-commutative98.7%
Simplified98.8%
associate-*l/98.9%
Applied egg-rr98.9%
Final simplification99.0%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (/ t y) y)))
(if (<= y -1.35e-98)
(+ x (/ (/ t_1 z) 3.0))
(if (<= y 1.25e-219)
(+ (/ (/ t z) (* y 3.0)) x)
(+ x (/ (* t_1 0.3333333333333333) z))))))
double code(double x, double y, double z, double t) {
double t_1 = (t / y) - y;
double tmp;
if (y <= -1.35e-98) {
tmp = x + ((t_1 / z) / 3.0);
} else if (y <= 1.25e-219) {
tmp = ((t / z) / (y * 3.0)) + x;
} else {
tmp = x + ((t_1 * 0.3333333333333333) / z);
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (t / y) - y
if (y <= (-1.35d-98)) then
tmp = x + ((t_1 / z) / 3.0d0)
else if (y <= 1.25d-219) then
tmp = ((t / z) / (y * 3.0d0)) + x
else
tmp = x + ((t_1 * 0.3333333333333333d0) / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (t / y) - y;
double tmp;
if (y <= -1.35e-98) {
tmp = x + ((t_1 / z) / 3.0);
} else if (y <= 1.25e-219) {
tmp = ((t / z) / (y * 3.0)) + x;
} else {
tmp = x + ((t_1 * 0.3333333333333333) / z);
}
return tmp;
}
def code(x, y, z, t): t_1 = (t / y) - y tmp = 0 if y <= -1.35e-98: tmp = x + ((t_1 / z) / 3.0) elif y <= 1.25e-219: tmp = ((t / z) / (y * 3.0)) + x else: tmp = x + ((t_1 * 0.3333333333333333) / z) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(t / y) - y) tmp = 0.0 if (y <= -1.35e-98) tmp = Float64(x + Float64(Float64(t_1 / z) / 3.0)); elseif (y <= 1.25e-219) tmp = Float64(Float64(Float64(t / z) / Float64(y * 3.0)) + x); else tmp = Float64(x + Float64(Float64(t_1 * 0.3333333333333333) / z)); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (t / y) - y; tmp = 0.0; if (y <= -1.35e-98) tmp = x + ((t_1 / z) / 3.0); elseif (y <= 1.25e-219) tmp = ((t / z) / (y * 3.0)) + x; else tmp = x + ((t_1 * 0.3333333333333333) / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision]}, If[LessEqual[y, -1.35e-98], N[(x + N[(N[(t$95$1 / z), $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.25e-219], N[(N[(N[(t / z), $MachinePrecision] / N[(y * 3.0), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision], N[(x + N[(N[(t$95$1 * 0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{y} - y\\
\mathbf{if}\;y \leq -1.35 \cdot 10^{-98}:\\
\;\;\;\;x + \frac{\frac{t\_1}{z}}{3}\\
\mathbf{elif}\;y \leq 1.25 \cdot 10^{-219}:\\
\;\;\;\;\frac{\frac{t}{z}}{y \cdot 3} + x\\
\mathbf{else}:\\
\;\;\;\;x + \frac{t\_1 \cdot 0.3333333333333333}{z}\\
\end{array}
\end{array}
if y < -1.3499999999999999e-98Initial program 95.8%
sub-neg95.8%
associate-+l+95.8%
+-commutative95.8%
remove-double-neg95.8%
distribute-frac-neg95.8%
distribute-neg-in95.8%
remove-double-neg95.8%
sub-neg95.8%
neg-mul-195.8%
times-frac98.7%
distribute-frac-neg98.7%
neg-mul-198.7%
*-commutative98.7%
associate-/l*98.7%
*-commutative98.7%
Simplified98.6%
Taylor expanded in z around 0 98.7%
*-commutative98.7%
metadata-eval98.7%
div-inv98.8%
Applied egg-rr98.8%
if -1.3499999999999999e-98 < y < 1.25e-219Initial program 79.6%
sub-neg79.6%
associate-+l+79.6%
+-commutative79.6%
remove-double-neg79.6%
distribute-frac-neg79.6%
distribute-neg-in79.6%
remove-double-neg79.6%
sub-neg79.6%
neg-mul-179.6%
times-frac85.9%
distribute-frac-neg85.9%
neg-mul-185.9%
*-commutative85.9%
associate-/l*85.9%
*-commutative85.9%
Simplified86.1%
Taylor expanded in y around 0 79.7%
+-commutative79.7%
Simplified79.7%
*-commutative79.7%
associate-/l/99.7%
metadata-eval99.7%
div-inv99.7%
associate-/r*99.8%
Applied egg-rr99.8%
if 1.25e-219 < y Initial program 96.7%
sub-neg96.7%
associate-+l+96.7%
+-commutative96.7%
remove-double-neg96.7%
distribute-frac-neg96.7%
distribute-neg-in96.7%
remove-double-neg96.7%
sub-neg96.7%
neg-mul-196.7%
times-frac98.8%
distribute-frac-neg98.8%
neg-mul-198.8%
*-commutative98.8%
associate-/l*98.7%
*-commutative98.7%
Simplified98.8%
associate-*l/98.9%
Applied egg-rr98.9%
Final simplification99.1%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (/ y (* z -3.0))))
(if (<= t 2e+102)
(+ (+ (/ (/ t z) (* y 3.0)) x) t_1)
(+ t_1 (+ x (/ t (* z (* y 3.0))))))))
double code(double x, double y, double z, double t) {
double t_1 = y / (z * -3.0);
double tmp;
if (t <= 2e+102) {
tmp = (((t / z) / (y * 3.0)) + x) + t_1;
} else {
tmp = t_1 + (x + (t / (z * (y * 3.0))));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = y / (z * (-3.0d0))
if (t <= 2d+102) then
tmp = (((t / z) / (y * 3.0d0)) + x) + t_1
else
tmp = t_1 + (x + (t / (z * (y * 3.0d0))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = y / (z * -3.0);
double tmp;
if (t <= 2e+102) {
tmp = (((t / z) / (y * 3.0)) + x) + t_1;
} else {
tmp = t_1 + (x + (t / (z * (y * 3.0))));
}
return tmp;
}
def code(x, y, z, t): t_1 = y / (z * -3.0) tmp = 0 if t <= 2e+102: tmp = (((t / z) / (y * 3.0)) + x) + t_1 else: tmp = t_1 + (x + (t / (z * (y * 3.0)))) return tmp
function code(x, y, z, t) t_1 = Float64(y / Float64(z * -3.0)) tmp = 0.0 if (t <= 2e+102) tmp = Float64(Float64(Float64(Float64(t / z) / Float64(y * 3.0)) + x) + t_1); else tmp = Float64(t_1 + Float64(x + Float64(t / Float64(z * Float64(y * 3.0))))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = y / (z * -3.0); tmp = 0.0; if (t <= 2e+102) tmp = (((t / z) / (y * 3.0)) + x) + t_1; else tmp = t_1 + (x + (t / (z * (y * 3.0)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(y / N[(z * -3.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, 2e+102], N[(N[(N[(N[(t / z), $MachinePrecision] / N[(y * 3.0), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision] + t$95$1), $MachinePrecision], N[(t$95$1 + N[(x + N[(t / N[(z * N[(y * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{y}{z \cdot -3}\\
\mathbf{if}\;t \leq 2 \cdot 10^{+102}:\\
\;\;\;\;\left(\frac{\frac{t}{z}}{y \cdot 3} + x\right) + t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_1 + \left(x + \frac{t}{z \cdot \left(y \cdot 3\right)}\right)\\
\end{array}
\end{array}
if t < 1.99999999999999995e102Initial program 91.4%
+-commutative91.4%
associate-+r-91.4%
sub-neg91.4%
associate-*l*91.4%
*-commutative91.4%
distribute-frac-neg291.4%
distribute-rgt-neg-in91.4%
metadata-eval91.4%
Simplified91.4%
associate-/r*98.9%
div-inv98.9%
Applied egg-rr98.9%
un-div-inv98.9%
Applied egg-rr98.9%
if 1.99999999999999995e102 < t Initial program 97.5%
+-commutative97.5%
associate-+r-97.5%
sub-neg97.5%
associate-*l*97.6%
*-commutative97.6%
distribute-frac-neg297.6%
distribute-rgt-neg-in97.6%
metadata-eval97.6%
Simplified97.6%
Final simplification98.7%
(FPCore (x y z t) :precision binary64 (if (<= t 1.8e-53) (+ x (/ (/ (- (/ t y) y) z) 3.0)) (+ (/ y (* z -3.0)) (+ x (/ t (* z (* y 3.0)))))))
double code(double x, double y, double z, double t) {
double tmp;
if (t <= 1.8e-53) {
tmp = x + ((((t / y) - y) / z) / 3.0);
} else {
tmp = (y / (z * -3.0)) + (x + (t / (z * (y * 3.0))));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (t <= 1.8d-53) then
tmp = x + ((((t / y) - y) / z) / 3.0d0)
else
tmp = (y / (z * (-3.0d0))) + (x + (t / (z * (y * 3.0d0))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (t <= 1.8e-53) {
tmp = x + ((((t / y) - y) / z) / 3.0);
} else {
tmp = (y / (z * -3.0)) + (x + (t / (z * (y * 3.0))));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if t <= 1.8e-53: tmp = x + ((((t / y) - y) / z) / 3.0) else: tmp = (y / (z * -3.0)) + (x + (t / (z * (y * 3.0)))) return tmp
function code(x, y, z, t) tmp = 0.0 if (t <= 1.8e-53) tmp = Float64(x + Float64(Float64(Float64(Float64(t / y) - y) / z) / 3.0)); else tmp = Float64(Float64(y / Float64(z * -3.0)) + Float64(x + Float64(t / Float64(z * Float64(y * 3.0))))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (t <= 1.8e-53) tmp = x + ((((t / y) - y) / z) / 3.0); else tmp = (y / (z * -3.0)) + (x + (t / (z * (y * 3.0)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[t, 1.8e-53], N[(x + N[(N[(N[(N[(t / y), $MachinePrecision] - y), $MachinePrecision] / z), $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision], N[(N[(y / N[(z * -3.0), $MachinePrecision]), $MachinePrecision] + N[(x + N[(t / N[(z * N[(y * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 1.8 \cdot 10^{-53}:\\
\;\;\;\;x + \frac{\frac{\frac{t}{y} - y}{z}}{3}\\
\mathbf{else}:\\
\;\;\;\;\frac{y}{z \cdot -3} + \left(x + \frac{t}{z \cdot \left(y \cdot 3\right)}\right)\\
\end{array}
\end{array}
if t < 1.7999999999999999e-53Initial program 90.4%
sub-neg90.4%
associate-+l+90.4%
+-commutative90.4%
remove-double-neg90.4%
distribute-frac-neg90.4%
distribute-neg-in90.4%
remove-double-neg90.4%
sub-neg90.4%
neg-mul-190.4%
times-frac97.8%
distribute-frac-neg97.8%
neg-mul-197.8%
*-commutative97.8%
associate-/l*97.8%
*-commutative97.8%
Simplified97.8%
Taylor expanded in z around 0 97.8%
*-commutative97.8%
metadata-eval97.8%
div-inv97.9%
Applied egg-rr97.9%
if 1.7999999999999999e-53 < t Initial program 98.2%
+-commutative98.2%
associate-+r-98.2%
sub-neg98.2%
associate-*l*98.4%
*-commutative98.4%
distribute-frac-neg298.4%
distribute-rgt-neg-in98.4%
metadata-eval98.4%
Simplified98.4%
Final simplification98.0%
(FPCore (x y z t)
:precision binary64
(if (<= y -65000000000.0)
(+ x (/ (* y -0.3333333333333333) z))
(if (<= y 25000.0)
(+ x (/ (* (/ t y) 0.3333333333333333) z))
(+ x (/ y (* z -3.0))))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -65000000000.0) {
tmp = x + ((y * -0.3333333333333333) / z);
} else if (y <= 25000.0) {
tmp = x + (((t / y) * 0.3333333333333333) / z);
} else {
tmp = x + (y / (z * -3.0));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (y <= (-65000000000.0d0)) then
tmp = x + ((y * (-0.3333333333333333d0)) / z)
else if (y <= 25000.0d0) then
tmp = x + (((t / y) * 0.3333333333333333d0) / z)
else
tmp = x + (y / (z * (-3.0d0)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -65000000000.0) {
tmp = x + ((y * -0.3333333333333333) / z);
} else if (y <= 25000.0) {
tmp = x + (((t / y) * 0.3333333333333333) / z);
} else {
tmp = x + (y / (z * -3.0));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -65000000000.0: tmp = x + ((y * -0.3333333333333333) / z) elif y <= 25000.0: tmp = x + (((t / y) * 0.3333333333333333) / z) else: tmp = x + (y / (z * -3.0)) return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -65000000000.0) tmp = Float64(x + Float64(Float64(y * -0.3333333333333333) / z)); elseif (y <= 25000.0) tmp = Float64(x + Float64(Float64(Float64(t / y) * 0.3333333333333333) / z)); else tmp = Float64(x + Float64(y / Float64(z * -3.0))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -65000000000.0) tmp = x + ((y * -0.3333333333333333) / z); elseif (y <= 25000.0) tmp = x + (((t / y) * 0.3333333333333333) / z); else tmp = x + (y / (z * -3.0)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -65000000000.0], N[(x + N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 25000.0], N[(x + N[(N[(N[(t / y), $MachinePrecision] * 0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], N[(x + N[(y / N[(z * -3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -65000000000:\\
\;\;\;\;x + \frac{y \cdot -0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 25000:\\
\;\;\;\;x + \frac{\frac{t}{y} \cdot 0.3333333333333333}{z}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{z \cdot -3}\\
\end{array}
\end{array}
if y < -6.5e10Initial program 97.1%
sub-neg97.1%
associate-+l+97.1%
+-commutative97.1%
remove-double-neg97.1%
distribute-frac-neg97.1%
distribute-neg-in97.1%
remove-double-neg97.1%
sub-neg97.1%
neg-mul-197.1%
times-frac99.7%
distribute-frac-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.7%
*-commutative99.7%
Simplified99.7%
Taylor expanded in t around 0 92.7%
+-commutative92.7%
Simplified92.7%
associate-*r/92.8%
Applied egg-rr92.8%
if -6.5e10 < y < 25000Initial program 86.8%
sub-neg86.8%
associate-+l+86.8%
+-commutative86.8%
remove-double-neg86.8%
distribute-frac-neg86.8%
distribute-neg-in86.8%
remove-double-neg86.8%
sub-neg86.8%
neg-mul-186.8%
times-frac92.0%
distribute-frac-neg92.0%
neg-mul-192.0%
*-commutative92.0%
associate-/l*92.0%
*-commutative92.0%
Simplified92.1%
Taylor expanded in y around 0 82.8%
+-commutative82.8%
Simplified82.8%
*-commutative82.8%
associate-/r*87.3%
associate-*l/87.3%
Applied egg-rr87.3%
if 25000 < y Initial program 99.8%
+-commutative99.8%
associate-+r-99.8%
sub-neg99.8%
associate-*l*99.8%
*-commutative99.8%
distribute-frac-neg299.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
associate-/r*90.8%
div-inv90.7%
Applied egg-rr90.7%
Taylor expanded in t around 0 92.9%
Final simplification89.9%
(FPCore (x y z t) :precision binary64 (if (<= y -180000000000.0) (+ x (/ (* y -0.3333333333333333) z)) (if (<= y 550000.0) (+ (/ (/ t z) (* y 3.0)) x) (+ x (/ y (* z -3.0))))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -180000000000.0) {
tmp = x + ((y * -0.3333333333333333) / z);
} else if (y <= 550000.0) {
tmp = ((t / z) / (y * 3.0)) + x;
} else {
tmp = x + (y / (z * -3.0));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (y <= (-180000000000.0d0)) then
tmp = x + ((y * (-0.3333333333333333d0)) / z)
else if (y <= 550000.0d0) then
tmp = ((t / z) / (y * 3.0d0)) + x
else
tmp = x + (y / (z * (-3.0d0)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -180000000000.0) {
tmp = x + ((y * -0.3333333333333333) / z);
} else if (y <= 550000.0) {
tmp = ((t / z) / (y * 3.0)) + x;
} else {
tmp = x + (y / (z * -3.0));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -180000000000.0: tmp = x + ((y * -0.3333333333333333) / z) elif y <= 550000.0: tmp = ((t / z) / (y * 3.0)) + x else: tmp = x + (y / (z * -3.0)) return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -180000000000.0) tmp = Float64(x + Float64(Float64(y * -0.3333333333333333) / z)); elseif (y <= 550000.0) tmp = Float64(Float64(Float64(t / z) / Float64(y * 3.0)) + x); else tmp = Float64(x + Float64(y / Float64(z * -3.0))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -180000000000.0) tmp = x + ((y * -0.3333333333333333) / z); elseif (y <= 550000.0) tmp = ((t / z) / (y * 3.0)) + x; else tmp = x + (y / (z * -3.0)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -180000000000.0], N[(x + N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 550000.0], N[(N[(N[(t / z), $MachinePrecision] / N[(y * 3.0), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision], N[(x + N[(y / N[(z * -3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -180000000000:\\
\;\;\;\;x + \frac{y \cdot -0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 550000:\\
\;\;\;\;\frac{\frac{t}{z}}{y \cdot 3} + x\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{z \cdot -3}\\
\end{array}
\end{array}
if y < -1.8e11Initial program 97.1%
sub-neg97.1%
associate-+l+97.1%
+-commutative97.1%
remove-double-neg97.1%
distribute-frac-neg97.1%
distribute-neg-in97.1%
remove-double-neg97.1%
sub-neg97.1%
neg-mul-197.1%
times-frac99.7%
distribute-frac-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.7%
*-commutative99.7%
Simplified99.7%
Taylor expanded in t around 0 92.7%
+-commutative92.7%
Simplified92.7%
associate-*r/92.8%
Applied egg-rr92.8%
if -1.8e11 < y < 5.5e5Initial program 86.8%
sub-neg86.8%
associate-+l+86.8%
+-commutative86.8%
remove-double-neg86.8%
distribute-frac-neg86.8%
distribute-neg-in86.8%
remove-double-neg86.8%
sub-neg86.8%
neg-mul-186.8%
times-frac92.0%
distribute-frac-neg92.0%
neg-mul-192.0%
*-commutative92.0%
associate-/l*92.0%
*-commutative92.0%
Simplified92.1%
Taylor expanded in y around 0 82.8%
+-commutative82.8%
Simplified82.8%
*-commutative82.8%
associate-/l/94.2%
metadata-eval94.2%
div-inv94.3%
associate-/r*94.3%
Applied egg-rr94.3%
if 5.5e5 < y Initial program 99.8%
+-commutative99.8%
associate-+r-99.8%
sub-neg99.8%
associate-*l*99.8%
*-commutative99.8%
distribute-frac-neg299.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
associate-/r*90.8%
div-inv90.7%
Applied egg-rr90.7%
Taylor expanded in t around 0 92.9%
Final simplification93.6%
(FPCore (x y z t) :precision binary64 (if (or (<= y -1.6e-35) (not (<= y 3600.0))) (+ x (* y (/ -0.3333333333333333 z))) (* 0.3333333333333333 (/ t (* z y)))))
double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -1.6e-35) || !(y <= 3600.0)) {
tmp = x + (y * (-0.3333333333333333 / z));
} else {
tmp = 0.3333333333333333 * (t / (z * y));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if ((y <= (-1.6d-35)) .or. (.not. (y <= 3600.0d0))) then
tmp = x + (y * ((-0.3333333333333333d0) / z))
else
tmp = 0.3333333333333333d0 * (t / (z * y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -1.6e-35) || !(y <= 3600.0)) {
tmp = x + (y * (-0.3333333333333333 / z));
} else {
tmp = 0.3333333333333333 * (t / (z * y));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (y <= -1.6e-35) or not (y <= 3600.0): tmp = x + (y * (-0.3333333333333333 / z)) else: tmp = 0.3333333333333333 * (t / (z * y)) return tmp
function code(x, y, z, t) tmp = 0.0 if ((y <= -1.6e-35) || !(y <= 3600.0)) tmp = Float64(x + Float64(y * Float64(-0.3333333333333333 / z))); else tmp = Float64(0.3333333333333333 * Float64(t / Float64(z * y))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((y <= -1.6e-35) || ~((y <= 3600.0))) tmp = x + (y * (-0.3333333333333333 / z)); else tmp = 0.3333333333333333 * (t / (z * y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[y, -1.6e-35], N[Not[LessEqual[y, 3600.0]], $MachinePrecision]], N[(x + N[(y * N[(-0.3333333333333333 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.3333333333333333 * N[(t / N[(z * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.6 \cdot 10^{-35} \lor \neg \left(y \leq 3600\right):\\
\;\;\;\;x + y \cdot \frac{-0.3333333333333333}{z}\\
\mathbf{else}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{t}{z \cdot y}\\
\end{array}
\end{array}
if y < -1.5999999999999999e-35 or 3600 < y Initial program 98.4%
sub-neg98.4%
associate-+l+98.4%
+-commutative98.4%
remove-double-neg98.4%
distribute-frac-neg98.4%
distribute-neg-in98.4%
remove-double-neg98.4%
sub-neg98.4%
neg-mul-198.4%
times-frac99.8%
distribute-frac-neg99.8%
neg-mul-199.8%
*-commutative99.8%
associate-/l*99.7%
*-commutative99.7%
Simplified99.7%
Taylor expanded in t around 0 91.1%
associate-*r/91.2%
*-commutative91.2%
associate-*r/91.1%
Simplified91.1%
if -1.5999999999999999e-35 < y < 3600Initial program 85.8%
sub-neg85.8%
associate-+l+85.8%
+-commutative85.8%
remove-double-neg85.8%
distribute-frac-neg85.8%
distribute-neg-in85.8%
remove-double-neg85.8%
sub-neg85.8%
neg-mul-185.8%
times-frac91.4%
distribute-frac-neg91.4%
neg-mul-191.4%
*-commutative91.4%
associate-/l*91.4%
*-commutative91.4%
Simplified91.5%
Taylor expanded in z around 0 91.4%
Taylor expanded in x around 0 73.6%
Taylor expanded in t around inf 66.3%
Final simplification79.1%
(FPCore (x y z t) :precision binary64 (if (or (<= y -3.1e-35) (not (<= y 450.0))) (+ x (/ y (* z -3.0))) (* 0.3333333333333333 (/ t (* z y)))))
double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -3.1e-35) || !(y <= 450.0)) {
tmp = x + (y / (z * -3.0));
} else {
tmp = 0.3333333333333333 * (t / (z * y));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if ((y <= (-3.1d-35)) .or. (.not. (y <= 450.0d0))) then
tmp = x + (y / (z * (-3.0d0)))
else
tmp = 0.3333333333333333d0 * (t / (z * y))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((y <= -3.1e-35) || !(y <= 450.0)) {
tmp = x + (y / (z * -3.0));
} else {
tmp = 0.3333333333333333 * (t / (z * y));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (y <= -3.1e-35) or not (y <= 450.0): tmp = x + (y / (z * -3.0)) else: tmp = 0.3333333333333333 * (t / (z * y)) return tmp
function code(x, y, z, t) tmp = 0.0 if ((y <= -3.1e-35) || !(y <= 450.0)) tmp = Float64(x + Float64(y / Float64(z * -3.0))); else tmp = Float64(0.3333333333333333 * Float64(t / Float64(z * y))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((y <= -3.1e-35) || ~((y <= 450.0))) tmp = x + (y / (z * -3.0)); else tmp = 0.3333333333333333 * (t / (z * y)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[Or[LessEqual[y, -3.1e-35], N[Not[LessEqual[y, 450.0]], $MachinePrecision]], N[(x + N[(y / N[(z * -3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.3333333333333333 * N[(t / N[(z * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -3.1 \cdot 10^{-35} \lor \neg \left(y \leq 450\right):\\
\;\;\;\;x + \frac{y}{z \cdot -3}\\
\mathbf{else}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{t}{z \cdot y}\\
\end{array}
\end{array}
if y < -3.10000000000000012e-35 or 450 < y Initial program 98.4%
+-commutative98.4%
associate-+r-98.4%
sub-neg98.4%
associate-*l*98.4%
*-commutative98.4%
distribute-frac-neg298.4%
distribute-rgt-neg-in98.4%
metadata-eval98.4%
Simplified98.4%
associate-/r*94.5%
div-inv94.5%
Applied egg-rr94.5%
Taylor expanded in t around 0 91.2%
if -3.10000000000000012e-35 < y < 450Initial program 85.8%
sub-neg85.8%
associate-+l+85.8%
+-commutative85.8%
remove-double-neg85.8%
distribute-frac-neg85.8%
distribute-neg-in85.8%
remove-double-neg85.8%
sub-neg85.8%
neg-mul-185.8%
times-frac91.4%
distribute-frac-neg91.4%
neg-mul-191.4%
*-commutative91.4%
associate-/l*91.4%
*-commutative91.4%
Simplified91.5%
Taylor expanded in z around 0 91.4%
Taylor expanded in x around 0 73.6%
Taylor expanded in t around inf 66.3%
Final simplification79.1%
(FPCore (x y z t) :precision binary64 (if (<= y -400000000000.0) (/ (/ y z) -3.0) (if (<= y 9.8e+34) (* 0.3333333333333333 (/ t (* z y))) (/ (/ y -3.0) z))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -400000000000.0) {
tmp = (y / z) / -3.0;
} else if (y <= 9.8e+34) {
tmp = 0.3333333333333333 * (t / (z * y));
} else {
tmp = (y / -3.0) / z;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (y <= (-400000000000.0d0)) then
tmp = (y / z) / (-3.0d0)
else if (y <= 9.8d+34) then
tmp = 0.3333333333333333d0 * (t / (z * y))
else
tmp = (y / (-3.0d0)) / z
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -400000000000.0) {
tmp = (y / z) / -3.0;
} else if (y <= 9.8e+34) {
tmp = 0.3333333333333333 * (t / (z * y));
} else {
tmp = (y / -3.0) / z;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -400000000000.0: tmp = (y / z) / -3.0 elif y <= 9.8e+34: tmp = 0.3333333333333333 * (t / (z * y)) else: tmp = (y / -3.0) / z return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -400000000000.0) tmp = Float64(Float64(y / z) / -3.0); elseif (y <= 9.8e+34) tmp = Float64(0.3333333333333333 * Float64(t / Float64(z * y))); else tmp = Float64(Float64(y / -3.0) / z); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -400000000000.0) tmp = (y / z) / -3.0; elseif (y <= 9.8e+34) tmp = 0.3333333333333333 * (t / (z * y)); else tmp = (y / -3.0) / z; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -400000000000.0], N[(N[(y / z), $MachinePrecision] / -3.0), $MachinePrecision], If[LessEqual[y, 9.8e+34], N[(0.3333333333333333 * N[(t / N[(z * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y / -3.0), $MachinePrecision] / z), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -400000000000:\\
\;\;\;\;\frac{\frac{y}{z}}{-3}\\
\mathbf{elif}\;y \leq 9.8 \cdot 10^{+34}:\\
\;\;\;\;0.3333333333333333 \cdot \frac{t}{z \cdot y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{y}{-3}}{z}\\
\end{array}
\end{array}
if y < -4e11Initial program 97.1%
sub-neg97.1%
associate-+l+97.1%
+-commutative97.1%
remove-double-neg97.1%
distribute-frac-neg97.1%
distribute-neg-in97.1%
remove-double-neg97.1%
sub-neg97.1%
neg-mul-197.1%
times-frac99.7%
distribute-frac-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.7%
*-commutative99.7%
Simplified99.7%
Taylor expanded in t around 0 92.7%
+-commutative92.7%
Simplified92.7%
Taylor expanded in y around inf 71.4%
associate-*r/71.5%
*-commutative71.5%
Simplified71.5%
associate-/l*71.3%
*-commutative71.3%
Applied egg-rr71.3%
metadata-eval71.3%
associate-/r*71.4%
*-commutative71.4%
associate-*l/71.4%
*-un-lft-identity71.4%
associate-/r*71.5%
Applied egg-rr71.5%
if -4e11 < y < 9.8000000000000005e34Initial program 87.5%
sub-neg87.5%
associate-+l+87.5%
+-commutative87.5%
remove-double-neg87.5%
distribute-frac-neg87.5%
distribute-neg-in87.5%
remove-double-neg87.5%
sub-neg87.5%
neg-mul-187.5%
times-frac92.4%
distribute-frac-neg92.4%
neg-mul-192.4%
*-commutative92.4%
associate-/l*92.4%
*-commutative92.4%
Simplified92.5%
Taylor expanded in z around 0 92.4%
Taylor expanded in x around 0 70.7%
Taylor expanded in t around inf 62.2%
if 9.8000000000000005e34 < y Initial program 99.8%
sub-neg99.8%
associate-+l+99.8%
+-commutative99.8%
remove-double-neg99.8%
distribute-frac-neg99.8%
distribute-neg-in99.8%
remove-double-neg99.8%
sub-neg99.8%
neg-mul-199.8%
times-frac99.9%
distribute-frac-neg99.9%
neg-mul-199.9%
*-commutative99.9%
associate-/l*99.7%
*-commutative99.7%
Simplified99.7%
Taylor expanded in t around 0 95.8%
+-commutative95.8%
Simplified95.8%
Taylor expanded in y around inf 73.0%
associate-*r/73.1%
*-commutative73.1%
Simplified73.1%
associate-/l*73.0%
*-commutative73.0%
Applied egg-rr73.0%
associate-*l/73.1%
associate-/l*73.0%
metadata-eval73.0%
times-frac73.1%
*-un-lft-identity73.1%
associate-/r*73.2%
Applied egg-rr73.2%
Final simplification66.7%
(FPCore (x y z t)
:precision binary64
(if (<= y -1.06e-34)
(+ x (/ (* y -0.3333333333333333) z))
(if (<= y 480.0)
(* 0.3333333333333333 (/ t (* z y)))
(+ x (/ y (* z -3.0))))))
double code(double x, double y, double z, double t) {
double tmp;
if (y <= -1.06e-34) {
tmp = x + ((y * -0.3333333333333333) / z);
} else if (y <= 480.0) {
tmp = 0.3333333333333333 * (t / (z * y));
} else {
tmp = x + (y / (z * -3.0));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (y <= (-1.06d-34)) then
tmp = x + ((y * (-0.3333333333333333d0)) / z)
else if (y <= 480.0d0) then
tmp = 0.3333333333333333d0 * (t / (z * y))
else
tmp = x + (y / (z * (-3.0d0)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (y <= -1.06e-34) {
tmp = x + ((y * -0.3333333333333333) / z);
} else if (y <= 480.0) {
tmp = 0.3333333333333333 * (t / (z * y));
} else {
tmp = x + (y / (z * -3.0));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if y <= -1.06e-34: tmp = x + ((y * -0.3333333333333333) / z) elif y <= 480.0: tmp = 0.3333333333333333 * (t / (z * y)) else: tmp = x + (y / (z * -3.0)) return tmp
function code(x, y, z, t) tmp = 0.0 if (y <= -1.06e-34) tmp = Float64(x + Float64(Float64(y * -0.3333333333333333) / z)); elseif (y <= 480.0) tmp = Float64(0.3333333333333333 * Float64(t / Float64(z * y))); else tmp = Float64(x + Float64(y / Float64(z * -3.0))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (y <= -1.06e-34) tmp = x + ((y * -0.3333333333333333) / z); elseif (y <= 480.0) tmp = 0.3333333333333333 * (t / (z * y)); else tmp = x + (y / (z * -3.0)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[y, -1.06e-34], N[(x + N[(N[(y * -0.3333333333333333), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 480.0], N[(0.3333333333333333 * N[(t / N[(z * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y / N[(z * -3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.06 \cdot 10^{-34}:\\
\;\;\;\;x + \frac{y \cdot -0.3333333333333333}{z}\\
\mathbf{elif}\;y \leq 480:\\
\;\;\;\;0.3333333333333333 \cdot \frac{t}{z \cdot y}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{z \cdot -3}\\
\end{array}
\end{array}
if y < -1.06000000000000006e-34Initial program 97.5%
sub-neg97.5%
associate-+l+97.5%
+-commutative97.5%
remove-double-neg97.5%
distribute-frac-neg97.5%
distribute-neg-in97.5%
remove-double-neg97.5%
sub-neg97.5%
neg-mul-197.5%
times-frac99.7%
distribute-frac-neg99.7%
neg-mul-199.7%
*-commutative99.7%
associate-/l*99.7%
*-commutative99.7%
Simplified99.7%
Taylor expanded in t around 0 89.9%
+-commutative89.9%
Simplified89.9%
associate-*r/90.0%
Applied egg-rr90.0%
if -1.06000000000000006e-34 < y < 480Initial program 85.8%
sub-neg85.8%
associate-+l+85.8%
+-commutative85.8%
remove-double-neg85.8%
distribute-frac-neg85.8%
distribute-neg-in85.8%
remove-double-neg85.8%
sub-neg85.8%
neg-mul-185.8%
times-frac91.4%
distribute-frac-neg91.4%
neg-mul-191.4%
*-commutative91.4%
associate-/l*91.4%
*-commutative91.4%
Simplified91.5%
Taylor expanded in z around 0 91.4%
Taylor expanded in x around 0 73.6%
Taylor expanded in t around inf 66.3%
if 480 < y Initial program 99.8%
+-commutative99.8%
associate-+r-99.8%
sub-neg99.8%
associate-*l*99.8%
*-commutative99.8%
distribute-frac-neg299.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
associate-/r*90.8%
div-inv90.7%
Applied egg-rr90.7%
Taylor expanded in t around 0 92.9%
Final simplification79.2%
(FPCore (x y z t) :precision binary64 (if (<= z -4.4e+183) x (if (<= z 2.2e-8) (* -0.3333333333333333 (/ y z)) x)))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -4.4e+183) {
tmp = x;
} else if (z <= 2.2e-8) {
tmp = -0.3333333333333333 * (y / z);
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (z <= (-4.4d+183)) then
tmp = x
else if (z <= 2.2d-8) then
tmp = (-0.3333333333333333d0) * (y / z)
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -4.4e+183) {
tmp = x;
} else if (z <= 2.2e-8) {
tmp = -0.3333333333333333 * (y / z);
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -4.4e+183: tmp = x elif z <= 2.2e-8: tmp = -0.3333333333333333 * (y / z) else: tmp = x return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -4.4e+183) tmp = x; elseif (z <= 2.2e-8) tmp = Float64(-0.3333333333333333 * Float64(y / z)); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -4.4e+183) tmp = x; elseif (z <= 2.2e-8) tmp = -0.3333333333333333 * (y / z); else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -4.4e+183], x, If[LessEqual[z, 2.2e-8], N[(-0.3333333333333333 * N[(y / z), $MachinePrecision]), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4.4 \cdot 10^{+183}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 2.2 \cdot 10^{-8}:\\
\;\;\;\;-0.3333333333333333 \cdot \frac{y}{z}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if z < -4.39999999999999981e183 or 2.1999999999999998e-8 < z Initial program 97.7%
sub-neg97.7%
associate-+l+97.7%
+-commutative97.7%
remove-double-neg97.7%
distribute-frac-neg97.7%
distribute-neg-in97.7%
remove-double-neg97.7%
sub-neg97.7%
neg-mul-197.7%
times-frac91.1%
distribute-frac-neg91.1%
neg-mul-191.1%
*-commutative91.1%
associate-/l*91.1%
*-commutative91.1%
Simplified91.1%
Taylor expanded in x around inf 51.6%
if -4.39999999999999981e183 < z < 2.1999999999999998e-8Initial program 89.6%
sub-neg89.6%
associate-+l+89.6%
+-commutative89.6%
remove-double-neg89.6%
distribute-frac-neg89.6%
distribute-neg-in89.6%
remove-double-neg89.6%
sub-neg89.6%
neg-mul-189.6%
times-frac98.1%
distribute-frac-neg98.1%
neg-mul-198.1%
*-commutative98.1%
associate-/l*98.1%
*-commutative98.1%
Simplified98.1%
Taylor expanded in t around 0 53.2%
+-commutative53.2%
Simplified53.2%
Taylor expanded in y around inf 45.1%
Final simplification47.3%
(FPCore (x y z t) :precision binary64 (if (<= z -4.4e+183) x (if (<= z 9.5e-10) (/ y (* z -3.0)) x)))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -4.4e+183) {
tmp = x;
} else if (z <= 9.5e-10) {
tmp = y / (z * -3.0);
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (z <= (-4.4d+183)) then
tmp = x
else if (z <= 9.5d-10) then
tmp = y / (z * (-3.0d0))
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -4.4e+183) {
tmp = x;
} else if (z <= 9.5e-10) {
tmp = y / (z * -3.0);
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -4.4e+183: tmp = x elif z <= 9.5e-10: tmp = y / (z * -3.0) else: tmp = x return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -4.4e+183) tmp = x; elseif (z <= 9.5e-10) tmp = Float64(y / Float64(z * -3.0)); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -4.4e+183) tmp = x; elseif (z <= 9.5e-10) tmp = y / (z * -3.0); else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -4.4e+183], x, If[LessEqual[z, 9.5e-10], N[(y / N[(z * -3.0), $MachinePrecision]), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4.4 \cdot 10^{+183}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 9.5 \cdot 10^{-10}:\\
\;\;\;\;\frac{y}{z \cdot -3}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if z < -4.39999999999999981e183 or 9.50000000000000028e-10 < z Initial program 97.7%
sub-neg97.7%
associate-+l+97.7%
+-commutative97.7%
remove-double-neg97.7%
distribute-frac-neg97.7%
distribute-neg-in97.7%
remove-double-neg97.7%
sub-neg97.7%
neg-mul-197.7%
times-frac91.1%
distribute-frac-neg91.1%
neg-mul-191.1%
*-commutative91.1%
associate-/l*91.1%
*-commutative91.1%
Simplified91.1%
Taylor expanded in x around inf 51.6%
if -4.39999999999999981e183 < z < 9.50000000000000028e-10Initial program 89.6%
sub-neg89.6%
associate-+l+89.6%
+-commutative89.6%
remove-double-neg89.6%
distribute-frac-neg89.6%
distribute-neg-in89.6%
remove-double-neg89.6%
sub-neg89.6%
neg-mul-189.6%
times-frac98.1%
distribute-frac-neg98.1%
neg-mul-198.1%
*-commutative98.1%
associate-/l*98.1%
*-commutative98.1%
Simplified98.1%
Taylor expanded in t around 0 53.2%
+-commutative53.2%
Simplified53.2%
Taylor expanded in y around inf 45.1%
associate-*r/45.1%
*-commutative45.1%
Simplified45.1%
associate-/l*45.1%
*-commutative45.1%
Applied egg-rr45.1%
Taylor expanded in z around 0 45.1%
metadata-eval45.1%
times-frac45.2%
neg-mul-145.2%
*-commutative45.2%
distribute-frac-neg45.2%
distribute-frac-neg245.2%
distribute-rgt-neg-in45.2%
metadata-eval45.2%
Simplified45.2%
Final simplification47.3%
(FPCore (x y z t) :precision binary64 (if (<= z -4.4e+183) x (if (<= z 2.25e-8) (/ (/ y -3.0) z) x)))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -4.4e+183) {
tmp = x;
} else if (z <= 2.25e-8) {
tmp = (y / -3.0) / z;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (z <= (-4.4d+183)) then
tmp = x
else if (z <= 2.25d-8) then
tmp = (y / (-3.0d0)) / z
else
tmp = x
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -4.4e+183) {
tmp = x;
} else if (z <= 2.25e-8) {
tmp = (y / -3.0) / z;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -4.4e+183: tmp = x elif z <= 2.25e-8: tmp = (y / -3.0) / z else: tmp = x return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -4.4e+183) tmp = x; elseif (z <= 2.25e-8) tmp = Float64(Float64(y / -3.0) / z); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -4.4e+183) tmp = x; elseif (z <= 2.25e-8) tmp = (y / -3.0) / z; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -4.4e+183], x, If[LessEqual[z, 2.25e-8], N[(N[(y / -3.0), $MachinePrecision] / z), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4.4 \cdot 10^{+183}:\\
\;\;\;\;x\\
\mathbf{elif}\;z \leq 2.25 \cdot 10^{-8}:\\
\;\;\;\;\frac{\frac{y}{-3}}{z}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if z < -4.39999999999999981e183 or 2.24999999999999996e-8 < z Initial program 97.7%
sub-neg97.7%
associate-+l+97.7%
+-commutative97.7%
remove-double-neg97.7%
distribute-frac-neg97.7%
distribute-neg-in97.7%
remove-double-neg97.7%
sub-neg97.7%
neg-mul-197.7%
times-frac91.1%
distribute-frac-neg91.1%
neg-mul-191.1%
*-commutative91.1%
associate-/l*91.1%
*-commutative91.1%
Simplified91.1%
Taylor expanded in x around inf 51.6%
if -4.39999999999999981e183 < z < 2.24999999999999996e-8Initial program 89.6%
sub-neg89.6%
associate-+l+89.6%
+-commutative89.6%
remove-double-neg89.6%
distribute-frac-neg89.6%
distribute-neg-in89.6%
remove-double-neg89.6%
sub-neg89.6%
neg-mul-189.6%
times-frac98.1%
distribute-frac-neg98.1%
neg-mul-198.1%
*-commutative98.1%
associate-/l*98.1%
*-commutative98.1%
Simplified98.1%
Taylor expanded in t around 0 53.2%
+-commutative53.2%
Simplified53.2%
Taylor expanded in y around inf 45.1%
associate-*r/45.1%
*-commutative45.1%
Simplified45.1%
associate-/l*45.1%
*-commutative45.1%
Applied egg-rr45.1%
associate-*l/45.1%
associate-/l*45.1%
metadata-eval45.1%
times-frac45.2%
*-un-lft-identity45.2%
associate-/r*45.2%
Applied egg-rr45.2%
Final simplification47.3%
(FPCore (x y z t) :precision binary64 x)
double code(double x, double y, double z, double t) {
return x;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = x
end function
public static double code(double x, double y, double z, double t) {
return x;
}
def code(x, y, z, t): return x
function code(x, y, z, t) return x end
function tmp = code(x, y, z, t) tmp = x; end
code[x_, y_, z_, t_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 92.3%
sub-neg92.3%
associate-+l+92.3%
+-commutative92.3%
remove-double-neg92.3%
distribute-frac-neg92.3%
distribute-neg-in92.3%
remove-double-neg92.3%
sub-neg92.3%
neg-mul-192.3%
times-frac95.7%
distribute-frac-neg95.7%
neg-mul-195.7%
*-commutative95.7%
associate-/l*95.7%
*-commutative95.7%
Simplified95.7%
Taylor expanded in x around inf 24.2%
Final simplification24.2%
(FPCore (x y z t) :precision binary64 (+ (- x (/ y (* z 3.0))) (/ (/ t (* z 3.0)) y)))
double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + ((t / (z * 3.0)) / y);
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x - (y / (z * 3.0d0))) + ((t / (z * 3.0d0)) / y)
end function
public static double code(double x, double y, double z, double t) {
return (x - (y / (z * 3.0))) + ((t / (z * 3.0)) / y);
}
def code(x, y, z, t): return (x - (y / (z * 3.0))) + ((t / (z * 3.0)) / y)
function code(x, y, z, t) return Float64(Float64(x - Float64(y / Float64(z * 3.0))) + Float64(Float64(t / Float64(z * 3.0)) / y)) end
function tmp = code(x, y, z, t) tmp = (x - (y / (z * 3.0))) + ((t / (z * 3.0)) / y); end
code[x_, y_, z_, t_] := N[(N[(x - N[(y / N[(z * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(t / N[(z * 3.0), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x - \frac{y}{z \cdot 3}\right) + \frac{\frac{t}{z \cdot 3}}{y}
\end{array}
herbie shell --seed 2024041
(FPCore (x y z t)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, H"
:precision binary64
:herbie-target
(+ (- x (/ y (* z 3.0))) (/ (/ t (* z 3.0)) y))
(+ (- x (/ y (* z 3.0))) (/ t (* (* z 3.0) y))))