
(FPCore (x y z t a b c i j k) :precision binary64 (- (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i)) (* (* j 27.0) k)))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k);
}
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
code = (((((((x * 18.0d0) * y) * z) * t) - ((a * 4.0d0) * t)) + (b * c)) - ((x * 4.0d0) * i)) - ((j * 27.0d0) * k)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k);
}
def code(x, y, z, t, a, b, c, i, j, k): return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k)
function code(x, y, z, t, a, b, c, i, j, k) return Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 18.0) * y) * z) * t) - Float64(Float64(a * 4.0) * t)) + Float64(b * c)) - Float64(Float64(x * 4.0) * i)) - Float64(Float64(j * 27.0) * k)) end
function tmp = code(x, y, z, t, a, b, c, i, j, k) tmp = (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k); end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := N[(N[(N[(N[(N[(N[(N[(N[(x * 18.0), $MachinePrecision] * y), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] - N[(N[(a * 4.0), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision] - N[(N[(x * 4.0), $MachinePrecision] * i), $MachinePrecision]), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t - \left(a \cdot 4\right) \cdot t\right) + b \cdot c\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 22 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b c i j k) :precision binary64 (- (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i)) (* (* j 27.0) k)))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k);
}
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
code = (((((((x * 18.0d0) * y) * z) * t) - ((a * 4.0d0) * t)) + (b * c)) - ((x * 4.0d0) * i)) - ((j * 27.0d0) * k)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k);
}
def code(x, y, z, t, a, b, c, i, j, k): return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k)
function code(x, y, z, t, a, b, c, i, j, k) return Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 18.0) * y) * z) * t) - Float64(Float64(a * 4.0) * t)) + Float64(b * c)) - Float64(Float64(x * 4.0) * i)) - Float64(Float64(j * 27.0) * k)) end
function tmp = code(x, y, z, t, a, b, c, i, j, k) tmp = (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k); end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := N[(N[(N[(N[(N[(N[(N[(N[(x * 18.0), $MachinePrecision] * y), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] - N[(N[(a * 4.0), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision] - N[(N[(x * 4.0), $MachinePrecision] * i), $MachinePrecision]), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t - \left(a \cdot 4\right) \cdot t\right) + b \cdot c\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k
\end{array}
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (fma z (* (* x 18.0) y) (* a -4.0))))
(if (<= t -3.4e+115)
(fma t_1 t (fma c b (- (fma k (* j 27.0) (* (* 4.0 x) i)))))
(if (<= t 1.16e-12)
(fma
-27.0
(* j k)
(fma (fma -4.0 i (* (* (* z t) y) 18.0)) x (fma (* a t) -4.0 (* b c))))
(fma (* -27.0 j) k (fma (* -4.0 x) i (fma t_1 t (* b c))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = fma(z, ((x * 18.0) * y), (a * -4.0));
double tmp;
if (t <= -3.4e+115) {
tmp = fma(t_1, t, fma(c, b, -fma(k, (j * 27.0), ((4.0 * x) * i))));
} else if (t <= 1.16e-12) {
tmp = fma(-27.0, (j * k), fma(fma(-4.0, i, (((z * t) * y) * 18.0)), x, fma((a * t), -4.0, (b * c))));
} else {
tmp = fma((-27.0 * j), k, fma((-4.0 * x), i, fma(t_1, t, (b * c))));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) t_1 = fma(z, Float64(Float64(x * 18.0) * y), Float64(a * -4.0)) tmp = 0.0 if (t <= -3.4e+115) tmp = fma(t_1, t, fma(c, b, Float64(-fma(k, Float64(j * 27.0), Float64(Float64(4.0 * x) * i))))); elseif (t <= 1.16e-12) tmp = fma(-27.0, Float64(j * k), fma(fma(-4.0, i, Float64(Float64(Float64(z * t) * y) * 18.0)), x, fma(Float64(a * t), -4.0, Float64(b * c)))); else tmp = fma(Float64(-27.0 * j), k, fma(Float64(-4.0 * x), i, fma(t_1, t, Float64(b * c)))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(z * N[(N[(x * 18.0), $MachinePrecision] * y), $MachinePrecision] + N[(a * -4.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -3.4e+115], N[(t$95$1 * t + N[(c * b + (-N[(k * N[(j * 27.0), $MachinePrecision] + N[(N[(4.0 * x), $MachinePrecision] * i), $MachinePrecision]), $MachinePrecision])), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.16e-12], N[(-27.0 * N[(j * k), $MachinePrecision] + N[(N[(-4.0 * i + N[(N[(N[(z * t), $MachinePrecision] * y), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * x + N[(N[(a * t), $MachinePrecision] * -4.0 + N[(b * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(-27.0 * j), $MachinePrecision] * k + N[(N[(-4.0 * x), $MachinePrecision] * i + N[(t$95$1 * t + N[(b * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(z, \left(x \cdot 18\right) \cdot y, a \cdot -4\right)\\
\mathbf{if}\;t \leq -3.4 \cdot 10^{+115}:\\
\;\;\;\;\mathsf{fma}\left(t\_1, t, \mathsf{fma}\left(c, b, -\mathsf{fma}\left(k, j \cdot 27, \left(4 \cdot x\right) \cdot i\right)\right)\right)\\
\mathbf{elif}\;t \leq 1.16 \cdot 10^{-12}:\\
\;\;\;\;\mathsf{fma}\left(-27, j \cdot k, \mathsf{fma}\left(\mathsf{fma}\left(-4, i, \left(\left(z \cdot t\right) \cdot y\right) \cdot 18\right), x, \mathsf{fma}\left(a \cdot t, -4, b \cdot c\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-27 \cdot j, k, \mathsf{fma}\left(-4 \cdot x, i, \mathsf{fma}\left(t\_1, t, b \cdot c\right)\right)\right)\\
\end{array}
\end{array}
if t < -3.4000000000000001e115Initial program 87.3%
lift--.f64N/A
lift--.f64N/A
associate--l-N/A
lift-+.f64N/A
associate--l+N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
distribute-rgt-out--N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites92.4%
if -3.4000000000000001e115 < t < 1.1599999999999999e-12Initial program 82.0%
Taylor expanded in x around 0
Applied rewrites89.3%
Applied rewrites96.4%
if 1.1599999999999999e-12 < t Initial program 92.1%
lift--.f64N/A
sub-negN/A
+-commutativeN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval96.0
lift--.f64N/A
sub-negN/A
+-commutativeN/A
Applied rewrites100.0%
Final simplification96.5%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1
(fma
(* -4.0 x)
i
(fma (fma z (* (* x 18.0) y) (* a -4.0)) t (* b c)))))
(if (<= t -3.4e+115)
(fma (- j) (* 27.0 k) t_1)
(if (<= t 1.06e-51)
(fma
-27.0
(* j k)
(fma (fma -4.0 i (* (* (* z t) y) 18.0)) x (fma (* a t) -4.0 (* b c))))
(fma (* -27.0 j) k t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = fma((-4.0 * x), i, fma(fma(z, ((x * 18.0) * y), (a * -4.0)), t, (b * c)));
double tmp;
if (t <= -3.4e+115) {
tmp = fma(-j, (27.0 * k), t_1);
} else if (t <= 1.06e-51) {
tmp = fma(-27.0, (j * k), fma(fma(-4.0, i, (((z * t) * y) * 18.0)), x, fma((a * t), -4.0, (b * c))));
} else {
tmp = fma((-27.0 * j), k, t_1);
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) t_1 = fma(Float64(-4.0 * x), i, fma(fma(z, Float64(Float64(x * 18.0) * y), Float64(a * -4.0)), t, Float64(b * c))) tmp = 0.0 if (t <= -3.4e+115) tmp = fma(Float64(-j), Float64(27.0 * k), t_1); elseif (t <= 1.06e-51) tmp = fma(-27.0, Float64(j * k), fma(fma(-4.0, i, Float64(Float64(Float64(z * t) * y) * 18.0)), x, fma(Float64(a * t), -4.0, Float64(b * c)))); else tmp = fma(Float64(-27.0 * j), k, t_1); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(-4.0 * x), $MachinePrecision] * i + N[(N[(z * N[(N[(x * 18.0), $MachinePrecision] * y), $MachinePrecision] + N[(a * -4.0), $MachinePrecision]), $MachinePrecision] * t + N[(b * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -3.4e+115], N[((-j) * N[(27.0 * k), $MachinePrecision] + t$95$1), $MachinePrecision], If[LessEqual[t, 1.06e-51], N[(-27.0 * N[(j * k), $MachinePrecision] + N[(N[(-4.0 * i + N[(N[(N[(z * t), $MachinePrecision] * y), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * x + N[(N[(a * t), $MachinePrecision] * -4.0 + N[(b * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(-27.0 * j), $MachinePrecision] * k + t$95$1), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(-4 \cdot x, i, \mathsf{fma}\left(\mathsf{fma}\left(z, \left(x \cdot 18\right) \cdot y, a \cdot -4\right), t, b \cdot c\right)\right)\\
\mathbf{if}\;t \leq -3.4 \cdot 10^{+115}:\\
\;\;\;\;\mathsf{fma}\left(-j, 27 \cdot k, t\_1\right)\\
\mathbf{elif}\;t \leq 1.06 \cdot 10^{-51}:\\
\;\;\;\;\mathsf{fma}\left(-27, j \cdot k, \mathsf{fma}\left(\mathsf{fma}\left(-4, i, \left(\left(z \cdot t\right) \cdot y\right) \cdot 18\right), x, \mathsf{fma}\left(a \cdot t, -4, b \cdot c\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-27 \cdot j, k, t\_1\right)\\
\end{array}
\end{array}
if t < -3.4000000000000001e115Initial program 87.3%
lift--.f64N/A
sub-negN/A
+-commutativeN/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
distribute-lft-neg-inN/A
lower-fma.f64N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6487.4
lift--.f64N/A
sub-negN/A
+-commutativeN/A
Applied rewrites89.9%
if -3.4000000000000001e115 < t < 1.0600000000000001e-51Initial program 81.7%
Taylor expanded in x around 0
Applied rewrites89.3%
Applied rewrites96.7%
if 1.0600000000000001e-51 < t Initial program 91.6%
lift--.f64N/A
sub-negN/A
+-commutativeN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval95.0
lift--.f64N/A
sub-negN/A
+-commutativeN/A
Applied rewrites98.5%
Final simplification96.1%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1
(fma
(* -27.0 j)
k
(fma
(* -4.0 x)
i
(fma (fma z (* (* x 18.0) y) (* a -4.0)) t (* b c))))))
(if (<= t -3.4e+115)
t_1
(if (<= t 1.06e-51)
(fma
-27.0
(* j k)
(fma (fma -4.0 i (* (* (* z t) y) 18.0)) x (fma (* a t) -4.0 (* b c))))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = fma((-27.0 * j), k, fma((-4.0 * x), i, fma(fma(z, ((x * 18.0) * y), (a * -4.0)), t, (b * c))));
double tmp;
if (t <= -3.4e+115) {
tmp = t_1;
} else if (t <= 1.06e-51) {
tmp = fma(-27.0, (j * k), fma(fma(-4.0, i, (((z * t) * y) * 18.0)), x, fma((a * t), -4.0, (b * c))));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) t_1 = fma(Float64(-27.0 * j), k, fma(Float64(-4.0 * x), i, fma(fma(z, Float64(Float64(x * 18.0) * y), Float64(a * -4.0)), t, Float64(b * c)))) tmp = 0.0 if (t <= -3.4e+115) tmp = t_1; elseif (t <= 1.06e-51) tmp = fma(-27.0, Float64(j * k), fma(fma(-4.0, i, Float64(Float64(Float64(z * t) * y) * 18.0)), x, fma(Float64(a * t), -4.0, Float64(b * c)))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(-27.0 * j), $MachinePrecision] * k + N[(N[(-4.0 * x), $MachinePrecision] * i + N[(N[(z * N[(N[(x * 18.0), $MachinePrecision] * y), $MachinePrecision] + N[(a * -4.0), $MachinePrecision]), $MachinePrecision] * t + N[(b * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -3.4e+115], t$95$1, If[LessEqual[t, 1.06e-51], N[(-27.0 * N[(j * k), $MachinePrecision] + N[(N[(-4.0 * i + N[(N[(N[(z * t), $MachinePrecision] * y), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * x + N[(N[(a * t), $MachinePrecision] * -4.0 + N[(b * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(-27 \cdot j, k, \mathsf{fma}\left(-4 \cdot x, i, \mathsf{fma}\left(\mathsf{fma}\left(z, \left(x \cdot 18\right) \cdot y, a \cdot -4\right), t, b \cdot c\right)\right)\right)\\
\mathbf{if}\;t \leq -3.4 \cdot 10^{+115}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 1.06 \cdot 10^{-51}:\\
\;\;\;\;\mathsf{fma}\left(-27, j \cdot k, \mathsf{fma}\left(\mathsf{fma}\left(-4, i, \left(\left(z \cdot t\right) \cdot y\right) \cdot 18\right), x, \mathsf{fma}\left(a \cdot t, -4, b \cdot c\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -3.4000000000000001e115 or 1.0600000000000001e-51 < t Initial program 89.8%
lift--.f64N/A
sub-negN/A
+-commutativeN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval91.9
lift--.f64N/A
sub-negN/A
+-commutativeN/A
Applied rewrites95.0%
if -3.4000000000000001e115 < t < 1.0600000000000001e-51Initial program 81.7%
Taylor expanded in x around 0
Applied rewrites89.3%
Applied rewrites96.7%
Final simplification96.1%
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= t -5e+14)
(fma
-27.0
(* j k)
(fma (* a t) -4.0 (fma c b (* (fma -4.0 i (* (* (* y z) t) 18.0)) x))))
(if (<= t 2.5e+88)
(fma
-27.0
(* j k)
(fma (fma -4.0 i (* (* (* z t) y) 18.0)) x (fma (* a t) -4.0 (* b c))))
(fma (* -27.0 j) k (fma (fma -4.0 a (* (* (* y z) x) 18.0)) t (* b c))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (t <= -5e+14) {
tmp = fma(-27.0, (j * k), fma((a * t), -4.0, fma(c, b, (fma(-4.0, i, (((y * z) * t) * 18.0)) * x))));
} else if (t <= 2.5e+88) {
tmp = fma(-27.0, (j * k), fma(fma(-4.0, i, (((z * t) * y) * 18.0)), x, fma((a * t), -4.0, (b * c))));
} else {
tmp = fma((-27.0 * j), k, fma(fma(-4.0, a, (((y * z) * x) * 18.0)), t, (b * c)));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (t <= -5e+14) tmp = fma(-27.0, Float64(j * k), fma(Float64(a * t), -4.0, fma(c, b, Float64(fma(-4.0, i, Float64(Float64(Float64(y * z) * t) * 18.0)) * x)))); elseif (t <= 2.5e+88) tmp = fma(-27.0, Float64(j * k), fma(fma(-4.0, i, Float64(Float64(Float64(z * t) * y) * 18.0)), x, fma(Float64(a * t), -4.0, Float64(b * c)))); else tmp = fma(Float64(-27.0 * j), k, fma(fma(-4.0, a, Float64(Float64(Float64(y * z) * x) * 18.0)), t, Float64(b * c))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[t, -5e+14], N[(-27.0 * N[(j * k), $MachinePrecision] + N[(N[(a * t), $MachinePrecision] * -4.0 + N[(c * b + N[(N[(-4.0 * i + N[(N[(N[(y * z), $MachinePrecision] * t), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.5e+88], N[(-27.0 * N[(j * k), $MachinePrecision] + N[(N[(-4.0 * i + N[(N[(N[(z * t), $MachinePrecision] * y), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * x + N[(N[(a * t), $MachinePrecision] * -4.0 + N[(b * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(-27.0 * j), $MachinePrecision] * k + N[(N[(-4.0 * a + N[(N[(N[(y * z), $MachinePrecision] * x), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * t + N[(b * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -5 \cdot 10^{+14}:\\
\;\;\;\;\mathsf{fma}\left(-27, j \cdot k, \mathsf{fma}\left(a \cdot t, -4, \mathsf{fma}\left(c, b, \mathsf{fma}\left(-4, i, \left(\left(y \cdot z\right) \cdot t\right) \cdot 18\right) \cdot x\right)\right)\right)\\
\mathbf{elif}\;t \leq 2.5 \cdot 10^{+88}:\\
\;\;\;\;\mathsf{fma}\left(-27, j \cdot k, \mathsf{fma}\left(\mathsf{fma}\left(-4, i, \left(\left(z \cdot t\right) \cdot y\right) \cdot 18\right), x, \mathsf{fma}\left(a \cdot t, -4, b \cdot c\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-27 \cdot j, k, \mathsf{fma}\left(\mathsf{fma}\left(-4, a, \left(\left(y \cdot z\right) \cdot x\right) \cdot 18\right), t, b \cdot c\right)\right)\\
\end{array}
\end{array}
if t < -5e14Initial program 86.9%
Taylor expanded in x around 0
Applied rewrites83.5%
Applied rewrites89.0%
if -5e14 < t < 2.49999999999999999e88Initial program 83.0%
Taylor expanded in x around 0
Applied rewrites89.7%
Applied rewrites96.9%
if 2.49999999999999999e88 < t Initial program 89.7%
Taylor expanded in i around 0
+-commutativeN/A
associate--r+N/A
metadata-evalN/A
cancel-sign-sub-invN/A
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
associate--l+N/A
Applied rewrites95.1%
Final simplification95.0%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (fma (fma i x (* a t)) -4.0 (fma c b (* (* j k) -27.0)))))
(if (<= (* 4.0 a) -1e+190)
t_1
(if (<= (* 4.0 a) 1e+48)
(fma -27.0 (* j k) (fma (fma -4.0 i (* (* (* y z) t) 18.0)) x (* b c)))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = fma(fma(i, x, (a * t)), -4.0, fma(c, b, ((j * k) * -27.0)));
double tmp;
if ((4.0 * a) <= -1e+190) {
tmp = t_1;
} else if ((4.0 * a) <= 1e+48) {
tmp = fma(-27.0, (j * k), fma(fma(-4.0, i, (((y * z) * t) * 18.0)), x, (b * c)));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) t_1 = fma(fma(i, x, Float64(a * t)), -4.0, fma(c, b, Float64(Float64(j * k) * -27.0))) tmp = 0.0 if (Float64(4.0 * a) <= -1e+190) tmp = t_1; elseif (Float64(4.0 * a) <= 1e+48) tmp = fma(-27.0, Float64(j * k), fma(fma(-4.0, i, Float64(Float64(Float64(y * z) * t) * 18.0)), x, Float64(b * c))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(i * x + N[(a * t), $MachinePrecision]), $MachinePrecision] * -4.0 + N[(c * b + N[(N[(j * k), $MachinePrecision] * -27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(4.0 * a), $MachinePrecision], -1e+190], t$95$1, If[LessEqual[N[(4.0 * a), $MachinePrecision], 1e+48], N[(-27.0 * N[(j * k), $MachinePrecision] + N[(N[(-4.0 * i + N[(N[(N[(y * z), $MachinePrecision] * t), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * x + N[(b * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(\mathsf{fma}\left(i, x, a \cdot t\right), -4, \mathsf{fma}\left(c, b, \left(j \cdot k\right) \cdot -27\right)\right)\\
\mathbf{if}\;4 \cdot a \leq -1 \cdot 10^{+190}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;4 \cdot a \leq 10^{+48}:\\
\;\;\;\;\mathsf{fma}\left(-27, j \cdot k, \mathsf{fma}\left(\mathsf{fma}\left(-4, i, \left(\left(y \cdot z\right) \cdot t\right) \cdot 18\right), x, b \cdot c\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (*.f64 a #s(literal 4 binary64)) < -1.0000000000000001e190 or 1.00000000000000004e48 < (*.f64 a #s(literal 4 binary64)) Initial program 83.4%
Taylor expanded in x around 0
Applied rewrites84.5%
Applied rewrites86.7%
Taylor expanded in y around 0
associate--r+N/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
cancel-sign-sub-invN/A
metadata-evalN/A
associate--r-N/A
cancel-sub-sign-invN/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
Applied rewrites86.3%
if -1.0000000000000001e190 < (*.f64 a #s(literal 4 binary64)) < 1.00000000000000004e48Initial program 85.6%
Taylor expanded in x around 0
Applied rewrites88.8%
Taylor expanded in a around 0
Applied rewrites84.9%
Final simplification85.4%
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= b 2.7e+76)
(fma
-27.0
(* j k)
(fma (* a t) -4.0 (fma c b (* (fma -4.0 i (* (* (* y z) t) 18.0)) x))))
(fma (* -27.0 k) j (* b c))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (b <= 2.7e+76) {
tmp = fma(-27.0, (j * k), fma((a * t), -4.0, fma(c, b, (fma(-4.0, i, (((y * z) * t) * 18.0)) * x))));
} else {
tmp = fma((-27.0 * k), j, (b * c));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (b <= 2.7e+76) tmp = fma(-27.0, Float64(j * k), fma(Float64(a * t), -4.0, fma(c, b, Float64(fma(-4.0, i, Float64(Float64(Float64(y * z) * t) * 18.0)) * x)))); else tmp = fma(Float64(-27.0 * k), j, Float64(b * c)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[b, 2.7e+76], N[(-27.0 * N[(j * k), $MachinePrecision] + N[(N[(a * t), $MachinePrecision] * -4.0 + N[(c * b + N[(N[(-4.0 * i + N[(N[(N[(y * z), $MachinePrecision] * t), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(-27.0 * k), $MachinePrecision] * j + N[(b * c), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq 2.7 \cdot 10^{+76}:\\
\;\;\;\;\mathsf{fma}\left(-27, j \cdot k, \mathsf{fma}\left(a \cdot t, -4, \mathsf{fma}\left(c, b, \mathsf{fma}\left(-4, i, \left(\left(y \cdot z\right) \cdot t\right) \cdot 18\right) \cdot x\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-27 \cdot k, j, b \cdot c\right)\\
\end{array}
\end{array}
if b < 2.6999999999999999e76Initial program 85.5%
Taylor expanded in x around 0
Applied rewrites89.5%
Applied rewrites90.5%
if 2.6999999999999999e76 < b Initial program 82.4%
Taylor expanded in x around 0
Applied rewrites78.4%
Applied rewrites82.3%
Taylor expanded in x around 0
associate--r+N/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6469.1
Applied rewrites69.1%
Taylor expanded in t around 0
Applied rewrites60.1%
Final simplification84.4%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1
(fma
(* -27.0 j)
k
(fma (fma -4.0 a (* (* (* y z) x) 18.0)) t (* b c)))))
(if (<= t -5.9e-118)
t_1
(if (<= t 2.65e+132)
(fma (fma i x (* a t)) -4.0 (fma c b (* (* j k) -27.0)))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = fma((-27.0 * j), k, fma(fma(-4.0, a, (((y * z) * x) * 18.0)), t, (b * c)));
double tmp;
if (t <= -5.9e-118) {
tmp = t_1;
} else if (t <= 2.65e+132) {
tmp = fma(fma(i, x, (a * t)), -4.0, fma(c, b, ((j * k) * -27.0)));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) t_1 = fma(Float64(-27.0 * j), k, fma(fma(-4.0, a, Float64(Float64(Float64(y * z) * x) * 18.0)), t, Float64(b * c))) tmp = 0.0 if (t <= -5.9e-118) tmp = t_1; elseif (t <= 2.65e+132) tmp = fma(fma(i, x, Float64(a * t)), -4.0, fma(c, b, Float64(Float64(j * k) * -27.0))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(-27.0 * j), $MachinePrecision] * k + N[(N[(-4.0 * a + N[(N[(N[(y * z), $MachinePrecision] * x), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * t + N[(b * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -5.9e-118], t$95$1, If[LessEqual[t, 2.65e+132], N[(N[(i * x + N[(a * t), $MachinePrecision]), $MachinePrecision] * -4.0 + N[(c * b + N[(N[(j * k), $MachinePrecision] * -27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(-27 \cdot j, k, \mathsf{fma}\left(\mathsf{fma}\left(-4, a, \left(\left(y \cdot z\right) \cdot x\right) \cdot 18\right), t, b \cdot c\right)\right)\\
\mathbf{if}\;t \leq -5.9 \cdot 10^{-118}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 2.65 \cdot 10^{+132}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(i, x, a \cdot t\right), -4, \mathsf{fma}\left(c, b, \left(j \cdot k\right) \cdot -27\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -5.9e-118 or 2.65e132 < t Initial program 88.5%
Taylor expanded in i around 0
+-commutativeN/A
associate--r+N/A
metadata-evalN/A
cancel-sign-sub-invN/A
+-commutativeN/A
associate--l+N/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
associate--l+N/A
Applied rewrites88.8%
if -5.9e-118 < t < 2.65e132Initial program 81.5%
Taylor expanded in x around 0
Applied rewrites87.4%
Applied rewrites88.9%
Taylor expanded in y around 0
associate--r+N/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
cancel-sign-sub-invN/A
metadata-evalN/A
associate--r-N/A
cancel-sub-sign-invN/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
Applied rewrites87.2%
Final simplification88.0%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (- (* (fma (* (* y z) x) 18.0 (* a -4.0)) t) (* (* j 27.0) k))))
(if (<= t -270.0)
t_1
(if (<= t 5.6e+132)
(fma (fma i x (* a t)) -4.0 (fma c b (* (* j k) -27.0)))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (fma(((y * z) * x), 18.0, (a * -4.0)) * t) - ((j * 27.0) * k);
double tmp;
if (t <= -270.0) {
tmp = t_1;
} else if (t <= 5.6e+132) {
tmp = fma(fma(i, x, (a * t)), -4.0, fma(c, b, ((j * k) * -27.0)));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(fma(Float64(Float64(y * z) * x), 18.0, Float64(a * -4.0)) * t) - Float64(Float64(j * 27.0) * k)) tmp = 0.0 if (t <= -270.0) tmp = t_1; elseif (t <= 5.6e+132) tmp = fma(fma(i, x, Float64(a * t)), -4.0, fma(c, b, Float64(Float64(j * k) * -27.0))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(N[(N[(N[(y * z), $MachinePrecision] * x), $MachinePrecision] * 18.0 + N[(a * -4.0), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -270.0], t$95$1, If[LessEqual[t, 5.6e+132], N[(N[(i * x + N[(a * t), $MachinePrecision]), $MachinePrecision] * -4.0 + N[(c * b + N[(N[(j * k), $MachinePrecision] * -27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(\left(y \cdot z\right) \cdot x, 18, a \cdot -4\right) \cdot t - \left(j \cdot 27\right) \cdot k\\
\mathbf{if}\;t \leq -270:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 5.6 \cdot 10^{+132}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(i, x, a \cdot t\right), -4, \mathsf{fma}\left(c, b, \left(j \cdot k\right) \cdot -27\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -270 or 5.5999999999999998e132 < t Initial program 90.0%
Taylor expanded in t around 0
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6435.1
Applied rewrites35.1%
Taylor expanded in b around 0
distribute-lft-outN/A
cancel-sub-sign-invN/A
associate-*r*N/A
metadata-evalN/A
distribute-lft-outN/A
lower-fma.f64N/A
lower-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
distribute-lft-outN/A
lower-*.f64N/A
lower-fma.f64N/A
lower-*.f6487.0
Applied rewrites87.0%
Taylor expanded in i around 0
Applied rewrites84.9%
if -270 < t < 5.5999999999999998e132Initial program 82.0%
Taylor expanded in x around 0
Applied rewrites88.1%
Applied rewrites88.7%
Taylor expanded in y around 0
associate--r+N/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
cancel-sign-sub-invN/A
metadata-evalN/A
associate--r-N/A
cancel-sub-sign-invN/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
Applied rewrites84.4%
Final simplification84.6%
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= t -6e+231)
(* (fma (* (* x y) z) 18.0 (* a -4.0)) t)
(if (<= t -3.35e-124)
(fma (* -4.0 t) a (fma b c (* (* -27.0 k) j)))
(if (<= t 0.0028)
(- (fma (* i x) -4.0 (* b c)) (* (* j 27.0) k))
(* (fma (* (* y z) x) 18.0 (* a -4.0)) t)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (t <= -6e+231) {
tmp = fma(((x * y) * z), 18.0, (a * -4.0)) * t;
} else if (t <= -3.35e-124) {
tmp = fma((-4.0 * t), a, fma(b, c, ((-27.0 * k) * j)));
} else if (t <= 0.0028) {
tmp = fma((i * x), -4.0, (b * c)) - ((j * 27.0) * k);
} else {
tmp = fma(((y * z) * x), 18.0, (a * -4.0)) * t;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (t <= -6e+231) tmp = Float64(fma(Float64(Float64(x * y) * z), 18.0, Float64(a * -4.0)) * t); elseif (t <= -3.35e-124) tmp = fma(Float64(-4.0 * t), a, fma(b, c, Float64(Float64(-27.0 * k) * j))); elseif (t <= 0.0028) tmp = Float64(fma(Float64(i * x), -4.0, Float64(b * c)) - Float64(Float64(j * 27.0) * k)); else tmp = Float64(fma(Float64(Float64(y * z) * x), 18.0, Float64(a * -4.0)) * t); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[t, -6e+231], N[(N[(N[(N[(x * y), $MachinePrecision] * z), $MachinePrecision] * 18.0 + N[(a * -4.0), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision], If[LessEqual[t, -3.35e-124], N[(N[(-4.0 * t), $MachinePrecision] * a + N[(b * c + N[(N[(-27.0 * k), $MachinePrecision] * j), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.0028], N[(N[(N[(i * x), $MachinePrecision] * -4.0 + N[(b * c), $MachinePrecision]), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(y * z), $MachinePrecision] * x), $MachinePrecision] * 18.0 + N[(a * -4.0), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -6 \cdot 10^{+231}:\\
\;\;\;\;\mathsf{fma}\left(\left(x \cdot y\right) \cdot z, 18, a \cdot -4\right) \cdot t\\
\mathbf{elif}\;t \leq -3.35 \cdot 10^{-124}:\\
\;\;\;\;\mathsf{fma}\left(-4 \cdot t, a, \mathsf{fma}\left(b, c, \left(-27 \cdot k\right) \cdot j\right)\right)\\
\mathbf{elif}\;t \leq 0.0028:\\
\;\;\;\;\mathsf{fma}\left(i \cdot x, -4, b \cdot c\right) - \left(j \cdot 27\right) \cdot k\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(y \cdot z\right) \cdot x, 18, a \cdot -4\right) \cdot t\\
\end{array}
\end{array}
if t < -6.0000000000000003e231Initial program 87.3%
Taylor expanded in x around 0
Applied rewrites81.9%
Taylor expanded in t around inf
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6491.8
Applied rewrites91.8%
if -6.0000000000000003e231 < t < -3.35e-124Initial program 86.5%
Taylor expanded in x around 0
Applied rewrites88.1%
Applied rewrites89.4%
Taylor expanded in x around 0
associate--r+N/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6467.2
Applied rewrites67.2%
Applied rewrites69.9%
if -3.35e-124 < t < 0.00279999999999999997Initial program 80.5%
Taylor expanded in t around 0
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6482.2
Applied rewrites82.2%
if 0.00279999999999999997 < t Initial program 91.8%
Taylor expanded in x around 0
Applied rewrites85.9%
Applied rewrites85.9%
Taylor expanded in t around inf
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6474.1
Applied rewrites74.1%
Final simplification77.6%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (fma (* -27.0 k) j (* (* a t) -4.0))))
(if (<= t -1.7e+167)
(* (fma -4.0 i (* (* (* y z) t) 18.0)) x)
(if (<= t -23000.0)
t_1
(if (<= t 2.4e-204)
(- (* b c) (* (* j 27.0) k))
(if (<= t 0.00075) (fma (* j k) -27.0 (* (* i x) -4.0)) t_1))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = fma((-27.0 * k), j, ((a * t) * -4.0));
double tmp;
if (t <= -1.7e+167) {
tmp = fma(-4.0, i, (((y * z) * t) * 18.0)) * x;
} else if (t <= -23000.0) {
tmp = t_1;
} else if (t <= 2.4e-204) {
tmp = (b * c) - ((j * 27.0) * k);
} else if (t <= 0.00075) {
tmp = fma((j * k), -27.0, ((i * x) * -4.0));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) t_1 = fma(Float64(-27.0 * k), j, Float64(Float64(a * t) * -4.0)) tmp = 0.0 if (t <= -1.7e+167) tmp = Float64(fma(-4.0, i, Float64(Float64(Float64(y * z) * t) * 18.0)) * x); elseif (t <= -23000.0) tmp = t_1; elseif (t <= 2.4e-204) tmp = Float64(Float64(b * c) - Float64(Float64(j * 27.0) * k)); elseif (t <= 0.00075) tmp = fma(Float64(j * k), -27.0, Float64(Float64(i * x) * -4.0)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(-27.0 * k), $MachinePrecision] * j + N[(N[(a * t), $MachinePrecision] * -4.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -1.7e+167], N[(N[(-4.0 * i + N[(N[(N[(y * z), $MachinePrecision] * t), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision], If[LessEqual[t, -23000.0], t$95$1, If[LessEqual[t, 2.4e-204], N[(N[(b * c), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.00075], N[(N[(j * k), $MachinePrecision] * -27.0 + N[(N[(i * x), $MachinePrecision] * -4.0), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(-27 \cdot k, j, \left(a \cdot t\right) \cdot -4\right)\\
\mathbf{if}\;t \leq -1.7 \cdot 10^{+167}:\\
\;\;\;\;\mathsf{fma}\left(-4, i, \left(\left(y \cdot z\right) \cdot t\right) \cdot 18\right) \cdot x\\
\mathbf{elif}\;t \leq -23000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 2.4 \cdot 10^{-204}:\\
\;\;\;\;b \cdot c - \left(j \cdot 27\right) \cdot k\\
\mathbf{elif}\;t \leq 0.00075:\\
\;\;\;\;\mathsf{fma}\left(j \cdot k, -27, \left(i \cdot x\right) \cdot -4\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -1.7e167Initial program 83.8%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6467.4
Applied rewrites67.4%
if -1.7e167 < t < -23000 or 7.5000000000000002e-4 < t Initial program 91.8%
Taylor expanded in t around 0
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6435.4
Applied rewrites35.4%
Taylor expanded in b around 0
distribute-lft-outN/A
cancel-sub-sign-invN/A
associate-*r*N/A
metadata-evalN/A
distribute-lft-outN/A
lower-fma.f64N/A
lower-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
distribute-lft-outN/A
lower-*.f64N/A
lower-fma.f64N/A
lower-*.f6486.7
Applied rewrites86.7%
Taylor expanded in x around 0
Applied rewrites58.8%
lift--.f64N/A
sub-negN/A
lift-*.f64N/A
distribute-lft-neg-outN/A
lift-*.f64N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
associate-*r*N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
Applied rewrites61.5%
if -23000 < t < 2.4e-204Initial program 81.4%
Taylor expanded in t around 0
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6482.8
Applied rewrites82.8%
Taylor expanded in x around 0
Applied rewrites70.7%
if 2.4e-204 < t < 7.5000000000000002e-4Initial program 82.2%
Taylor expanded in t around 0
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6471.5
Applied rewrites71.5%
Taylor expanded in x around inf
Applied rewrites55.4%
lift--.f64N/A
sub-negN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
lift-*.f64N/A
*-commutativeN/A
distribute-lft-neg-inN/A
metadata-evalN/A
lift-*.f64N/A
lift-*.f64N/A
+-commutativeN/A
Applied rewrites55.5%
Final simplification64.4%
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= t -2.15e+257)
(* (fma (* (* x y) z) 18.0 (* a -4.0)) t)
(if (<= t 1.08e+134)
(fma (fma i x (* a t)) -4.0 (fma c b (* (* j k) -27.0)))
(* (fma (* (* y z) x) 18.0 (* a -4.0)) t))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (t <= -2.15e+257) {
tmp = fma(((x * y) * z), 18.0, (a * -4.0)) * t;
} else if (t <= 1.08e+134) {
tmp = fma(fma(i, x, (a * t)), -4.0, fma(c, b, ((j * k) * -27.0)));
} else {
tmp = fma(((y * z) * x), 18.0, (a * -4.0)) * t;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (t <= -2.15e+257) tmp = Float64(fma(Float64(Float64(x * y) * z), 18.0, Float64(a * -4.0)) * t); elseif (t <= 1.08e+134) tmp = fma(fma(i, x, Float64(a * t)), -4.0, fma(c, b, Float64(Float64(j * k) * -27.0))); else tmp = Float64(fma(Float64(Float64(y * z) * x), 18.0, Float64(a * -4.0)) * t); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[t, -2.15e+257], N[(N[(N[(N[(x * y), $MachinePrecision] * z), $MachinePrecision] * 18.0 + N[(a * -4.0), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision], If[LessEqual[t, 1.08e+134], N[(N[(i * x + N[(a * t), $MachinePrecision]), $MachinePrecision] * -4.0 + N[(c * b + N[(N[(j * k), $MachinePrecision] * -27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(y * z), $MachinePrecision] * x), $MachinePrecision] * 18.0 + N[(a * -4.0), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -2.15 \cdot 10^{+257}:\\
\;\;\;\;\mathsf{fma}\left(\left(x \cdot y\right) \cdot z, 18, a \cdot -4\right) \cdot t\\
\mathbf{elif}\;t \leq 1.08 \cdot 10^{+134}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(i, x, a \cdot t\right), -4, \mathsf{fma}\left(c, b, \left(j \cdot k\right) \cdot -27\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(y \cdot z\right) \cdot x, 18, a \cdot -4\right) \cdot t\\
\end{array}
\end{array}
if t < -2.1499999999999999e257Initial program 84.4%
Taylor expanded in x around 0
Applied rewrites77.7%
Taylor expanded in t around inf
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6497.0
Applied rewrites97.0%
if -2.1499999999999999e257 < t < 1.0800000000000001e134Initial program 83.4%
Taylor expanded in x around 0
Applied rewrites87.7%
Applied rewrites89.1%
Taylor expanded in y around 0
associate--r+N/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
cancel-sign-sub-invN/A
metadata-evalN/A
associate--r-N/A
cancel-sub-sign-invN/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
Applied rewrites81.9%
if 1.0800000000000001e134 < t Initial program 94.1%
Taylor expanded in x around 0
Applied rewrites88.5%
Applied rewrites88.5%
Taylor expanded in t around inf
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6483.0
Applied rewrites83.0%
Final simplification82.9%
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= t -8500.0)
(* (fma (* (* x y) z) 18.0 (* a -4.0)) t)
(if (<= t 2.4e-204)
(- (* b c) (* (* j 27.0) k))
(if (<= t 0.001)
(fma (* j k) -27.0 (* (* i x) -4.0))
(* (fma (* (* y z) x) 18.0 (* a -4.0)) t)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (t <= -8500.0) {
tmp = fma(((x * y) * z), 18.0, (a * -4.0)) * t;
} else if (t <= 2.4e-204) {
tmp = (b * c) - ((j * 27.0) * k);
} else if (t <= 0.001) {
tmp = fma((j * k), -27.0, ((i * x) * -4.0));
} else {
tmp = fma(((y * z) * x), 18.0, (a * -4.0)) * t;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (t <= -8500.0) tmp = Float64(fma(Float64(Float64(x * y) * z), 18.0, Float64(a * -4.0)) * t); elseif (t <= 2.4e-204) tmp = Float64(Float64(b * c) - Float64(Float64(j * 27.0) * k)); elseif (t <= 0.001) tmp = fma(Float64(j * k), -27.0, Float64(Float64(i * x) * -4.0)); else tmp = Float64(fma(Float64(Float64(y * z) * x), 18.0, Float64(a * -4.0)) * t); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[t, -8500.0], N[(N[(N[(N[(x * y), $MachinePrecision] * z), $MachinePrecision] * 18.0 + N[(a * -4.0), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision], If[LessEqual[t, 2.4e-204], N[(N[(b * c), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.001], N[(N[(j * k), $MachinePrecision] * -27.0 + N[(N[(i * x), $MachinePrecision] * -4.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(y * z), $MachinePrecision] * x), $MachinePrecision] * 18.0 + N[(a * -4.0), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -8500:\\
\;\;\;\;\mathsf{fma}\left(\left(x \cdot y\right) \cdot z, 18, a \cdot -4\right) \cdot t\\
\mathbf{elif}\;t \leq 2.4 \cdot 10^{-204}:\\
\;\;\;\;b \cdot c - \left(j \cdot 27\right) \cdot k\\
\mathbf{elif}\;t \leq 0.001:\\
\;\;\;\;\mathsf{fma}\left(j \cdot k, -27, \left(i \cdot x\right) \cdot -4\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(y \cdot z\right) \cdot x, 18, a \cdot -4\right) \cdot t\\
\end{array}
\end{array}
if t < -8500Initial program 87.3%
Taylor expanded in x around 0
Applied rewrites84.1%
Taylor expanded in t around inf
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6467.8
Applied rewrites67.8%
if -8500 < t < 2.4e-204Initial program 81.4%
Taylor expanded in t around 0
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6482.8
Applied rewrites82.8%
Taylor expanded in x around 0
Applied rewrites70.7%
if 2.4e-204 < t < 1e-3Initial program 82.2%
Taylor expanded in t around 0
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6471.5
Applied rewrites71.5%
Taylor expanded in x around inf
Applied rewrites55.4%
lift--.f64N/A
sub-negN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
lift-*.f64N/A
*-commutativeN/A
distribute-lft-neg-inN/A
metadata-evalN/A
lift-*.f64N/A
lift-*.f64N/A
+-commutativeN/A
Applied rewrites55.5%
if 1e-3 < t Initial program 91.8%
Taylor expanded in x around 0
Applied rewrites85.9%
Applied rewrites85.9%
Taylor expanded in t around inf
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6474.1
Applied rewrites74.1%
Final simplification67.5%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* (fma (* (* x y) z) 18.0 (* a -4.0)) t)))
(if (<= t -8500.0)
t_1
(if (<= t 2.4e-204)
(- (* b c) (* (* j 27.0) k))
(if (<= t 0.001) (fma (* j k) -27.0 (* (* i x) -4.0)) t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = fma(((x * y) * z), 18.0, (a * -4.0)) * t;
double tmp;
if (t <= -8500.0) {
tmp = t_1;
} else if (t <= 2.4e-204) {
tmp = (b * c) - ((j * 27.0) * k);
} else if (t <= 0.001) {
tmp = fma((j * k), -27.0, ((i * x) * -4.0));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(fma(Float64(Float64(x * y) * z), 18.0, Float64(a * -4.0)) * t) tmp = 0.0 if (t <= -8500.0) tmp = t_1; elseif (t <= 2.4e-204) tmp = Float64(Float64(b * c) - Float64(Float64(j * 27.0) * k)); elseif (t <= 0.001) tmp = fma(Float64(j * k), -27.0, Float64(Float64(i * x) * -4.0)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(N[(N[(x * y), $MachinePrecision] * z), $MachinePrecision] * 18.0 + N[(a * -4.0), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision]}, If[LessEqual[t, -8500.0], t$95$1, If[LessEqual[t, 2.4e-204], N[(N[(b * c), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.001], N[(N[(j * k), $MachinePrecision] * -27.0 + N[(N[(i * x), $MachinePrecision] * -4.0), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(\left(x \cdot y\right) \cdot z, 18, a \cdot -4\right) \cdot t\\
\mathbf{if}\;t \leq -8500:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 2.4 \cdot 10^{-204}:\\
\;\;\;\;b \cdot c - \left(j \cdot 27\right) \cdot k\\
\mathbf{elif}\;t \leq 0.001:\\
\;\;\;\;\mathsf{fma}\left(j \cdot k, -27, \left(i \cdot x\right) \cdot -4\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -8500 or 1e-3 < t Initial program 89.4%
Taylor expanded in x around 0
Applied rewrites84.9%
Taylor expanded in t around inf
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6470.7
Applied rewrites70.7%
if -8500 < t < 2.4e-204Initial program 81.4%
Taylor expanded in t around 0
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6482.8
Applied rewrites82.8%
Taylor expanded in x around 0
Applied rewrites70.7%
if 2.4e-204 < t < 1e-3Initial program 82.2%
Taylor expanded in t around 0
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6471.5
Applied rewrites71.5%
Taylor expanded in x around inf
Applied rewrites55.4%
lift--.f64N/A
sub-negN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
lift-*.f64N/A
*-commutativeN/A
distribute-lft-neg-inN/A
metadata-evalN/A
lift-*.f64N/A
lift-*.f64N/A
+-commutativeN/A
Applied rewrites55.5%
Final simplification67.5%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* (* j 27.0) k)))
(if (<= t_1 -1e+83)
(* (* j k) -27.0)
(if (<= t_1 5e+152) (* (* i x) -4.0) (* (* -27.0 j) k)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (j * 27.0) * k;
double tmp;
if (t_1 <= -1e+83) {
tmp = (j * k) * -27.0;
} else if (t_1 <= 5e+152) {
tmp = (i * x) * -4.0;
} else {
tmp = (-27.0 * j) * k;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = (j * 27.0d0) * k
if (t_1 <= (-1d+83)) then
tmp = (j * k) * (-27.0d0)
else if (t_1 <= 5d+152) then
tmp = (i * x) * (-4.0d0)
else
tmp = ((-27.0d0) * j) * k
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 i, double j, double k) {
double t_1 = (j * 27.0) * k;
double tmp;
if (t_1 <= -1e+83) {
tmp = (j * k) * -27.0;
} else if (t_1 <= 5e+152) {
tmp = (i * x) * -4.0;
} else {
tmp = (-27.0 * j) * k;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i, j, k): t_1 = (j * 27.0) * k tmp = 0 if t_1 <= -1e+83: tmp = (j * k) * -27.0 elif t_1 <= 5e+152: tmp = (i * x) * -4.0 else: tmp = (-27.0 * j) * k return tmp
function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(j * 27.0) * k) tmp = 0.0 if (t_1 <= -1e+83) tmp = Float64(Float64(j * k) * -27.0); elseif (t_1 <= 5e+152) tmp = Float64(Float64(i * x) * -4.0); else tmp = Float64(Float64(-27.0 * j) * k); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k) t_1 = (j * 27.0) * k; tmp = 0.0; if (t_1 <= -1e+83) tmp = (j * k) * -27.0; elseif (t_1 <= 5e+152) tmp = (i * x) * -4.0; else tmp = (-27.0 * j) * k; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]}, If[LessEqual[t$95$1, -1e+83], N[(N[(j * k), $MachinePrecision] * -27.0), $MachinePrecision], If[LessEqual[t$95$1, 5e+152], N[(N[(i * x), $MachinePrecision] * -4.0), $MachinePrecision], N[(N[(-27.0 * j), $MachinePrecision] * k), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(j \cdot 27\right) \cdot k\\
\mathbf{if}\;t\_1 \leq -1 \cdot 10^{+83}:\\
\;\;\;\;\left(j \cdot k\right) \cdot -27\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{+152}:\\
\;\;\;\;\left(i \cdot x\right) \cdot -4\\
\mathbf{else}:\\
\;\;\;\;\left(-27 \cdot j\right) \cdot k\\
\end{array}
\end{array}
if (*.f64 (*.f64 j #s(literal 27 binary64)) k) < -1.00000000000000003e83Initial program 81.6%
Taylor expanded in j around inf
lower-*.f64N/A
*-commutativeN/A
lower-*.f6455.1
Applied rewrites55.1%
if -1.00000000000000003e83 < (*.f64 (*.f64 j #s(literal 27 binary64)) k) < 5e152Initial program 85.8%
Taylor expanded in i around inf
*-commutativeN/A
lower-*.f64N/A
lower-*.f6424.5
Applied rewrites24.5%
if 5e152 < (*.f64 (*.f64 j #s(literal 27 binary64)) k) Initial program 85.6%
Taylor expanded in j around inf
lower-*.f64N/A
*-commutativeN/A
lower-*.f6472.3
Applied rewrites72.3%
Applied rewrites74.7%
Final simplification39.2%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (fma (* -27.0 k) j (* (* a t) -4.0))))
(if (<= t -23000.0)
t_1
(if (<= t 2.4e-204)
(- (* b c) (* (* j 27.0) k))
(if (<= t 0.00075) (fma (* j k) -27.0 (* (* i x) -4.0)) t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = fma((-27.0 * k), j, ((a * t) * -4.0));
double tmp;
if (t <= -23000.0) {
tmp = t_1;
} else if (t <= 2.4e-204) {
tmp = (b * c) - ((j * 27.0) * k);
} else if (t <= 0.00075) {
tmp = fma((j * k), -27.0, ((i * x) * -4.0));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) t_1 = fma(Float64(-27.0 * k), j, Float64(Float64(a * t) * -4.0)) tmp = 0.0 if (t <= -23000.0) tmp = t_1; elseif (t <= 2.4e-204) tmp = Float64(Float64(b * c) - Float64(Float64(j * 27.0) * k)); elseif (t <= 0.00075) tmp = fma(Float64(j * k), -27.0, Float64(Float64(i * x) * -4.0)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(-27.0 * k), $MachinePrecision] * j + N[(N[(a * t), $MachinePrecision] * -4.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -23000.0], t$95$1, If[LessEqual[t, 2.4e-204], N[(N[(b * c), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.00075], N[(N[(j * k), $MachinePrecision] * -27.0 + N[(N[(i * x), $MachinePrecision] * -4.0), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(-27 \cdot k, j, \left(a \cdot t\right) \cdot -4\right)\\
\mathbf{if}\;t \leq -23000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 2.4 \cdot 10^{-204}:\\
\;\;\;\;b \cdot c - \left(j \cdot 27\right) \cdot k\\
\mathbf{elif}\;t \leq 0.00075:\\
\;\;\;\;\mathsf{fma}\left(j \cdot k, -27, \left(i \cdot x\right) \cdot -4\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -23000 or 7.5000000000000002e-4 < t Initial program 89.4%
Taylor expanded in t around 0
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6436.9
Applied rewrites36.9%
Taylor expanded in b around 0
distribute-lft-outN/A
cancel-sub-sign-invN/A
associate-*r*N/A
metadata-evalN/A
distribute-lft-outN/A
lower-fma.f64N/A
lower-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
distribute-lft-outN/A
lower-*.f64N/A
lower-fma.f64N/A
lower-*.f6486.8
Applied rewrites86.8%
Taylor expanded in x around 0
Applied rewrites56.1%
lift--.f64N/A
sub-negN/A
lift-*.f64N/A
distribute-lft-neg-outN/A
lift-*.f64N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
associate-*r*N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
Applied rewrites57.1%
if -23000 < t < 2.4e-204Initial program 81.4%
Taylor expanded in t around 0
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6482.8
Applied rewrites82.8%
Taylor expanded in x around 0
Applied rewrites70.7%
if 2.4e-204 < t < 7.5000000000000002e-4Initial program 82.2%
Taylor expanded in t around 0
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6471.5
Applied rewrites71.5%
Taylor expanded in x around inf
Applied rewrites55.4%
lift--.f64N/A
sub-negN/A
lift-*.f64N/A
distribute-lft-neg-inN/A
lift-*.f64N/A
*-commutativeN/A
distribute-lft-neg-inN/A
metadata-evalN/A
lift-*.f64N/A
lift-*.f64N/A
+-commutativeN/A
Applied rewrites55.5%
Final simplification61.9%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* (fma -4.0 i (* (* (* y z) t) 18.0)) x)))
(if (<= x -2.8e+46)
t_1
(if (<= x 1.2e+39) (fma (* -27.0 k) j (fma (* a t) -4.0 (* b c))) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = fma(-4.0, i, (((y * z) * t) * 18.0)) * x;
double tmp;
if (x <= -2.8e+46) {
tmp = t_1;
} else if (x <= 1.2e+39) {
tmp = fma((-27.0 * k), j, fma((a * t), -4.0, (b * c)));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(fma(-4.0, i, Float64(Float64(Float64(y * z) * t) * 18.0)) * x) tmp = 0.0 if (x <= -2.8e+46) tmp = t_1; elseif (x <= 1.2e+39) tmp = fma(Float64(-27.0 * k), j, fma(Float64(a * t), -4.0, Float64(b * c))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(-4.0 * i + N[(N[(N[(y * z), $MachinePrecision] * t), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]}, If[LessEqual[x, -2.8e+46], t$95$1, If[LessEqual[x, 1.2e+39], N[(N[(-27.0 * k), $MachinePrecision] * j + N[(N[(a * t), $MachinePrecision] * -4.0 + N[(b * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(-4, i, \left(\left(y \cdot z\right) \cdot t\right) \cdot 18\right) \cdot x\\
\mathbf{if}\;x \leq -2.8 \cdot 10^{+46}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;x \leq 1.2 \cdot 10^{+39}:\\
\;\;\;\;\mathsf{fma}\left(-27 \cdot k, j, \mathsf{fma}\left(a \cdot t, -4, b \cdot c\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if x < -2.80000000000000018e46 or 1.2e39 < x Initial program 73.5%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6463.6
Applied rewrites63.6%
if -2.80000000000000018e46 < x < 1.2e39Initial program 93.2%
Taylor expanded in x around 0
Applied rewrites85.4%
Applied rewrites86.1%
Taylor expanded in x around 0
associate--r+N/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6479.4
Applied rewrites79.4%
Final simplification72.7%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* (fma -4.0 i (* (* (* y z) t) 18.0)) x)))
(if (<= x -2.8e+46)
t_1
(if (<= x 1.2e+39) (fma -27.0 (* j k) (fma -4.0 (* a t) (* b c))) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = fma(-4.0, i, (((y * z) * t) * 18.0)) * x;
double tmp;
if (x <= -2.8e+46) {
tmp = t_1;
} else if (x <= 1.2e+39) {
tmp = fma(-27.0, (j * k), fma(-4.0, (a * t), (b * c)));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(fma(-4.0, i, Float64(Float64(Float64(y * z) * t) * 18.0)) * x) tmp = 0.0 if (x <= -2.8e+46) tmp = t_1; elseif (x <= 1.2e+39) tmp = fma(-27.0, Float64(j * k), fma(-4.0, Float64(a * t), Float64(b * c))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(-4.0 * i + N[(N[(N[(y * z), $MachinePrecision] * t), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]}, If[LessEqual[x, -2.8e+46], t$95$1, If[LessEqual[x, 1.2e+39], N[(-27.0 * N[(j * k), $MachinePrecision] + N[(-4.0 * N[(a * t), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(-4, i, \left(\left(y \cdot z\right) \cdot t\right) \cdot 18\right) \cdot x\\
\mathbf{if}\;x \leq -2.8 \cdot 10^{+46}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;x \leq 1.2 \cdot 10^{+39}:\\
\;\;\;\;\mathsf{fma}\left(-27, j \cdot k, \mathsf{fma}\left(-4, a \cdot t, b \cdot c\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if x < -2.80000000000000018e46 or 1.2e39 < x Initial program 73.5%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6463.6
Applied rewrites63.6%
if -2.80000000000000018e46 < x < 1.2e39Initial program 93.2%
Taylor expanded in x around 0
Applied rewrites85.4%
Taylor expanded in x around 0
associate--r+N/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6478.1
Applied rewrites78.1%
Final simplification71.9%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (fma (* -27.0 k) j (* (* a t) -4.0))))
(if (<= t -23000.0)
t_1
(if (<= t 0.0038) (fma (* -27.0 k) j (* b c)) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = fma((-27.0 * k), j, ((a * t) * -4.0));
double tmp;
if (t <= -23000.0) {
tmp = t_1;
} else if (t <= 0.0038) {
tmp = fma((-27.0 * k), j, (b * c));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) t_1 = fma(Float64(-27.0 * k), j, Float64(Float64(a * t) * -4.0)) tmp = 0.0 if (t <= -23000.0) tmp = t_1; elseif (t <= 0.0038) tmp = fma(Float64(-27.0 * k), j, Float64(b * c)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(-27.0 * k), $MachinePrecision] * j + N[(N[(a * t), $MachinePrecision] * -4.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -23000.0], t$95$1, If[LessEqual[t, 0.0038], N[(N[(-27.0 * k), $MachinePrecision] * j + N[(b * c), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(-27 \cdot k, j, \left(a \cdot t\right) \cdot -4\right)\\
\mathbf{if}\;t \leq -23000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 0.0038:\\
\;\;\;\;\mathsf{fma}\left(-27 \cdot k, j, b \cdot c\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -23000 or 0.00379999999999999999 < t Initial program 89.3%
Taylor expanded in t around 0
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6437.2
Applied rewrites37.2%
Taylor expanded in b around 0
distribute-lft-outN/A
cancel-sub-sign-invN/A
associate-*r*N/A
metadata-evalN/A
distribute-lft-outN/A
lower-fma.f64N/A
lower-*.f64N/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
distribute-lft-outN/A
lower-*.f64N/A
lower-fma.f64N/A
lower-*.f6486.7
Applied rewrites86.7%
Taylor expanded in x around 0
Applied rewrites56.6%
lift--.f64N/A
sub-negN/A
lift-*.f64N/A
distribute-lft-neg-outN/A
lift-*.f64N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
associate-*r*N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
Applied rewrites57.6%
if -23000 < t < 0.00379999999999999999Initial program 81.8%
Taylor expanded in x around 0
Applied rewrites89.0%
Applied rewrites89.7%
Taylor expanded in x around 0
associate--r+N/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6467.3
Applied rewrites67.3%
Taylor expanded in t around 0
Applied rewrites61.6%
Final simplification60.0%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* (* i x) -4.0)))
(if (<= i -3.8e+111)
t_1
(if (<= i 1.65e+261) (fma (* -27.0 k) j (* b c)) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (i * x) * -4.0;
double tmp;
if (i <= -3.8e+111) {
tmp = t_1;
} else if (i <= 1.65e+261) {
tmp = fma((-27.0 * k), j, (b * c));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(i * x) * -4.0) tmp = 0.0 if (i <= -3.8e+111) tmp = t_1; elseif (i <= 1.65e+261) tmp = fma(Float64(-27.0 * k), j, Float64(b * c)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(i * x), $MachinePrecision] * -4.0), $MachinePrecision]}, If[LessEqual[i, -3.8e+111], t$95$1, If[LessEqual[i, 1.65e+261], N[(N[(-27.0 * k), $MachinePrecision] * j + N[(b * c), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(i \cdot x\right) \cdot -4\\
\mathbf{if}\;i \leq -3.8 \cdot 10^{+111}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;i \leq 1.65 \cdot 10^{+261}:\\
\;\;\;\;\mathsf{fma}\left(-27 \cdot k, j, b \cdot c\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if i < -3.79999999999999976e111 or 1.65e261 < i Initial program 84.2%
Taylor expanded in i around inf
*-commutativeN/A
lower-*.f64N/A
lower-*.f6461.9
Applied rewrites61.9%
if -3.79999999999999976e111 < i < 1.65e261Initial program 85.0%
Taylor expanded in x around 0
Applied rewrites86.9%
Applied rewrites88.3%
Taylor expanded in x around 0
associate--r+N/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6468.5
Applied rewrites68.5%
Taylor expanded in t around 0
Applied rewrites51.4%
(FPCore (x y z t a b c i j k) :precision binary64 (* (* -27.0 k) j))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (-27.0 * k) * j;
}
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
code = ((-27.0d0) * k) * j
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (-27.0 * k) * j;
}
def code(x, y, z, t, a, b, c, i, j, k): return (-27.0 * k) * j
function code(x, y, z, t, a, b, c, i, j, k) return Float64(Float64(-27.0 * k) * j) end
function tmp = code(x, y, z, t, a, b, c, i, j, k) tmp = (-27.0 * k) * j; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := N[(N[(-27.0 * k), $MachinePrecision] * j), $MachinePrecision]
\begin{array}{l}
\\
\left(-27 \cdot k\right) \cdot j
\end{array}
Initial program 84.8%
Taylor expanded in j around inf
lower-*.f64N/A
*-commutativeN/A
lower-*.f6427.7
Applied rewrites27.7%
Applied rewrites27.7%
(FPCore (x y z t a b c i j k) :precision binary64 (* (* -27.0 j) k))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (-27.0 * j) * k;
}
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
code = ((-27.0d0) * j) * k
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (-27.0 * j) * k;
}
def code(x, y, z, t, a, b, c, i, j, k): return (-27.0 * j) * k
function code(x, y, z, t, a, b, c, i, j, k) return Float64(Float64(-27.0 * j) * k) end
function tmp = code(x, y, z, t, a, b, c, i, j, k) tmp = (-27.0 * j) * k; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := N[(N[(-27.0 * j), $MachinePrecision] * k), $MachinePrecision]
\begin{array}{l}
\\
\left(-27 \cdot j\right) \cdot k
\end{array}
Initial program 84.8%
Taylor expanded in j around inf
lower-*.f64N/A
*-commutativeN/A
lower-*.f6427.7
Applied rewrites27.7%
Applied rewrites28.1%
(FPCore (x y z t a b c i j k) :precision binary64 (* (* j k) -27.0))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (j * k) * -27.0;
}
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
code = (j * k) * (-27.0d0)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (j * k) * -27.0;
}
def code(x, y, z, t, a, b, c, i, j, k): return (j * k) * -27.0
function code(x, y, z, t, a, b, c, i, j, k) return Float64(Float64(j * k) * -27.0) end
function tmp = code(x, y, z, t, a, b, c, i, j, k) tmp = (j * k) * -27.0; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := N[(N[(j * k), $MachinePrecision] * -27.0), $MachinePrecision]
\begin{array}{l}
\\
\left(j \cdot k\right) \cdot -27
\end{array}
Initial program 84.8%
Taylor expanded in j around inf
lower-*.f64N/A
*-commutativeN/A
lower-*.f6427.7
Applied rewrites27.7%
Final simplification27.7%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* (+ (* a t) (* i x)) 4.0))
(t_2
(-
(- (* (* 18.0 t) (* (* x y) z)) t_1)
(- (* (* k j) 27.0) (* c b)))))
(if (< t -1.6210815397541398e-69)
t_2
(if (< t 165.68027943805222)
(+ (- (* (* 18.0 y) (* x (* z t))) t_1) (- (* c b) (* 27.0 (* k j))))
t_2))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = ((a * t) + (i * x)) * 4.0;
double t_2 = (((18.0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0) - (c * b));
double tmp;
if (t < -1.6210815397541398e-69) {
tmp = t_2;
} else if (t < 165.68027943805222) {
tmp = (((18.0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0 * (k * j)));
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i, j, k)
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), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = ((a * t) + (i * x)) * 4.0d0
t_2 = (((18.0d0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0d0) - (c * b))
if (t < (-1.6210815397541398d-69)) then
tmp = t_2
else if (t < 165.68027943805222d0) then
tmp = (((18.0d0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0d0 * (k * j)))
else
tmp = 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 i, double j, double k) {
double t_1 = ((a * t) + (i * x)) * 4.0;
double t_2 = (((18.0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0) - (c * b));
double tmp;
if (t < -1.6210815397541398e-69) {
tmp = t_2;
} else if (t < 165.68027943805222) {
tmp = (((18.0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0 * (k * j)));
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i, j, k): t_1 = ((a * t) + (i * x)) * 4.0 t_2 = (((18.0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0) - (c * b)) tmp = 0 if t < -1.6210815397541398e-69: tmp = t_2 elif t < 165.68027943805222: tmp = (((18.0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0 * (k * j))) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(Float64(a * t) + Float64(i * x)) * 4.0) t_2 = Float64(Float64(Float64(Float64(18.0 * t) * Float64(Float64(x * y) * z)) - t_1) - Float64(Float64(Float64(k * j) * 27.0) - Float64(c * b))) tmp = 0.0 if (t < -1.6210815397541398e-69) tmp = t_2; elseif (t < 165.68027943805222) tmp = Float64(Float64(Float64(Float64(18.0 * y) * Float64(x * Float64(z * t))) - t_1) + Float64(Float64(c * b) - Float64(27.0 * Float64(k * j)))); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k) t_1 = ((a * t) + (i * x)) * 4.0; t_2 = (((18.0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0) - (c * b)); tmp = 0.0; if (t < -1.6210815397541398e-69) tmp = t_2; elseif (t < 165.68027943805222) tmp = (((18.0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0 * (k * j))); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(N[(a * t), $MachinePrecision] + N[(i * x), $MachinePrecision]), $MachinePrecision] * 4.0), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(N[(18.0 * t), $MachinePrecision] * N[(N[(x * y), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision] - N[(N[(N[(k * j), $MachinePrecision] * 27.0), $MachinePrecision] - N[(c * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Less[t, -1.6210815397541398e-69], t$95$2, If[Less[t, 165.68027943805222], N[(N[(N[(N[(18.0 * y), $MachinePrecision] * N[(x * N[(z * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision] + N[(N[(c * b), $MachinePrecision] - N[(27.0 * N[(k * j), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$2]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(a \cdot t + i \cdot x\right) \cdot 4\\
t_2 := \left(\left(18 \cdot t\right) \cdot \left(\left(x \cdot y\right) \cdot z\right) - t\_1\right) - \left(\left(k \cdot j\right) \cdot 27 - c \cdot b\right)\\
\mathbf{if}\;t < -1.6210815397541398 \cdot 10^{-69}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t < 165.68027943805222:\\
\;\;\;\;\left(\left(18 \cdot y\right) \cdot \left(x \cdot \left(z \cdot t\right)\right) - t\_1\right) + \left(c \cdot b - 27 \cdot \left(k \cdot j\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
herbie shell --seed 2024298
(FPCore (x y z t a b c i j k)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, E"
:precision binary64
:alt
(! :herbie-platform default (if (< t -8105407698770699/5000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (- (- (* (* 18 t) (* (* x y) z)) (* (+ (* a t) (* i x)) 4)) (- (* (* k j) 27) (* c b))) (if (< t 8284013971902611/50000000000000) (+ (- (* (* 18 y) (* x (* z t))) (* (+ (* a t) (* i x)) 4)) (- (* c b) (* 27 (* k j)))) (- (- (* (* 18 t) (* (* x y) z)) (* (+ (* a t) (* i x)) 4)) (- (* (* k j) 27) (* c b))))))
(- (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i)) (* (* j 27.0) k)))