
(FPCore (x y z t a b c) :precision binary64 (/ (+ (- (* (* x 9.0) y) (* (* (* z 4.0) t) a)) b) (* z c)))
double code(double x, double y, double z, double t, double a, double b, double c) {
return ((((x * 9.0) * y) - (((z * 4.0) * t) * a)) + b) / (z * c);
}
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
code = ((((x * 9.0d0) * y) - (((z * 4.0d0) * t) * a)) + b) / (z * c)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return ((((x * 9.0) * y) - (((z * 4.0) * t) * a)) + b) / (z * c);
}
def code(x, y, z, t, a, b, c): return ((((x * 9.0) * y) - (((z * 4.0) * t) * a)) + b) / (z * c)
function code(x, y, z, t, a, b, c) return Float64(Float64(Float64(Float64(Float64(x * 9.0) * y) - Float64(Float64(Float64(z * 4.0) * t) * a)) + b) / Float64(z * c)) end
function tmp = code(x, y, z, t, a, b, c) tmp = ((((x * 9.0) * y) - (((z * 4.0) * t) * a)) + b) / (z * c); end
code[x_, y_, z_, t_, a_, b_, c_] := N[(N[(N[(N[(N[(x * 9.0), $MachinePrecision] * y), $MachinePrecision] - N[(N[(N[(z * 4.0), $MachinePrecision] * t), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision] / N[(z * c), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(\left(x \cdot 9\right) \cdot y - \left(\left(z \cdot 4\right) \cdot t\right) \cdot a\right) + b}{z \cdot c}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 19 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b c) :precision binary64 (/ (+ (- (* (* x 9.0) y) (* (* (* z 4.0) t) a)) b) (* z c)))
double code(double x, double y, double z, double t, double a, double b, double c) {
return ((((x * 9.0) * y) - (((z * 4.0) * t) * a)) + b) / (z * c);
}
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
code = ((((x * 9.0d0) * y) - (((z * 4.0d0) * t) * a)) + b) / (z * c)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return ((((x * 9.0) * y) - (((z * 4.0) * t) * a)) + b) / (z * c);
}
def code(x, y, z, t, a, b, c): return ((((x * 9.0) * y) - (((z * 4.0) * t) * a)) + b) / (z * c)
function code(x, y, z, t, a, b, c) return Float64(Float64(Float64(Float64(Float64(x * 9.0) * y) - Float64(Float64(Float64(z * 4.0) * t) * a)) + b) / Float64(z * c)) end
function tmp = code(x, y, z, t, a, b, c) tmp = ((((x * 9.0) * y) - (((z * 4.0) * t) * a)) + b) / (z * c); end
code[x_, y_, z_, t_, a_, b_, c_] := N[(N[(N[(N[(N[(x * 9.0), $MachinePrecision] * y), $MachinePrecision] - N[(N[(N[(z * 4.0), $MachinePrecision] * t), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision] / N[(z * c), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(\left(x \cdot 9\right) \cdot y - \left(\left(z \cdot 4\right) \cdot t\right) \cdot a\right) + b}{z \cdot c}
\end{array}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c)
:precision binary64
(if (<= z -1.25e-82)
(/ (+ (/ (+ b (* x (* 9.0 y))) z) (* t (* a -4.0))) c)
(if (<= z 4.4e+24)
(/ (+ b (- (* y (* x 9.0)) (* t (* a (* z 4.0))))) (* z c))
(/ (* y (- (/ (- (/ b z) (* 4.0 (* t a))) y) (/ (* x -9.0) z))) c))))assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (z <= -1.25e-82) {
tmp = (((b + (x * (9.0 * y))) / z) + (t * (a * -4.0))) / c;
} else if (z <= 4.4e+24) {
tmp = (b + ((y * (x * 9.0)) - (t * (a * (z * 4.0))))) / (z * c);
} else {
tmp = (y * ((((b / z) - (4.0 * (t * a))) / y) - ((x * -9.0) / z))) / c;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if (z <= (-1.25d-82)) then
tmp = (((b + (x * (9.0d0 * y))) / z) + (t * (a * (-4.0d0)))) / c
else if (z <= 4.4d+24) then
tmp = (b + ((y * (x * 9.0d0)) - (t * (a * (z * 4.0d0))))) / (z * c)
else
tmp = (y * ((((b / z) - (4.0d0 * (t * a))) / y) - ((x * (-9.0d0)) / z))) / c
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (z <= -1.25e-82) {
tmp = (((b + (x * (9.0 * y))) / z) + (t * (a * -4.0))) / c;
} else if (z <= 4.4e+24) {
tmp = (b + ((y * (x * 9.0)) - (t * (a * (z * 4.0))))) / (z * c);
} else {
tmp = (y * ((((b / z) - (4.0 * (t * a))) / y) - ((x * -9.0) / z))) / c;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): tmp = 0 if z <= -1.25e-82: tmp = (((b + (x * (9.0 * y))) / z) + (t * (a * -4.0))) / c elif z <= 4.4e+24: tmp = (b + ((y * (x * 9.0)) - (t * (a * (z * 4.0))))) / (z * c) else: tmp = (y * ((((b / z) - (4.0 * (t * a))) / y) - ((x * -9.0) / z))) / c return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) tmp = 0.0 if (z <= -1.25e-82) tmp = Float64(Float64(Float64(Float64(b + Float64(x * Float64(9.0 * y))) / z) + Float64(t * Float64(a * -4.0))) / c); elseif (z <= 4.4e+24) tmp = Float64(Float64(b + Float64(Float64(y * Float64(x * 9.0)) - Float64(t * Float64(a * Float64(z * 4.0))))) / Float64(z * c)); else tmp = Float64(Float64(y * Float64(Float64(Float64(Float64(b / z) - Float64(4.0 * Float64(t * a))) / y) - Float64(Float64(x * -9.0) / z))) / c); end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
tmp = 0.0;
if (z <= -1.25e-82)
tmp = (((b + (x * (9.0 * y))) / z) + (t * (a * -4.0))) / c;
elseif (z <= 4.4e+24)
tmp = (b + ((y * (x * 9.0)) - (t * (a * (z * 4.0))))) / (z * c);
else
tmp = (y * ((((b / z) - (4.0 * (t * a))) / y) - ((x * -9.0) / z))) / c;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[z, -1.25e-82], N[(N[(N[(N[(b + N[(x * N[(9.0 * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision] + N[(t * N[(a * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], If[LessEqual[z, 4.4e+24], N[(N[(b + N[(N[(y * N[(x * 9.0), $MachinePrecision]), $MachinePrecision] - N[(t * N[(a * N[(z * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(z * c), $MachinePrecision]), $MachinePrecision], N[(N[(y * N[(N[(N[(N[(b / z), $MachinePrecision] - N[(4.0 * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision] - N[(N[(x * -9.0), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.25 \cdot 10^{-82}:\\
\;\;\;\;\frac{\frac{b + x \cdot \left(9 \cdot y\right)}{z} + t \cdot \left(a \cdot -4\right)}{c}\\
\mathbf{elif}\;z \leq 4.4 \cdot 10^{+24}:\\
\;\;\;\;\frac{b + \left(y \cdot \left(x \cdot 9\right) - t \cdot \left(a \cdot \left(z \cdot 4\right)\right)\right)}{z \cdot c}\\
\mathbf{else}:\\
\;\;\;\;\frac{y \cdot \left(\frac{\frac{b}{z} - 4 \cdot \left(t \cdot a\right)}{y} - \frac{x \cdot -9}{z}\right)}{c}\\
\end{array}
\end{array}
if z < -1.25e-82Initial program 73.1%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified90.3%
if -1.25e-82 < z < 4.40000000000000003e24Initial program 94.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6494.9%
Applied egg-rr94.9%
if 4.40000000000000003e24 < z Initial program 68.2%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified91.5%
Taylor expanded in y around -inf
mul-1-negN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
neg-sub0N/A
--lowering--.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r/N/A
mul-1-negN/A
distribute-neg-inN/A
distribute-lft-neg-inN/A
metadata-evalN/A
sub-negN/A
/-lowering-/.f64N/A
Simplified94.1%
Final simplification93.3%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (* t (* a -4.0))))
(if (<= z -9e-83)
(/ (+ (/ (+ b (* x (* 9.0 y))) z) t_1) c)
(if (<= z 1.9e+20)
(/ (+ b (- (* y (* x 9.0)) (* t (* a (* z 4.0))))) (* z c))
(/ (+ t_1 (* y (+ (/ (* x 9.0) z) (/ b (* z y))))) c)))))assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = t * (a * -4.0);
double tmp;
if (z <= -9e-83) {
tmp = (((b + (x * (9.0 * y))) / z) + t_1) / c;
} else if (z <= 1.9e+20) {
tmp = (b + ((y * (x * 9.0)) - (t * (a * (z * 4.0))))) / (z * c);
} else {
tmp = (t_1 + (y * (((x * 9.0) / z) + (b / (z * y))))) / c;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: tmp
t_1 = t * (a * (-4.0d0))
if (z <= (-9d-83)) then
tmp = (((b + (x * (9.0d0 * y))) / z) + t_1) / c
else if (z <= 1.9d+20) then
tmp = (b + ((y * (x * 9.0d0)) - (t * (a * (z * 4.0d0))))) / (z * c)
else
tmp = (t_1 + (y * (((x * 9.0d0) / z) + (b / (z * y))))) / c
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = t * (a * -4.0);
double tmp;
if (z <= -9e-83) {
tmp = (((b + (x * (9.0 * y))) / z) + t_1) / c;
} else if (z <= 1.9e+20) {
tmp = (b + ((y * (x * 9.0)) - (t * (a * (z * 4.0))))) / (z * c);
} else {
tmp = (t_1 + (y * (((x * 9.0) / z) + (b / (z * y))))) / c;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): t_1 = t * (a * -4.0) tmp = 0 if z <= -9e-83: tmp = (((b + (x * (9.0 * y))) / z) + t_1) / c elif z <= 1.9e+20: tmp = (b + ((y * (x * 9.0)) - (t * (a * (z * 4.0))))) / (z * c) else: tmp = (t_1 + (y * (((x * 9.0) / z) + (b / (z * y))))) / c return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) t_1 = Float64(t * Float64(a * -4.0)) tmp = 0.0 if (z <= -9e-83) tmp = Float64(Float64(Float64(Float64(b + Float64(x * Float64(9.0 * y))) / z) + t_1) / c); elseif (z <= 1.9e+20) tmp = Float64(Float64(b + Float64(Float64(y * Float64(x * 9.0)) - Float64(t * Float64(a * Float64(z * 4.0))))) / Float64(z * c)); else tmp = Float64(Float64(t_1 + Float64(y * Float64(Float64(Float64(x * 9.0) / z) + Float64(b / Float64(z * y))))) / c); end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
t_1 = t * (a * -4.0);
tmp = 0.0;
if (z <= -9e-83)
tmp = (((b + (x * (9.0 * y))) / z) + t_1) / c;
elseif (z <= 1.9e+20)
tmp = (b + ((y * (x * 9.0)) - (t * (a * (z * 4.0))))) / (z * c);
else
tmp = (t_1 + (y * (((x * 9.0) / z) + (b / (z * y))))) / c;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(t * N[(a * -4.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -9e-83], N[(N[(N[(N[(b + N[(x * N[(9.0 * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision] + t$95$1), $MachinePrecision] / c), $MachinePrecision], If[LessEqual[z, 1.9e+20], N[(N[(b + N[(N[(y * N[(x * 9.0), $MachinePrecision]), $MachinePrecision] - N[(t * N[(a * N[(z * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(z * c), $MachinePrecision]), $MachinePrecision], N[(N[(t$95$1 + N[(y * N[(N[(N[(x * 9.0), $MachinePrecision] / z), $MachinePrecision] + N[(b / N[(z * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision]]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
t_1 := t \cdot \left(a \cdot -4\right)\\
\mathbf{if}\;z \leq -9 \cdot 10^{-83}:\\
\;\;\;\;\frac{\frac{b + x \cdot \left(9 \cdot y\right)}{z} + t\_1}{c}\\
\mathbf{elif}\;z \leq 1.9 \cdot 10^{+20}:\\
\;\;\;\;\frac{b + \left(y \cdot \left(x \cdot 9\right) - t \cdot \left(a \cdot \left(z \cdot 4\right)\right)\right)}{z \cdot c}\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_1 + y \cdot \left(\frac{x \cdot 9}{z} + \frac{b}{z \cdot y}\right)}{c}\\
\end{array}
\end{array}
if z < -8.99999999999999995e-83Initial program 73.1%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified90.3%
if -8.99999999999999995e-83 < z < 1.9e20Initial program 94.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6494.9%
Applied egg-rr94.9%
if 1.9e20 < z Initial program 68.2%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified91.5%
Taylor expanded in y around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6495.7%
Simplified95.7%
Final simplification93.6%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (* t (* a -4.0))))
(if (<= b -1.65e+41)
(/ (+ t_1 (/ b z)) c)
(if (<= b -8.2e-77)
(/ (+ (* x (* 9.0 y)) (* a (* -4.0 (* z t)))) (* z c))
(if (<= b 1e-76)
(/ (+ t_1 (* 9.0 (/ (* x y) z))) c)
(/ (/ (+ b (* y (* x 9.0))) z) c))))))assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = t * (a * -4.0);
double tmp;
if (b <= -1.65e+41) {
tmp = (t_1 + (b / z)) / c;
} else if (b <= -8.2e-77) {
tmp = ((x * (9.0 * y)) + (a * (-4.0 * (z * t)))) / (z * c);
} else if (b <= 1e-76) {
tmp = (t_1 + (9.0 * ((x * y) / z))) / c;
} else {
tmp = ((b + (y * (x * 9.0))) / z) / c;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: tmp
t_1 = t * (a * (-4.0d0))
if (b <= (-1.65d+41)) then
tmp = (t_1 + (b / z)) / c
else if (b <= (-8.2d-77)) then
tmp = ((x * (9.0d0 * y)) + (a * ((-4.0d0) * (z * t)))) / (z * c)
else if (b <= 1d-76) then
tmp = (t_1 + (9.0d0 * ((x * y) / z))) / c
else
tmp = ((b + (y * (x * 9.0d0))) / z) / c
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = t * (a * -4.0);
double tmp;
if (b <= -1.65e+41) {
tmp = (t_1 + (b / z)) / c;
} else if (b <= -8.2e-77) {
tmp = ((x * (9.0 * y)) + (a * (-4.0 * (z * t)))) / (z * c);
} else if (b <= 1e-76) {
tmp = (t_1 + (9.0 * ((x * y) / z))) / c;
} else {
tmp = ((b + (y * (x * 9.0))) / z) / c;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): t_1 = t * (a * -4.0) tmp = 0 if b <= -1.65e+41: tmp = (t_1 + (b / z)) / c elif b <= -8.2e-77: tmp = ((x * (9.0 * y)) + (a * (-4.0 * (z * t)))) / (z * c) elif b <= 1e-76: tmp = (t_1 + (9.0 * ((x * y) / z))) / c else: tmp = ((b + (y * (x * 9.0))) / z) / c return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) t_1 = Float64(t * Float64(a * -4.0)) tmp = 0.0 if (b <= -1.65e+41) tmp = Float64(Float64(t_1 + Float64(b / z)) / c); elseif (b <= -8.2e-77) tmp = Float64(Float64(Float64(x * Float64(9.0 * y)) + Float64(a * Float64(-4.0 * Float64(z * t)))) / Float64(z * c)); elseif (b <= 1e-76) tmp = Float64(Float64(t_1 + Float64(9.0 * Float64(Float64(x * y) / z))) / c); else tmp = Float64(Float64(Float64(b + Float64(y * Float64(x * 9.0))) / z) / c); end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
t_1 = t * (a * -4.0);
tmp = 0.0;
if (b <= -1.65e+41)
tmp = (t_1 + (b / z)) / c;
elseif (b <= -8.2e-77)
tmp = ((x * (9.0 * y)) + (a * (-4.0 * (z * t)))) / (z * c);
elseif (b <= 1e-76)
tmp = (t_1 + (9.0 * ((x * y) / z))) / c;
else
tmp = ((b + (y * (x * 9.0))) / z) / c;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(t * N[(a * -4.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1.65e+41], N[(N[(t$95$1 + N[(b / z), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], If[LessEqual[b, -8.2e-77], N[(N[(N[(x * N[(9.0 * y), $MachinePrecision]), $MachinePrecision] + N[(a * N[(-4.0 * N[(z * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(z * c), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1e-76], N[(N[(t$95$1 + N[(9.0 * N[(N[(x * y), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], N[(N[(N[(b + N[(y * N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision] / c), $MachinePrecision]]]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
t_1 := t \cdot \left(a \cdot -4\right)\\
\mathbf{if}\;b \leq -1.65 \cdot 10^{+41}:\\
\;\;\;\;\frac{t\_1 + \frac{b}{z}}{c}\\
\mathbf{elif}\;b \leq -8.2 \cdot 10^{-77}:\\
\;\;\;\;\frac{x \cdot \left(9 \cdot y\right) + a \cdot \left(-4 \cdot \left(z \cdot t\right)\right)}{z \cdot c}\\
\mathbf{elif}\;b \leq 10^{-76}:\\
\;\;\;\;\frac{t\_1 + 9 \cdot \frac{x \cdot y}{z}}{c}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{b + y \cdot \left(x \cdot 9\right)}{z}}{c}\\
\end{array}
\end{array}
if b < -1.65e41Initial program 85.7%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified81.3%
Taylor expanded in x around 0
/-lowering-/.f6473.9%
Simplified73.9%
if -1.65e41 < b < -8.19999999999999925e-77Initial program 87.9%
Taylor expanded in b around 0
/-lowering-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6471.3%
Simplified71.3%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f6471.4%
Applied egg-rr71.4%
if -8.19999999999999925e-77 < b < 9.99999999999999927e-77Initial program 82.0%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified93.2%
Taylor expanded in x around inf
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6487.7%
Simplified87.7%
if 9.99999999999999927e-77 < b Initial program 77.3%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified86.8%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6486.7%
Applied egg-rr86.7%
Taylor expanded in z around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6473.6%
Simplified73.6%
*-commutativeN/A
div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6473.8%
Applied egg-rr73.8%
Final simplification79.3%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (* t (* a -4.0))))
(if (<= b -1.55e+41)
(/ (+ t_1 (/ b z)) c)
(if (<= b -8.2e-76)
(/ (+ (* a (* -4.0 (* z t))) (* 9.0 (* x y))) (* z c))
(if (<= b 1e-76)
(/ (+ t_1 (* 9.0 (/ (* x y) z))) c)
(/ (/ (+ b (* y (* x 9.0))) z) c))))))assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = t * (a * -4.0);
double tmp;
if (b <= -1.55e+41) {
tmp = (t_1 + (b / z)) / c;
} else if (b <= -8.2e-76) {
tmp = ((a * (-4.0 * (z * t))) + (9.0 * (x * y))) / (z * c);
} else if (b <= 1e-76) {
tmp = (t_1 + (9.0 * ((x * y) / z))) / c;
} else {
tmp = ((b + (y * (x * 9.0))) / z) / c;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: tmp
t_1 = t * (a * (-4.0d0))
if (b <= (-1.55d+41)) then
tmp = (t_1 + (b / z)) / c
else if (b <= (-8.2d-76)) then
tmp = ((a * ((-4.0d0) * (z * t))) + (9.0d0 * (x * y))) / (z * c)
else if (b <= 1d-76) then
tmp = (t_1 + (9.0d0 * ((x * y) / z))) / c
else
tmp = ((b + (y * (x * 9.0d0))) / z) / c
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = t * (a * -4.0);
double tmp;
if (b <= -1.55e+41) {
tmp = (t_1 + (b / z)) / c;
} else if (b <= -8.2e-76) {
tmp = ((a * (-4.0 * (z * t))) + (9.0 * (x * y))) / (z * c);
} else if (b <= 1e-76) {
tmp = (t_1 + (9.0 * ((x * y) / z))) / c;
} else {
tmp = ((b + (y * (x * 9.0))) / z) / c;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): t_1 = t * (a * -4.0) tmp = 0 if b <= -1.55e+41: tmp = (t_1 + (b / z)) / c elif b <= -8.2e-76: tmp = ((a * (-4.0 * (z * t))) + (9.0 * (x * y))) / (z * c) elif b <= 1e-76: tmp = (t_1 + (9.0 * ((x * y) / z))) / c else: tmp = ((b + (y * (x * 9.0))) / z) / c return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) t_1 = Float64(t * Float64(a * -4.0)) tmp = 0.0 if (b <= -1.55e+41) tmp = Float64(Float64(t_1 + Float64(b / z)) / c); elseif (b <= -8.2e-76) tmp = Float64(Float64(Float64(a * Float64(-4.0 * Float64(z * t))) + Float64(9.0 * Float64(x * y))) / Float64(z * c)); elseif (b <= 1e-76) tmp = Float64(Float64(t_1 + Float64(9.0 * Float64(Float64(x * y) / z))) / c); else tmp = Float64(Float64(Float64(b + Float64(y * Float64(x * 9.0))) / z) / c); end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
t_1 = t * (a * -4.0);
tmp = 0.0;
if (b <= -1.55e+41)
tmp = (t_1 + (b / z)) / c;
elseif (b <= -8.2e-76)
tmp = ((a * (-4.0 * (z * t))) + (9.0 * (x * y))) / (z * c);
elseif (b <= 1e-76)
tmp = (t_1 + (9.0 * ((x * y) / z))) / c;
else
tmp = ((b + (y * (x * 9.0))) / z) / c;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(t * N[(a * -4.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1.55e+41], N[(N[(t$95$1 + N[(b / z), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], If[LessEqual[b, -8.2e-76], N[(N[(N[(a * N[(-4.0 * N[(z * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(9.0 * N[(x * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(z * c), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1e-76], N[(N[(t$95$1 + N[(9.0 * N[(N[(x * y), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], N[(N[(N[(b + N[(y * N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision] / c), $MachinePrecision]]]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
t_1 := t \cdot \left(a \cdot -4\right)\\
\mathbf{if}\;b \leq -1.55 \cdot 10^{+41}:\\
\;\;\;\;\frac{t\_1 + \frac{b}{z}}{c}\\
\mathbf{elif}\;b \leq -8.2 \cdot 10^{-76}:\\
\;\;\;\;\frac{a \cdot \left(-4 \cdot \left(z \cdot t\right)\right) + 9 \cdot \left(x \cdot y\right)}{z \cdot c}\\
\mathbf{elif}\;b \leq 10^{-76}:\\
\;\;\;\;\frac{t\_1 + 9 \cdot \frac{x \cdot y}{z}}{c}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{b + y \cdot \left(x \cdot 9\right)}{z}}{c}\\
\end{array}
\end{array}
if b < -1.55e41Initial program 85.7%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified81.3%
Taylor expanded in x around 0
/-lowering-/.f6473.9%
Simplified73.9%
if -1.55e41 < b < -8.1999999999999996e-76Initial program 87.9%
Taylor expanded in b around 0
/-lowering-/.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6471.3%
Simplified71.3%
if -8.1999999999999996e-76 < b < 9.99999999999999927e-77Initial program 82.0%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified93.2%
Taylor expanded in x around inf
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6487.7%
Simplified87.7%
if 9.99999999999999927e-77 < b Initial program 77.3%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified86.8%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6486.7%
Applied egg-rr86.7%
Taylor expanded in z around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6473.6%
Simplified73.6%
*-commutativeN/A
div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6473.8%
Applied egg-rr73.8%
Final simplification79.2%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (/ (+ (/ (+ b (* x (* 9.0 y))) z) (* t (* a -4.0))) c)))
(if (<= z -1.25e-82)
t_1
(if (<= z 1e-8)
(/ (+ b (- (* y (* x 9.0)) (* t (* a (* z 4.0))))) (* z c))
t_1))))assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (((b + (x * (9.0 * y))) / z) + (t * (a * -4.0))) / c;
double tmp;
if (z <= -1.25e-82) {
tmp = t_1;
} else if (z <= 1e-8) {
tmp = (b + ((y * (x * 9.0)) - (t * (a * (z * 4.0))))) / (z * c);
} else {
tmp = t_1;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: tmp
t_1 = (((b + (x * (9.0d0 * y))) / z) + (t * (a * (-4.0d0)))) / c
if (z <= (-1.25d-82)) then
tmp = t_1
else if (z <= 1d-8) then
tmp = (b + ((y * (x * 9.0d0)) - (t * (a * (z * 4.0d0))))) / (z * c)
else
tmp = t_1
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (((b + (x * (9.0 * y))) / z) + (t * (a * -4.0))) / c;
double tmp;
if (z <= -1.25e-82) {
tmp = t_1;
} else if (z <= 1e-8) {
tmp = (b + ((y * (x * 9.0)) - (t * (a * (z * 4.0))))) / (z * c);
} else {
tmp = t_1;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): t_1 = (((b + (x * (9.0 * y))) / z) + (t * (a * -4.0))) / c tmp = 0 if z <= -1.25e-82: tmp = t_1 elif z <= 1e-8: tmp = (b + ((y * (x * 9.0)) - (t * (a * (z * 4.0))))) / (z * c) else: tmp = t_1 return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(Float64(Float64(b + Float64(x * Float64(9.0 * y))) / z) + Float64(t * Float64(a * -4.0))) / c) tmp = 0.0 if (z <= -1.25e-82) tmp = t_1; elseif (z <= 1e-8) tmp = Float64(Float64(b + Float64(Float64(y * Float64(x * 9.0)) - Float64(t * Float64(a * Float64(z * 4.0))))) / Float64(z * c)); else tmp = t_1; end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
t_1 = (((b + (x * (9.0 * y))) / z) + (t * (a * -4.0))) / c;
tmp = 0.0;
if (z <= -1.25e-82)
tmp = t_1;
elseif (z <= 1e-8)
tmp = (b + ((y * (x * 9.0)) - (t * (a * (z * 4.0))))) / (z * c);
else
tmp = t_1;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(N[(N[(b + N[(x * N[(9.0 * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision] + N[(t * N[(a * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision]}, If[LessEqual[z, -1.25e-82], t$95$1, If[LessEqual[z, 1e-8], N[(N[(b + N[(N[(y * N[(x * 9.0), $MachinePrecision]), $MachinePrecision] - N[(t * N[(a * N[(z * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(z * c), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
t_1 := \frac{\frac{b + x \cdot \left(9 \cdot y\right)}{z} + t \cdot \left(a \cdot -4\right)}{c}\\
\mathbf{if}\;z \leq -1.25 \cdot 10^{-82}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 10^{-8}:\\
\;\;\;\;\frac{b + \left(y \cdot \left(x \cdot 9\right) - t \cdot \left(a \cdot \left(z \cdot 4\right)\right)\right)}{z \cdot c}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.25e-82 or 1e-8 < z Initial program 71.9%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified90.6%
if -1.25e-82 < z < 1e-8Initial program 95.4%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6495.4%
Applied egg-rr95.4%
Final simplification92.7%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (* x (* 9.0 y))) (t_2 (/ (+ (/ (+ b t_1) z) (* t (* a -4.0))) c)))
(if (<= z -1.25e-82)
t_2
(if (<= z 6.8e-93) (/ (+ b (- t_1 (* a (* t (* z 4.0))))) (* z c)) t_2))))assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = x * (9.0 * y);
double t_2 = (((b + t_1) / z) + (t * (a * -4.0))) / c;
double tmp;
if (z <= -1.25e-82) {
tmp = t_2;
} else if (z <= 6.8e-93) {
tmp = (b + (t_1 - (a * (t * (z * 4.0))))) / (z * c);
} else {
tmp = t_2;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = x * (9.0d0 * y)
t_2 = (((b + t_1) / z) + (t * (a * (-4.0d0)))) / c
if (z <= (-1.25d-82)) then
tmp = t_2
else if (z <= 6.8d-93) then
tmp = (b + (t_1 - (a * (t * (z * 4.0d0))))) / (z * c)
else
tmp = t_2
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = x * (9.0 * y);
double t_2 = (((b + t_1) / z) + (t * (a * -4.0))) / c;
double tmp;
if (z <= -1.25e-82) {
tmp = t_2;
} else if (z <= 6.8e-93) {
tmp = (b + (t_1 - (a * (t * (z * 4.0))))) / (z * c);
} else {
tmp = t_2;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): t_1 = x * (9.0 * y) t_2 = (((b + t_1) / z) + (t * (a * -4.0))) / c tmp = 0 if z <= -1.25e-82: tmp = t_2 elif z <= 6.8e-93: tmp = (b + (t_1 - (a * (t * (z * 4.0))))) / (z * c) else: tmp = t_2 return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) t_1 = Float64(x * Float64(9.0 * y)) t_2 = Float64(Float64(Float64(Float64(b + t_1) / z) + Float64(t * Float64(a * -4.0))) / c) tmp = 0.0 if (z <= -1.25e-82) tmp = t_2; elseif (z <= 6.8e-93) tmp = Float64(Float64(b + Float64(t_1 - Float64(a * Float64(t * Float64(z * 4.0))))) / Float64(z * c)); else tmp = t_2; end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
t_1 = x * (9.0 * y);
t_2 = (((b + t_1) / z) + (t * (a * -4.0))) / c;
tmp = 0.0;
if (z <= -1.25e-82)
tmp = t_2;
elseif (z <= 6.8e-93)
tmp = (b + (t_1 - (a * (t * (z * 4.0))))) / (z * c);
else
tmp = t_2;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(x * N[(9.0 * y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(N[(b + t$95$1), $MachinePrecision] / z), $MachinePrecision] + N[(t * N[(a * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision]}, If[LessEqual[z, -1.25e-82], t$95$2, If[LessEqual[z, 6.8e-93], N[(N[(b + N[(t$95$1 - N[(a * N[(t * N[(z * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(z * c), $MachinePrecision]), $MachinePrecision], t$95$2]]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
t_1 := x \cdot \left(9 \cdot y\right)\\
t_2 := \frac{\frac{b + t\_1}{z} + t \cdot \left(a \cdot -4\right)}{c}\\
\mathbf{if}\;z \leq -1.25 \cdot 10^{-82}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;z \leq 6.8 \cdot 10^{-93}:\\
\;\;\;\;\frac{b + \left(t\_1 - a \cdot \left(t \cdot \left(z \cdot 4\right)\right)\right)}{z \cdot c}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if z < -1.25e-82 or 6.80000000000000002e-93 < z Initial program 74.2%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified91.0%
if -1.25e-82 < z < 6.80000000000000002e-93Initial program 95.7%
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6495.6%
Applied egg-rr95.6%
Final simplification92.7%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (* t (* a -4.0))))
(if (<= b -2e-53)
(/ (+ t_1 (/ b z)) c)
(if (<= b 8.2e-77)
(/ (+ t_1 (* 9.0 (/ (* x y) z))) c)
(/ (/ (+ b (* y (* x 9.0))) z) c)))))assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = t * (a * -4.0);
double tmp;
if (b <= -2e-53) {
tmp = (t_1 + (b / z)) / c;
} else if (b <= 8.2e-77) {
tmp = (t_1 + (9.0 * ((x * y) / z))) / c;
} else {
tmp = ((b + (y * (x * 9.0))) / z) / c;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: tmp
t_1 = t * (a * (-4.0d0))
if (b <= (-2d-53)) then
tmp = (t_1 + (b / z)) / c
else if (b <= 8.2d-77) then
tmp = (t_1 + (9.0d0 * ((x * y) / z))) / c
else
tmp = ((b + (y * (x * 9.0d0))) / z) / c
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = t * (a * -4.0);
double tmp;
if (b <= -2e-53) {
tmp = (t_1 + (b / z)) / c;
} else if (b <= 8.2e-77) {
tmp = (t_1 + (9.0 * ((x * y) / z))) / c;
} else {
tmp = ((b + (y * (x * 9.0))) / z) / c;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): t_1 = t * (a * -4.0) tmp = 0 if b <= -2e-53: tmp = (t_1 + (b / z)) / c elif b <= 8.2e-77: tmp = (t_1 + (9.0 * ((x * y) / z))) / c else: tmp = ((b + (y * (x * 9.0))) / z) / c return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) t_1 = Float64(t * Float64(a * -4.0)) tmp = 0.0 if (b <= -2e-53) tmp = Float64(Float64(t_1 + Float64(b / z)) / c); elseif (b <= 8.2e-77) tmp = Float64(Float64(t_1 + Float64(9.0 * Float64(Float64(x * y) / z))) / c); else tmp = Float64(Float64(Float64(b + Float64(y * Float64(x * 9.0))) / z) / c); end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
t_1 = t * (a * -4.0);
tmp = 0.0;
if (b <= -2e-53)
tmp = (t_1 + (b / z)) / c;
elseif (b <= 8.2e-77)
tmp = (t_1 + (9.0 * ((x * y) / z))) / c;
else
tmp = ((b + (y * (x * 9.0))) / z) / c;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(t * N[(a * -4.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -2e-53], N[(N[(t$95$1 + N[(b / z), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], If[LessEqual[b, 8.2e-77], N[(N[(t$95$1 + N[(9.0 * N[(N[(x * y), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], N[(N[(N[(b + N[(y * N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision] / c), $MachinePrecision]]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
t_1 := t \cdot \left(a \cdot -4\right)\\
\mathbf{if}\;b \leq -2 \cdot 10^{-53}:\\
\;\;\;\;\frac{t\_1 + \frac{b}{z}}{c}\\
\mathbf{elif}\;b \leq 8.2 \cdot 10^{-77}:\\
\;\;\;\;\frac{t\_1 + 9 \cdot \frac{x \cdot y}{z}}{c}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{b + y \cdot \left(x \cdot 9\right)}{z}}{c}\\
\end{array}
\end{array}
if b < -2.00000000000000006e-53Initial program 85.1%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified81.6%
Taylor expanded in x around 0
/-lowering-/.f6471.3%
Simplified71.3%
if -2.00000000000000006e-53 < b < 8.19999999999999925e-77Initial program 83.1%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified91.1%
Taylor expanded in x around inf
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6485.2%
Simplified85.2%
if 8.19999999999999925e-77 < b Initial program 77.3%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified86.8%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6486.7%
Applied egg-rr86.7%
Taylor expanded in z around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6473.6%
Simplified73.6%
*-commutativeN/A
div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6473.8%
Applied egg-rr73.8%
Final simplification78.0%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c)
:precision binary64
(if (<= b -0.0105)
(/ b (* z c))
(if (<= b -2e-133)
(/ (* y (* x 9.0)) (* z c))
(if (<= b 1e-76) (/ (* a (* t -4.0)) c) (/ (/ b c) z)))))assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -0.0105) {
tmp = b / (z * c);
} else if (b <= -2e-133) {
tmp = (y * (x * 9.0)) / (z * c);
} else if (b <= 1e-76) {
tmp = (a * (t * -4.0)) / c;
} else {
tmp = (b / c) / z;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if (b <= (-0.0105d0)) then
tmp = b / (z * c)
else if (b <= (-2d-133)) then
tmp = (y * (x * 9.0d0)) / (z * c)
else if (b <= 1d-76) then
tmp = (a * (t * (-4.0d0))) / c
else
tmp = (b / c) / z
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -0.0105) {
tmp = b / (z * c);
} else if (b <= -2e-133) {
tmp = (y * (x * 9.0)) / (z * c);
} else if (b <= 1e-76) {
tmp = (a * (t * -4.0)) / c;
} else {
tmp = (b / c) / z;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): tmp = 0 if b <= -0.0105: tmp = b / (z * c) elif b <= -2e-133: tmp = (y * (x * 9.0)) / (z * c) elif b <= 1e-76: tmp = (a * (t * -4.0)) / c else: tmp = (b / c) / z return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) tmp = 0.0 if (b <= -0.0105) tmp = Float64(b / Float64(z * c)); elseif (b <= -2e-133) tmp = Float64(Float64(y * Float64(x * 9.0)) / Float64(z * c)); elseif (b <= 1e-76) tmp = Float64(Float64(a * Float64(t * -4.0)) / c); else tmp = Float64(Float64(b / c) / z); end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
tmp = 0.0;
if (b <= -0.0105)
tmp = b / (z * c);
elseif (b <= -2e-133)
tmp = (y * (x * 9.0)) / (z * c);
elseif (b <= 1e-76)
tmp = (a * (t * -4.0)) / c;
else
tmp = (b / c) / z;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[b, -0.0105], N[(b / N[(z * c), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -2e-133], N[(N[(y * N[(x * 9.0), $MachinePrecision]), $MachinePrecision] / N[(z * c), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1e-76], N[(N[(a * N[(t * -4.0), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], N[(N[(b / c), $MachinePrecision] / z), $MachinePrecision]]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
\mathbf{if}\;b \leq -0.0105:\\
\;\;\;\;\frac{b}{z \cdot c}\\
\mathbf{elif}\;b \leq -2 \cdot 10^{-133}:\\
\;\;\;\;\frac{y \cdot \left(x \cdot 9\right)}{z \cdot c}\\
\mathbf{elif}\;b \leq 10^{-76}:\\
\;\;\;\;\frac{a \cdot \left(t \cdot -4\right)}{c}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{b}{c}}{z}\\
\end{array}
\end{array}
if b < -0.0105000000000000007Initial program 86.0%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6461.0%
Simplified61.0%
if -0.0105000000000000007 < b < -2.0000000000000001e-133Initial program 90.4%
div-invN/A
associate-*r/N/A
*-commutativeN/A
times-fracN/A
*-lowering-*.f64N/A
Applied egg-rr84.3%
Taylor expanded in x around inf
*-lowering-*.f64N/A
*-lowering-*.f6452.2%
Simplified52.2%
un-div-invN/A
associate-/r*N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6458.3%
Applied egg-rr58.3%
if -2.0000000000000001e-133 < b < 9.99999999999999927e-77Initial program 80.2%
Taylor expanded in z around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6455.9%
Simplified55.9%
if 9.99999999999999927e-77 < b Initial program 77.3%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6457.7%
Simplified57.7%
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6462.8%
Applied egg-rr62.8%
Final simplification59.4%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c)
:precision binary64
(if (<= b -0.0077)
(/ b (* z c))
(if (<= b -1.05e-133)
(/ (* 9.0 (* x y)) (* z c))
(if (<= b 8.2e-77) (/ (* a (* t -4.0)) c) (/ (/ b c) z)))))assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -0.0077) {
tmp = b / (z * c);
} else if (b <= -1.05e-133) {
tmp = (9.0 * (x * y)) / (z * c);
} else if (b <= 8.2e-77) {
tmp = (a * (t * -4.0)) / c;
} else {
tmp = (b / c) / z;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if (b <= (-0.0077d0)) then
tmp = b / (z * c)
else if (b <= (-1.05d-133)) then
tmp = (9.0d0 * (x * y)) / (z * c)
else if (b <= 8.2d-77) then
tmp = (a * (t * (-4.0d0))) / c
else
tmp = (b / c) / z
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -0.0077) {
tmp = b / (z * c);
} else if (b <= -1.05e-133) {
tmp = (9.0 * (x * y)) / (z * c);
} else if (b <= 8.2e-77) {
tmp = (a * (t * -4.0)) / c;
} else {
tmp = (b / c) / z;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): tmp = 0 if b <= -0.0077: tmp = b / (z * c) elif b <= -1.05e-133: tmp = (9.0 * (x * y)) / (z * c) elif b <= 8.2e-77: tmp = (a * (t * -4.0)) / c else: tmp = (b / c) / z return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) tmp = 0.0 if (b <= -0.0077) tmp = Float64(b / Float64(z * c)); elseif (b <= -1.05e-133) tmp = Float64(Float64(9.0 * Float64(x * y)) / Float64(z * c)); elseif (b <= 8.2e-77) tmp = Float64(Float64(a * Float64(t * -4.0)) / c); else tmp = Float64(Float64(b / c) / z); end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
tmp = 0.0;
if (b <= -0.0077)
tmp = b / (z * c);
elseif (b <= -1.05e-133)
tmp = (9.0 * (x * y)) / (z * c);
elseif (b <= 8.2e-77)
tmp = (a * (t * -4.0)) / c;
else
tmp = (b / c) / z;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[b, -0.0077], N[(b / N[(z * c), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -1.05e-133], N[(N[(9.0 * N[(x * y), $MachinePrecision]), $MachinePrecision] / N[(z * c), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 8.2e-77], N[(N[(a * N[(t * -4.0), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], N[(N[(b / c), $MachinePrecision] / z), $MachinePrecision]]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
\mathbf{if}\;b \leq -0.0077:\\
\;\;\;\;\frac{b}{z \cdot c}\\
\mathbf{elif}\;b \leq -1.05 \cdot 10^{-133}:\\
\;\;\;\;\frac{9 \cdot \left(x \cdot y\right)}{z \cdot c}\\
\mathbf{elif}\;b \leq 8.2 \cdot 10^{-77}:\\
\;\;\;\;\frac{a \cdot \left(t \cdot -4\right)}{c}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{b}{c}}{z}\\
\end{array}
\end{array}
if b < -0.0077000000000000002Initial program 86.0%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6461.0%
Simplified61.0%
if -0.0077000000000000002 < b < -1.05e-133Initial program 90.4%
Taylor expanded in x around inf
*-lowering-*.f64N/A
*-lowering-*.f6458.2%
Simplified58.2%
if -1.05e-133 < b < 8.19999999999999925e-77Initial program 80.2%
Taylor expanded in z around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6455.9%
Simplified55.9%
if 8.19999999999999925e-77 < b Initial program 77.3%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6457.7%
Simplified57.7%
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6462.8%
Applied egg-rr62.8%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c)
:precision binary64
(if (<= b -7.2e-53)
(/ b (* z c))
(if (<= b -4.8e-137)
(* (/ (* 9.0 y) z) (/ x c))
(if (<= b 9.8e-77) (/ (* a (* t -4.0)) c) (/ (/ b c) z)))))assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -7.2e-53) {
tmp = b / (z * c);
} else if (b <= -4.8e-137) {
tmp = ((9.0 * y) / z) * (x / c);
} else if (b <= 9.8e-77) {
tmp = (a * (t * -4.0)) / c;
} else {
tmp = (b / c) / z;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if (b <= (-7.2d-53)) then
tmp = b / (z * c)
else if (b <= (-4.8d-137)) then
tmp = ((9.0d0 * y) / z) * (x / c)
else if (b <= 9.8d-77) then
tmp = (a * (t * (-4.0d0))) / c
else
tmp = (b / c) / z
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -7.2e-53) {
tmp = b / (z * c);
} else if (b <= -4.8e-137) {
tmp = ((9.0 * y) / z) * (x / c);
} else if (b <= 9.8e-77) {
tmp = (a * (t * -4.0)) / c;
} else {
tmp = (b / c) / z;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): tmp = 0 if b <= -7.2e-53: tmp = b / (z * c) elif b <= -4.8e-137: tmp = ((9.0 * y) / z) * (x / c) elif b <= 9.8e-77: tmp = (a * (t * -4.0)) / c else: tmp = (b / c) / z return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) tmp = 0.0 if (b <= -7.2e-53) tmp = Float64(b / Float64(z * c)); elseif (b <= -4.8e-137) tmp = Float64(Float64(Float64(9.0 * y) / z) * Float64(x / c)); elseif (b <= 9.8e-77) tmp = Float64(Float64(a * Float64(t * -4.0)) / c); else tmp = Float64(Float64(b / c) / z); end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
tmp = 0.0;
if (b <= -7.2e-53)
tmp = b / (z * c);
elseif (b <= -4.8e-137)
tmp = ((9.0 * y) / z) * (x / c);
elseif (b <= 9.8e-77)
tmp = (a * (t * -4.0)) / c;
else
tmp = (b / c) / z;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[b, -7.2e-53], N[(b / N[(z * c), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -4.8e-137], N[(N[(N[(9.0 * y), $MachinePrecision] / z), $MachinePrecision] * N[(x / c), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 9.8e-77], N[(N[(a * N[(t * -4.0), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], N[(N[(b / c), $MachinePrecision] / z), $MachinePrecision]]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
\mathbf{if}\;b \leq -7.2 \cdot 10^{-53}:\\
\;\;\;\;\frac{b}{z \cdot c}\\
\mathbf{elif}\;b \leq -4.8 \cdot 10^{-137}:\\
\;\;\;\;\frac{9 \cdot y}{z} \cdot \frac{x}{c}\\
\mathbf{elif}\;b \leq 9.8 \cdot 10^{-77}:\\
\;\;\;\;\frac{a \cdot \left(t \cdot -4\right)}{c}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{b}{c}}{z}\\
\end{array}
\end{array}
if b < -7.1999999999999998e-53Initial program 85.1%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6456.5%
Simplified56.5%
if -7.1999999999999998e-53 < b < -4.8000000000000001e-137Initial program 91.3%
div-invN/A
associate-*r/N/A
*-commutativeN/A
times-fracN/A
*-lowering-*.f64N/A
Applied egg-rr75.0%
Taylor expanded in x around inf
*-lowering-*.f64N/A
*-lowering-*.f6456.9%
Simplified56.9%
un-div-invN/A
associate-/l/N/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6452.7%
Applied egg-rr52.7%
if -4.8000000000000001e-137 < b < 9.7999999999999994e-77Initial program 81.0%
Taylor expanded in z around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6456.5%
Simplified56.5%
if 9.7999999999999994e-77 < b Initial program 77.3%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6457.7%
Simplified57.7%
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6462.8%
Applied egg-rr62.8%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c)
:precision binary64
(if (<= b -7.5e-53)
(/ b (* z c))
(if (<= b -7.8e-134)
(* 9.0 (/ (/ (* x y) z) c))
(if (<= b 4.8e-77) (/ (* a (* t -4.0)) c) (/ (/ b c) z)))))assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -7.5e-53) {
tmp = b / (z * c);
} else if (b <= -7.8e-134) {
tmp = 9.0 * (((x * y) / z) / c);
} else if (b <= 4.8e-77) {
tmp = (a * (t * -4.0)) / c;
} else {
tmp = (b / c) / z;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if (b <= (-7.5d-53)) then
tmp = b / (z * c)
else if (b <= (-7.8d-134)) then
tmp = 9.0d0 * (((x * y) / z) / c)
else if (b <= 4.8d-77) then
tmp = (a * (t * (-4.0d0))) / c
else
tmp = (b / c) / z
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -7.5e-53) {
tmp = b / (z * c);
} else if (b <= -7.8e-134) {
tmp = 9.0 * (((x * y) / z) / c);
} else if (b <= 4.8e-77) {
tmp = (a * (t * -4.0)) / c;
} else {
tmp = (b / c) / z;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): tmp = 0 if b <= -7.5e-53: tmp = b / (z * c) elif b <= -7.8e-134: tmp = 9.0 * (((x * y) / z) / c) elif b <= 4.8e-77: tmp = (a * (t * -4.0)) / c else: tmp = (b / c) / z return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) tmp = 0.0 if (b <= -7.5e-53) tmp = Float64(b / Float64(z * c)); elseif (b <= -7.8e-134) tmp = Float64(9.0 * Float64(Float64(Float64(x * y) / z) / c)); elseif (b <= 4.8e-77) tmp = Float64(Float64(a * Float64(t * -4.0)) / c); else tmp = Float64(Float64(b / c) / z); end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
tmp = 0.0;
if (b <= -7.5e-53)
tmp = b / (z * c);
elseif (b <= -7.8e-134)
tmp = 9.0 * (((x * y) / z) / c);
elseif (b <= 4.8e-77)
tmp = (a * (t * -4.0)) / c;
else
tmp = (b / c) / z;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[b, -7.5e-53], N[(b / N[(z * c), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -7.8e-134], N[(9.0 * N[(N[(N[(x * y), $MachinePrecision] / z), $MachinePrecision] / c), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 4.8e-77], N[(N[(a * N[(t * -4.0), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], N[(N[(b / c), $MachinePrecision] / z), $MachinePrecision]]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
\mathbf{if}\;b \leq -7.5 \cdot 10^{-53}:\\
\;\;\;\;\frac{b}{z \cdot c}\\
\mathbf{elif}\;b \leq -7.8 \cdot 10^{-134}:\\
\;\;\;\;9 \cdot \frac{\frac{x \cdot y}{z}}{c}\\
\mathbf{elif}\;b \leq 4.8 \cdot 10^{-77}:\\
\;\;\;\;\frac{a \cdot \left(t \cdot -4\right)}{c}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{b}{c}}{z}\\
\end{array}
\end{array}
if b < -7.5000000000000001e-53Initial program 85.1%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6456.5%
Simplified56.5%
if -7.5000000000000001e-53 < b < -7.8000000000000002e-134Initial program 95.2%
Taylor expanded in x around inf
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6460.5%
Simplified60.5%
if -7.8000000000000002e-134 < b < 4.7999999999999998e-77Initial program 80.2%
Taylor expanded in z around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6455.9%
Simplified55.9%
if 4.7999999999999998e-77 < b Initial program 77.3%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6457.7%
Simplified57.7%
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6462.8%
Applied egg-rr62.8%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (/ (+ (* t (* a -4.0)) (/ b z)) c)))
(if (<= a -1.3e-192)
t_1
(if (<= a 1.26e+94) (/ (/ (+ b (* y (* x 9.0))) z) c) t_1))))assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((t * (a * -4.0)) + (b / z)) / c;
double tmp;
if (a <= -1.3e-192) {
tmp = t_1;
} else if (a <= 1.26e+94) {
tmp = ((b + (y * (x * 9.0))) / z) / c;
} else {
tmp = t_1;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: tmp
t_1 = ((t * (a * (-4.0d0))) + (b / z)) / c
if (a <= (-1.3d-192)) then
tmp = t_1
else if (a <= 1.26d+94) then
tmp = ((b + (y * (x * 9.0d0))) / z) / c
else
tmp = t_1
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((t * (a * -4.0)) + (b / z)) / c;
double tmp;
if (a <= -1.3e-192) {
tmp = t_1;
} else if (a <= 1.26e+94) {
tmp = ((b + (y * (x * 9.0))) / z) / c;
} else {
tmp = t_1;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): t_1 = ((t * (a * -4.0)) + (b / z)) / c tmp = 0 if a <= -1.3e-192: tmp = t_1 elif a <= 1.26e+94: tmp = ((b + (y * (x * 9.0))) / z) / c else: tmp = t_1 return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(Float64(t * Float64(a * -4.0)) + Float64(b / z)) / c) tmp = 0.0 if (a <= -1.3e-192) tmp = t_1; elseif (a <= 1.26e+94) tmp = Float64(Float64(Float64(b + Float64(y * Float64(x * 9.0))) / z) / c); else tmp = t_1; end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
t_1 = ((t * (a * -4.0)) + (b / z)) / c;
tmp = 0.0;
if (a <= -1.3e-192)
tmp = t_1;
elseif (a <= 1.26e+94)
tmp = ((b + (y * (x * 9.0))) / z) / c;
else
tmp = t_1;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(N[(t * N[(a * -4.0), $MachinePrecision]), $MachinePrecision] + N[(b / z), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision]}, If[LessEqual[a, -1.3e-192], t$95$1, If[LessEqual[a, 1.26e+94], N[(N[(N[(b + N[(y * N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision] / c), $MachinePrecision], t$95$1]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
t_1 := \frac{t \cdot \left(a \cdot -4\right) + \frac{b}{z}}{c}\\
\mathbf{if}\;a \leq -1.3 \cdot 10^{-192}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 1.26 \cdot 10^{+94}:\\
\;\;\;\;\frac{\frac{b + y \cdot \left(x \cdot 9\right)}{z}}{c}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -1.3000000000000001e-192 or 1.25999999999999997e94 < a Initial program 80.4%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified82.6%
Taylor expanded in x around 0
/-lowering-/.f6472.5%
Simplified72.5%
if -1.3000000000000001e-192 < a < 1.25999999999999997e94Initial program 84.1%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified93.1%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6492.9%
Applied egg-rr92.9%
Taylor expanded in z around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6479.9%
Simplified79.9%
*-commutativeN/A
div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6480.0%
Applied egg-rr80.0%
Final simplification75.8%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (/ (+ (* t (* a -4.0)) (/ b z)) c)))
(if (<= a -1.3e-192)
t_1
(if (<= a 1.56e+94) (/ (/ (+ b (* 9.0 (* x y))) z) c) t_1))))assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((t * (a * -4.0)) + (b / z)) / c;
double tmp;
if (a <= -1.3e-192) {
tmp = t_1;
} else if (a <= 1.56e+94) {
tmp = ((b + (9.0 * (x * y))) / z) / c;
} else {
tmp = t_1;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: tmp
t_1 = ((t * (a * (-4.0d0))) + (b / z)) / c
if (a <= (-1.3d-192)) then
tmp = t_1
else if (a <= 1.56d+94) then
tmp = ((b + (9.0d0 * (x * y))) / z) / c
else
tmp = t_1
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((t * (a * -4.0)) + (b / z)) / c;
double tmp;
if (a <= -1.3e-192) {
tmp = t_1;
} else if (a <= 1.56e+94) {
tmp = ((b + (9.0 * (x * y))) / z) / c;
} else {
tmp = t_1;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): t_1 = ((t * (a * -4.0)) + (b / z)) / c tmp = 0 if a <= -1.3e-192: tmp = t_1 elif a <= 1.56e+94: tmp = ((b + (9.0 * (x * y))) / z) / c else: tmp = t_1 return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(Float64(t * Float64(a * -4.0)) + Float64(b / z)) / c) tmp = 0.0 if (a <= -1.3e-192) tmp = t_1; elseif (a <= 1.56e+94) tmp = Float64(Float64(Float64(b + Float64(9.0 * Float64(x * y))) / z) / c); else tmp = t_1; end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
t_1 = ((t * (a * -4.0)) + (b / z)) / c;
tmp = 0.0;
if (a <= -1.3e-192)
tmp = t_1;
elseif (a <= 1.56e+94)
tmp = ((b + (9.0 * (x * y))) / z) / c;
else
tmp = t_1;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(N[(t * N[(a * -4.0), $MachinePrecision]), $MachinePrecision] + N[(b / z), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision]}, If[LessEqual[a, -1.3e-192], t$95$1, If[LessEqual[a, 1.56e+94], N[(N[(N[(b + N[(9.0 * N[(x * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision] / c), $MachinePrecision], t$95$1]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
t_1 := \frac{t \cdot \left(a \cdot -4\right) + \frac{b}{z}}{c}\\
\mathbf{if}\;a \leq -1.3 \cdot 10^{-192}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 1.56 \cdot 10^{+94}:\\
\;\;\;\;\frac{\frac{b + 9 \cdot \left(x \cdot y\right)}{z}}{c}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -1.3000000000000001e-192 or 1.5600000000000001e94 < a Initial program 80.4%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified82.6%
Taylor expanded in x around 0
/-lowering-/.f6472.5%
Simplified72.5%
if -1.3000000000000001e-192 < a < 1.5600000000000001e94Initial program 84.1%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified93.1%
Taylor expanded in z around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6480.0%
Simplified80.0%
Final simplification75.7%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. (FPCore (x y z t a b c) :precision binary64 (if (<= b -1.55e+45) (/ b (* z c)) (if (<= b 1e-76) (* a (/ (* t -4.0) c)) (/ (/ b c) z))))
assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -1.55e+45) {
tmp = b / (z * c);
} else if (b <= 1e-76) {
tmp = a * ((t * -4.0) / c);
} else {
tmp = (b / c) / z;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if (b <= (-1.55d+45)) then
tmp = b / (z * c)
else if (b <= 1d-76) then
tmp = a * ((t * (-4.0d0)) / c)
else
tmp = (b / c) / z
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -1.55e+45) {
tmp = b / (z * c);
} else if (b <= 1e-76) {
tmp = a * ((t * -4.0) / c);
} else {
tmp = (b / c) / z;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): tmp = 0 if b <= -1.55e+45: tmp = b / (z * c) elif b <= 1e-76: tmp = a * ((t * -4.0) / c) else: tmp = (b / c) / z return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) tmp = 0.0 if (b <= -1.55e+45) tmp = Float64(b / Float64(z * c)); elseif (b <= 1e-76) tmp = Float64(a * Float64(Float64(t * -4.0) / c)); else tmp = Float64(Float64(b / c) / z); end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
tmp = 0.0;
if (b <= -1.55e+45)
tmp = b / (z * c);
elseif (b <= 1e-76)
tmp = a * ((t * -4.0) / c);
else
tmp = (b / c) / z;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[b, -1.55e+45], N[(b / N[(z * c), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1e-76], N[(a * N[(N[(t * -4.0), $MachinePrecision] / c), $MachinePrecision]), $MachinePrecision], N[(N[(b / c), $MachinePrecision] / z), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.55 \cdot 10^{+45}:\\
\;\;\;\;\frac{b}{z \cdot c}\\
\mathbf{elif}\;b \leq 10^{-76}:\\
\;\;\;\;a \cdot \frac{t \cdot -4}{c}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{b}{c}}{z}\\
\end{array}
\end{array}
if b < -1.54999999999999994e45Initial program 87.1%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6466.5%
Simplified66.5%
if -1.54999999999999994e45 < b < 9.99999999999999927e-77Initial program 82.6%
Taylor expanded in z around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6449.3%
Simplified49.3%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6448.2%
Applied egg-rr48.2%
if 9.99999999999999927e-77 < b Initial program 77.3%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6457.7%
Simplified57.7%
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6462.8%
Applied egg-rr62.8%
Final simplification56.2%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. (FPCore (x y z t a b c) :precision binary64 (if (<= b -4e+45) (/ b (* z c)) (if (<= b 9e-77) (* (* t a) (/ -4.0 c)) (/ (/ b c) z))))
assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -4e+45) {
tmp = b / (z * c);
} else if (b <= 9e-77) {
tmp = (t * a) * (-4.0 / c);
} else {
tmp = (b / c) / z;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if (b <= (-4d+45)) then
tmp = b / (z * c)
else if (b <= 9d-77) then
tmp = (t * a) * ((-4.0d0) / c)
else
tmp = (b / c) / z
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -4e+45) {
tmp = b / (z * c);
} else if (b <= 9e-77) {
tmp = (t * a) * (-4.0 / c);
} else {
tmp = (b / c) / z;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): tmp = 0 if b <= -4e+45: tmp = b / (z * c) elif b <= 9e-77: tmp = (t * a) * (-4.0 / c) else: tmp = (b / c) / z return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) tmp = 0.0 if (b <= -4e+45) tmp = Float64(b / Float64(z * c)); elseif (b <= 9e-77) tmp = Float64(Float64(t * a) * Float64(-4.0 / c)); else tmp = Float64(Float64(b / c) / z); end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
tmp = 0.0;
if (b <= -4e+45)
tmp = b / (z * c);
elseif (b <= 9e-77)
tmp = (t * a) * (-4.0 / c);
else
tmp = (b / c) / z;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[b, -4e+45], N[(b / N[(z * c), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 9e-77], N[(N[(t * a), $MachinePrecision] * N[(-4.0 / c), $MachinePrecision]), $MachinePrecision], N[(N[(b / c), $MachinePrecision] / z), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
\mathbf{if}\;b \leq -4 \cdot 10^{+45}:\\
\;\;\;\;\frac{b}{z \cdot c}\\
\mathbf{elif}\;b \leq 9 \cdot 10^{-77}:\\
\;\;\;\;\left(t \cdot a\right) \cdot \frac{-4}{c}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{b}{c}}{z}\\
\end{array}
\end{array}
if b < -3.9999999999999997e45Initial program 87.1%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6466.5%
Simplified66.5%
if -3.9999999999999997e45 < b < 9.0000000000000001e-77Initial program 82.6%
Taylor expanded in z around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6449.3%
Simplified49.3%
associate-*r*N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6449.3%
Applied egg-rr49.3%
if 9.0000000000000001e-77 < b Initial program 77.3%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6457.7%
Simplified57.7%
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6462.8%
Applied egg-rr62.8%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. (FPCore (x y z t a b c) :precision binary64 (if (<= b -2.6e+45) (/ b (* z c)) (if (<= b 1e-76) (* t (/ (* a -4.0) c)) (/ (/ b c) z))))
assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -2.6e+45) {
tmp = b / (z * c);
} else if (b <= 1e-76) {
tmp = t * ((a * -4.0) / c);
} else {
tmp = (b / c) / z;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if (b <= (-2.6d+45)) then
tmp = b / (z * c)
else if (b <= 1d-76) then
tmp = t * ((a * (-4.0d0)) / c)
else
tmp = (b / c) / z
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -2.6e+45) {
tmp = b / (z * c);
} else if (b <= 1e-76) {
tmp = t * ((a * -4.0) / c);
} else {
tmp = (b / c) / z;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): tmp = 0 if b <= -2.6e+45: tmp = b / (z * c) elif b <= 1e-76: tmp = t * ((a * -4.0) / c) else: tmp = (b / c) / z return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) tmp = 0.0 if (b <= -2.6e+45) tmp = Float64(b / Float64(z * c)); elseif (b <= 1e-76) tmp = Float64(t * Float64(Float64(a * -4.0) / c)); else tmp = Float64(Float64(b / c) / z); end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
tmp = 0.0;
if (b <= -2.6e+45)
tmp = b / (z * c);
elseif (b <= 1e-76)
tmp = t * ((a * -4.0) / c);
else
tmp = (b / c) / z;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[b, -2.6e+45], N[(b / N[(z * c), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1e-76], N[(t * N[(N[(a * -4.0), $MachinePrecision] / c), $MachinePrecision]), $MachinePrecision], N[(N[(b / c), $MachinePrecision] / z), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
\mathbf{if}\;b \leq -2.6 \cdot 10^{+45}:\\
\;\;\;\;\frac{b}{z \cdot c}\\
\mathbf{elif}\;b \leq 10^{-76}:\\
\;\;\;\;t \cdot \frac{a \cdot -4}{c}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{b}{c}}{z}\\
\end{array}
\end{array}
if b < -2.60000000000000007e45Initial program 87.1%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6466.5%
Simplified66.5%
if -2.60000000000000007e45 < b < 9.99999999999999927e-77Initial program 82.6%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified90.1%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
associate-*r*N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6490.8%
Applied egg-rr90.8%
Taylor expanded in z around inf
associate-*r/N/A
associate-*r*N/A
associate-*l/N/A
associate-*r/N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f6450.1%
Simplified50.1%
if 9.99999999999999927e-77 < b Initial program 77.3%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6457.7%
Simplified57.7%
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
/-lowering-/.f6462.8%
Applied egg-rr62.8%
Final simplification57.1%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. (FPCore (x y z t a b c) :precision binary64 (/ (+ (/ (+ b (* x (* 9.0 y))) z) (* t (* a -4.0))) c))
assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
return (((b + (x * (9.0 * y))) / z) + (t * (a * -4.0))) / c;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
code = (((b + (x * (9.0d0 * y))) / z) + (t * (a * (-4.0d0)))) / c
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return (((b + (x * (9.0 * y))) / z) + (t * (a * -4.0))) / c;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): return (((b + (x * (9.0 * y))) / z) + (t * (a * -4.0))) / c
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) return Float64(Float64(Float64(Float64(b + Float64(x * Float64(9.0 * y))) / z) + Float64(t * Float64(a * -4.0))) / c) end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp = code(x, y, z, t, a, b, c)
tmp = (((b + (x * (9.0 * y))) / z) + (t * (a * -4.0))) / c;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_] := N[(N[(N[(N[(b + N[(x * N[(9.0 * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision] + N[(t * N[(a * -4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\frac{\frac{b + x \cdot \left(9 \cdot y\right)}{z} + t \cdot \left(a \cdot -4\right)}{c}
\end{array}
Initial program 82.0%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified87.2%
Final simplification87.2%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. (FPCore (x y z t a b c) :precision binary64 (if (<= y 1.05e+249) (/ (+ (* t (* a -4.0)) (/ b z)) c) (/ (/ y (/ z (* x 9.0))) c)))
assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (y <= 1.05e+249) {
tmp = ((t * (a * -4.0)) + (b / z)) / c;
} else {
tmp = (y / (z / (x * 9.0))) / c;
}
return tmp;
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if (y <= 1.05d+249) then
tmp = ((t * (a * (-4.0d0))) + (b / z)) / c
else
tmp = (y / (z / (x * 9.0d0))) / c
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (y <= 1.05e+249) {
tmp = ((t * (a * -4.0)) + (b / z)) / c;
} else {
tmp = (y / (z / (x * 9.0))) / c;
}
return tmp;
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): tmp = 0 if y <= 1.05e+249: tmp = ((t * (a * -4.0)) + (b / z)) / c else: tmp = (y / (z / (x * 9.0))) / c return tmp
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) tmp = 0.0 if (y <= 1.05e+249) tmp = Float64(Float64(Float64(t * Float64(a * -4.0)) + Float64(b / z)) / c); else tmp = Float64(Float64(y / Float64(z / Float64(x * 9.0))) / c); end return tmp end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp_2 = code(x, y, z, t, a, b, c)
tmp = 0.0;
if (y <= 1.05e+249)
tmp = ((t * (a * -4.0)) + (b / z)) / c;
else
tmp = (y / (z / (x * 9.0))) / c;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[y, 1.05e+249], N[(N[(N[(t * N[(a * -4.0), $MachinePrecision]), $MachinePrecision] + N[(b / z), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], N[(N[(y / N[(z / N[(x * 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\begin{array}{l}
\mathbf{if}\;y \leq 1.05 \cdot 10^{+249}:\\
\;\;\;\;\frac{t \cdot \left(a \cdot -4\right) + \frac{b}{z}}{c}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{y}{\frac{z}{x \cdot 9}}}{c}\\
\end{array}
\end{array}
if y < 1.0499999999999999e249Initial program 82.0%
associate-/r*N/A
/-lowering-/.f64N/A
Simplified87.4%
Taylor expanded in x around 0
/-lowering-/.f6468.1%
Simplified68.1%
if 1.0499999999999999e249 < y Initial program 82.9%
div-invN/A
associate-*r/N/A
*-commutativeN/A
times-fracN/A
*-lowering-*.f64N/A
Applied egg-rr91.5%
Taylor expanded in x around inf
*-lowering-*.f64N/A
*-lowering-*.f6474.9%
Simplified74.9%
associate-*l/N/A
un-div-invN/A
associate-*r*N/A
associate-*l/N/A
*-commutativeN/A
/-lowering-/.f64N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6474.4%
Applied egg-rr74.4%
Final simplification68.4%
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. (FPCore (x y z t a b c) :precision binary64 (/ b (* z c)))
assert(x < y && y < z && z < t && t < a && a < b && b < c);
assert(x < y && y < z && z < t && t < a && a < b && b < c);
double code(double x, double y, double z, double t, double a, double b, double c) {
return b / (z * c);
}
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function.
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
code = b / (z * c)
end function
assert x < y && y < z && z < t && t < a && a < b && b < c;
assert x < y && y < z && z < t && t < a && a < b && b < c;
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return b / (z * c);
}
[x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) [x, y, z, t, a, b, c] = sort([x, y, z, t, a, b, c]) def code(x, y, z, t, a, b, c): return b / (z * c)
x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) x, y, z, t, a, b, c = sort([x, y, z, t, a, b, c]) function code(x, y, z, t, a, b, c) return Float64(b / Float64(z * c)) end
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
x, y, z, t, a, b, c = num2cell(sort([x, y, z, t, a, b, c])){:}
function tmp = code(x, y, z, t, a, b, c)
tmp = b / (z * c);
end
NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. NOTE: x, y, z, t, a, b, and c should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_] := N[(b / N[(z * c), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\\\
[x, y, z, t, a, b, c] = \mathsf{sort}([x, y, z, t, a, b, c])\\
\\
\frac{b}{z \cdot c}
\end{array}
Initial program 82.0%
Taylor expanded in b around inf
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6439.3%
Simplified39.3%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (/ b (* c z)))
(t_2 (* 4.0 (/ (* a t) c)))
(t_3 (* (* x 9.0) y))
(t_4 (+ (- t_3 (* (* (* z 4.0) t) a)) b))
(t_5 (/ t_4 (* z c)))
(t_6 (/ (+ (- t_3 (* (* z 4.0) (* t a))) b) (* z c))))
(if (< t_5 -1.100156740804105e-171)
t_6
(if (< t_5 0.0)
(/ (/ t_4 z) c)
(if (< t_5 1.1708877911747488e-53)
t_6
(if (< t_5 2.876823679546137e+130)
(- (+ (* (* 9.0 (/ y c)) (/ x z)) t_1) t_2)
(if (< t_5 1.3838515042456319e+158)
t_6
(- (+ (* 9.0 (* (/ y (* c z)) x)) t_1) t_2))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = b / (c * z);
double t_2 = 4.0 * ((a * t) / c);
double t_3 = (x * 9.0) * y;
double t_4 = (t_3 - (((z * 4.0) * t) * a)) + b;
double t_5 = t_4 / (z * c);
double t_6 = ((t_3 - ((z * 4.0) * (t * a))) + b) / (z * c);
double tmp;
if (t_5 < -1.100156740804105e-171) {
tmp = t_6;
} else if (t_5 < 0.0) {
tmp = (t_4 / z) / c;
} else if (t_5 < 1.1708877911747488e-53) {
tmp = t_6;
} else if (t_5 < 2.876823679546137e+130) {
tmp = (((9.0 * (y / c)) * (x / z)) + t_1) - t_2;
} else if (t_5 < 1.3838515042456319e+158) {
tmp = t_6;
} else {
tmp = ((9.0 * ((y / (c * z)) * x)) + t_1) - t_2;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
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), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: t_4
real(8) :: t_5
real(8) :: t_6
real(8) :: tmp
t_1 = b / (c * z)
t_2 = 4.0d0 * ((a * t) / c)
t_3 = (x * 9.0d0) * y
t_4 = (t_3 - (((z * 4.0d0) * t) * a)) + b
t_5 = t_4 / (z * c)
t_6 = ((t_3 - ((z * 4.0d0) * (t * a))) + b) / (z * c)
if (t_5 < (-1.100156740804105d-171)) then
tmp = t_6
else if (t_5 < 0.0d0) then
tmp = (t_4 / z) / c
else if (t_5 < 1.1708877911747488d-53) then
tmp = t_6
else if (t_5 < 2.876823679546137d+130) then
tmp = (((9.0d0 * (y / c)) * (x / z)) + t_1) - t_2
else if (t_5 < 1.3838515042456319d+158) then
tmp = t_6
else
tmp = ((9.0d0 * ((y / (c * z)) * x)) + t_1) - t_2
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = b / (c * z);
double t_2 = 4.0 * ((a * t) / c);
double t_3 = (x * 9.0) * y;
double t_4 = (t_3 - (((z * 4.0) * t) * a)) + b;
double t_5 = t_4 / (z * c);
double t_6 = ((t_3 - ((z * 4.0) * (t * a))) + b) / (z * c);
double tmp;
if (t_5 < -1.100156740804105e-171) {
tmp = t_6;
} else if (t_5 < 0.0) {
tmp = (t_4 / z) / c;
} else if (t_5 < 1.1708877911747488e-53) {
tmp = t_6;
} else if (t_5 < 2.876823679546137e+130) {
tmp = (((9.0 * (y / c)) * (x / z)) + t_1) - t_2;
} else if (t_5 < 1.3838515042456319e+158) {
tmp = t_6;
} else {
tmp = ((9.0 * ((y / (c * z)) * x)) + t_1) - t_2;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = b / (c * z) t_2 = 4.0 * ((a * t) / c) t_3 = (x * 9.0) * y t_4 = (t_3 - (((z * 4.0) * t) * a)) + b t_5 = t_4 / (z * c) t_6 = ((t_3 - ((z * 4.0) * (t * a))) + b) / (z * c) tmp = 0 if t_5 < -1.100156740804105e-171: tmp = t_6 elif t_5 < 0.0: tmp = (t_4 / z) / c elif t_5 < 1.1708877911747488e-53: tmp = t_6 elif t_5 < 2.876823679546137e+130: tmp = (((9.0 * (y / c)) * (x / z)) + t_1) - t_2 elif t_5 < 1.3838515042456319e+158: tmp = t_6 else: tmp = ((9.0 * ((y / (c * z)) * x)) + t_1) - t_2 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(b / Float64(c * z)) t_2 = Float64(4.0 * Float64(Float64(a * t) / c)) t_3 = Float64(Float64(x * 9.0) * y) t_4 = Float64(Float64(t_3 - Float64(Float64(Float64(z * 4.0) * t) * a)) + b) t_5 = Float64(t_4 / Float64(z * c)) t_6 = Float64(Float64(Float64(t_3 - Float64(Float64(z * 4.0) * Float64(t * a))) + b) / Float64(z * c)) tmp = 0.0 if (t_5 < -1.100156740804105e-171) tmp = t_6; elseif (t_5 < 0.0) tmp = Float64(Float64(t_4 / z) / c); elseif (t_5 < 1.1708877911747488e-53) tmp = t_6; elseif (t_5 < 2.876823679546137e+130) tmp = Float64(Float64(Float64(Float64(9.0 * Float64(y / c)) * Float64(x / z)) + t_1) - t_2); elseif (t_5 < 1.3838515042456319e+158) tmp = t_6; else tmp = Float64(Float64(Float64(9.0 * Float64(Float64(y / Float64(c * z)) * x)) + t_1) - t_2); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = b / (c * z); t_2 = 4.0 * ((a * t) / c); t_3 = (x * 9.0) * y; t_4 = (t_3 - (((z * 4.0) * t) * a)) + b; t_5 = t_4 / (z * c); t_6 = ((t_3 - ((z * 4.0) * (t * a))) + b) / (z * c); tmp = 0.0; if (t_5 < -1.100156740804105e-171) tmp = t_6; elseif (t_5 < 0.0) tmp = (t_4 / z) / c; elseif (t_5 < 1.1708877911747488e-53) tmp = t_6; elseif (t_5 < 2.876823679546137e+130) tmp = (((9.0 * (y / c)) * (x / z)) + t_1) - t_2; elseif (t_5 < 1.3838515042456319e+158) tmp = t_6; else tmp = ((9.0 * ((y / (c * z)) * x)) + t_1) - t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(b / N[(c * z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(4.0 * N[(N[(a * t), $MachinePrecision] / c), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(x * 9.0), $MachinePrecision] * y), $MachinePrecision]}, Block[{t$95$4 = N[(N[(t$95$3 - N[(N[(N[(z * 4.0), $MachinePrecision] * t), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]}, Block[{t$95$5 = N[(t$95$4 / N[(z * c), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$6 = N[(N[(N[(t$95$3 - N[(N[(z * 4.0), $MachinePrecision] * N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision] / N[(z * c), $MachinePrecision]), $MachinePrecision]}, If[Less[t$95$5, -1.100156740804105e-171], t$95$6, If[Less[t$95$5, 0.0], N[(N[(t$95$4 / z), $MachinePrecision] / c), $MachinePrecision], If[Less[t$95$5, 1.1708877911747488e-53], t$95$6, If[Less[t$95$5, 2.876823679546137e+130], N[(N[(N[(N[(9.0 * N[(y / c), $MachinePrecision]), $MachinePrecision] * N[(x / z), $MachinePrecision]), $MachinePrecision] + t$95$1), $MachinePrecision] - t$95$2), $MachinePrecision], If[Less[t$95$5, 1.3838515042456319e+158], t$95$6, N[(N[(N[(9.0 * N[(N[(y / N[(c * z), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]), $MachinePrecision] + t$95$1), $MachinePrecision] - t$95$2), $MachinePrecision]]]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{b}{c \cdot z}\\
t_2 := 4 \cdot \frac{a \cdot t}{c}\\
t_3 := \left(x \cdot 9\right) \cdot y\\
t_4 := \left(t\_3 - \left(\left(z \cdot 4\right) \cdot t\right) \cdot a\right) + b\\
t_5 := \frac{t\_4}{z \cdot c}\\
t_6 := \frac{\left(t\_3 - \left(z \cdot 4\right) \cdot \left(t \cdot a\right)\right) + b}{z \cdot c}\\
\mathbf{if}\;t\_5 < -1.100156740804105 \cdot 10^{-171}:\\
\;\;\;\;t\_6\\
\mathbf{elif}\;t\_5 < 0:\\
\;\;\;\;\frac{\frac{t\_4}{z}}{c}\\
\mathbf{elif}\;t\_5 < 1.1708877911747488 \cdot 10^{-53}:\\
\;\;\;\;t\_6\\
\mathbf{elif}\;t\_5 < 2.876823679546137 \cdot 10^{+130}:\\
\;\;\;\;\left(\left(9 \cdot \frac{y}{c}\right) \cdot \frac{x}{z} + t\_1\right) - t\_2\\
\mathbf{elif}\;t\_5 < 1.3838515042456319 \cdot 10^{+158}:\\
\;\;\;\;t\_6\\
\mathbf{else}:\\
\;\;\;\;\left(9 \cdot \left(\frac{y}{c \cdot z} \cdot x\right) + t\_1\right) - t\_2\\
\end{array}
\end{array}
herbie shell --seed 2024158
(FPCore (x y z t a b c)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, J"
:precision binary64
:alt
(! :herbie-platform default (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) -220031348160821/200000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (/ (+ (- (* (* x 9) y) (* (* z 4) (* t a))) b) (* z c)) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) 0) (/ (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) z) c) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) 365902434742109/31250000000000000000000000000000000000000000000000000000000000000000) (/ (+ (- (* (* x 9) y) (* (* z 4) (* t a))) b) (* z c)) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) 28768236795461370000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (- (+ (* (* 9 (/ y c)) (/ x z)) (/ b (* c z))) (* 4 (/ (* a t) c))) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) 138385150424563190000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (/ (+ (- (* (* x 9) y) (* (* z 4) (* t a))) b) (* z c)) (- (+ (* 9 (* (/ y (* c z)) x)) (/ b (* c z))) (* 4 (/ (* a t) c)))))))))
(/ (+ (- (* (* x 9.0) y) (* (* (* z 4.0) t) a)) b) (* z c)))