\left(\left(\left(\left(\left(x \cdot y - z \cdot t\right) \cdot \left(a \cdot b - c \cdot i\right) - \left(x \cdot j - z \cdot k\right) \cdot \left(y0 \cdot b - y1 \cdot i\right)\right) + \left(x \cdot y2 - z \cdot y3\right) \cdot \left(y0 \cdot c - y1 \cdot a\right)\right) + \left(t \cdot j - y \cdot k\right) \cdot \left(y4 \cdot b - y5 \cdot i\right)\right) - \left(t \cdot y2 - y \cdot y3\right) \cdot \left(y4 \cdot c - y5 \cdot a\right)\right) + \left(k \cdot y2 - j \cdot y3\right) \cdot \left(y4 \cdot y1 - y5 \cdot y0\right)
\begin{array}{l}
t_1 := x \cdot \left(\left(a \cdot \left(y \cdot b\right) + \left(c \cdot \left(y0 \cdot y2\right) + y1 \cdot \left(i \cdot j\right)\right)\right) - \left(c \cdot \left(y \cdot i\right) + \left(y0 \cdot \left(b \cdot j\right) + y1 \cdot \left(a \cdot y2\right)\right)\right)\right)\\
\mathbf{if}\;x \leq -4.528176101828275 \cdot 10^{+220}:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;\begin{array}{l}
t_2 := y2 \cdot t - y \cdot y3\\
t_3 := x \cdot j - k \cdot z\\
t_4 := y1 \cdot i - b \cdot y0\\
t_5 := a \cdot y5 - c \cdot y4\\
t_6 := x \cdot y - t \cdot z\\
t_7 := y2 \cdot k - j \cdot y3\\
t_8 := a \cdot b - c \cdot i\\
t_9 := y1 \cdot y4 - y0 \cdot y5\\
\mathbf{if}\;x \leq -9.863007287958035 \cdot 10^{+138}:\\
\;\;\;\;\mathsf{fma}\left(t_7, t_9, \mathsf{fma}\left(t_2, t_5, \mathsf{fma}\left(t_6, t_8, \mathsf{fma}\left(t_3, t_4, y3 \cdot \left(y1 \cdot \left(a \cdot z\right) - c \cdot \left(y0 \cdot z\right)\right)\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\begin{array}{l}
t_10 := y0 \cdot \left(y3 \cdot z\right)\\
t_11 := y1 \cdot \left(x \cdot y2\right)\\
t_12 := \left(b \cdot j\right) \cdot t\\
t_13 := y \cdot \left(x \cdot b\right)\\
t_14 := y0 \cdot \left(y2 \cdot y5\right)\\
t_15 := y4 \cdot \left(y2 \cdot y1\right)\\
t_16 := y0 \cdot \left(b \cdot z\right)\\
t_17 := k \cdot t_16\\
t_18 := y1 \cdot \left(j \cdot y3\right)\\
t_19 := y3 \cdot \left(j \cdot y5\right)\\
t_20 := a \cdot \left(y3 \cdot y5\right)\\
t_21 := t \cdot \left(b \cdot z\right)\\
t_22 := a \cdot t_21\\
t_23 := y1 \cdot \left(y3 \cdot z\right)\\
t_24 := y4 \cdot \left(y \cdot y3\right)\\
t_25 := y4 \cdot \left(y2 \cdot t\right)\\
t_26 := t \cdot \left(y2 \cdot y5\right)\\
t_27 := y0 \cdot \left(x \cdot y2\right)\\
t_28 := \left(y \cdot b\right) \cdot y4\\
t_29 := \mathsf{fma}\left(a, t_26, \mathsf{fma}\left(y4, t_12, \mathsf{fma}\left(k, t_15, \mathsf{fma}\left(y0, t_19, \mathsf{fma}\left(a, t_23, \mathsf{fma}\left(c, t_24, \mathsf{fma}\left(c, t_27, \mathsf{fma}\left(a, t_13, t_17\right)\right)\right)\right)\right)\right)\right)\right) - \mathsf{fma}\left(y, t_20, \mathsf{fma}\left(c, t_10, \mathsf{fma}\left(c, t_25, \mathsf{fma}\left(k, t_28, \mathsf{fma}\left(a, t_11, \mathsf{fma}\left(k, t_14, \mathsf{fma}\left(y4, t_18, \mathsf{fma}\left(y0, j \cdot \left(x \cdot b\right), t_22\right)\right)\right)\right)\right)\right)\right)\right)\\
\mathbf{if}\;x \leq -1.1434668545282551 \cdot 10^{+132}:\\
\;\;\;\;t_29\\
\mathbf{else}:\\
\;\;\;\;\begin{array}{l}
t_30 := \mathsf{fma}\left(t_7, t_9, \mathsf{fma}\left(t_2, t_5, \mathsf{fma}\left(t_6, t_8, \mathsf{fma}\left(t_3, t_4, \mathsf{fma}\left(x \cdot y2 - y3 \cdot z, c \cdot y0 - a \cdot y1, \left(j \cdot t - y \cdot k\right) \cdot \left(b \cdot y4 - i \cdot y5\right)\right)\right)\right)\right)\right)\\
\mathbf{if}\;x \leq -1.2934549413283296 \cdot 10^{-21}:\\
\;\;\;\;t_30\\
\mathbf{else}:\\
\;\;\;\;\begin{array}{l}
t_31 := y \cdot \left(i \cdot y5\right)\\
t_32 := k \cdot t_31\\
t_33 := y \cdot \left(x \cdot i\right)\\
t_34 := b \cdot \left(x \cdot j\right)\\
t_35 := k \cdot \left(y1 \cdot \left(i \cdot z\right)\right)\\
t_36 := i \cdot \left(t \cdot z\right)\\
t_37 := c \cdot t_36\\
t_38 := t \cdot \left(j \cdot y5\right)\\
t_39 := k \cdot t_28\\
\mathbf{if}\;x \leq -2.289257968400257 \cdot 10^{-89}:\\
\;\;\;\;\left(a \cdot t_26 + \left(y4 \cdot t_12 + \left(k \cdot t_15 + \left(t_37 + \left(c \cdot t_27 + \left(a \cdot t_13 + \left(t_17 + \left(y1 \cdot \left(i \cdot \left(x \cdot j\right)\right) + t_32\right)\right)\right)\right)\right)\right)\right)\right) - \left(c \cdot t_33 + \left(t_22 + \left(c \cdot t_25 + \left(i \cdot t_38 + \left(t_39 + \left(a \cdot t_11 + \left(k \cdot t_14 + \left(y0 \cdot t_34 + t_35\right)\right)\right)\right)\right)\right)\right)\right)\\
\mathbf{elif}\;x \leq -4.620735289744297 \cdot 10^{-142}:\\
\;\;\;\;\mathsf{fma}\left(t_7, t_9, \mathsf{fma}\left(t_2, t_5, \mathsf{fma}\left(t_6, t_8, \mathsf{fma}\left(t_3, t_4, \mathsf{fma}\left(c, t_27, \mathsf{fma}\left(y4, t_12, t_32\right)\right) - \mathsf{fma}\left(c, t_10, \mathsf{fma}\left(i, t_38, t_39\right)\right)\right)\right)\right)\right)\\
\mathbf{elif}\;x \leq -4.2139772218421345 \cdot 10^{-164}:\\
\;\;\;\;\mathsf{fma}\left(a, t_26, \mathsf{fma}\left(k, t_15, \mathsf{fma}\left(a, t_23, \mathsf{fma}\left(c, t_24, \mathsf{fma}\left(c, t_36, \mathsf{fma}\left(c, t_27, \mathsf{fma}\left(a, t_13, k \cdot \left(t_16 + t_31\right)\right)\right)\right)\right)\right)\right)\right) - \mathsf{fma}\left(c, t_10, \mathsf{fma}\left(a, y \cdot \left(y3 \cdot y5\right), \mathsf{fma}\left(c, t_33, \mathsf{fma}\left(c, t_25, \mathsf{fma}\left(k, t_28, \mathsf{fma}\left(a, t_11, \mathsf{fma}\left(k, t_14, \mathsf{fma}\left(a, t_21, t_35\right)\right)\right)\right)\right)\right)\right)\right)\\
\mathbf{elif}\;x \leq -9.622849633415526 \cdot 10^{-185}:\\
\;\;\;\;t_29\\
\mathbf{elif}\;x \leq -2.4025835670476306 \cdot 10^{-276}:\\
\;\;\;\;t_30\\
\mathbf{else}:\\
\;\;\;\;\begin{array}{l}
t_40 := y1 \cdot \left(x \cdot j\right)\\
\mathbf{if}\;x \leq 2.1939841572173253 \cdot 10^{-171}:\\
\;\;\;\;\mathsf{fma}\left(t_7, t_9, \mathsf{fma}\left(t_2, t_5, \mathsf{fma}\left(c, t_27, \mathsf{fma}\left(k, t_16, \mathsf{fma}\left(y4, t_12, \mathsf{fma}\left(i, t_40, \mathsf{fma}\left(a, t_23, t_37\right)\right)\right)\right)\right) - \mathsf{fma}\left(c, t_10, \mathsf{fma}\left(a, t_11, \mathsf{fma}\left(y0, t_34, \mathsf{fma}\left(a, t_21, \mathsf{fma}\left(i, t_38, t_35\right)\right)\right)\right)\right)\right)\right)\\
\mathbf{elif}\;x \leq 4.7340441146328203 \cdot 10^{+126}:\\
\;\;\;\;\mathsf{fma}\left(a, t_26, \mathsf{fma}\left(y4, t_12, \mathsf{fma}\left(k, t_15, \mathsf{fma}\left(y0, t_19, \mathsf{fma}\left(a, t_23, \mathsf{fma}\left(c, t_24, \mathsf{fma}\left(c, t_36, \mathsf{fma}\left(c, t_27, \mathsf{fma}\left(y, a \cdot \left(x \cdot b\right), \mathsf{fma}\left(k, t_16, \mathsf{fma}\left(i, t_40, t_32\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right) - \mathsf{fma}\left(c, t_10, \mathsf{fma}\left(y, t_20, \mathsf{fma}\left(c, t_33, \mathsf{fma}\left(a, t_21, \mathsf{fma}\left(c, t_25, \mathsf{fma}\left(i, t_38, \mathsf{fma}\left(k, t_28, \mathsf{fma}\left(k, t_14, \mathsf{fma}\left(y1, a \cdot \left(x \cdot y2\right), \mathsf{fma}\left(y4, t_18, \mathsf{fma}\left(y0, t_34, t_35\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}\\
\end{array}\\
\end{array}\\
\end{array}\\
\end{array}\\
\end{array}
(FPCore (x y z t a b c i j k y0 y1 y2 y3 y4 y5)
:precision binary64
(+
(-
(+
(+
(-
(* (- (* x y) (* z t)) (- (* a b) (* c i)))
(* (- (* x j) (* z k)) (- (* y0 b) (* y1 i))))
(* (- (* x y2) (* z y3)) (- (* y0 c) (* y1 a))))
(* (- (* t j) (* y k)) (- (* y4 b) (* y5 i))))
(* (- (* t y2) (* y y3)) (- (* y4 c) (* y5 a))))
(* (- (* k y2) (* j y3)) (- (* y4 y1) (* y5 y0)))))(FPCore (x y z t a b c i j k y0 y1 y2 y3 y4 y5)
:precision binary64
(let* ((t_1
(*
x
(-
(+ (* a (* y b)) (+ (* c (* y0 y2)) (* y1 (* i j))))
(+ (* c (* y i)) (+ (* y0 (* b j)) (* y1 (* a y2))))))))
(if (<= x -4.528176101828275e+220)
t_1
(let* ((t_2 (- (* y2 t) (* y y3)))
(t_3 (- (* x j) (* k z)))
(t_4 (- (* y1 i) (* b y0)))
(t_5 (- (* a y5) (* c y4)))
(t_6 (- (* x y) (* t z)))
(t_7 (- (* y2 k) (* j y3)))
(t_8 (- (* a b) (* c i)))
(t_9 (- (* y1 y4) (* y0 y5))))
(if (<= x -9.863007287958035e+138)
(fma
t_7
t_9
(fma
t_2
t_5
(fma
t_6
t_8
(fma t_3 t_4 (* y3 (- (* y1 (* a z)) (* c (* y0 z))))))))
(let* ((t_10 (* y0 (* y3 z)))
(t_11 (* y1 (* x y2)))
(t_12 (* (* b j) t))
(t_13 (* y (* x b)))
(t_14 (* y0 (* y2 y5)))
(t_15 (* y4 (* y2 y1)))
(t_16 (* y0 (* b z)))
(t_17 (* k t_16))
(t_18 (* y1 (* j y3)))
(t_19 (* y3 (* j y5)))
(t_20 (* a (* y3 y5)))
(t_21 (* t (* b z)))
(t_22 (* a t_21))
(t_23 (* y1 (* y3 z)))
(t_24 (* y4 (* y y3)))
(t_25 (* y4 (* y2 t)))
(t_26 (* t (* y2 y5)))
(t_27 (* y0 (* x y2)))
(t_28 (* (* y b) y4))
(t_29
(-
(fma
a
t_26
(fma
y4
t_12
(fma
k
t_15
(fma
y0
t_19
(fma
a
t_23
(fma c t_24 (fma c t_27 (fma a t_13 t_17))))))))
(fma
y
t_20
(fma
c
t_10
(fma
c
t_25
(fma
k
t_28
(fma
a
t_11
(fma
k
t_14
(fma y4 t_18 (fma y0 (* j (* x b)) t_22)))))))))))
(if (<= x -1.1434668545282551e+132)
t_29
(let* ((t_30
(fma
t_7
t_9
(fma
t_2
t_5
(fma
t_6
t_8
(fma
t_3
t_4
(fma
(- (* x y2) (* y3 z))
(- (* c y0) (* a y1))
(* (- (* j t) (* y k)) (- (* b y4) (* i y5))))))))))
(if (<= x -1.2934549413283296e-21)
t_30
(let* ((t_31 (* y (* i y5)))
(t_32 (* k t_31))
(t_33 (* y (* x i)))
(t_34 (* b (* x j)))
(t_35 (* k (* y1 (* i z))))
(t_36 (* i (* t z)))
(t_37 (* c t_36))
(t_38 (* t (* j y5)))
(t_39 (* k t_28)))
(if (<= x -2.289257968400257e-89)
(-
(+
(* a t_26)
(+
(* y4 t_12)
(+
(* k t_15)
(+
t_37
(+
(* c t_27)
(+
(* a t_13)
(+ t_17 (+ (* y1 (* i (* x j))) t_32))))))))
(+
(* c t_33)
(+
t_22
(+
(* c t_25)
(+
(* i t_38)
(+
t_39
(+
(* a t_11)
(+ (* k t_14) (+ (* y0 t_34) t_35)))))))))
(if (<= x -4.620735289744297e-142)
(fma
t_7
t_9
(fma
t_2
t_5
(fma
t_6
t_8
(fma
t_3
t_4
(-
(fma c t_27 (fma y4 t_12 t_32))
(fma c t_10 (fma i t_38 t_39)))))))
(if (<= x -4.2139772218421345e-164)
(-
(fma
a
t_26
(fma
k
t_15
(fma
a
t_23
(fma
c
t_24
(fma
c
t_36
(fma
c
t_27
(fma a t_13 (* k (+ t_16 t_31)))))))))
(fma
c
t_10
(fma
a
(* y (* y3 y5))
(fma
c
t_33
(fma
c
t_25
(fma
k
t_28
(fma
a
t_11
(fma k t_14 (fma a t_21 t_35)))))))))
(if (<= x -9.622849633415526e-185)
t_29
(if (<= x -2.4025835670476306e-276)
t_30
(let* ((t_40 (* y1 (* x j))))
(if (<= x 2.1939841572173253e-171)
(fma
t_7
t_9
(fma
t_2
t_5
(-
(fma
c
t_27
(fma
k
t_16
(fma
y4
t_12
(fma i t_40 (fma a t_23 t_37)))))
(fma
c
t_10
(fma
a
t_11
(fma
y0
t_34
(fma a t_21 (fma i t_38 t_35))))))))
(if (<= x 4.7340441146328203e+126)
(-
(fma
a
t_26
(fma
y4
t_12
(fma
k
t_15
(fma
y0
t_19
(fma
a
t_23
(fma
c
t_24
(fma
c
t_36
(fma
c
t_27
(fma
y
(* a (* x b))
(fma
k
t_16
(fma i t_40 t_32)))))))))))
(fma
c
t_10
(fma
y
t_20
(fma
c
t_33
(fma
a
t_21
(fma
c
t_25
(fma
i
t_38
(fma
k
t_28
(fma
k
t_14
(fma
y1
(* a (* x y2))
(fma
y4
t_18
(fma y0 t_34 t_35))))))))))))
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 y0, double y1, double y2, double y3, double y4, double y5) {
return (((((((x * y) - (z * t)) * ((a * b) - (c * i))) - (((x * j) - (z * k)) * ((y0 * b) - (y1 * i)))) + (((x * y2) - (z * y3)) * ((y0 * c) - (y1 * a)))) + (((t * j) - (y * k)) * ((y4 * b) - (y5 * i)))) - (((t * y2) - (y * y3)) * ((y4 * c) - (y5 * a)))) + (((k * y2) - (j * y3)) * ((y4 * y1) - (y5 * y0)));
}
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k, double y0, double y1, double y2, double y3, double y4, double y5) {
double t_1 = x * (((a * (y * b)) + ((c * (y0 * y2)) + (y1 * (i * j)))) - ((c * (y * i)) + ((y0 * (b * j)) + (y1 * (a * y2)))));
double tmp;
if (x <= -4.528176101828275e+220) {
tmp = t_1;
} else {
double t_2 = (y2 * t) - (y * y3);
double t_3 = (x * j) - (k * z);
double t_4 = (y1 * i) - (b * y0);
double t_5 = (a * y5) - (c * y4);
double t_6 = (x * y) - (t * z);
double t_7 = (y2 * k) - (j * y3);
double t_8 = (a * b) - (c * i);
double t_9 = (y1 * y4) - (y0 * y5);
double tmp_1;
if (x <= -9.863007287958035e+138) {
tmp_1 = fma(t_7, t_9, fma(t_2, t_5, fma(t_6, t_8, fma(t_3, t_4, (y3 * ((y1 * (a * z)) - (c * (y0 * z))))))));
} else {
double t_10 = y0 * (y3 * z);
double t_11 = y1 * (x * y2);
double t_12 = (b * j) * t;
double t_13 = y * (x * b);
double t_14 = y0 * (y2 * y5);
double t_15 = y4 * (y2 * y1);
double t_16 = y0 * (b * z);
double t_17 = k * t_16;
double t_18 = y1 * (j * y3);
double t_19 = y3 * (j * y5);
double t_20 = a * (y3 * y5);
double t_21 = t * (b * z);
double t_22 = a * t_21;
double t_23 = y1 * (y3 * z);
double t_24 = y4 * (y * y3);
double t_25 = y4 * (y2 * t);
double t_26 = t * (y2 * y5);
double t_27 = y0 * (x * y2);
double t_28 = (y * b) * y4;
double t_29 = fma(a, t_26, fma(y4, t_12, fma(k, t_15, fma(y0, t_19, fma(a, t_23, fma(c, t_24, fma(c, t_27, fma(a, t_13, t_17)))))))) - fma(y, t_20, fma(c, t_10, fma(c, t_25, fma(k, t_28, fma(a, t_11, fma(k, t_14, fma(y4, t_18, fma(y0, (j * (x * b)), t_22))))))));
double tmp_2;
if (x <= -1.1434668545282551e+132) {
tmp_2 = t_29;
} else {
double t_30 = fma(t_7, t_9, fma(t_2, t_5, fma(t_6, t_8, fma(t_3, t_4, fma(((x * y2) - (y3 * z)), ((c * y0) - (a * y1)), (((j * t) - (y * k)) * ((b * y4) - (i * y5))))))));
double tmp_3;
if (x <= -1.2934549413283296e-21) {
tmp_3 = t_30;
} else {
double t_31 = y * (i * y5);
double t_32 = k * t_31;
double t_33 = y * (x * i);
double t_34 = b * (x * j);
double t_35 = k * (y1 * (i * z));
double t_36 = i * (t * z);
double t_37 = c * t_36;
double t_38 = t * (j * y5);
double t_39 = k * t_28;
double tmp_4;
if (x <= -2.289257968400257e-89) {
tmp_4 = ((a * t_26) + ((y4 * t_12) + ((k * t_15) + (t_37 + ((c * t_27) + ((a * t_13) + (t_17 + ((y1 * (i * (x * j))) + t_32)))))))) - ((c * t_33) + (t_22 + ((c * t_25) + ((i * t_38) + (t_39 + ((a * t_11) + ((k * t_14) + ((y0 * t_34) + t_35))))))));
} else if (x <= -4.620735289744297e-142) {
tmp_4 = fma(t_7, t_9, fma(t_2, t_5, fma(t_6, t_8, fma(t_3, t_4, (fma(c, t_27, fma(y4, t_12, t_32)) - fma(c, t_10, fma(i, t_38, t_39)))))));
} else if (x <= -4.2139772218421345e-164) {
tmp_4 = fma(a, t_26, fma(k, t_15, fma(a, t_23, fma(c, t_24, fma(c, t_36, fma(c, t_27, fma(a, t_13, (k * (t_16 + t_31))))))))) - fma(c, t_10, fma(a, (y * (y3 * y5)), fma(c, t_33, fma(c, t_25, fma(k, t_28, fma(a, t_11, fma(k, t_14, fma(a, t_21, t_35))))))));
} else if (x <= -9.622849633415526e-185) {
tmp_4 = t_29;
} else if (x <= -2.4025835670476306e-276) {
tmp_4 = t_30;
} else {
double t_40 = y1 * (x * j);
double tmp_5;
if (x <= 2.1939841572173253e-171) {
tmp_5 = fma(t_7, t_9, fma(t_2, t_5, (fma(c, t_27, fma(k, t_16, fma(y4, t_12, fma(i, t_40, fma(a, t_23, t_37))))) - fma(c, t_10, fma(a, t_11, fma(y0, t_34, fma(a, t_21, fma(i, t_38, t_35))))))));
} else if (x <= 4.7340441146328203e+126) {
tmp_5 = fma(a, t_26, fma(y4, t_12, fma(k, t_15, fma(y0, t_19, fma(a, t_23, fma(c, t_24, fma(c, t_36, fma(c, t_27, fma(y, (a * (x * b)), fma(k, t_16, fma(i, t_40, t_32))))))))))) - fma(c, t_10, fma(y, t_20, fma(c, t_33, fma(a, t_21, fma(c, t_25, fma(i, t_38, fma(k, t_28, fma(k, t_14, fma(y1, (a * (x * y2)), fma(y4, t_18, fma(y0, t_34, t_35)))))))))));
} else {
tmp_5 = t_1;
}
tmp_4 = tmp_5;
}
tmp_3 = tmp_4;
}
tmp_2 = tmp_3;
}
tmp_1 = tmp_2;
}
tmp = tmp_1;
}
return tmp;
}
Bits error versus x
Bits error versus y
Bits error versus z
Bits error versus t
Bits error versus a
Bits error versus b
Bits error versus c
Bits error versus i
Bits error versus j
Bits error versus k
Bits error versus y0
Bits error versus y1
Bits error versus y2
Bits error versus y3
Bits error versus y4
Bits error versus y5
| Original | 27.1 |
|---|---|
| Target | 30.7 |
| Herbie | 29.9 |
if x < -4.52817610182827526e220 or 4.73404411463282034e126 < x Initial program 33.2
Simplified33.2
Taylor expanded in x around inf 28.6
if -4.52817610182827526e220 < x < -9.8630072879580355e138Initial program 30.3
Simplified30.3
Taylor expanded in y3 around inf 37.3
if -9.8630072879580355e138 < x < -1.14346685452825512e132 or -4.21397722184213455e-164 < x < -9.6228496334155264e-185Initial program 22.6
Simplified22.6
Taylor expanded in i around 0 36.3
Simplified36.3
if -1.14346685452825512e132 < x < -1.2934549413283296e-21 or -9.6228496334155264e-185 < x < -2.4025835670476306e-276Initial program 26.4
Simplified26.4
if -1.2934549413283296e-21 < x < -2.2892579684002569e-89Initial program 25.1
Simplified25.1
Taylor expanded in y3 around 0 33.4
if -2.2892579684002569e-89 < x < -4.62073528974429701e-142Initial program 24.0
Simplified24.0
Taylor expanded in a around 0 30.6
Simplified30.6
if -4.62073528974429701e-142 < x < -4.21397722184213455e-164Initial program 28.5
Simplified28.4
Taylor expanded in j around 0 31.2
Simplified31.2
if -2.4025835670476306e-276 < x < 2.19398415721732528e-171Initial program 28.2
Simplified28.2
Taylor expanded in y around 0 31.4
Simplified31.4
if 2.19398415721732528e-171 < x < 4.73404411463282034e126Initial program 25.5
Simplified25.5
Taylor expanded in k around 0 29.3
Simplified29.3
Final simplification29.9
herbie shell --seed 2021307
(FPCore (x y z t a b c i j k y0 y1 y2 y3 y4 y5)
:name "Linear.Matrix:det44 from linear-1.19.1.3"
:precision binary64
:herbie-target
(if (< y4 -7.206256231996481e+60) (- (- (* (- (* b a) (* i c)) (- (* y x) (* t z))) (- (* (- (* j x) (* k z)) (- (* y0 b) (* i y1))) (* (- (* j t) (* k y)) (- (* y4 b) (* y5 i))))) (- (/ (- (* y2 t) (* y3 y)) (/ 1.0 (- (* y4 c) (* y5 a)))) (* (- (* y2 k) (* y3 j)) (- (* y4 y1) (* y5 y0))))) (if (< y4 -3.364603505246317e-66) (+ (- (- (- (* (* t c) (* i z)) (* (* a t) (* b z))) (* (* y c) (* i x))) (* (- (* b y0) (* i y1)) (- (* j x) (* k z)))) (- (* (- (* y0 c) (* a y1)) (- (* x y2) (* z y3))) (- (* (- (* t y2) (* y y3)) (- (* y4 c) (* a y5))) (* (- (* y1 y4) (* y5 y0)) (- (* k y2) (* j y3)))))) (if (< y4 -1.2000065055686116e-105) (+ (+ (- (* (- (* j t) (* k y)) (- (* y4 b) (* y5 i))) (* (* y3 y) (- (* y5 a) (* y4 c)))) (+ (* (* y5 a) (* t y2)) (* (- (* k y2) (* j y3)) (- (* y4 y1) (* y5 y0))))) (- (* (- (* x y2) (* z y3)) (- (* c y0) (* a y1))) (- (* (- (* b y0) (* i y1)) (- (* j x) (* k z))) (* (- (* y x) (* z t)) (- (* b a) (* i c)))))) (if (< y4 6.718963124057495e-279) (+ (- (- (- (* (* k y) (* y5 i)) (* (* y b) (* y4 k))) (* (* y5 t) (* i j))) (- (* (- (* y2 t) (* y3 y)) (- (* y4 c) (* y5 a))) (* (- (* y2 k) (* y3 j)) (- (* y4 y1) (* y5 y0))))) (- (* (- (* b a) (* i c)) (- (* y x) (* t z))) (- (* (- (* j x) (* k z)) (- (* y0 b) (* i y1))) (* (- (* y2 x) (* y3 z)) (- (* c y0) (* y1 a)))))) (if (< y4 4.77962681403792e-222) (+ (+ (- (* (- (* j t) (* k y)) (- (* y4 b) (* y5 i))) (* (* y3 y) (- (* y5 a) (* y4 c)))) (+ (* (* y5 a) (* t y2)) (* (- (* k y2) (* j y3)) (- (* y4 y1) (* y5 y0))))) (- (* (- (* x y2) (* z y3)) (- (* c y0) (* a y1))) (- (* (- (* b y0) (* i y1)) (- (* j x) (* k z))) (* (- (* y x) (* z t)) (- (* b a) (* i c)))))) (if (< y4 2.2852241541266835e-175) (+ (- (- (- (* (* k y) (* y5 i)) (* (* y b) (* y4 k))) (* (* y5 t) (* i j))) (- (* (- (* y2 t) (* y3 y)) (- (* y4 c) (* y5 a))) (* (- (* y2 k) (* y3 j)) (- (* y4 y1) (* y5 y0))))) (- (* (- (* b a) (* i c)) (- (* y x) (* t z))) (- (* (- (* j x) (* k z)) (- (* y0 b) (* i y1))) (* (- (* y2 x) (* y3 z)) (- (* c y0) (* y1 a)))))) (+ (- (+ (+ (- (* (- (* x y) (* z t)) (- (* a b) (* c i))) (- (* k (* i (* z y1))) (+ (* j (* i (* x y1))) (* y0 (* k (* z b)))))) (- (* z (* y3 (* a y1))) (+ (* y2 (* x (* a y1))) (* y0 (* z (* c y3)))))) (* (- (* t j) (* y k)) (- (* y4 b) (* y5 i)))) (* (- (* t y2) (* y y3)) (- (* y4 c) (* y5 a)))) (* (- (* k y2) (* j y3)) (- (* y4 y1) (* y5 y0))))))))))
(+ (- (+ (+ (- (* (- (* x y) (* z t)) (- (* a b) (* c i))) (* (- (* x j) (* z k)) (- (* y0 b) (* y1 i)))) (* (- (* x y2) (* z y3)) (- (* y0 c) (* y1 a)))) (* (- (* t j) (* y k)) (- (* y4 b) (* y5 i)))) (* (- (* t y2) (* y y3)) (- (* y4 c) (* y5 a)))) (* (- (* k y2) (* j y3)) (- (* y4 y1) (* y5 y0)))))