
(FPCore (x y z t a b c i j) :precision binary64 (+ (- (* x (- (* y z) (* t a))) (* b (- (* c z) (* t i)))) (* j (- (* c a) (* y i)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
return ((x * ((y * z) - (t * a))) - (b * ((c * z) - (t * i)))) + (j * ((c * a) - (y * i)));
}
real(8) function code(x, y, z, t, a, b, c, i, j)
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
code = ((x * ((y * z) - (t * a))) - (b * ((c * z) - (t * i)))) + (j * ((c * a) - (y * i)))
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
return ((x * ((y * z) - (t * a))) - (b * ((c * z) - (t * i)))) + (j * ((c * a) - (y * i)));
}
def code(x, y, z, t, a, b, c, i, j): return ((x * ((y * z) - (t * a))) - (b * ((c * z) - (t * i)))) + (j * ((c * a) - (y * i)))
function code(x, y, z, t, a, b, c, i, j) return Float64(Float64(Float64(x * Float64(Float64(y * z) - Float64(t * a))) - Float64(b * Float64(Float64(c * z) - Float64(t * i)))) + Float64(j * Float64(Float64(c * a) - Float64(y * i)))) end
function tmp = code(x, y, z, t, a, b, c, i, j) tmp = ((x * ((y * z) - (t * a))) - (b * ((c * z) - (t * i)))) + (j * ((c * a) - (y * i))); end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := N[(N[(N[(x * N[(N[(y * z), $MachinePrecision] - N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(b * N[(N[(c * z), $MachinePrecision] - N[(t * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(j * N[(N[(c * a), $MachinePrecision] - N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - t \cdot i\right)\right) + j \cdot \left(c \cdot a - y \cdot i\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 23 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b c i j) :precision binary64 (+ (- (* x (- (* y z) (* t a))) (* b (- (* c z) (* t i)))) (* j (- (* c a) (* y i)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
return ((x * ((y * z) - (t * a))) - (b * ((c * z) - (t * i)))) + (j * ((c * a) - (y * i)));
}
real(8) function code(x, y, z, t, a, b, c, i, j)
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
code = ((x * ((y * z) - (t * a))) - (b * ((c * z) - (t * i)))) + (j * ((c * a) - (y * i)))
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
return ((x * ((y * z) - (t * a))) - (b * ((c * z) - (t * i)))) + (j * ((c * a) - (y * i)));
}
def code(x, y, z, t, a, b, c, i, j): return ((x * ((y * z) - (t * a))) - (b * ((c * z) - (t * i)))) + (j * ((c * a) - (y * i)))
function code(x, y, z, t, a, b, c, i, j) return Float64(Float64(Float64(x * Float64(Float64(y * z) - Float64(t * a))) - Float64(b * Float64(Float64(c * z) - Float64(t * i)))) + Float64(j * Float64(Float64(c * a) - Float64(y * i)))) end
function tmp = code(x, y, z, t, a, b, c, i, j) tmp = ((x * ((y * z) - (t * a))) - (b * ((c * z) - (t * i)))) + (j * ((c * a) - (y * i))); end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := N[(N[(N[(x * N[(N[(y * z), $MachinePrecision] - N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(b * N[(N[(c * z), $MachinePrecision] - N[(t * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(j * N[(N[(c * a), $MachinePrecision] - N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - t \cdot i\right)\right) + j \cdot \left(c \cdot a - y \cdot i\right)
\end{array}
(FPCore (x y z t a b c i j)
:precision binary64
(let* ((t_1 (fma c (- (* a j) (* z b)) (* x (fma t (- a) (* y z))))))
(if (<= c -4.6e+71)
t_1
(if (<= c 1.1e+115)
(fma
a
(fma j c (- (* x t)))
(fma b (- (* t i) (* c z)) (* y (fma j (- i) (* x z)))))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = fma(c, ((a * j) - (z * b)), (x * fma(t, -a, (y * z))));
double tmp;
if (c <= -4.6e+71) {
tmp = t_1;
} else if (c <= 1.1e+115) {
tmp = fma(a, fma(j, c, -(x * t)), fma(b, ((t * i) - (c * z)), (y * fma(j, -i, (x * z)))));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j) t_1 = fma(c, Float64(Float64(a * j) - Float64(z * b)), Float64(x * fma(t, Float64(-a), Float64(y * z)))) tmp = 0.0 if (c <= -4.6e+71) tmp = t_1; elseif (c <= 1.1e+115) tmp = fma(a, fma(j, c, Float64(-Float64(x * t))), fma(b, Float64(Float64(t * i) - Float64(c * z)), Float64(y * fma(j, Float64(-i), Float64(x * z))))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(c * N[(N[(a * j), $MachinePrecision] - N[(z * b), $MachinePrecision]), $MachinePrecision] + N[(x * N[(t * (-a) + N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[c, -4.6e+71], t$95$1, If[LessEqual[c, 1.1e+115], N[(a * N[(j * c + (-N[(x * t), $MachinePrecision])), $MachinePrecision] + N[(b * N[(N[(t * i), $MachinePrecision] - N[(c * z), $MachinePrecision]), $MachinePrecision] + N[(y * N[(j * (-i) + N[(x * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(c, a \cdot j - z \cdot b, x \cdot \mathsf{fma}\left(t, -a, y \cdot z\right)\right)\\
\mathbf{if}\;c \leq -4.6 \cdot 10^{+71}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;c \leq 1.1 \cdot 10^{+115}:\\
\;\;\;\;\mathsf{fma}\left(a, \mathsf{fma}\left(j, c, -x \cdot t\right), \mathsf{fma}\left(b, t \cdot i - c \cdot z, y \cdot \mathsf{fma}\left(j, -i, x \cdot z\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if c < -4.6000000000000005e71 or 1.1e115 < c Initial program 60.1%
Taylor expanded in i around 0
sub-negN/A
+-commutativeN/A
associate-+l+N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
distribute-lft-neg-inN/A
mul-1-negN/A
distribute-rgt-inN/A
mul-1-negN/A
sub-negN/A
+-commutativeN/A
lower-fma.f64N/A
Applied rewrites85.2%
if -4.6000000000000005e71 < c < 1.1e115Initial program 76.9%
Taylor expanded in y around 0
Applied rewrites85.0%
Final simplification85.1%
(FPCore (x y z t a b c i j)
:precision binary64
(let* ((t_1 (* i (fma j (- y) (* t b)))) (t_2 (* c (- (* a j) (* z b)))))
(if (<= c -4.1e+101)
t_2
(if (<= c -1.25e-203)
(* y (fma j (- i) (* x z)))
(if (<= c 2.7e-155)
t_1
(if (<= c 2.45e-49)
(* i (* t (- b (/ (* x a) i))))
(if (<= c 3.8e+109) t_1 t_2)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = i * fma(j, -y, (t * b));
double t_2 = c * ((a * j) - (z * b));
double tmp;
if (c <= -4.1e+101) {
tmp = t_2;
} else if (c <= -1.25e-203) {
tmp = y * fma(j, -i, (x * z));
} else if (c <= 2.7e-155) {
tmp = t_1;
} else if (c <= 2.45e-49) {
tmp = i * (t * (b - ((x * a) / i)));
} else if (c <= 3.8e+109) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j) t_1 = Float64(i * fma(j, Float64(-y), Float64(t * b))) t_2 = Float64(c * Float64(Float64(a * j) - Float64(z * b))) tmp = 0.0 if (c <= -4.1e+101) tmp = t_2; elseif (c <= -1.25e-203) tmp = Float64(y * fma(j, Float64(-i), Float64(x * z))); elseif (c <= 2.7e-155) tmp = t_1; elseif (c <= 2.45e-49) tmp = Float64(i * Float64(t * Float64(b - Float64(Float64(x * a) / i)))); elseif (c <= 3.8e+109) tmp = t_1; else tmp = t_2; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(i * N[(j * (-y) + N[(t * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(c * N[(N[(a * j), $MachinePrecision] - N[(z * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[c, -4.1e+101], t$95$2, If[LessEqual[c, -1.25e-203], N[(y * N[(j * (-i) + N[(x * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 2.7e-155], t$95$1, If[LessEqual[c, 2.45e-49], N[(i * N[(t * N[(b - N[(N[(x * a), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 3.8e+109], t$95$1, t$95$2]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := i \cdot \mathsf{fma}\left(j, -y, t \cdot b\right)\\
t_2 := c \cdot \left(a \cdot j - z \cdot b\right)\\
\mathbf{if}\;c \leq -4.1 \cdot 10^{+101}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;c \leq -1.25 \cdot 10^{-203}:\\
\;\;\;\;y \cdot \mathsf{fma}\left(j, -i, x \cdot z\right)\\
\mathbf{elif}\;c \leq 2.7 \cdot 10^{-155}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;c \leq 2.45 \cdot 10^{-49}:\\
\;\;\;\;i \cdot \left(t \cdot \left(b - \frac{x \cdot a}{i}\right)\right)\\
\mathbf{elif}\;c \leq 3.8 \cdot 10^{+109}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if c < -4.1e101 or 3.80000000000000039e109 < c Initial program 55.7%
Taylor expanded in c around inf
lower-*.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6476.7
Applied rewrites76.7%
if -4.1e101 < c < -1.25e-203Initial program 79.8%
Taylor expanded in y around inf
lower-*.f64N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6453.4
Applied rewrites53.4%
if -1.25e-203 < c < 2.69999999999999981e-155 or 2.4500000000000001e-49 < c < 3.80000000000000039e109Initial program 75.8%
Taylor expanded in i around inf
lower-*.f64N/A
cancel-sign-sub-invN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
metadata-evalN/A
*-lft-identityN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-*.f6464.4
Applied rewrites64.4%
if 2.69999999999999981e-155 < c < 2.4500000000000001e-49Initial program 82.7%
Taylor expanded in i around -inf
Applied rewrites70.7%
Taylor expanded in t around inf
Applied rewrites70.6%
Final simplification65.8%
(FPCore (x y z t a b c i j)
:precision binary64
(let* ((t_1 (* i (fma j (- y) (* t b)))) (t_2 (* c (- (* a j) (* z b)))))
(if (<= c -4.1e+101)
t_2
(if (<= c -1.25e-203)
(* y (fma j (- i) (* x z)))
(if (<= c 7e-199)
t_1
(if (<= c 1.56e-49)
(* a (* t (fma c (/ j t) (- x))))
(if (<= c 3.8e+109) t_1 t_2)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = i * fma(j, -y, (t * b));
double t_2 = c * ((a * j) - (z * b));
double tmp;
if (c <= -4.1e+101) {
tmp = t_2;
} else if (c <= -1.25e-203) {
tmp = y * fma(j, -i, (x * z));
} else if (c <= 7e-199) {
tmp = t_1;
} else if (c <= 1.56e-49) {
tmp = a * (t * fma(c, (j / t), -x));
} else if (c <= 3.8e+109) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j) t_1 = Float64(i * fma(j, Float64(-y), Float64(t * b))) t_2 = Float64(c * Float64(Float64(a * j) - Float64(z * b))) tmp = 0.0 if (c <= -4.1e+101) tmp = t_2; elseif (c <= -1.25e-203) tmp = Float64(y * fma(j, Float64(-i), Float64(x * z))); elseif (c <= 7e-199) tmp = t_1; elseif (c <= 1.56e-49) tmp = Float64(a * Float64(t * fma(c, Float64(j / t), Float64(-x)))); elseif (c <= 3.8e+109) tmp = t_1; else tmp = t_2; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(i * N[(j * (-y) + N[(t * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(c * N[(N[(a * j), $MachinePrecision] - N[(z * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[c, -4.1e+101], t$95$2, If[LessEqual[c, -1.25e-203], N[(y * N[(j * (-i) + N[(x * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 7e-199], t$95$1, If[LessEqual[c, 1.56e-49], N[(a * N[(t * N[(c * N[(j / t), $MachinePrecision] + (-x)), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 3.8e+109], t$95$1, t$95$2]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := i \cdot \mathsf{fma}\left(j, -y, t \cdot b\right)\\
t_2 := c \cdot \left(a \cdot j - z \cdot b\right)\\
\mathbf{if}\;c \leq -4.1 \cdot 10^{+101}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;c \leq -1.25 \cdot 10^{-203}:\\
\;\;\;\;y \cdot \mathsf{fma}\left(j, -i, x \cdot z\right)\\
\mathbf{elif}\;c \leq 7 \cdot 10^{-199}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;c \leq 1.56 \cdot 10^{-49}:\\
\;\;\;\;a \cdot \left(t \cdot \mathsf{fma}\left(c, \frac{j}{t}, -x\right)\right)\\
\mathbf{elif}\;c \leq 3.8 \cdot 10^{+109}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if c < -4.1e101 or 3.80000000000000039e109 < c Initial program 55.7%
Taylor expanded in c around inf
lower-*.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6476.7
Applied rewrites76.7%
if -4.1e101 < c < -1.25e-203Initial program 79.8%
Taylor expanded in y around inf
lower-*.f64N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6453.4
Applied rewrites53.4%
if -1.25e-203 < c < 6.9999999999999998e-199 or 1.56000000000000008e-49 < c < 3.80000000000000039e109Initial program 76.2%
Taylor expanded in i around inf
lower-*.f64N/A
cancel-sign-sub-invN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
metadata-evalN/A
*-lft-identityN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-*.f6466.1
Applied rewrites66.1%
if 6.9999999999999998e-199 < c < 1.56000000000000008e-49Initial program 80.2%
Taylor expanded in a around inf
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6458.0
Applied rewrites58.0%
Taylor expanded in t around inf
Applied rewrites61.5%
Final simplification65.4%
(FPCore (x y z t a b c i j)
:precision binary64
(let* ((t_1 (fma c (- (* a j) (* z b)) (* x (fma t (- a) (* y z))))))
(if (<= c -6.2e-132)
t_1
(if (<= c 4.3e+77)
(fma j (fma a c (* i (- y))) (* t (fma a (- x) (* b i))))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = fma(c, ((a * j) - (z * b)), (x * fma(t, -a, (y * z))));
double tmp;
if (c <= -6.2e-132) {
tmp = t_1;
} else if (c <= 4.3e+77) {
tmp = fma(j, fma(a, c, (i * -y)), (t * fma(a, -x, (b * i))));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j) t_1 = fma(c, Float64(Float64(a * j) - Float64(z * b)), Float64(x * fma(t, Float64(-a), Float64(y * z)))) tmp = 0.0 if (c <= -6.2e-132) tmp = t_1; elseif (c <= 4.3e+77) tmp = fma(j, fma(a, c, Float64(i * Float64(-y))), Float64(t * fma(a, Float64(-x), Float64(b * i)))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(c * N[(N[(a * j), $MachinePrecision] - N[(z * b), $MachinePrecision]), $MachinePrecision] + N[(x * N[(t * (-a) + N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[c, -6.2e-132], t$95$1, If[LessEqual[c, 4.3e+77], N[(j * N[(a * c + N[(i * (-y)), $MachinePrecision]), $MachinePrecision] + N[(t * N[(a * (-x) + N[(b * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(c, a \cdot j - z \cdot b, x \cdot \mathsf{fma}\left(t, -a, y \cdot z\right)\right)\\
\mathbf{if}\;c \leq -6.2 \cdot 10^{-132}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;c \leq 4.3 \cdot 10^{+77}:\\
\;\;\;\;\mathsf{fma}\left(j, \mathsf{fma}\left(a, c, i \cdot \left(-y\right)\right), t \cdot \mathsf{fma}\left(a, -x, b \cdot i\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if c < -6.20000000000000016e-132 or 4.29999999999999991e77 < c Initial program 64.5%
Taylor expanded in i around 0
sub-negN/A
+-commutativeN/A
associate-+l+N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
distribute-lft-neg-inN/A
mul-1-negN/A
distribute-rgt-inN/A
mul-1-negN/A
sub-negN/A
+-commutativeN/A
lower-fma.f64N/A
Applied rewrites77.2%
if -6.20000000000000016e-132 < c < 4.29999999999999991e77Initial program 79.2%
Taylor expanded in z around 0
+-commutativeN/A
associate--l+N/A
sub-negN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
mul-1-negN/A
remove-double-negN/A
associate-*r*N/A
distribute-rgt-inN/A
*-lft-identityN/A
metadata-evalN/A
cancel-sign-sub-invN/A
Applied rewrites71.5%
Final simplification74.6%
(FPCore (x y z t a b c i j)
:precision binary64
(if (<= i -2e+68)
(+ (* b (* t i)) (* j (- (* c a) (* y i))))
(if (<= i 2.4e+78)
(fma c (- (* a j) (* z b)) (* x (fma t (- a) (* y z))))
(* i (fma j (- y) (* t b))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double tmp;
if (i <= -2e+68) {
tmp = (b * (t * i)) + (j * ((c * a) - (y * i)));
} else if (i <= 2.4e+78) {
tmp = fma(c, ((a * j) - (z * b)), (x * fma(t, -a, (y * z))));
} else {
tmp = i * fma(j, -y, (t * b));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j) tmp = 0.0 if (i <= -2e+68) tmp = Float64(Float64(b * Float64(t * i)) + Float64(j * Float64(Float64(c * a) - Float64(y * i)))); elseif (i <= 2.4e+78) tmp = fma(c, Float64(Float64(a * j) - Float64(z * b)), Float64(x * fma(t, Float64(-a), Float64(y * z)))); else tmp = Float64(i * fma(j, Float64(-y), Float64(t * b))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := If[LessEqual[i, -2e+68], N[(N[(b * N[(t * i), $MachinePrecision]), $MachinePrecision] + N[(j * N[(N[(c * a), $MachinePrecision] - N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[i, 2.4e+78], N[(c * N[(N[(a * j), $MachinePrecision] - N[(z * b), $MachinePrecision]), $MachinePrecision] + N[(x * N[(t * (-a) + N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(i * N[(j * (-y) + N[(t * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -2 \cdot 10^{+68}:\\
\;\;\;\;b \cdot \left(t \cdot i\right) + j \cdot \left(c \cdot a - y \cdot i\right)\\
\mathbf{elif}\;i \leq 2.4 \cdot 10^{+78}:\\
\;\;\;\;\mathsf{fma}\left(c, a \cdot j - z \cdot b, x \cdot \mathsf{fma}\left(t, -a, y \cdot z\right)\right)\\
\mathbf{else}:\\
\;\;\;\;i \cdot \mathsf{fma}\left(j, -y, t \cdot b\right)\\
\end{array}
\end{array}
if i < -1.99999999999999991e68Initial program 72.2%
lift-*.f64N/A
lift--.f64N/A
flip--N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
clear-numN/A
flip--N/A
lift--.f64N/A
lower-/.f6472.2
Applied rewrites72.2%
Taylor expanded in i around inf
lower-*.f64N/A
*-commutativeN/A
lower-*.f6474.7
Applied rewrites74.7%
if -1.99999999999999991e68 < i < 2.3999999999999999e78Initial program 73.8%
Taylor expanded in i around 0
sub-negN/A
+-commutativeN/A
associate-+l+N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
distribute-lft-neg-inN/A
mul-1-negN/A
distribute-rgt-inN/A
mul-1-negN/A
sub-negN/A
+-commutativeN/A
lower-fma.f64N/A
Applied rewrites73.3%
if 2.3999999999999999e78 < i Initial program 63.8%
Taylor expanded in i around inf
lower-*.f64N/A
cancel-sign-sub-invN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
metadata-evalN/A
*-lft-identityN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-*.f6472.3
Applied rewrites72.3%
Final simplification73.4%
(FPCore (x y z t a b c i j)
:precision binary64
(let* ((t_1 (* z (- (* c b)))))
(if (<= c -6.6e+103)
t_1
(if (<= c -5.2e-204)
(* z (* x y))
(if (<= c 6.4e-199)
(* b (* t i))
(if (<= c 4e-38)
(* a (- (* x t)))
(if (<= c 5.4e+70)
(* j (* i (- y)))
(if (<= c 1.22e+268) t_1 (* j (* c a))))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = z * -(c * b);
double tmp;
if (c <= -6.6e+103) {
tmp = t_1;
} else if (c <= -5.2e-204) {
tmp = z * (x * y);
} else if (c <= 6.4e-199) {
tmp = b * (t * i);
} else if (c <= 4e-38) {
tmp = a * -(x * t);
} else if (c <= 5.4e+70) {
tmp = j * (i * -y);
} else if (c <= 1.22e+268) {
tmp = t_1;
} else {
tmp = j * (c * a);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i, j)
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) :: t_1
real(8) :: tmp
t_1 = z * -(c * b)
if (c <= (-6.6d+103)) then
tmp = t_1
else if (c <= (-5.2d-204)) then
tmp = z * (x * y)
else if (c <= 6.4d-199) then
tmp = b * (t * i)
else if (c <= 4d-38) then
tmp = a * -(x * t)
else if (c <= 5.4d+70) then
tmp = j * (i * -y)
else if (c <= 1.22d+268) then
tmp = t_1
else
tmp = j * (c * a)
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 t_1 = z * -(c * b);
double tmp;
if (c <= -6.6e+103) {
tmp = t_1;
} else if (c <= -5.2e-204) {
tmp = z * (x * y);
} else if (c <= 6.4e-199) {
tmp = b * (t * i);
} else if (c <= 4e-38) {
tmp = a * -(x * t);
} else if (c <= 5.4e+70) {
tmp = j * (i * -y);
} else if (c <= 1.22e+268) {
tmp = t_1;
} else {
tmp = j * (c * a);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i, j): t_1 = z * -(c * b) tmp = 0 if c <= -6.6e+103: tmp = t_1 elif c <= -5.2e-204: tmp = z * (x * y) elif c <= 6.4e-199: tmp = b * (t * i) elif c <= 4e-38: tmp = a * -(x * t) elif c <= 5.4e+70: tmp = j * (i * -y) elif c <= 1.22e+268: tmp = t_1 else: tmp = j * (c * a) return tmp
function code(x, y, z, t, a, b, c, i, j) t_1 = Float64(z * Float64(-Float64(c * b))) tmp = 0.0 if (c <= -6.6e+103) tmp = t_1; elseif (c <= -5.2e-204) tmp = Float64(z * Float64(x * y)); elseif (c <= 6.4e-199) tmp = Float64(b * Float64(t * i)); elseif (c <= 4e-38) tmp = Float64(a * Float64(-Float64(x * t))); elseif (c <= 5.4e+70) tmp = Float64(j * Float64(i * Float64(-y))); elseif (c <= 1.22e+268) tmp = t_1; else tmp = Float64(j * Float64(c * a)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i, j) t_1 = z * -(c * b); tmp = 0.0; if (c <= -6.6e+103) tmp = t_1; elseif (c <= -5.2e-204) tmp = z * (x * y); elseif (c <= 6.4e-199) tmp = b * (t * i); elseif (c <= 4e-38) tmp = a * -(x * t); elseif (c <= 5.4e+70) tmp = j * (i * -y); elseif (c <= 1.22e+268) tmp = t_1; else tmp = j * (c * a); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(z * (-N[(c * b), $MachinePrecision])), $MachinePrecision]}, If[LessEqual[c, -6.6e+103], t$95$1, If[LessEqual[c, -5.2e-204], N[(z * N[(x * y), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 6.4e-199], N[(b * N[(t * i), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 4e-38], N[(a * (-N[(x * t), $MachinePrecision])), $MachinePrecision], If[LessEqual[c, 5.4e+70], N[(j * N[(i * (-y)), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 1.22e+268], t$95$1, N[(j * N[(c * a), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := z \cdot \left(-c \cdot b\right)\\
\mathbf{if}\;c \leq -6.6 \cdot 10^{+103}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;c \leq -5.2 \cdot 10^{-204}:\\
\;\;\;\;z \cdot \left(x \cdot y\right)\\
\mathbf{elif}\;c \leq 6.4 \cdot 10^{-199}:\\
\;\;\;\;b \cdot \left(t \cdot i\right)\\
\mathbf{elif}\;c \leq 4 \cdot 10^{-38}:\\
\;\;\;\;a \cdot \left(-x \cdot t\right)\\
\mathbf{elif}\;c \leq 5.4 \cdot 10^{+70}:\\
\;\;\;\;j \cdot \left(i \cdot \left(-y\right)\right)\\
\mathbf{elif}\;c \leq 1.22 \cdot 10^{+268}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;j \cdot \left(c \cdot a\right)\\
\end{array}
\end{array}
if c < -6.60000000000000017e103 or 5.3999999999999999e70 < c < 1.2200000000000001e268Initial program 53.8%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6460.6
Applied rewrites60.6%
Taylor expanded in c around inf
Applied rewrites50.5%
if -6.60000000000000017e103 < c < -5.19999999999999965e-204Initial program 79.8%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6446.1
Applied rewrites46.1%
Taylor expanded in c around 0
Applied rewrites38.8%
if -5.19999999999999965e-204 < c < 6.3999999999999999e-199Initial program 84.4%
Taylor expanded in t around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
mul-1-negN/A
remove-double-negN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6446.2
Applied rewrites46.2%
Taylor expanded in a around 0
Applied rewrites36.4%
if 6.3999999999999999e-199 < c < 3.9999999999999998e-38Initial program 78.3%
Taylor expanded in a around inf
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6454.5
Applied rewrites54.5%
Taylor expanded in j around 0
Applied rewrites45.7%
if 3.9999999999999998e-38 < c < 5.3999999999999999e70Initial program 65.0%
Taylor expanded in j around inf
lower-*.f64N/A
sub-negN/A
lower-fma.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6465.3
Applied rewrites65.3%
Taylor expanded in a around 0
Applied rewrites54.3%
if 1.2200000000000001e268 < c Initial program 75.0%
Taylor expanded in j around inf
lower-*.f64N/A
sub-negN/A
lower-fma.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6487.3
Applied rewrites87.3%
Taylor expanded in a around inf
Applied rewrites99.8%
Final simplification45.7%
(FPCore (x y z t a b c i j)
:precision binary64
(let* ((t_1 (* j (* i (- y)))))
(if (<= i -2.5e+24)
t_1
(if (<= i -1.85e-124)
(* y (* x z))
(if (<= i 3.15e-86)
(* b (* c (- z)))
(if (<= i 5.8e+78)
(* t (- (* x a)))
(if (<= i 1.35e+274) (* b (* t i)) t_1)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = j * (i * -y);
double tmp;
if (i <= -2.5e+24) {
tmp = t_1;
} else if (i <= -1.85e-124) {
tmp = y * (x * z);
} else if (i <= 3.15e-86) {
tmp = b * (c * -z);
} else if (i <= 5.8e+78) {
tmp = t * -(x * a);
} else if (i <= 1.35e+274) {
tmp = b * (t * i);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i, j)
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) :: t_1
real(8) :: tmp
t_1 = j * (i * -y)
if (i <= (-2.5d+24)) then
tmp = t_1
else if (i <= (-1.85d-124)) then
tmp = y * (x * z)
else if (i <= 3.15d-86) then
tmp = b * (c * -z)
else if (i <= 5.8d+78) then
tmp = t * -(x * a)
else if (i <= 1.35d+274) then
tmp = b * (t * i)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = j * (i * -y);
double tmp;
if (i <= -2.5e+24) {
tmp = t_1;
} else if (i <= -1.85e-124) {
tmp = y * (x * z);
} else if (i <= 3.15e-86) {
tmp = b * (c * -z);
} else if (i <= 5.8e+78) {
tmp = t * -(x * a);
} else if (i <= 1.35e+274) {
tmp = b * (t * i);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i, j): t_1 = j * (i * -y) tmp = 0 if i <= -2.5e+24: tmp = t_1 elif i <= -1.85e-124: tmp = y * (x * z) elif i <= 3.15e-86: tmp = b * (c * -z) elif i <= 5.8e+78: tmp = t * -(x * a) elif i <= 1.35e+274: tmp = b * (t * i) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i, j) t_1 = Float64(j * Float64(i * Float64(-y))) tmp = 0.0 if (i <= -2.5e+24) tmp = t_1; elseif (i <= -1.85e-124) tmp = Float64(y * Float64(x * z)); elseif (i <= 3.15e-86) tmp = Float64(b * Float64(c * Float64(-z))); elseif (i <= 5.8e+78) tmp = Float64(t * Float64(-Float64(x * a))); elseif (i <= 1.35e+274) tmp = Float64(b * Float64(t * i)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i, j) t_1 = j * (i * -y); tmp = 0.0; if (i <= -2.5e+24) tmp = t_1; elseif (i <= -1.85e-124) tmp = y * (x * z); elseif (i <= 3.15e-86) tmp = b * (c * -z); elseif (i <= 5.8e+78) tmp = t * -(x * a); elseif (i <= 1.35e+274) tmp = b * (t * i); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(j * N[(i * (-y)), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[i, -2.5e+24], t$95$1, If[LessEqual[i, -1.85e-124], N[(y * N[(x * z), $MachinePrecision]), $MachinePrecision], If[LessEqual[i, 3.15e-86], N[(b * N[(c * (-z)), $MachinePrecision]), $MachinePrecision], If[LessEqual[i, 5.8e+78], N[(t * (-N[(x * a), $MachinePrecision])), $MachinePrecision], If[LessEqual[i, 1.35e+274], N[(b * N[(t * i), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := j \cdot \left(i \cdot \left(-y\right)\right)\\
\mathbf{if}\;i \leq -2.5 \cdot 10^{+24}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;i \leq -1.85 \cdot 10^{-124}:\\
\;\;\;\;y \cdot \left(x \cdot z\right)\\
\mathbf{elif}\;i \leq 3.15 \cdot 10^{-86}:\\
\;\;\;\;b \cdot \left(c \cdot \left(-z\right)\right)\\
\mathbf{elif}\;i \leq 5.8 \cdot 10^{+78}:\\
\;\;\;\;t \cdot \left(-x \cdot a\right)\\
\mathbf{elif}\;i \leq 1.35 \cdot 10^{+274}:\\
\;\;\;\;b \cdot \left(t \cdot i\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if i < -2.50000000000000023e24 or 1.3500000000000001e274 < i Initial program 66.9%
Taylor expanded in j around inf
lower-*.f64N/A
sub-negN/A
lower-fma.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6455.2
Applied rewrites55.2%
Taylor expanded in a around 0
Applied rewrites46.2%
if -2.50000000000000023e24 < i < -1.84999999999999995e-124Initial program 74.9%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6453.7
Applied rewrites53.7%
Taylor expanded in c around 0
Applied rewrites35.5%
Applied rewrites46.0%
if -1.84999999999999995e-124 < i < 3.15e-86Initial program 75.3%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6452.8
Applied rewrites52.8%
Taylor expanded in c around inf
Applied rewrites35.2%
if 3.15e-86 < i < 5.80000000000000034e78Initial program 76.3%
Taylor expanded in t around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
mul-1-negN/A
remove-double-negN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6461.6
Applied rewrites61.6%
Taylor expanded in a around inf
Applied rewrites47.5%
if 5.80000000000000034e78 < i < 1.3500000000000001e274Initial program 64.1%
Taylor expanded in t around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
mul-1-negN/A
remove-double-negN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6461.3
Applied rewrites61.3%
Taylor expanded in a around 0
Applied rewrites56.2%
Final simplification43.9%
(FPCore (x y z t a b c i j)
:precision binary64
(let* ((t_1 (* j (* i (- y)))))
(if (<= i -2.5e+24)
t_1
(if (<= i -1.85e-124)
(* y (* x z))
(if (<= i 3.15e-86)
(* b (* c (- z)))
(if (<= i 5.8e+78)
(* a (- (* x t)))
(if (<= i 1.35e+274) (* b (* t i)) t_1)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = j * (i * -y);
double tmp;
if (i <= -2.5e+24) {
tmp = t_1;
} else if (i <= -1.85e-124) {
tmp = y * (x * z);
} else if (i <= 3.15e-86) {
tmp = b * (c * -z);
} else if (i <= 5.8e+78) {
tmp = a * -(x * t);
} else if (i <= 1.35e+274) {
tmp = b * (t * i);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i, j)
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) :: t_1
real(8) :: tmp
t_1 = j * (i * -y)
if (i <= (-2.5d+24)) then
tmp = t_1
else if (i <= (-1.85d-124)) then
tmp = y * (x * z)
else if (i <= 3.15d-86) then
tmp = b * (c * -z)
else if (i <= 5.8d+78) then
tmp = a * -(x * t)
else if (i <= 1.35d+274) then
tmp = b * (t * i)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = j * (i * -y);
double tmp;
if (i <= -2.5e+24) {
tmp = t_1;
} else if (i <= -1.85e-124) {
tmp = y * (x * z);
} else if (i <= 3.15e-86) {
tmp = b * (c * -z);
} else if (i <= 5.8e+78) {
tmp = a * -(x * t);
} else if (i <= 1.35e+274) {
tmp = b * (t * i);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i, j): t_1 = j * (i * -y) tmp = 0 if i <= -2.5e+24: tmp = t_1 elif i <= -1.85e-124: tmp = y * (x * z) elif i <= 3.15e-86: tmp = b * (c * -z) elif i <= 5.8e+78: tmp = a * -(x * t) elif i <= 1.35e+274: tmp = b * (t * i) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i, j) t_1 = Float64(j * Float64(i * Float64(-y))) tmp = 0.0 if (i <= -2.5e+24) tmp = t_1; elseif (i <= -1.85e-124) tmp = Float64(y * Float64(x * z)); elseif (i <= 3.15e-86) tmp = Float64(b * Float64(c * Float64(-z))); elseif (i <= 5.8e+78) tmp = Float64(a * Float64(-Float64(x * t))); elseif (i <= 1.35e+274) tmp = Float64(b * Float64(t * i)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i, j) t_1 = j * (i * -y); tmp = 0.0; if (i <= -2.5e+24) tmp = t_1; elseif (i <= -1.85e-124) tmp = y * (x * z); elseif (i <= 3.15e-86) tmp = b * (c * -z); elseif (i <= 5.8e+78) tmp = a * -(x * t); elseif (i <= 1.35e+274) tmp = b * (t * i); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(j * N[(i * (-y)), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[i, -2.5e+24], t$95$1, If[LessEqual[i, -1.85e-124], N[(y * N[(x * z), $MachinePrecision]), $MachinePrecision], If[LessEqual[i, 3.15e-86], N[(b * N[(c * (-z)), $MachinePrecision]), $MachinePrecision], If[LessEqual[i, 5.8e+78], N[(a * (-N[(x * t), $MachinePrecision])), $MachinePrecision], If[LessEqual[i, 1.35e+274], N[(b * N[(t * i), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := j \cdot \left(i \cdot \left(-y\right)\right)\\
\mathbf{if}\;i \leq -2.5 \cdot 10^{+24}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;i \leq -1.85 \cdot 10^{-124}:\\
\;\;\;\;y \cdot \left(x \cdot z\right)\\
\mathbf{elif}\;i \leq 3.15 \cdot 10^{-86}:\\
\;\;\;\;b \cdot \left(c \cdot \left(-z\right)\right)\\
\mathbf{elif}\;i \leq 5.8 \cdot 10^{+78}:\\
\;\;\;\;a \cdot \left(-x \cdot t\right)\\
\mathbf{elif}\;i \leq 1.35 \cdot 10^{+274}:\\
\;\;\;\;b \cdot \left(t \cdot i\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if i < -2.50000000000000023e24 or 1.3500000000000001e274 < i Initial program 66.9%
Taylor expanded in j around inf
lower-*.f64N/A
sub-negN/A
lower-fma.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6455.2
Applied rewrites55.2%
Taylor expanded in a around 0
Applied rewrites46.2%
if -2.50000000000000023e24 < i < -1.84999999999999995e-124Initial program 74.9%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6453.7
Applied rewrites53.7%
Taylor expanded in c around 0
Applied rewrites35.5%
Applied rewrites46.0%
if -1.84999999999999995e-124 < i < 3.15e-86Initial program 75.3%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6452.8
Applied rewrites52.8%
Taylor expanded in c around inf
Applied rewrites35.2%
if 3.15e-86 < i < 5.80000000000000034e78Initial program 76.3%
Taylor expanded in a around inf
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6454.4
Applied rewrites54.4%
Taylor expanded in j around 0
Applied rewrites47.4%
if 5.80000000000000034e78 < i < 1.3500000000000001e274Initial program 64.1%
Taylor expanded in t around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
mul-1-negN/A
remove-double-negN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6461.3
Applied rewrites61.3%
Taylor expanded in a around 0
Applied rewrites56.2%
Final simplification43.9%
(FPCore (x y z t a b c i j)
:precision binary64
(if (<= b -3.9e+17)
(* b (fma c (- z) (* t i)))
(if (<= b 6.4e+137)
(fma a (fma j c (- (* x t))) (* x (* y z)))
(* b (- (* t i) (* c z))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double tmp;
if (b <= -3.9e+17) {
tmp = b * fma(c, -z, (t * i));
} else if (b <= 6.4e+137) {
tmp = fma(a, fma(j, c, -(x * t)), (x * (y * z)));
} else {
tmp = b * ((t * i) - (c * z));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j) tmp = 0.0 if (b <= -3.9e+17) tmp = Float64(b * fma(c, Float64(-z), Float64(t * i))); elseif (b <= 6.4e+137) tmp = fma(a, fma(j, c, Float64(-Float64(x * t))), Float64(x * Float64(y * z))); else tmp = Float64(b * Float64(Float64(t * i) - Float64(c * z))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := If[LessEqual[b, -3.9e+17], N[(b * N[(c * (-z) + N[(t * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 6.4e+137], N[(a * N[(j * c + (-N[(x * t), $MachinePrecision])), $MachinePrecision] + N[(x * N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(b * N[(N[(t * i), $MachinePrecision] - N[(c * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -3.9 \cdot 10^{+17}:\\
\;\;\;\;b \cdot \mathsf{fma}\left(c, -z, t \cdot i\right)\\
\mathbf{elif}\;b \leq 6.4 \cdot 10^{+137}:\\
\;\;\;\;\mathsf{fma}\left(a, \mathsf{fma}\left(j, c, -x \cdot t\right), x \cdot \left(y \cdot z\right)\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot \left(t \cdot i - c \cdot z\right)\\
\end{array}
\end{array}
if b < -3.9e17Initial program 75.4%
Taylor expanded in y around 0
Applied rewrites75.4%
Taylor expanded in b around inf
lower-*.f64N/A
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-*.f6462.2
Applied rewrites62.2%
if -3.9e17 < b < 6.40000000000000038e137Initial program 68.3%
Taylor expanded in y around 0
Applied rewrites73.1%
Taylor expanded in x around inf
Applied rewrites62.8%
if 6.40000000000000038e137 < b Initial program 75.0%
Taylor expanded in b around inf
cancel-sign-sub-invN/A
remove-double-negN/A
distribute-lft-neg-inN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
sub-negN/A
distribute-neg-inN/A
remove-double-negN/A
+-commutativeN/A
sub-negN/A
lower--.f64N/A
lower-*.f64N/A
lower-*.f6481.6
Applied rewrites81.6%
Final simplification65.5%
(FPCore (x y z t a b c i j)
:precision binary64
(if (<= b -1.05e+27)
(* t (* b i))
(if (<= b -4.6e-57)
(* j (* c a))
(if (<= b 6.8e-273)
(* i (* y (- j)))
(if (<= b 2.3e-52)
(* a (- (* x t)))
(if (<= b 3.6e+148) (* x (* y z)) (* i (* t b))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double tmp;
if (b <= -1.05e+27) {
tmp = t * (b * i);
} else if (b <= -4.6e-57) {
tmp = j * (c * a);
} else if (b <= 6.8e-273) {
tmp = i * (y * -j);
} else if (b <= 2.3e-52) {
tmp = a * -(x * t);
} else if (b <= 3.6e+148) {
tmp = x * (y * z);
} else {
tmp = i * (t * b);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i, j)
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) :: tmp
if (b <= (-1.05d+27)) then
tmp = t * (b * i)
else if (b <= (-4.6d-57)) then
tmp = j * (c * a)
else if (b <= 6.8d-273) then
tmp = i * (y * -j)
else if (b <= 2.3d-52) then
tmp = a * -(x * t)
else if (b <= 3.6d+148) then
tmp = x * (y * z)
else
tmp = i * (t * b)
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 tmp;
if (b <= -1.05e+27) {
tmp = t * (b * i);
} else if (b <= -4.6e-57) {
tmp = j * (c * a);
} else if (b <= 6.8e-273) {
tmp = i * (y * -j);
} else if (b <= 2.3e-52) {
tmp = a * -(x * t);
} else if (b <= 3.6e+148) {
tmp = x * (y * z);
} else {
tmp = i * (t * b);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i, j): tmp = 0 if b <= -1.05e+27: tmp = t * (b * i) elif b <= -4.6e-57: tmp = j * (c * a) elif b <= 6.8e-273: tmp = i * (y * -j) elif b <= 2.3e-52: tmp = a * -(x * t) elif b <= 3.6e+148: tmp = x * (y * z) else: tmp = i * (t * b) return tmp
function code(x, y, z, t, a, b, c, i, j) tmp = 0.0 if (b <= -1.05e+27) tmp = Float64(t * Float64(b * i)); elseif (b <= -4.6e-57) tmp = Float64(j * Float64(c * a)); elseif (b <= 6.8e-273) tmp = Float64(i * Float64(y * Float64(-j))); elseif (b <= 2.3e-52) tmp = Float64(a * Float64(-Float64(x * t))); elseif (b <= 3.6e+148) tmp = Float64(x * Float64(y * z)); else tmp = Float64(i * Float64(t * b)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i, j) tmp = 0.0; if (b <= -1.05e+27) tmp = t * (b * i); elseif (b <= -4.6e-57) tmp = j * (c * a); elseif (b <= 6.8e-273) tmp = i * (y * -j); elseif (b <= 2.3e-52) tmp = a * -(x * t); elseif (b <= 3.6e+148) tmp = x * (y * z); else tmp = i * (t * b); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := If[LessEqual[b, -1.05e+27], N[(t * N[(b * i), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -4.6e-57], N[(j * N[(c * a), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 6.8e-273], N[(i * N[(y * (-j)), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.3e-52], N[(a * (-N[(x * t), $MachinePrecision])), $MachinePrecision], If[LessEqual[b, 3.6e+148], N[(x * N[(y * z), $MachinePrecision]), $MachinePrecision], N[(i * N[(t * b), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.05 \cdot 10^{+27}:\\
\;\;\;\;t \cdot \left(b \cdot i\right)\\
\mathbf{elif}\;b \leq -4.6 \cdot 10^{-57}:\\
\;\;\;\;j \cdot \left(c \cdot a\right)\\
\mathbf{elif}\;b \leq 6.8 \cdot 10^{-273}:\\
\;\;\;\;i \cdot \left(y \cdot \left(-j\right)\right)\\
\mathbf{elif}\;b \leq 2.3 \cdot 10^{-52}:\\
\;\;\;\;a \cdot \left(-x \cdot t\right)\\
\mathbf{elif}\;b \leq 3.6 \cdot 10^{+148}:\\
\;\;\;\;x \cdot \left(y \cdot z\right)\\
\mathbf{else}:\\
\;\;\;\;i \cdot \left(t \cdot b\right)\\
\end{array}
\end{array}
if b < -1.04999999999999997e27Initial program 75.1%
Taylor expanded in t around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
mul-1-negN/A
remove-double-negN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6448.3
Applied rewrites48.3%
Taylor expanded in a around 0
Applied rewrites41.8%
if -1.04999999999999997e27 < b < -4.6e-57Initial program 71.4%
Taylor expanded in j around inf
lower-*.f64N/A
sub-negN/A
lower-fma.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6443.9
Applied rewrites43.9%
Taylor expanded in a around inf
Applied rewrites34.4%
if -4.6e-57 < b < 6.79999999999999982e-273Initial program 69.7%
Taylor expanded in j around inf
lower-*.f64N/A
sub-negN/A
lower-fma.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6461.6
Applied rewrites61.6%
Taylor expanded in a around 0
Applied rewrites45.3%
if 6.79999999999999982e-273 < b < 2.29999999999999994e-52Initial program 65.1%
Taylor expanded in a around inf
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6452.8
Applied rewrites52.8%
Taylor expanded in j around 0
Applied rewrites41.0%
if 2.29999999999999994e-52 < b < 3.60000000000000006e148Initial program 75.2%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6463.2
Applied rewrites63.2%
Taylor expanded in c around 0
Applied rewrites41.5%
if 3.60000000000000006e148 < b Initial program 71.4%
Taylor expanded in i around -inf
Applied rewrites66.6%
Taylor expanded in b around inf
Applied rewrites83.8%
Taylor expanded in t around inf
Applied rewrites55.5%
Final simplification43.6%
(FPCore (x y z t a b c i j)
:precision binary64
(let* ((t_1 (* c (- (* a j) (* z b)))))
(if (<= c -4.1e+101)
t_1
(if (<= c -1.25e-203)
(* y (fma j (- i) (* x z)))
(if (<= c 3.8e+109) (* i (fma j (- y) (* t b))) t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = c * ((a * j) - (z * b));
double tmp;
if (c <= -4.1e+101) {
tmp = t_1;
} else if (c <= -1.25e-203) {
tmp = y * fma(j, -i, (x * z));
} else if (c <= 3.8e+109) {
tmp = i * fma(j, -y, (t * b));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j) t_1 = Float64(c * Float64(Float64(a * j) - Float64(z * b))) tmp = 0.0 if (c <= -4.1e+101) tmp = t_1; elseif (c <= -1.25e-203) tmp = Float64(y * fma(j, Float64(-i), Float64(x * z))); elseif (c <= 3.8e+109) tmp = Float64(i * fma(j, Float64(-y), Float64(t * b))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(c * N[(N[(a * j), $MachinePrecision] - N[(z * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[c, -4.1e+101], t$95$1, If[LessEqual[c, -1.25e-203], N[(y * N[(j * (-i) + N[(x * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 3.8e+109], N[(i * N[(j * (-y) + N[(t * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := c \cdot \left(a \cdot j - z \cdot b\right)\\
\mathbf{if}\;c \leq -4.1 \cdot 10^{+101}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;c \leq -1.25 \cdot 10^{-203}:\\
\;\;\;\;y \cdot \mathsf{fma}\left(j, -i, x \cdot z\right)\\
\mathbf{elif}\;c \leq 3.8 \cdot 10^{+109}:\\
\;\;\;\;i \cdot \mathsf{fma}\left(j, -y, t \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if c < -4.1e101 or 3.80000000000000039e109 < c Initial program 55.7%
Taylor expanded in c around inf
lower-*.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6476.7
Applied rewrites76.7%
if -4.1e101 < c < -1.25e-203Initial program 79.8%
Taylor expanded in y around inf
lower-*.f64N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6453.4
Applied rewrites53.4%
if -1.25e-203 < c < 3.80000000000000039e109Initial program 77.3%
Taylor expanded in i around inf
lower-*.f64N/A
cancel-sign-sub-invN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
metadata-evalN/A
*-lft-identityN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-*.f6457.1
Applied rewrites57.1%
Final simplification62.1%
(FPCore (x y z t a b c i j)
:precision binary64
(let* ((t_1 (* c (- (* a j) (* z b)))))
(if (<= c -2.5e+102)
t_1
(if (<= c -7.5e-206)
(* x (fma t (- a) (* y z)))
(if (<= c 3.8e+109) (* i (fma j (- y) (* t b))) t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = c * ((a * j) - (z * b));
double tmp;
if (c <= -2.5e+102) {
tmp = t_1;
} else if (c <= -7.5e-206) {
tmp = x * fma(t, -a, (y * z));
} else if (c <= 3.8e+109) {
tmp = i * fma(j, -y, (t * b));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j) t_1 = Float64(c * Float64(Float64(a * j) - Float64(z * b))) tmp = 0.0 if (c <= -2.5e+102) tmp = t_1; elseif (c <= -7.5e-206) tmp = Float64(x * fma(t, Float64(-a), Float64(y * z))); elseif (c <= 3.8e+109) tmp = Float64(i * fma(j, Float64(-y), Float64(t * b))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(c * N[(N[(a * j), $MachinePrecision] - N[(z * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[c, -2.5e+102], t$95$1, If[LessEqual[c, -7.5e-206], N[(x * N[(t * (-a) + N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 3.8e+109], N[(i * N[(j * (-y) + N[(t * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := c \cdot \left(a \cdot j - z \cdot b\right)\\
\mathbf{if}\;c \leq -2.5 \cdot 10^{+102}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;c \leq -7.5 \cdot 10^{-206}:\\
\;\;\;\;x \cdot \mathsf{fma}\left(t, -a, y \cdot z\right)\\
\mathbf{elif}\;c \leq 3.8 \cdot 10^{+109}:\\
\;\;\;\;i \cdot \mathsf{fma}\left(j, -y, t \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if c < -2.5e102 or 3.80000000000000039e109 < c Initial program 55.7%
Taylor expanded in c around inf
lower-*.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6476.7
Applied rewrites76.7%
if -2.5e102 < c < -7.5e-206Initial program 80.1%
Taylor expanded in x around inf
lower-*.f64N/A
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-*.f6451.9
Applied rewrites51.9%
if -7.5e-206 < c < 3.80000000000000039e109Initial program 77.0%
Taylor expanded in i around inf
lower-*.f64N/A
cancel-sign-sub-invN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
metadata-evalN/A
*-lft-identityN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-*.f6457.6
Applied rewrites57.6%
Final simplification61.9%
(FPCore (x y z t a b c i j)
:precision binary64
(if (<= b -3.4e+17)
(* b (fma c (- z) (* t i)))
(if (<= b 3.2e-238)
(* j (fma a c (* i (- y))))
(if (<= b 1.9e+52)
(* a (fma j c (- (* x t))))
(* b (- (* t i) (* c z)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double tmp;
if (b <= -3.4e+17) {
tmp = b * fma(c, -z, (t * i));
} else if (b <= 3.2e-238) {
tmp = j * fma(a, c, (i * -y));
} else if (b <= 1.9e+52) {
tmp = a * fma(j, c, -(x * t));
} else {
tmp = b * ((t * i) - (c * z));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j) tmp = 0.0 if (b <= -3.4e+17) tmp = Float64(b * fma(c, Float64(-z), Float64(t * i))); elseif (b <= 3.2e-238) tmp = Float64(j * fma(a, c, Float64(i * Float64(-y)))); elseif (b <= 1.9e+52) tmp = Float64(a * fma(j, c, Float64(-Float64(x * t)))); else tmp = Float64(b * Float64(Float64(t * i) - Float64(c * z))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := If[LessEqual[b, -3.4e+17], N[(b * N[(c * (-z) + N[(t * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 3.2e-238], N[(j * N[(a * c + N[(i * (-y)), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.9e+52], N[(a * N[(j * c + (-N[(x * t), $MachinePrecision])), $MachinePrecision]), $MachinePrecision], N[(b * N[(N[(t * i), $MachinePrecision] - N[(c * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -3.4 \cdot 10^{+17}:\\
\;\;\;\;b \cdot \mathsf{fma}\left(c, -z, t \cdot i\right)\\
\mathbf{elif}\;b \leq 3.2 \cdot 10^{-238}:\\
\;\;\;\;j \cdot \mathsf{fma}\left(a, c, i \cdot \left(-y\right)\right)\\
\mathbf{elif}\;b \leq 1.9 \cdot 10^{+52}:\\
\;\;\;\;a \cdot \mathsf{fma}\left(j, c, -x \cdot t\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot \left(t \cdot i - c \cdot z\right)\\
\end{array}
\end{array}
if b < -3.4e17Initial program 75.4%
Taylor expanded in y around 0
Applied rewrites75.4%
Taylor expanded in b around inf
lower-*.f64N/A
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-*.f6462.2
Applied rewrites62.2%
if -3.4e17 < b < 3.2000000000000002e-238Initial program 69.2%
Taylor expanded in j around inf
lower-*.f64N/A
sub-negN/A
lower-fma.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6457.9
Applied rewrites57.9%
if 3.2000000000000002e-238 < b < 1.9e52Initial program 64.4%
Taylor expanded in a around inf
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6452.6
Applied rewrites52.6%
if 1.9e52 < b Initial program 75.5%
Taylor expanded in b around inf
cancel-sign-sub-invN/A
remove-double-negN/A
distribute-lft-neg-inN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
sub-negN/A
distribute-neg-inN/A
remove-double-negN/A
+-commutativeN/A
sub-negN/A
lower--.f64N/A
lower-*.f64N/A
lower-*.f6469.6
Applied rewrites69.6%
Final simplification60.6%
(FPCore (x y z t a b c i j)
:precision binary64
(let* ((t_1 (* c (- (* a j) (* z b)))))
(if (<= c -2.6e+100)
t_1
(if (<= c -4.2e+85)
(* z (* x y))
(if (<= c 3.8e+109) (* i (fma j (- y) (* t b))) t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = c * ((a * j) - (z * b));
double tmp;
if (c <= -2.6e+100) {
tmp = t_1;
} else if (c <= -4.2e+85) {
tmp = z * (x * y);
} else if (c <= 3.8e+109) {
tmp = i * fma(j, -y, (t * b));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j) t_1 = Float64(c * Float64(Float64(a * j) - Float64(z * b))) tmp = 0.0 if (c <= -2.6e+100) tmp = t_1; elseif (c <= -4.2e+85) tmp = Float64(z * Float64(x * y)); elseif (c <= 3.8e+109) tmp = Float64(i * fma(j, Float64(-y), Float64(t * b))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(c * N[(N[(a * j), $MachinePrecision] - N[(z * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[c, -2.6e+100], t$95$1, If[LessEqual[c, -4.2e+85], N[(z * N[(x * y), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 3.8e+109], N[(i * N[(j * (-y) + N[(t * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := c \cdot \left(a \cdot j - z \cdot b\right)\\
\mathbf{if}\;c \leq -2.6 \cdot 10^{+100}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;c \leq -4.2 \cdot 10^{+85}:\\
\;\;\;\;z \cdot \left(x \cdot y\right)\\
\mathbf{elif}\;c \leq 3.8 \cdot 10^{+109}:\\
\;\;\;\;i \cdot \mathsf{fma}\left(j, -y, t \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if c < -2.6000000000000002e100 or 3.80000000000000039e109 < c Initial program 55.7%
Taylor expanded in c around inf
lower-*.f64N/A
lower--.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6476.7
Applied rewrites76.7%
if -2.6000000000000002e100 < c < -4.2000000000000002e85Initial program 83.1%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6485.5
Applied rewrites85.5%
Taylor expanded in c around 0
Applied rewrites85.5%
if -4.2000000000000002e85 < c < 3.80000000000000039e109Initial program 78.1%
Taylor expanded in i around inf
lower-*.f64N/A
cancel-sign-sub-invN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
metadata-evalN/A
*-lft-identityN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-*.f6452.3
Applied rewrites52.3%
Final simplification60.5%
(FPCore (x y z t a b c i j)
:precision binary64
(let* ((t_1 (* a (fma j c (- (* x t))))))
(if (<= a -1.25e+25)
t_1
(if (<= a 2.5e-206)
(* z (* x y))
(if (<= a 0.0037) (* b (* c (- z))) t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = a * fma(j, c, -(x * t));
double tmp;
if (a <= -1.25e+25) {
tmp = t_1;
} else if (a <= 2.5e-206) {
tmp = z * (x * y);
} else if (a <= 0.0037) {
tmp = b * (c * -z);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j) t_1 = Float64(a * fma(j, c, Float64(-Float64(x * t)))) tmp = 0.0 if (a <= -1.25e+25) tmp = t_1; elseif (a <= 2.5e-206) tmp = Float64(z * Float64(x * y)); elseif (a <= 0.0037) tmp = Float64(b * Float64(c * Float64(-z))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(a * N[(j * c + (-N[(x * t), $MachinePrecision])), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -1.25e+25], t$95$1, If[LessEqual[a, 2.5e-206], N[(z * N[(x * y), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 0.0037], N[(b * N[(c * (-z)), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := a \cdot \mathsf{fma}\left(j, c, -x \cdot t\right)\\
\mathbf{if}\;a \leq -1.25 \cdot 10^{+25}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 2.5 \cdot 10^{-206}:\\
\;\;\;\;z \cdot \left(x \cdot y\right)\\
\mathbf{elif}\;a \leq 0.0037:\\
\;\;\;\;b \cdot \left(c \cdot \left(-z\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -1.25000000000000006e25 or 0.0037000000000000002 < a Initial program 65.7%
Taylor expanded in a around inf
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6460.4
Applied rewrites60.4%
if -1.25000000000000006e25 < a < 2.5e-206Initial program 78.5%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6449.3
Applied rewrites49.3%
Taylor expanded in c around 0
Applied rewrites34.5%
if 2.5e-206 < a < 0.0037000000000000002Initial program 72.0%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6439.2
Applied rewrites39.2%
Taylor expanded in c around inf
Applied rewrites35.1%
Final simplification46.8%
(FPCore (x y z t a b c i j)
:precision binary64
(if (<= b -1.05e+27)
(* t (* b i))
(if (<= b -4.6e-57)
(* j (* c a))
(if (<= b 1.08e-252)
(* i (* y (- j)))
(if (<= b 2.2e+133) (* y (* x z)) (* i (* t b)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double tmp;
if (b <= -1.05e+27) {
tmp = t * (b * i);
} else if (b <= -4.6e-57) {
tmp = j * (c * a);
} else if (b <= 1.08e-252) {
tmp = i * (y * -j);
} else if (b <= 2.2e+133) {
tmp = y * (x * z);
} else {
tmp = i * (t * b);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i, j)
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) :: tmp
if (b <= (-1.05d+27)) then
tmp = t * (b * i)
else if (b <= (-4.6d-57)) then
tmp = j * (c * a)
else if (b <= 1.08d-252) then
tmp = i * (y * -j)
else if (b <= 2.2d+133) then
tmp = y * (x * z)
else
tmp = i * (t * b)
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 tmp;
if (b <= -1.05e+27) {
tmp = t * (b * i);
} else if (b <= -4.6e-57) {
tmp = j * (c * a);
} else if (b <= 1.08e-252) {
tmp = i * (y * -j);
} else if (b <= 2.2e+133) {
tmp = y * (x * z);
} else {
tmp = i * (t * b);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i, j): tmp = 0 if b <= -1.05e+27: tmp = t * (b * i) elif b <= -4.6e-57: tmp = j * (c * a) elif b <= 1.08e-252: tmp = i * (y * -j) elif b <= 2.2e+133: tmp = y * (x * z) else: tmp = i * (t * b) return tmp
function code(x, y, z, t, a, b, c, i, j) tmp = 0.0 if (b <= -1.05e+27) tmp = Float64(t * Float64(b * i)); elseif (b <= -4.6e-57) tmp = Float64(j * Float64(c * a)); elseif (b <= 1.08e-252) tmp = Float64(i * Float64(y * Float64(-j))); elseif (b <= 2.2e+133) tmp = Float64(y * Float64(x * z)); else tmp = Float64(i * Float64(t * b)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i, j) tmp = 0.0; if (b <= -1.05e+27) tmp = t * (b * i); elseif (b <= -4.6e-57) tmp = j * (c * a); elseif (b <= 1.08e-252) tmp = i * (y * -j); elseif (b <= 2.2e+133) tmp = y * (x * z); else tmp = i * (t * b); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := If[LessEqual[b, -1.05e+27], N[(t * N[(b * i), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -4.6e-57], N[(j * N[(c * a), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.08e-252], N[(i * N[(y * (-j)), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.2e+133], N[(y * N[(x * z), $MachinePrecision]), $MachinePrecision], N[(i * N[(t * b), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.05 \cdot 10^{+27}:\\
\;\;\;\;t \cdot \left(b \cdot i\right)\\
\mathbf{elif}\;b \leq -4.6 \cdot 10^{-57}:\\
\;\;\;\;j \cdot \left(c \cdot a\right)\\
\mathbf{elif}\;b \leq 1.08 \cdot 10^{-252}:\\
\;\;\;\;i \cdot \left(y \cdot \left(-j\right)\right)\\
\mathbf{elif}\;b \leq 2.2 \cdot 10^{+133}:\\
\;\;\;\;y \cdot \left(x \cdot z\right)\\
\mathbf{else}:\\
\;\;\;\;i \cdot \left(t \cdot b\right)\\
\end{array}
\end{array}
if b < -1.04999999999999997e27Initial program 75.1%
Taylor expanded in t around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
mul-1-negN/A
remove-double-negN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6448.3
Applied rewrites48.3%
Taylor expanded in a around 0
Applied rewrites41.8%
if -1.04999999999999997e27 < b < -4.6e-57Initial program 71.4%
Taylor expanded in j around inf
lower-*.f64N/A
sub-negN/A
lower-fma.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6443.9
Applied rewrites43.9%
Taylor expanded in a around inf
Applied rewrites34.4%
if -4.6e-57 < b < 1.08e-252Initial program 70.2%
Taylor expanded in j around inf
lower-*.f64N/A
sub-negN/A
lower-fma.f64N/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6461.2
Applied rewrites61.2%
Taylor expanded in a around 0
Applied rewrites44.2%
if 1.08e-252 < b < 2.2e133Initial program 66.1%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6447.2
Applied rewrites47.2%
Taylor expanded in c around 0
Applied rewrites31.6%
Applied rewrites35.9%
if 2.2e133 < b Initial program 75.6%
Taylor expanded in i around -inf
Applied rewrites64.5%
Taylor expanded in b around inf
Applied rewrites76.8%
Taylor expanded in t around inf
Applied rewrites50.1%
Final simplification41.5%
(FPCore (x y z t a b c i j) :precision binary64 (if (<= b -3.9e+17) (* b (fma c (- z) (* t i))) (if (<= b 1.9e+52) (* a (fma j c (- (* x t)))) (* b (- (* t i) (* c z))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double tmp;
if (b <= -3.9e+17) {
tmp = b * fma(c, -z, (t * i));
} else if (b <= 1.9e+52) {
tmp = a * fma(j, c, -(x * t));
} else {
tmp = b * ((t * i) - (c * z));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j) tmp = 0.0 if (b <= -3.9e+17) tmp = Float64(b * fma(c, Float64(-z), Float64(t * i))); elseif (b <= 1.9e+52) tmp = Float64(a * fma(j, c, Float64(-Float64(x * t)))); else tmp = Float64(b * Float64(Float64(t * i) - Float64(c * z))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := If[LessEqual[b, -3.9e+17], N[(b * N[(c * (-z) + N[(t * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.9e+52], N[(a * N[(j * c + (-N[(x * t), $MachinePrecision])), $MachinePrecision]), $MachinePrecision], N[(b * N[(N[(t * i), $MachinePrecision] - N[(c * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -3.9 \cdot 10^{+17}:\\
\;\;\;\;b \cdot \mathsf{fma}\left(c, -z, t \cdot i\right)\\
\mathbf{elif}\;b \leq 1.9 \cdot 10^{+52}:\\
\;\;\;\;a \cdot \mathsf{fma}\left(j, c, -x \cdot t\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot \left(t \cdot i - c \cdot z\right)\\
\end{array}
\end{array}
if b < -3.9e17Initial program 75.4%
Taylor expanded in y around 0
Applied rewrites75.4%
Taylor expanded in b around inf
lower-*.f64N/A
sub-negN/A
+-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
lower-*.f6462.2
Applied rewrites62.2%
if -3.9e17 < b < 1.9e52Initial program 67.4%
Taylor expanded in a around inf
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6447.3
Applied rewrites47.3%
if 1.9e52 < b Initial program 75.5%
Taylor expanded in b around inf
cancel-sign-sub-invN/A
remove-double-negN/A
distribute-lft-neg-inN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
sub-negN/A
distribute-neg-inN/A
remove-double-negN/A
+-commutativeN/A
sub-negN/A
lower--.f64N/A
lower-*.f64N/A
lower-*.f6469.6
Applied rewrites69.6%
Final simplification56.2%
(FPCore (x y z t a b c i j)
:precision binary64
(let* ((t_1 (* b (- (* t i) (* c z)))))
(if (<= b -3.9e+17)
t_1
(if (<= b 1.9e+52) (* a (fma j c (- (* x t)))) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = b * ((t * i) - (c * z));
double tmp;
if (b <= -3.9e+17) {
tmp = t_1;
} else if (b <= 1.9e+52) {
tmp = a * fma(j, c, -(x * t));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j) t_1 = Float64(b * Float64(Float64(t * i) - Float64(c * z))) tmp = 0.0 if (b <= -3.9e+17) tmp = t_1; elseif (b <= 1.9e+52) tmp = Float64(a * fma(j, c, Float64(-Float64(x * t)))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(b * N[(N[(t * i), $MachinePrecision] - N[(c * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -3.9e+17], t$95$1, If[LessEqual[b, 1.9e+52], N[(a * N[(j * c + (-N[(x * t), $MachinePrecision])), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := b \cdot \left(t \cdot i - c \cdot z\right)\\
\mathbf{if}\;b \leq -3.9 \cdot 10^{+17}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 1.9 \cdot 10^{+52}:\\
\;\;\;\;a \cdot \mathsf{fma}\left(j, c, -x \cdot t\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -3.9e17 or 1.9e52 < b Initial program 75.4%
Taylor expanded in b around inf
cancel-sign-sub-invN/A
remove-double-negN/A
distribute-lft-neg-inN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
sub-negN/A
distribute-neg-inN/A
remove-double-negN/A
+-commutativeN/A
sub-negN/A
lower--.f64N/A
lower-*.f64N/A
lower-*.f6464.5
Applied rewrites64.5%
if -3.9e17 < b < 1.9e52Initial program 67.4%
Taylor expanded in a around inf
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6447.3
Applied rewrites47.3%
Final simplification55.8%
(FPCore (x y z t a b c i j) :precision binary64 (if (<= b -1.75e+41) (* t (* b i)) (if (<= b 2.2e+133) (* y (* x z)) (* i (* t b)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double tmp;
if (b <= -1.75e+41) {
tmp = t * (b * i);
} else if (b <= 2.2e+133) {
tmp = y * (x * z);
} else {
tmp = i * (t * b);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i, j)
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) :: tmp
if (b <= (-1.75d+41)) then
tmp = t * (b * i)
else if (b <= 2.2d+133) then
tmp = y * (x * z)
else
tmp = i * (t * b)
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 tmp;
if (b <= -1.75e+41) {
tmp = t * (b * i);
} else if (b <= 2.2e+133) {
tmp = y * (x * z);
} else {
tmp = i * (t * b);
}
return tmp;
}
def code(x, y, z, t, a, b, c, i, j): tmp = 0 if b <= -1.75e+41: tmp = t * (b * i) elif b <= 2.2e+133: tmp = y * (x * z) else: tmp = i * (t * b) return tmp
function code(x, y, z, t, a, b, c, i, j) tmp = 0.0 if (b <= -1.75e+41) tmp = Float64(t * Float64(b * i)); elseif (b <= 2.2e+133) tmp = Float64(y * Float64(x * z)); else tmp = Float64(i * Float64(t * b)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i, j) tmp = 0.0; if (b <= -1.75e+41) tmp = t * (b * i); elseif (b <= 2.2e+133) tmp = y * (x * z); else tmp = i * (t * b); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := If[LessEqual[b, -1.75e+41], N[(t * N[(b * i), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.2e+133], N[(y * N[(x * z), $MachinePrecision]), $MachinePrecision], N[(i * N[(t * b), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.75 \cdot 10^{+41}:\\
\;\;\;\;t \cdot \left(b \cdot i\right)\\
\mathbf{elif}\;b \leq 2.2 \cdot 10^{+133}:\\
\;\;\;\;y \cdot \left(x \cdot z\right)\\
\mathbf{else}:\\
\;\;\;\;i \cdot \left(t \cdot b\right)\\
\end{array}
\end{array}
if b < -1.75e41Initial program 77.0%
Taylor expanded in t around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
mul-1-negN/A
remove-double-negN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6448.9
Applied rewrites48.9%
Taylor expanded in a around 0
Applied rewrites42.1%
if -1.75e41 < b < 2.2e133Initial program 67.8%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6439.1
Applied rewrites39.1%
Taylor expanded in c around 0
Applied rewrites27.0%
Applied rewrites28.8%
if 2.2e133 < b Initial program 75.6%
Taylor expanded in i around -inf
Applied rewrites64.5%
Taylor expanded in b around inf
Applied rewrites76.8%
Taylor expanded in t around inf
Applied rewrites50.1%
Final simplification35.6%
(FPCore (x y z t a b c i j) :precision binary64 (let* ((t_1 (* i (* t b)))) (if (<= b -1.75e+41) t_1 (if (<= b 2.2e+133) (* y (* x z)) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = i * (t * b);
double tmp;
if (b <= -1.75e+41) {
tmp = t_1;
} else if (b <= 2.2e+133) {
tmp = y * (x * z);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i, j)
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) :: t_1
real(8) :: tmp
t_1 = i * (t * b)
if (b <= (-1.75d+41)) then
tmp = t_1
else if (b <= 2.2d+133) then
tmp = y * (x * z)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = i * (t * b);
double tmp;
if (b <= -1.75e+41) {
tmp = t_1;
} else if (b <= 2.2e+133) {
tmp = y * (x * z);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i, j): t_1 = i * (t * b) tmp = 0 if b <= -1.75e+41: tmp = t_1 elif b <= 2.2e+133: tmp = y * (x * z) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i, j) t_1 = Float64(i * Float64(t * b)) tmp = 0.0 if (b <= -1.75e+41) tmp = t_1; elseif (b <= 2.2e+133) tmp = Float64(y * Float64(x * z)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i, j) t_1 = i * (t * b); tmp = 0.0; if (b <= -1.75e+41) tmp = t_1; elseif (b <= 2.2e+133) tmp = y * (x * z); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(i * N[(t * b), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1.75e+41], t$95$1, If[LessEqual[b, 2.2e+133], N[(y * N[(x * z), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := i \cdot \left(t \cdot b\right)\\
\mathbf{if}\;b \leq -1.75 \cdot 10^{+41}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 2.2 \cdot 10^{+133}:\\
\;\;\;\;y \cdot \left(x \cdot z\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.75e41 or 2.2e133 < b Initial program 76.5%
Taylor expanded in i around -inf
Applied rewrites66.9%
Taylor expanded in b around inf
Applied rewrites65.1%
Taylor expanded in t around inf
Applied rewrites43.5%
if -1.75e41 < b < 2.2e133Initial program 67.8%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6439.1
Applied rewrites39.1%
Taylor expanded in c around 0
Applied rewrites27.0%
Applied rewrites28.8%
Final simplification34.9%
(FPCore (x y z t a b c i j) :precision binary64 (let* ((t_1 (* b (* t i)))) (if (<= b -1.75e+41) t_1 (if (<= b 2.2e+133) (* y (* x z)) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = b * (t * i);
double tmp;
if (b <= -1.75e+41) {
tmp = t_1;
} else if (b <= 2.2e+133) {
tmp = y * (x * z);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i, j)
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) :: t_1
real(8) :: tmp
t_1 = b * (t * i)
if (b <= (-1.75d+41)) then
tmp = t_1
else if (b <= 2.2d+133) then
tmp = y * (x * z)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = b * (t * i);
double tmp;
if (b <= -1.75e+41) {
tmp = t_1;
} else if (b <= 2.2e+133) {
tmp = y * (x * z);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i, j): t_1 = b * (t * i) tmp = 0 if b <= -1.75e+41: tmp = t_1 elif b <= 2.2e+133: tmp = y * (x * z) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i, j) t_1 = Float64(b * Float64(t * i)) tmp = 0.0 if (b <= -1.75e+41) tmp = t_1; elseif (b <= 2.2e+133) tmp = Float64(y * Float64(x * z)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i, j) t_1 = b * (t * i); tmp = 0.0; if (b <= -1.75e+41) tmp = t_1; elseif (b <= 2.2e+133) tmp = y * (x * z); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(b * N[(t * i), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1.75e+41], t$95$1, If[LessEqual[b, 2.2e+133], N[(y * N[(x * z), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := b \cdot \left(t \cdot i\right)\\
\mathbf{if}\;b \leq -1.75 \cdot 10^{+41}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 2.2 \cdot 10^{+133}:\\
\;\;\;\;y \cdot \left(x \cdot z\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.75e41 or 2.2e133 < b Initial program 76.5%
Taylor expanded in t around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
mul-1-negN/A
remove-double-negN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6448.7
Applied rewrites48.7%
Taylor expanded in a around 0
Applied rewrites42.5%
if -1.75e41 < b < 2.2e133Initial program 67.8%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6439.1
Applied rewrites39.1%
Taylor expanded in c around 0
Applied rewrites27.0%
Applied rewrites28.8%
Final simplification34.5%
(FPCore (x y z t a b c i j) :precision binary64 (let* ((t_1 (* x (* y z)))) (if (<= z -1.15e+52) t_1 (if (<= z 2.8e+113) (* b (* t i)) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = x * (y * z);
double tmp;
if (z <= -1.15e+52) {
tmp = t_1;
} else if (z <= 2.8e+113) {
tmp = b * (t * i);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i, j)
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) :: t_1
real(8) :: tmp
t_1 = x * (y * z)
if (z <= (-1.15d+52)) then
tmp = t_1
else if (z <= 2.8d+113) then
tmp = b * (t * i)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = x * (y * z);
double tmp;
if (z <= -1.15e+52) {
tmp = t_1;
} else if (z <= 2.8e+113) {
tmp = b * (t * i);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i, j): t_1 = x * (y * z) tmp = 0 if z <= -1.15e+52: tmp = t_1 elif z <= 2.8e+113: tmp = b * (t * i) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i, j) t_1 = Float64(x * Float64(y * z)) tmp = 0.0 if (z <= -1.15e+52) tmp = t_1; elseif (z <= 2.8e+113) tmp = Float64(b * Float64(t * i)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i, j) t_1 = x * (y * z); tmp = 0.0; if (z <= -1.15e+52) tmp = t_1; elseif (z <= 2.8e+113) tmp = b * (t * i); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(x * N[(y * z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.15e+52], t$95$1, If[LessEqual[z, 2.8e+113], N[(b * N[(t * i), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(y \cdot z\right)\\
\mathbf{if}\;z \leq -1.15 \cdot 10^{+52}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 2.8 \cdot 10^{+113}:\\
\;\;\;\;b \cdot \left(t \cdot i\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.15e52 or 2.79999999999999998e113 < z Initial program 58.3%
Taylor expanded in z around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
neg-mul-1N/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6465.5
Applied rewrites65.5%
Taylor expanded in c around 0
Applied rewrites41.9%
if -1.15e52 < z < 2.79999999999999998e113Initial program 79.1%
Taylor expanded in t around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
mul-1-negN/A
remove-double-negN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6446.6
Applied rewrites46.6%
Taylor expanded in a around 0
Applied rewrites29.2%
(FPCore (x y z t a b c i j) :precision binary64 (* b (* t i)))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
return b * (t * i);
}
real(8) function code(x, y, z, t, a, b, c, i, j)
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
code = b * (t * i)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
return b * (t * i);
}
def code(x, y, z, t, a, b, c, i, j): return b * (t * i)
function code(x, y, z, t, a, b, c, i, j) return Float64(b * Float64(t * i)) end
function tmp = code(x, y, z, t, a, b, c, i, j) tmp = b * (t * i); end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := N[(b * N[(t * i), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
b \cdot \left(t \cdot i\right)
\end{array}
Initial program 71.4%
Taylor expanded in t around inf
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
mul-1-negN/A
remove-double-negN/A
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f6440.8
Applied rewrites40.8%
Taylor expanded in a around 0
Applied rewrites23.0%
(FPCore (x y z t a b c i j)
:precision binary64
(let* ((t_1 (* j (- (* c a) (* y i))))
(t_2
(+
(-
(* x (- (* y z) (* t a)))
(/
(* b (- (pow (* c z) 2.0) (pow (* t i) 2.0)))
(+ (* c z) (* t i))))
t_1)))
(if (< x -1.469694296777705e-64)
t_2
(if (< x 3.2113527362226803e-147)
(- (* (- (* b i) (* x a)) t) (- (* z (* c b)) t_1))
t_2))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = j * ((c * a) - (y * i));
double t_2 = ((x * ((y * z) - (t * a))) - ((b * (pow((c * z), 2.0) - pow((t * i), 2.0))) / ((c * z) + (t * i)))) + t_1;
double tmp;
if (x < -1.469694296777705e-64) {
tmp = t_2;
} else if (x < 3.2113527362226803e-147) {
tmp = (((b * i) - (x * a)) * t) - ((z * (c * b)) - t_1);
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i, j)
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) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = j * ((c * a) - (y * i))
t_2 = ((x * ((y * z) - (t * a))) - ((b * (((c * z) ** 2.0d0) - ((t * i) ** 2.0d0))) / ((c * z) + (t * i)))) + t_1
if (x < (-1.469694296777705d-64)) then
tmp = t_2
else if (x < 3.2113527362226803d-147) then
tmp = (((b * i) - (x * a)) * t) - ((z * (c * b)) - t_1)
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 t_1 = j * ((c * a) - (y * i));
double t_2 = ((x * ((y * z) - (t * a))) - ((b * (Math.pow((c * z), 2.0) - Math.pow((t * i), 2.0))) / ((c * z) + (t * i)))) + t_1;
double tmp;
if (x < -1.469694296777705e-64) {
tmp = t_2;
} else if (x < 3.2113527362226803e-147) {
tmp = (((b * i) - (x * a)) * t) - ((z * (c * b)) - t_1);
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i, j): t_1 = j * ((c * a) - (y * i)) t_2 = ((x * ((y * z) - (t * a))) - ((b * (math.pow((c * z), 2.0) - math.pow((t * i), 2.0))) / ((c * z) + (t * i)))) + t_1 tmp = 0 if x < -1.469694296777705e-64: tmp = t_2 elif x < 3.2113527362226803e-147: tmp = (((b * i) - (x * a)) * t) - ((z * (c * b)) - t_1) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b, c, i, j) t_1 = Float64(j * Float64(Float64(c * a) - Float64(y * i))) t_2 = Float64(Float64(Float64(x * Float64(Float64(y * z) - Float64(t * a))) - Float64(Float64(b * Float64((Float64(c * z) ^ 2.0) - (Float64(t * i) ^ 2.0))) / Float64(Float64(c * z) + Float64(t * i)))) + t_1) tmp = 0.0 if (x < -1.469694296777705e-64) tmp = t_2; elseif (x < 3.2113527362226803e-147) tmp = Float64(Float64(Float64(Float64(b * i) - Float64(x * a)) * t) - Float64(Float64(z * Float64(c * b)) - t_1)); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i, j) t_1 = j * ((c * a) - (y * i)); t_2 = ((x * ((y * z) - (t * a))) - ((b * (((c * z) ^ 2.0) - ((t * i) ^ 2.0))) / ((c * z) + (t * i)))) + t_1; tmp = 0.0; if (x < -1.469694296777705e-64) tmp = t_2; elseif (x < 3.2113527362226803e-147) tmp = (((b * i) - (x * a)) * t) - ((z * (c * b)) - t_1); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(j * N[(N[(c * a), $MachinePrecision] - N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(x * N[(N[(y * z), $MachinePrecision] - N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(b * N[(N[Power[N[(c * z), $MachinePrecision], 2.0], $MachinePrecision] - N[Power[N[(t * i), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(c * z), $MachinePrecision] + N[(t * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t$95$1), $MachinePrecision]}, If[Less[x, -1.469694296777705e-64], t$95$2, If[Less[x, 3.2113527362226803e-147], N[(N[(N[(N[(b * i), $MachinePrecision] - N[(x * a), $MachinePrecision]), $MachinePrecision] * t), $MachinePrecision] - N[(N[(z * N[(c * b), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision]), $MachinePrecision], t$95$2]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := j \cdot \left(c \cdot a - y \cdot i\right)\\
t_2 := \left(x \cdot \left(y \cdot z - t \cdot a\right) - \frac{b \cdot \left({\left(c \cdot z\right)}^{2} - {\left(t \cdot i\right)}^{2}\right)}{c \cdot z + t \cdot i}\right) + t\_1\\
\mathbf{if}\;x < -1.469694296777705 \cdot 10^{-64}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;x < 3.2113527362226803 \cdot 10^{-147}:\\
\;\;\;\;\left(b \cdot i - x \cdot a\right) \cdot t - \left(z \cdot \left(c \cdot b\right) - t\_1\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
herbie shell --seed 2024220
(FPCore (x y z t a b c i j)
:name "Data.Colour.Matrix:determinant from colour-2.3.3, A"
:precision binary64
:alt
(! :herbie-platform default (if (< x -293938859355541/2000000000000000000000000000000000000000000000000000000000000000000000000000000) (+ (- (* x (- (* y z) (* t a))) (/ (* b (- (pow (* c z) 2) (pow (* t i) 2))) (+ (* c z) (* t i)))) (* j (- (* c a) (* y i)))) (if (< x 32113527362226803/10000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (- (* (- (* b i) (* x a)) t) (- (* z (* c b)) (* j (- (* c a) (* y i))))) (+ (- (* x (- (* y z) (* t a))) (/ (* b (- (pow (* c z) 2) (pow (* t i) 2))) (+ (* c z) (* t i)))) (* j (- (* c a) (* y i)))))))
(+ (- (* x (- (* y z) (* t a))) (* b (- (* c z) (* t i)))) (* j (- (* c a) (* y i)))))