(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 (- (* a y5) (* c y4)))
(t_2 (- (* t j) (* y k)))
(t_3 (fma x y2 (* z (- y3))))
(t_4 (fma a b (- (* c i))))
(t_5 (* y4 (- (* y y3) (* y2 t))))
(t_6 (- (* y x) (* z t)))
(t_7 (- (* x y2) (* z y3)))
(t_8 (- (* x j) (* z k)))
(t_9 (- (* y2 t) (* y y3)))
(t_10 (- (* c y0) (* a y1)))
(t_11 (* y4 (* y3 y1)))
(t_12 (- (* a b) (* c i)))
(t_13 (* x t_4))
(t_14 (- (* b y4) (* i y5)))
(t_15 (- (* y4 y1) (* y5 y0)))
(t_16 (* y2 t_15))
(t_17 (- (* i y1) (* b y0)))
(t_18
(fma
y
t_13
(fma
t_1
t_9
(fma
k
t_16
(-
(* t_10 (* x y2))
(* z (+ (* t_4 t) (fma t_10 y3 (* k t_17)))))))))
(t_19 (- (* y3 t_15) (fma t t_14 (* x t_17))))
(t_20 (* j t_19))
(t_21 (* z t_17))
(t_22 (- (* y2 k) (* y3 j)))
(t_23
(+
(* y0 (* b (- (* z k) (* x j))))
(+
(* t_6 t_12)
(+
(* t_14 t_2)
(+
(* y0 (* y5 (- (* y3 j) (* y2 k))))
(+
(* t_1 t_9)
(+
(* y1 (+ (* a (- (* z y3) (* x y2))) (+ (* y4 t_22) (* i t_8))))
(* c (* y0 t_7)))))))))
(t_24
(fma
t_22
t_15
(fma
t_9
t_1
(fma
t_7
t_10
(+
(* i (* y5 (- (* y k) (* t j))))
(+
(* b (+ (* a t_6) (- (* y4 t_2) (* y0 t_8))))
(+ (* i (* y1 t_8)) (* c (* i (- (* z t) (* y x))))))))))))
(if (<= z -0.01)
(- t_18 (fma j t_19 (* k (* y t_14))))
(if (<= z -1.6314152204079688e-97)
(-
(fma
c
(* y0 t_3)
(fma
k
t_16
(+
(* c (- t_5 (* i t_6)))
(* a (- (fma t_9 y5 (* b t_6)) (* y1 t_3))))))
(fma k (fma y t_14 t_21) t_20))
(if (<= z -4.556730256669597e-143)
(-
(-
(fma
y0
(-
(fma c t_3 (* k (* y5 (- y2))))
(fma j (fma (- y3) y5 (* x b)) (* (* z b) (- k))))
(fma k (* y2 (* y4 y1)) (fma t_1 t_9 (* t_4 t_6))))
(* y1 (* a t_3)))
(fma
k
(fma y t_14 (* y1 (* z i)))
(* j (- t_11 (fma y1 (* x i) (* t t_14))))))
(if (<= z -5.3061788076370425e-193)
t_24
(if (<= z -6.363787525205459e-307)
t_23
(if (<= z 5.204480699894442e-279)
t_24
(if (<= z 2.9764907766571937e-144)
(fma
t_12
(fma x y (* z (- t)))
(fma
t_17
t_8
(fma
y4
(* b t_2)
(fma
t_7
t_10
(+
(+
(* y4 (* y1 t_22))
(* y5 (- (* a t_9) (fma y0 t_22 (* i t_2)))))
(* c t_5))))))
(if (<= z 4.284117068250017e-124)
t_23
(if (<= z 5e-52)
(-
(-
(fma
k
t_16
(fma
t_1
(* y2 t)
(fma
t_10
t_3
(*
y
(+
(+ t_13 (* y3 (- (* c y4) (* a y5))))
(* k (- (* i y5) (* b y4))))))))
(* t (* z t_4)))
(fma k t_21 t_20))
(if (<= z 1e+72)
t_24
(+
t_18
(+
(* y5 (+ (* k (* y i)) (* j (- (* y3 y0) (* i t)))))
(-
(*
j
(-
(- (* y4 (* b t)) (* x (- (* b y0) (* i y1))))
t_11))
(* k (* y4 (* y b))))))))))))))))))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 = (a * y5) - (c * y4);
double t_2 = (t * j) - (y * k);
double t_3 = fma(x, y2, (z * -y3));
double t_4 = fma(a, b, -(c * i));
double t_5 = y4 * ((y * y3) - (y2 * t));
double t_6 = (y * x) - (z * t);
double t_7 = (x * y2) - (z * y3);
double t_8 = (x * j) - (z * k);
double t_9 = (y2 * t) - (y * y3);
double t_10 = (c * y0) - (a * y1);
double t_11 = y4 * (y3 * y1);
double t_12 = (a * b) - (c * i);
double t_13 = x * t_4;
double t_14 = (b * y4) - (i * y5);
double t_15 = (y4 * y1) - (y5 * y0);
double t_16 = y2 * t_15;
double t_17 = (i * y1) - (b * y0);
double t_18 = fma(y, t_13, fma(t_1, t_9, fma(k, t_16, ((t_10 * (x * y2)) - (z * ((t_4 * t) + fma(t_10, y3, (k * t_17))))))));
double t_19 = (y3 * t_15) - fma(t, t_14, (x * t_17));
double t_20 = j * t_19;
double t_21 = z * t_17;
double t_22 = (y2 * k) - (y3 * j);
double t_23 = (y0 * (b * ((z * k) - (x * j)))) + ((t_6 * t_12) + ((t_14 * t_2) + ((y0 * (y5 * ((y3 * j) - (y2 * k)))) + ((t_1 * t_9) + ((y1 * ((a * ((z * y3) - (x * y2))) + ((y4 * t_22) + (i * t_8)))) + (c * (y0 * t_7)))))));
double t_24 = fma(t_22, t_15, fma(t_9, t_1, fma(t_7, t_10, ((i * (y5 * ((y * k) - (t * j)))) + ((b * ((a * t_6) + ((y4 * t_2) - (y0 * t_8)))) + ((i * (y1 * t_8)) + (c * (i * ((z * t) - (y * x))))))))));
double tmp;
if (z <= -0.01) {
tmp = t_18 - fma(j, t_19, (k * (y * t_14)));
} else if (z <= -1.6314152204079688e-97) {
tmp = fma(c, (y0 * t_3), fma(k, t_16, ((c * (t_5 - (i * t_6))) + (a * (fma(t_9, y5, (b * t_6)) - (y1 * t_3)))))) - fma(k, fma(y, t_14, t_21), t_20);
} else if (z <= -4.556730256669597e-143) {
tmp = (fma(y0, (fma(c, t_3, (k * (y5 * -y2))) - fma(j, fma(-y3, y5, (x * b)), ((z * b) * -k))), fma(k, (y2 * (y4 * y1)), fma(t_1, t_9, (t_4 * t_6)))) - (y1 * (a * t_3))) - fma(k, fma(y, t_14, (y1 * (z * i))), (j * (t_11 - fma(y1, (x * i), (t * t_14)))));
} else if (z <= -5.3061788076370425e-193) {
tmp = t_24;
} else if (z <= -6.363787525205459e-307) {
tmp = t_23;
} else if (z <= 5.204480699894442e-279) {
tmp = t_24;
} else if (z <= 2.9764907766571937e-144) {
tmp = fma(t_12, fma(x, y, (z * -t)), fma(t_17, t_8, fma(y4, (b * t_2), fma(t_7, t_10, (((y4 * (y1 * t_22)) + (y5 * ((a * t_9) - fma(y0, t_22, (i * t_2))))) + (c * t_5))))));
} else if (z <= 4.284117068250017e-124) {
tmp = t_23;
} else if (z <= 5e-52) {
tmp = (fma(k, t_16, fma(t_1, (y2 * t), fma(t_10, t_3, (y * ((t_13 + (y3 * ((c * y4) - (a * y5)))) + (k * ((i * y5) - (b * y4)))))))) - (t * (z * t_4))) - fma(k, t_21, t_20);
} else if (z <= 1e+72) {
tmp = t_24;
} else {
tmp = t_18 + ((y5 * ((k * (y * i)) + (j * ((y3 * y0) - (i * t))))) + ((j * (((y4 * (b * t)) - (x * ((b * y0) - (i * y1)))) - t_11)) - (k * (y4 * (y * b)))));
}
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(a * y5) - Float64(c * y4)) t_2 = Float64(Float64(t * j) - Float64(y * k)) t_3 = fma(x, y2, Float64(z * Float64(-y3))) t_4 = fma(a, b, Float64(-Float64(c * i))) t_5 = Float64(y4 * Float64(Float64(y * y3) - Float64(y2 * t))) t_6 = Float64(Float64(y * x) - Float64(z * t)) t_7 = Float64(Float64(x * y2) - Float64(z * y3)) t_8 = Float64(Float64(x * j) - Float64(z * k)) t_9 = Float64(Float64(y2 * t) - Float64(y * y3)) t_10 = Float64(Float64(c * y0) - Float64(a * y1)) t_11 = Float64(y4 * Float64(y3 * y1)) t_12 = Float64(Float64(a * b) - Float64(c * i)) t_13 = Float64(x * t_4) t_14 = Float64(Float64(b * y4) - Float64(i * y5)) t_15 = Float64(Float64(y4 * y1) - Float64(y5 * y0)) t_16 = Float64(y2 * t_15) t_17 = Float64(Float64(i * y1) - Float64(b * y0)) t_18 = fma(y, t_13, fma(t_1, t_9, fma(k, t_16, Float64(Float64(t_10 * Float64(x * y2)) - Float64(z * Float64(Float64(t_4 * t) + fma(t_10, y3, Float64(k * t_17)))))))) t_19 = Float64(Float64(y3 * t_15) - fma(t, t_14, Float64(x * t_17))) t_20 = Float64(j * t_19) t_21 = Float64(z * t_17) t_22 = Float64(Float64(y2 * k) - Float64(y3 * j)) t_23 = Float64(Float64(y0 * Float64(b * Float64(Float64(z * k) - Float64(x * j)))) + Float64(Float64(t_6 * t_12) + Float64(Float64(t_14 * t_2) + Float64(Float64(y0 * Float64(y5 * Float64(Float64(y3 * j) - Float64(y2 * k)))) + Float64(Float64(t_1 * t_9) + Float64(Float64(y1 * Float64(Float64(a * Float64(Float64(z * y3) - Float64(x * y2))) + Float64(Float64(y4 * t_22) + Float64(i * t_8)))) + Float64(c * Float64(y0 * t_7)))))))) t_24 = fma(t_22, t_15, fma(t_9, t_1, fma(t_7, t_10, Float64(Float64(i * Float64(y5 * Float64(Float64(y * k) - Float64(t * j)))) + Float64(Float64(b * Float64(Float64(a * t_6) + Float64(Float64(y4 * t_2) - Float64(y0 * t_8)))) + Float64(Float64(i * Float64(y1 * t_8)) + Float64(c * Float64(i * Float64(Float64(z * t) - Float64(y * x)))))))))) tmp = 0.0 if (z <= -0.01) tmp = Float64(t_18 - fma(j, t_19, Float64(k * Float64(y * t_14)))); elseif (z <= -1.6314152204079688e-97) tmp = Float64(fma(c, Float64(y0 * t_3), fma(k, t_16, Float64(Float64(c * Float64(t_5 - Float64(i * t_6))) + Float64(a * Float64(fma(t_9, y5, Float64(b * t_6)) - Float64(y1 * t_3)))))) - fma(k, fma(y, t_14, t_21), t_20)); elseif (z <= -4.556730256669597e-143) tmp = Float64(Float64(fma(y0, Float64(fma(c, t_3, Float64(k * Float64(y5 * Float64(-y2)))) - fma(j, fma(Float64(-y3), y5, Float64(x * b)), Float64(Float64(z * b) * Float64(-k)))), fma(k, Float64(y2 * Float64(y4 * y1)), fma(t_1, t_9, Float64(t_4 * t_6)))) - Float64(y1 * Float64(a * t_3))) - fma(k, fma(y, t_14, Float64(y1 * Float64(z * i))), Float64(j * Float64(t_11 - fma(y1, Float64(x * i), Float64(t * t_14)))))); elseif (z <= -5.3061788076370425e-193) tmp = t_24; elseif (z <= -6.363787525205459e-307) tmp = t_23; elseif (z <= 5.204480699894442e-279) tmp = t_24; elseif (z <= 2.9764907766571937e-144) tmp = fma(t_12, fma(x, y, Float64(z * Float64(-t))), fma(t_17, t_8, fma(y4, Float64(b * t_2), fma(t_7, t_10, Float64(Float64(Float64(y4 * Float64(y1 * t_22)) + Float64(y5 * Float64(Float64(a * t_9) - fma(y0, t_22, Float64(i * t_2))))) + Float64(c * t_5)))))); elseif (z <= 4.284117068250017e-124) tmp = t_23; elseif (z <= 5e-52) tmp = Float64(Float64(fma(k, t_16, fma(t_1, Float64(y2 * t), fma(t_10, t_3, Float64(y * Float64(Float64(t_13 + Float64(y3 * Float64(Float64(c * y4) - Float64(a * y5)))) + Float64(k * Float64(Float64(i * y5) - Float64(b * y4)))))))) - Float64(t * Float64(z * t_4))) - fma(k, t_21, t_20)); elseif (z <= 1e+72) tmp = t_24; else tmp = Float64(t_18 + Float64(Float64(y5 * Float64(Float64(k * Float64(y * i)) + Float64(j * Float64(Float64(y3 * y0) - Float64(i * t))))) + Float64(Float64(j * Float64(Float64(Float64(y4 * Float64(b * t)) - Float64(x * Float64(Float64(b * y0) - Float64(i * y1)))) - t_11)) - Float64(k * Float64(y4 * Float64(y * b)))))); 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[(a * y5), $MachinePrecision] - N[(c * y4), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(t * j), $MachinePrecision] - N[(y * k), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(x * y2 + N[(z * (-y3)), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(a * b + (-N[(c * i), $MachinePrecision])), $MachinePrecision]}, Block[{t$95$5 = N[(y4 * N[(N[(y * y3), $MachinePrecision] - N[(y2 * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$6 = N[(N[(y * x), $MachinePrecision] - N[(z * t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$7 = N[(N[(x * y2), $MachinePrecision] - N[(z * y3), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$8 = N[(N[(x * j), $MachinePrecision] - N[(z * k), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$9 = N[(N[(y2 * t), $MachinePrecision] - N[(y * y3), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$10 = N[(N[(c * y0), $MachinePrecision] - N[(a * y1), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$11 = N[(y4 * N[(y3 * y1), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$12 = N[(N[(a * b), $MachinePrecision] - N[(c * i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$13 = N[(x * t$95$4), $MachinePrecision]}, Block[{t$95$14 = N[(N[(b * y4), $MachinePrecision] - N[(i * y5), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$15 = N[(N[(y4 * y1), $MachinePrecision] - N[(y5 * y0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$16 = N[(y2 * t$95$15), $MachinePrecision]}, Block[{t$95$17 = N[(N[(i * y1), $MachinePrecision] - N[(b * y0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$18 = N[(y * t$95$13 + N[(t$95$1 * t$95$9 + N[(k * t$95$16 + N[(N[(t$95$10 * N[(x * y2), $MachinePrecision]), $MachinePrecision] - N[(z * N[(N[(t$95$4 * t), $MachinePrecision] + N[(t$95$10 * y3 + N[(k * t$95$17), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$19 = N[(N[(y3 * t$95$15), $MachinePrecision] - N[(t * t$95$14 + N[(x * t$95$17), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$20 = N[(j * t$95$19), $MachinePrecision]}, Block[{t$95$21 = N[(z * t$95$17), $MachinePrecision]}, Block[{t$95$22 = N[(N[(y2 * k), $MachinePrecision] - N[(y3 * j), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$23 = N[(N[(y0 * N[(b * N[(N[(z * k), $MachinePrecision] - N[(x * j), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(t$95$6 * t$95$12), $MachinePrecision] + N[(N[(t$95$14 * t$95$2), $MachinePrecision] + N[(N[(y0 * N[(y5 * N[(N[(y3 * j), $MachinePrecision] - N[(y2 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(t$95$1 * t$95$9), $MachinePrecision] + N[(N[(y1 * N[(N[(a * N[(N[(z * y3), $MachinePrecision] - N[(x * y2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(y4 * t$95$22), $MachinePrecision] + N[(i * t$95$8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(c * N[(y0 * t$95$7), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$24 = N[(t$95$22 * t$95$15 + N[(t$95$9 * t$95$1 + N[(t$95$7 * t$95$10 + N[(N[(i * N[(y5 * N[(N[(y * k), $MachinePrecision] - N[(t * j), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(b * N[(N[(a * t$95$6), $MachinePrecision] + N[(N[(y4 * t$95$2), $MachinePrecision] - N[(y0 * t$95$8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(i * N[(y1 * t$95$8), $MachinePrecision]), $MachinePrecision] + N[(c * N[(i * N[(N[(z * t), $MachinePrecision] - N[(y * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -0.01], N[(t$95$18 - N[(j * t$95$19 + N[(k * N[(y * t$95$14), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -1.6314152204079688e-97], N[(N[(c * N[(y0 * t$95$3), $MachinePrecision] + N[(k * t$95$16 + N[(N[(c * N[(t$95$5 - N[(i * t$95$6), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a * N[(N[(t$95$9 * y5 + N[(b * t$95$6), $MachinePrecision]), $MachinePrecision] - N[(y1 * t$95$3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(k * N[(y * t$95$14 + t$95$21), $MachinePrecision] + t$95$20), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -4.556730256669597e-143], N[(N[(N[(y0 * N[(N[(c * t$95$3 + N[(k * N[(y5 * (-y2)), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(j * N[((-y3) * y5 + N[(x * b), $MachinePrecision]), $MachinePrecision] + N[(N[(z * b), $MachinePrecision] * (-k)), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(k * N[(y2 * N[(y4 * y1), $MachinePrecision]), $MachinePrecision] + N[(t$95$1 * t$95$9 + N[(t$95$4 * t$95$6), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(y1 * N[(a * t$95$3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(k * N[(y * t$95$14 + N[(y1 * N[(z * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(j * N[(t$95$11 - N[(y1 * N[(x * i), $MachinePrecision] + N[(t * t$95$14), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -5.3061788076370425e-193], t$95$24, If[LessEqual[z, -6.363787525205459e-307], t$95$23, If[LessEqual[z, 5.204480699894442e-279], t$95$24, If[LessEqual[z, 2.9764907766571937e-144], N[(t$95$12 * N[(x * y + N[(z * (-t)), $MachinePrecision]), $MachinePrecision] + N[(t$95$17 * t$95$8 + N[(y4 * N[(b * t$95$2), $MachinePrecision] + N[(t$95$7 * t$95$10 + N[(N[(N[(y4 * N[(y1 * t$95$22), $MachinePrecision]), $MachinePrecision] + N[(y5 * N[(N[(a * t$95$9), $MachinePrecision] - N[(y0 * t$95$22 + N[(i * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(c * t$95$5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 4.284117068250017e-124], t$95$23, If[LessEqual[z, 5e-52], N[(N[(N[(k * t$95$16 + N[(t$95$1 * N[(y2 * t), $MachinePrecision] + N[(t$95$10 * t$95$3 + N[(y * N[(N[(t$95$13 + N[(y3 * N[(N[(c * y4), $MachinePrecision] - N[(a * y5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(k * N[(N[(i * y5), $MachinePrecision] - N[(b * y4), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(t * N[(z * t$95$4), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(k * t$95$21 + t$95$20), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 1e+72], t$95$24, N[(t$95$18 + N[(N[(y5 * N[(N[(k * N[(y * i), $MachinePrecision]), $MachinePrecision] + N[(j * N[(N[(y3 * y0), $MachinePrecision] - N[(i * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(j * N[(N[(N[(y4 * N[(b * t), $MachinePrecision]), $MachinePrecision] - N[(x * N[(N[(b * y0), $MachinePrecision] - N[(i * y1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t$95$11), $MachinePrecision]), $MachinePrecision] - N[(k * N[(y4 * N[(y * b), $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 := a \cdot y5 - c \cdot y4\\
t_2 := t \cdot j - y \cdot k\\
t_3 := \mathsf{fma}\left(x, y2, z \cdot \left(-y3\right)\right)\\
t_4 := \mathsf{fma}\left(a, b, -c \cdot i\right)\\
t_5 := y4 \cdot \left(y \cdot y3 - y2 \cdot t\right)\\
t_6 := y \cdot x - z \cdot t\\
t_7 := x \cdot y2 - z \cdot y3\\
t_8 := x \cdot j - z \cdot k\\
t_9 := y2 \cdot t - y \cdot y3\\
t_10 := c \cdot y0 - a \cdot y1\\
t_11 := y4 \cdot \left(y3 \cdot y1\right)\\
t_12 := a \cdot b - c \cdot i\\
t_13 := x \cdot t_4\\
t_14 := b \cdot y4 - i \cdot y5\\
t_15 := y4 \cdot y1 - y5 \cdot y0\\
t_16 := y2 \cdot t_15\\
t_17 := i \cdot y1 - b \cdot y0\\
t_18 := \mathsf{fma}\left(y, t_13, \mathsf{fma}\left(t_1, t_9, \mathsf{fma}\left(k, t_16, t_10 \cdot \left(x \cdot y2\right) - z \cdot \left(t_4 \cdot t + \mathsf{fma}\left(t_10, y3, k \cdot t_17\right)\right)\right)\right)\right)\\
t_19 := y3 \cdot t_15 - \mathsf{fma}\left(t, t_14, x \cdot t_17\right)\\
t_20 := j \cdot t_19\\
t_21 := z \cdot t_17\\
t_22 := y2 \cdot k - y3 \cdot j\\
t_23 := y0 \cdot \left(b \cdot \left(z \cdot k - x \cdot j\right)\right) + \left(t_6 \cdot t_12 + \left(t_14 \cdot t_2 + \left(y0 \cdot \left(y5 \cdot \left(y3 \cdot j - y2 \cdot k\right)\right) + \left(t_1 \cdot t_9 + \left(y1 \cdot \left(a \cdot \left(z \cdot y3 - x \cdot y2\right) + \left(y4 \cdot t_22 + i \cdot t_8\right)\right) + c \cdot \left(y0 \cdot t_7\right)\right)\right)\right)\right)\right)\\
t_24 := \mathsf{fma}\left(t_22, t_15, \mathsf{fma}\left(t_9, t_1, \mathsf{fma}\left(t_7, t_10, i \cdot \left(y5 \cdot \left(y \cdot k - t \cdot j\right)\right) + \left(b \cdot \left(a \cdot t_6 + \left(y4 \cdot t_2 - y0 \cdot t_8\right)\right) + \left(i \cdot \left(y1 \cdot t_8\right) + c \cdot \left(i \cdot \left(z \cdot t - y \cdot x\right)\right)\right)\right)\right)\right)\right)\\
\mathbf{if}\;z \leq -0.01:\\
\;\;\;\;t_18 - \mathsf{fma}\left(j, t_19, k \cdot \left(y \cdot t_14\right)\right)\\
\mathbf{elif}\;z \leq -1.6314152204079688 \cdot 10^{-97}:\\
\;\;\;\;\mathsf{fma}\left(c, y0 \cdot t_3, \mathsf{fma}\left(k, t_16, c \cdot \left(t_5 - i \cdot t_6\right) + a \cdot \left(\mathsf{fma}\left(t_9, y5, b \cdot t_6\right) - y1 \cdot t_3\right)\right)\right) - \mathsf{fma}\left(k, \mathsf{fma}\left(y, t_14, t_21\right), t_20\right)\\
\mathbf{elif}\;z \leq -4.556730256669597 \cdot 10^{-143}:\\
\;\;\;\;\left(\mathsf{fma}\left(y0, \mathsf{fma}\left(c, t_3, k \cdot \left(y5 \cdot \left(-y2\right)\right)\right) - \mathsf{fma}\left(j, \mathsf{fma}\left(-y3, y5, x \cdot b\right), \left(z \cdot b\right) \cdot \left(-k\right)\right), \mathsf{fma}\left(k, y2 \cdot \left(y4 \cdot y1\right), \mathsf{fma}\left(t_1, t_9, t_4 \cdot t_6\right)\right)\right) - y1 \cdot \left(a \cdot t_3\right)\right) - \mathsf{fma}\left(k, \mathsf{fma}\left(y, t_14, y1 \cdot \left(z \cdot i\right)\right), j \cdot \left(t_11 - \mathsf{fma}\left(y1, x \cdot i, t \cdot t_14\right)\right)\right)\\
\mathbf{elif}\;z \leq -5.3061788076370425 \cdot 10^{-193}:\\
\;\;\;\;t_24\\
\mathbf{elif}\;z \leq -6.363787525205459 \cdot 10^{-307}:\\
\;\;\;\;t_23\\
\mathbf{elif}\;z \leq 5.204480699894442 \cdot 10^{-279}:\\
\;\;\;\;t_24\\
\mathbf{elif}\;z \leq 2.9764907766571937 \cdot 10^{-144}:\\
\;\;\;\;\mathsf{fma}\left(t_12, \mathsf{fma}\left(x, y, z \cdot \left(-t\right)\right), \mathsf{fma}\left(t_17, t_8, \mathsf{fma}\left(y4, b \cdot t_2, \mathsf{fma}\left(t_7, t_10, \left(y4 \cdot \left(y1 \cdot t_22\right) + y5 \cdot \left(a \cdot t_9 - \mathsf{fma}\left(y0, t_22, i \cdot t_2\right)\right)\right) + c \cdot t_5\right)\right)\right)\right)\\
\mathbf{elif}\;z \leq 4.284117068250017 \cdot 10^{-124}:\\
\;\;\;\;t_23\\
\mathbf{elif}\;z \leq 5 \cdot 10^{-52}:\\
\;\;\;\;\left(\mathsf{fma}\left(k, t_16, \mathsf{fma}\left(t_1, y2 \cdot t, \mathsf{fma}\left(t_10, t_3, y \cdot \left(\left(t_13 + y3 \cdot \left(c \cdot y4 - a \cdot y5\right)\right) + k \cdot \left(i \cdot y5 - b \cdot y4\right)\right)\right)\right)\right) - t \cdot \left(z \cdot t_4\right)\right) - \mathsf{fma}\left(k, t_21, t_20\right)\\
\mathbf{elif}\;z \leq 10^{+72}:\\
\;\;\;\;t_24\\
\mathbf{else}:\\
\;\;\;\;t_18 + \left(y5 \cdot \left(k \cdot \left(y \cdot i\right) + j \cdot \left(y3 \cdot y0 - i \cdot t\right)\right) + \left(j \cdot \left(\left(y4 \cdot \left(b \cdot t\right) - x \cdot \left(b \cdot y0 - i \cdot y1\right)\right) - t_11\right) - k \cdot \left(y4 \cdot \left(y \cdot b\right)\right)\right)\right)\\
\end{array}
| Original | 27.2 |
|---|---|
| Target | 31.3 |
| Herbie | 24.6 |
if z < -0.0100000000000000002Initial program 29.7
Simplified29.7
Taylor expanded in j around -inf 29.1
Simplified29.1
Taylor expanded in z around -inf 21.3
Simplified21.3
if -0.0100000000000000002 < z < -1.63141522040796883e-97Initial program 23.0
Simplified23.0
Taylor expanded in j around -inf 22.9
Simplified22.9
Taylor expanded in a around 0 22.2
Simplified22.2
if -1.63141522040796883e-97 < z < -4.556730256669597e-143Initial program 26.6
Simplified26.6
Taylor expanded in j around -inf 26.0
Simplified26.0
Taylor expanded in y0 around inf 28.7
Simplified29.7
if -4.556730256669597e-143 < z < -5.3061788076370425e-193 or -6.36378752520545925e-307 < z < 5.20448069989444205e-279 or 5e-52 < z < 9.99999999999999944e71Initial program 25.4
Simplified25.4
Taylor expanded in b around 0 26.7
if -5.3061788076370425e-193 < z < -6.36378752520545925e-307 or 2.9764907766571937e-144 < z < 4.2841170682500172e-124Initial program 28.1
Simplified28.1
Taylor expanded in y1 around 0 24.7
if 5.20448069989444205e-279 < z < 2.9764907766571937e-144Initial program 27.3
Simplified27.3
Taylor expanded in y5 around -inf 29.3
Simplified29.3
if 4.2841170682500172e-124 < z < 5e-52Initial program 24.1
Simplified24.1
Taylor expanded in j around -inf 23.4
Simplified23.4
Taylor expanded in y around 0 24.5
Simplified24.5
if 9.99999999999999944e71 < z Initial program 29.9
Simplified29.9
Taylor expanded in j around -inf 30.6
Simplified30.6
Taylor expanded in z around -inf 18.0
Simplified18.0
Taylor expanded in y5 around -inf 20.3
Final simplification24.6
herbie shell --seed 2022210
(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)))))