(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 y2) (* z y3)) (- (* c y0) (* a y1))))
(t_2 (* c (* y4 (* y y3))))
(t_3 (- (* t y2) (* y y3)))
(t_4 (* c (* y0 (* x y2))))
(t_5 (- (* x y) (* z t)))
(t_6 (- (* a b) (* c i)))
(t_7 (* c (* i (* z t))))
(t_8 (* k (* y4 (* y b))))
(t_9 (* c (* y (* x i))))
(t_10 (* (- (* b y0) (* i y1)) (- (* z k) (* x j))))
(t_11 (* c (* y4 (* t y2))))
(t_12 (- (* a y5) (* c y4)))
(t_13 (* k (* y0 (* z b))))
(t_14 (* y4 (* t (* b j))))
(t_15 (* i (* t (* j y5))))
(t_16 (- (* t j) (* y k)))
(t_17 (- (* k y2) (* j y3)))
(t_18 (* c (* y0 (* z y3))))
(t_19 (* y1 (* i (* x j))))
(t_20 (* k (* y1 (* z i))))
(t_21 (- (* b y4) (* i y5)))
(t_22 (- (* y1 y4) (* y0 y5)))
(t_23
(+
(+ (+ (+ (+ (* t_5 t_6) t_10) t_1) (* t_16 t_21)) (* t_3 t_12))
(* t_17 t_22)))
(t_24 (* y0 (* j (* x b)))))
(if (<= t_23 -2e+297)
(fma
t_17
t_22
(-
(+ t_4 (+ (+ t_14 (+ t_19 (+ t_7 (+ t_2 (* k (* y (* i y5))))))) t_13))
(+ t_18 (+ t_24 (+ t_9 (+ t_11 (+ t_15 (+ t_20 t_8))))))))
(if (<= t_23 2e+300)
(fma
t_17
t_22
(fma t_3 t_12 (+ t_1 (fma t_16 t_21 (fma t_5 t_6 t_10)))))
(fma
t_17
t_22
(-
(+
t_13
(+
(* a (* t (* y2 y5)))
(+
t_14
(+
(* a (* y1 (* z y3)))
(+
t_2
(+
t_7
(+
t_4
(+ (* y (* a (* x b))) (+ t_19 (* k (* i (* y y5))))))))))))
(+
(* y (* a (* y3 y5)))
(+
t_18
(+
t_9
(+
(* a (* t (* z b)))
(+
t_11
(+
t_15
(+ t_8 (+ (* y1 (* a (* x y2))) (+ t_24 t_20)))))))))))))))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 * y2) - (z * y3)) * ((c * y0) - (a * y1));
double t_2 = c * (y4 * (y * y3));
double t_3 = (t * y2) - (y * y3);
double t_4 = c * (y0 * (x * y2));
double t_5 = (x * y) - (z * t);
double t_6 = (a * b) - (c * i);
double t_7 = c * (i * (z * t));
double t_8 = k * (y4 * (y * b));
double t_9 = c * (y * (x * i));
double t_10 = ((b * y0) - (i * y1)) * ((z * k) - (x * j));
double t_11 = c * (y4 * (t * y2));
double t_12 = (a * y5) - (c * y4);
double t_13 = k * (y0 * (z * b));
double t_14 = y4 * (t * (b * j));
double t_15 = i * (t * (j * y5));
double t_16 = (t * j) - (y * k);
double t_17 = (k * y2) - (j * y3);
double t_18 = c * (y0 * (z * y3));
double t_19 = y1 * (i * (x * j));
double t_20 = k * (y1 * (z * i));
double t_21 = (b * y4) - (i * y5);
double t_22 = (y1 * y4) - (y0 * y5);
double t_23 = (((((t_5 * t_6) + t_10) + t_1) + (t_16 * t_21)) + (t_3 * t_12)) + (t_17 * t_22);
double t_24 = y0 * (j * (x * b));
double tmp;
if (t_23 <= -2e+297) {
tmp = fma(t_17, t_22, ((t_4 + ((t_14 + (t_19 + (t_7 + (t_2 + (k * (y * (i * y5))))))) + t_13)) - (t_18 + (t_24 + (t_9 + (t_11 + (t_15 + (t_20 + t_8))))))));
} else if (t_23 <= 2e+300) {
tmp = fma(t_17, t_22, fma(t_3, t_12, (t_1 + fma(t_16, t_21, fma(t_5, t_6, t_10)))));
} else {
tmp = fma(t_17, t_22, ((t_13 + ((a * (t * (y2 * y5))) + (t_14 + ((a * (y1 * (z * y3))) + (t_2 + (t_7 + (t_4 + ((y * (a * (x * b))) + (t_19 + (k * (i * (y * y5)))))))))))) - ((y * (a * (y3 * y5))) + (t_18 + (t_9 + ((a * (t * (z * b))) + (t_11 + (t_15 + (t_8 + ((y1 * (a * (x * y2))) + (t_24 + t_20)))))))))));
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j, k, y0, y1, y2, y3, y4, y5) return Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * y) - Float64(z * t)) * Float64(Float64(a * b) - Float64(c * i))) - Float64(Float64(Float64(x * j) - Float64(z * k)) * Float64(Float64(y0 * b) - Float64(y1 * i)))) + Float64(Float64(Float64(x * y2) - Float64(z * y3)) * Float64(Float64(y0 * c) - Float64(y1 * a)))) + Float64(Float64(Float64(t * j) - Float64(y * k)) * Float64(Float64(y4 * b) - Float64(y5 * i)))) - Float64(Float64(Float64(t * y2) - Float64(y * y3)) * Float64(Float64(y4 * c) - Float64(y5 * a)))) + Float64(Float64(Float64(k * y2) - Float64(j * y3)) * Float64(Float64(y4 * y1) - Float64(y5 * y0)))) end
function code(x, y, z, t, a, b, c, i, j, k, y0, y1, y2, y3, y4, y5) t_1 = Float64(Float64(Float64(x * y2) - Float64(z * y3)) * Float64(Float64(c * y0) - Float64(a * y1))) t_2 = Float64(c * Float64(y4 * Float64(y * y3))) t_3 = Float64(Float64(t * y2) - Float64(y * y3)) t_4 = Float64(c * Float64(y0 * Float64(x * y2))) t_5 = Float64(Float64(x * y) - Float64(z * t)) t_6 = Float64(Float64(a * b) - Float64(c * i)) t_7 = Float64(c * Float64(i * Float64(z * t))) t_8 = Float64(k * Float64(y4 * Float64(y * b))) t_9 = Float64(c * Float64(y * Float64(x * i))) t_10 = Float64(Float64(Float64(b * y0) - Float64(i * y1)) * Float64(Float64(z * k) - Float64(x * j))) t_11 = Float64(c * Float64(y4 * Float64(t * y2))) t_12 = Float64(Float64(a * y5) - Float64(c * y4)) t_13 = Float64(k * Float64(y0 * Float64(z * b))) t_14 = Float64(y4 * Float64(t * Float64(b * j))) t_15 = Float64(i * Float64(t * Float64(j * y5))) t_16 = Float64(Float64(t * j) - Float64(y * k)) t_17 = Float64(Float64(k * y2) - Float64(j * y3)) t_18 = Float64(c * Float64(y0 * Float64(z * y3))) t_19 = Float64(y1 * Float64(i * Float64(x * j))) t_20 = Float64(k * Float64(y1 * Float64(z * i))) t_21 = Float64(Float64(b * y4) - Float64(i * y5)) t_22 = Float64(Float64(y1 * y4) - Float64(y0 * y5)) t_23 = Float64(Float64(Float64(Float64(Float64(Float64(t_5 * t_6) + t_10) + t_1) + Float64(t_16 * t_21)) + Float64(t_3 * t_12)) + Float64(t_17 * t_22)) t_24 = Float64(y0 * Float64(j * Float64(x * b))) tmp = 0.0 if (t_23 <= -2e+297) tmp = fma(t_17, t_22, Float64(Float64(t_4 + Float64(Float64(t_14 + Float64(t_19 + Float64(t_7 + Float64(t_2 + Float64(k * Float64(y * Float64(i * y5))))))) + t_13)) - Float64(t_18 + Float64(t_24 + Float64(t_9 + Float64(t_11 + Float64(t_15 + Float64(t_20 + t_8)))))))); elseif (t_23 <= 2e+300) tmp = fma(t_17, t_22, fma(t_3, t_12, Float64(t_1 + fma(t_16, t_21, fma(t_5, t_6, t_10))))); else tmp = fma(t_17, t_22, Float64(Float64(t_13 + Float64(Float64(a * Float64(t * Float64(y2 * y5))) + Float64(t_14 + Float64(Float64(a * Float64(y1 * Float64(z * y3))) + Float64(t_2 + Float64(t_7 + Float64(t_4 + Float64(Float64(y * Float64(a * Float64(x * b))) + Float64(t_19 + Float64(k * Float64(i * Float64(y * y5)))))))))))) - Float64(Float64(y * Float64(a * Float64(y3 * y5))) + Float64(t_18 + Float64(t_9 + Float64(Float64(a * Float64(t * Float64(z * b))) + Float64(t_11 + Float64(t_15 + Float64(t_8 + Float64(Float64(y1 * Float64(a * Float64(x * y2))) + Float64(t_24 + t_20))))))))))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_, y0_, y1_, y2_, y3_, y4_, y5_] := N[(N[(N[(N[(N[(N[(N[(N[(x * y), $MachinePrecision] - N[(z * t), $MachinePrecision]), $MachinePrecision] * N[(N[(a * b), $MachinePrecision] - N[(c * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(x * j), $MachinePrecision] - N[(z * k), $MachinePrecision]), $MachinePrecision] * N[(N[(y0 * b), $MachinePrecision] - N[(y1 * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(x * y2), $MachinePrecision] - N[(z * y3), $MachinePrecision]), $MachinePrecision] * N[(N[(y0 * c), $MachinePrecision] - N[(y1 * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(t * j), $MachinePrecision] - N[(y * k), $MachinePrecision]), $MachinePrecision] * N[(N[(y4 * b), $MachinePrecision] - N[(y5 * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(t * y2), $MachinePrecision] - N[(y * y3), $MachinePrecision]), $MachinePrecision] * N[(N[(y4 * c), $MachinePrecision] - N[(y5 * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(k * y2), $MachinePrecision] - N[(j * y3), $MachinePrecision]), $MachinePrecision] * N[(N[(y4 * y1), $MachinePrecision] - N[(y5 * y0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_, y0_, y1_, y2_, y3_, y4_, y5_] := Block[{t$95$1 = N[(N[(N[(x * y2), $MachinePrecision] - N[(z * y3), $MachinePrecision]), $MachinePrecision] * N[(N[(c * y0), $MachinePrecision] - N[(a * y1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(c * N[(y4 * N[(y * y3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(t * y2), $MachinePrecision] - N[(y * y3), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(c * N[(y0 * N[(x * y2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$5 = N[(N[(x * y), $MachinePrecision] - N[(z * t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$6 = N[(N[(a * b), $MachinePrecision] - N[(c * i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$7 = N[(c * N[(i * N[(z * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$8 = N[(k * N[(y4 * N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$9 = N[(c * N[(y * N[(x * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$10 = N[(N[(N[(b * y0), $MachinePrecision] - N[(i * y1), $MachinePrecision]), $MachinePrecision] * N[(N[(z * k), $MachinePrecision] - N[(x * j), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$11 = N[(c * N[(y4 * N[(t * y2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$12 = N[(N[(a * y5), $MachinePrecision] - N[(c * y4), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$13 = N[(k * N[(y0 * N[(z * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$14 = N[(y4 * N[(t * N[(b * j), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$15 = N[(i * N[(t * N[(j * y5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$16 = N[(N[(t * j), $MachinePrecision] - N[(y * k), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$17 = N[(N[(k * y2), $MachinePrecision] - N[(j * y3), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$18 = N[(c * N[(y0 * N[(z * y3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$19 = N[(y1 * N[(i * N[(x * j), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$20 = N[(k * N[(y1 * N[(z * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$21 = N[(N[(b * y4), $MachinePrecision] - N[(i * y5), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$22 = N[(N[(y1 * y4), $MachinePrecision] - N[(y0 * y5), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$23 = N[(N[(N[(N[(N[(N[(t$95$5 * t$95$6), $MachinePrecision] + t$95$10), $MachinePrecision] + t$95$1), $MachinePrecision] + N[(t$95$16 * t$95$21), $MachinePrecision]), $MachinePrecision] + N[(t$95$3 * t$95$12), $MachinePrecision]), $MachinePrecision] + N[(t$95$17 * t$95$22), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$24 = N[(y0 * N[(j * N[(x * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$23, -2e+297], N[(t$95$17 * t$95$22 + N[(N[(t$95$4 + N[(N[(t$95$14 + N[(t$95$19 + N[(t$95$7 + N[(t$95$2 + N[(k * N[(y * N[(i * y5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t$95$13), $MachinePrecision]), $MachinePrecision] - N[(t$95$18 + N[(t$95$24 + N[(t$95$9 + N[(t$95$11 + N[(t$95$15 + N[(t$95$20 + t$95$8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$23, 2e+300], N[(t$95$17 * t$95$22 + N[(t$95$3 * t$95$12 + N[(t$95$1 + N[(t$95$16 * t$95$21 + N[(t$95$5 * t$95$6 + t$95$10), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$17 * t$95$22 + N[(N[(t$95$13 + N[(N[(a * N[(t * N[(y2 * y5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$14 + N[(N[(a * N[(y1 * N[(z * y3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$2 + N[(t$95$7 + N[(t$95$4 + N[(N[(y * N[(a * N[(x * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$19 + N[(k * N[(i * N[(y * y5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(y * N[(a * N[(y3 * y5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$18 + N[(t$95$9 + N[(N[(a * N[(t * N[(z * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$11 + N[(t$95$15 + N[(t$95$8 + N[(N[(y1 * N[(a * N[(x * y2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$24 + t$95$20), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]]]]]]]]]]]]]]]]]]]]
\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 := \left(x \cdot y2 - z \cdot y3\right) \cdot \left(c \cdot y0 - a \cdot y1\right)\\
t_2 := c \cdot \left(y4 \cdot \left(y \cdot y3\right)\right)\\
t_3 := t \cdot y2 - y \cdot y3\\
t_4 := c \cdot \left(y0 \cdot \left(x \cdot y2\right)\right)\\
t_5 := x \cdot y - z \cdot t\\
t_6 := a \cdot b - c \cdot i\\
t_7 := c \cdot \left(i \cdot \left(z \cdot t\right)\right)\\
t_8 := k \cdot \left(y4 \cdot \left(y \cdot b\right)\right)\\
t_9 := c \cdot \left(y \cdot \left(x \cdot i\right)\right)\\
t_10 := \left(b \cdot y0 - i \cdot y1\right) \cdot \left(z \cdot k - x \cdot j\right)\\
t_11 := c \cdot \left(y4 \cdot \left(t \cdot y2\right)\right)\\
t_12 := a \cdot y5 - c \cdot y4\\
t_13 := k \cdot \left(y0 \cdot \left(z \cdot b\right)\right)\\
t_14 := y4 \cdot \left(t \cdot \left(b \cdot j\right)\right)\\
t_15 := i \cdot \left(t \cdot \left(j \cdot y5\right)\right)\\
t_16 := t \cdot j - y \cdot k\\
t_17 := k \cdot y2 - j \cdot y3\\
t_18 := c \cdot \left(y0 \cdot \left(z \cdot y3\right)\right)\\
t_19 := y1 \cdot \left(i \cdot \left(x \cdot j\right)\right)\\
t_20 := k \cdot \left(y1 \cdot \left(z \cdot i\right)\right)\\
t_21 := b \cdot y4 - i \cdot y5\\
t_22 := y1 \cdot y4 - y0 \cdot y5\\
t_23 := \left(\left(\left(\left(t_5 \cdot t_6 + t_10\right) + t_1\right) + t_16 \cdot t_21\right) + t_3 \cdot t_12\right) + t_17 \cdot t_22\\
t_24 := y0 \cdot \left(j \cdot \left(x \cdot b\right)\right)\\
\mathbf{if}\;t_23 \leq -2 \cdot 10^{+297}:\\
\;\;\;\;\mathsf{fma}\left(t_17, t_22, \left(t_4 + \left(\left(t_14 + \left(t_19 + \left(t_7 + \left(t_2 + k \cdot \left(y \cdot \left(i \cdot y5\right)\right)\right)\right)\right)\right) + t_13\right)\right) - \left(t_18 + \left(t_24 + \left(t_9 + \left(t_11 + \left(t_15 + \left(t_20 + t_8\right)\right)\right)\right)\right)\right)\right)\\
\mathbf{elif}\;t_23 \leq 2 \cdot 10^{+300}:\\
\;\;\;\;\mathsf{fma}\left(t_17, t_22, \mathsf{fma}\left(t_3, t_12, t_1 + \mathsf{fma}\left(t_16, t_21, \mathsf{fma}\left(t_5, t_6, t_10\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t_17, t_22, \left(t_13 + \left(a \cdot \left(t \cdot \left(y2 \cdot y5\right)\right) + \left(t_14 + \left(a \cdot \left(y1 \cdot \left(z \cdot y3\right)\right) + \left(t_2 + \left(t_7 + \left(t_4 + \left(y \cdot \left(a \cdot \left(x \cdot b\right)\right) + \left(t_19 + k \cdot \left(i \cdot \left(y \cdot y5\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\right) - \left(y \cdot \left(a \cdot \left(y3 \cdot y5\right)\right) + \left(t_18 + \left(t_9 + \left(a \cdot \left(t \cdot \left(z \cdot b\right)\right) + \left(t_11 + \left(t_15 + \left(t_8 + \left(y1 \cdot \left(a \cdot \left(x \cdot y2\right)\right) + \left(t_24 + t_20\right)\right)\right)\right)\right)\right)\right)\right)\right)\right)\\
\end{array}




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 | 31.2 |
| Herbie | 20.9 |
if (+.f64 (-.f64 (+.f64 (+.f64 (-.f64 (*.f64 (-.f64 (*.f64 x y) (*.f64 z t)) (-.f64 (*.f64 a b) (*.f64 c i))) (*.f64 (-.f64 (*.f64 x j) (*.f64 z k)) (-.f64 (*.f64 y0 b) (*.f64 y1 i)))) (*.f64 (-.f64 (*.f64 x y2) (*.f64 z y3)) (-.f64 (*.f64 y0 c) (*.f64 y1 a)))) (*.f64 (-.f64 (*.f64 t j) (*.f64 y k)) (-.f64 (*.f64 y4 b) (*.f64 y5 i)))) (*.f64 (-.f64 (*.f64 t y2) (*.f64 y y3)) (-.f64 (*.f64 y4 c) (*.f64 y5 a)))) (*.f64 (-.f64 (*.f64 k y2) (*.f64 j y3)) (-.f64 (*.f64 y4 y1) (*.f64 y5 y0)))) < -2e297Initial program 58.4
Simplified58.4
Taylor expanded in a around 0 46.4
if -2e297 < (+.f64 (-.f64 (+.f64 (+.f64 (-.f64 (*.f64 (-.f64 (*.f64 x y) (*.f64 z t)) (-.f64 (*.f64 a b) (*.f64 c i))) (*.f64 (-.f64 (*.f64 x j) (*.f64 z k)) (-.f64 (*.f64 y0 b) (*.f64 y1 i)))) (*.f64 (-.f64 (*.f64 x y2) (*.f64 z y3)) (-.f64 (*.f64 y0 c) (*.f64 y1 a)))) (*.f64 (-.f64 (*.f64 t j) (*.f64 y k)) (-.f64 (*.f64 y4 b) (*.f64 y5 i)))) (*.f64 (-.f64 (*.f64 t y2) (*.f64 y y3)) (-.f64 (*.f64 y4 c) (*.f64 y5 a)))) (*.f64 (-.f64 (*.f64 k y2) (*.f64 j y3)) (-.f64 (*.f64 y4 y1) (*.f64 y5 y0)))) < 2.0000000000000001e300Initial program 0.7
Simplified0.7
Applied egg-rr0.7
if 2.0000000000000001e300 < (+.f64 (-.f64 (+.f64 (+.f64 (-.f64 (*.f64 (-.f64 (*.f64 x y) (*.f64 z t)) (-.f64 (*.f64 a b) (*.f64 c i))) (*.f64 (-.f64 (*.f64 x j) (*.f64 z k)) (-.f64 (*.f64 y0 b) (*.f64 y1 i)))) (*.f64 (-.f64 (*.f64 x y2) (*.f64 z y3)) (-.f64 (*.f64 y0 c) (*.f64 y1 a)))) (*.f64 (-.f64 (*.f64 t j) (*.f64 y k)) (-.f64 (*.f64 y4 b) (*.f64 y5 i)))) (*.f64 (-.f64 (*.f64 t y2) (*.f64 y y3)) (-.f64 (*.f64 y4 c) (*.f64 y5 a)))) (*.f64 (-.f64 (*.f64 k y2) (*.f64 j y3)) (-.f64 (*.f64 y4 y1) (*.f64 y5 y0)))) Initial program 61.4
Simplified61.4
Taylor expanded in t around 0 46.2
Final simplification20.9
herbie shell --seed 2022166
(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)))))