\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}
\mathbf{if}\;i \le -1.885291032353528 \cdot 10^{+143}:\\
\;\;\;\;\left(\left(\left(\left(y0 \cdot c - a \cdot y1\right) \cdot \left(x \cdot y2 - y3 \cdot z\right) + \left(\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - c \cdot i\right) - \left(i \cdot \left(z \cdot \left(k \cdot y1\right)\right) - \left(i \cdot \left(\left(x \cdot y1\right) \cdot j\right) + \left(\left(b \cdot y0\right) \cdot z\right) \cdot k\right)\right)\right)\right) + \left(b \cdot y4 - y5 \cdot i\right) \cdot \left(j \cdot t - y \cdot k\right)\right) - \left(c \cdot y4 - a \cdot y5\right) \cdot \left(t \cdot y2 - y3 \cdot y\right)\right) + \left(y2 \cdot k - y3 \cdot j\right) \cdot \left(y1 \cdot y4 - y5 \cdot y0\right)\\
\mathbf{elif}\;i \le -7.571795248650015 \cdot 10^{-30}:\\
\;\;\;\;\left(y2 \cdot k - y3 \cdot j\right) \cdot \left(y1 \cdot y4 - y5 \cdot y0\right) + \left(\left(\left(b \cdot y4 - y5 \cdot i\right) \cdot \left(j \cdot t - y \cdot k\right) + \left(\left(\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - c \cdot i\right) - \left(x \cdot j - k \cdot z\right) \cdot \left(b \cdot y0 - i \cdot y1\right)\right) + \left(y0 \cdot c - a \cdot y1\right) \cdot \left(x \cdot y2 - y3 \cdot z\right)\right)\right) - \left(a \cdot \left(y3 \cdot \left(y \cdot y5\right)\right) - \left(a \cdot \left(y2 \cdot \left(y5 \cdot t\right)\right) + \left(y4 \cdot \left(y \cdot c\right)\right) \cdot y3\right)\right)\right)\\
\mathbf{elif}\;i \le -4.0262567889326796 \cdot 10^{-265}:\\
\;\;\;\;\left(\left(j \cdot \left(y5 \cdot y0\right)\right) \cdot y3 - \left(\left(\left(j \cdot y4\right) \cdot y3\right) \cdot y1 + \left(\left(y5 \cdot y0\right) \cdot y2\right) \cdot k\right)\right) + \left(\left(\left(b \cdot y4 - y5 \cdot i\right) \cdot \left(j \cdot t - y \cdot k\right) + \left(\left(\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - c \cdot i\right) - \left(x \cdot j - k \cdot z\right) \cdot \left(b \cdot y0 - i \cdot y1\right)\right) + \left(y0 \cdot c - a \cdot y1\right) \cdot \left(x \cdot y2 - y3 \cdot z\right)\right)\right) - \left(c \cdot y4 - a \cdot y5\right) \cdot \left(t \cdot y2 - y3 \cdot y\right)\right)\\
\mathbf{elif}\;i \le 7.006953314056233 \cdot 10^{-244}:\\
\;\;\;\;\left(\left(\left(\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - c \cdot i\right) - \left(x \cdot j - k \cdot z\right) \cdot \left(b \cdot y0 - i \cdot y1\right)\right) + \left(y0 \cdot c - a \cdot y1\right) \cdot \left(x \cdot y2 - y3 \cdot z\right)\right) - \left(c \cdot y4 - a \cdot y5\right) \cdot \left(t \cdot y2 - y3 \cdot y\right)\right) + \left(y2 \cdot k - y3 \cdot j\right) \cdot \left(y1 \cdot y4 - y5 \cdot y0\right)\\
\mathbf{else}:\\
\;\;\;\;\left(y2 \cdot k - y3 \cdot j\right) \cdot \left(y1 \cdot y4 - y5 \cdot y0\right) + \left(\left(\left(b \cdot y4 - y5 \cdot i\right) \cdot \left(j \cdot t - y \cdot k\right) + \left(\left(\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - c \cdot i\right) - \left(x \cdot j - k \cdot z\right) \cdot \left(b \cdot y0 - i \cdot y1\right)\right) + \left(y0 \cdot c - a \cdot y1\right) \cdot \left(x \cdot y2 - y3 \cdot z\right)\right)\right) - \sqrt[3]{c \cdot y4 - a \cdot y5} \cdot \left(\left(\sqrt[3]{c \cdot y4 - a \cdot y5} \cdot \sqrt[3]{c \cdot y4 - a \cdot y5}\right) \cdot \left(t \cdot y2 - y3 \cdot y\right)\right)\right)\\
\end{array}double f(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 r6456116 = x;
double r6456117 = y;
double r6456118 = r6456116 * r6456117;
double r6456119 = z;
double r6456120 = t;
double r6456121 = r6456119 * r6456120;
double r6456122 = r6456118 - r6456121;
double r6456123 = a;
double r6456124 = b;
double r6456125 = r6456123 * r6456124;
double r6456126 = c;
double r6456127 = i;
double r6456128 = r6456126 * r6456127;
double r6456129 = r6456125 - r6456128;
double r6456130 = r6456122 * r6456129;
double r6456131 = j;
double r6456132 = r6456116 * r6456131;
double r6456133 = k;
double r6456134 = r6456119 * r6456133;
double r6456135 = r6456132 - r6456134;
double r6456136 = y0;
double r6456137 = r6456136 * r6456124;
double r6456138 = y1;
double r6456139 = r6456138 * r6456127;
double r6456140 = r6456137 - r6456139;
double r6456141 = r6456135 * r6456140;
double r6456142 = r6456130 - r6456141;
double r6456143 = y2;
double r6456144 = r6456116 * r6456143;
double r6456145 = y3;
double r6456146 = r6456119 * r6456145;
double r6456147 = r6456144 - r6456146;
double r6456148 = r6456136 * r6456126;
double r6456149 = r6456138 * r6456123;
double r6456150 = r6456148 - r6456149;
double r6456151 = r6456147 * r6456150;
double r6456152 = r6456142 + r6456151;
double r6456153 = r6456120 * r6456131;
double r6456154 = r6456117 * r6456133;
double r6456155 = r6456153 - r6456154;
double r6456156 = y4;
double r6456157 = r6456156 * r6456124;
double r6456158 = y5;
double r6456159 = r6456158 * r6456127;
double r6456160 = r6456157 - r6456159;
double r6456161 = r6456155 * r6456160;
double r6456162 = r6456152 + r6456161;
double r6456163 = r6456120 * r6456143;
double r6456164 = r6456117 * r6456145;
double r6456165 = r6456163 - r6456164;
double r6456166 = r6456156 * r6456126;
double r6456167 = r6456158 * r6456123;
double r6456168 = r6456166 - r6456167;
double r6456169 = r6456165 * r6456168;
double r6456170 = r6456162 - r6456169;
double r6456171 = r6456133 * r6456143;
double r6456172 = r6456131 * r6456145;
double r6456173 = r6456171 - r6456172;
double r6456174 = r6456156 * r6456138;
double r6456175 = r6456158 * r6456136;
double r6456176 = r6456174 - r6456175;
double r6456177 = r6456173 * r6456176;
double r6456178 = r6456170 + r6456177;
return r6456178;
}
double f(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 r6456179 = i;
double r6456180 = -1.885291032353528e+143;
bool r6456181 = r6456179 <= r6456180;
double r6456182 = y0;
double r6456183 = c;
double r6456184 = r6456182 * r6456183;
double r6456185 = a;
double r6456186 = y1;
double r6456187 = r6456185 * r6456186;
double r6456188 = r6456184 - r6456187;
double r6456189 = x;
double r6456190 = y2;
double r6456191 = r6456189 * r6456190;
double r6456192 = y3;
double r6456193 = z;
double r6456194 = r6456192 * r6456193;
double r6456195 = r6456191 - r6456194;
double r6456196 = r6456188 * r6456195;
double r6456197 = y;
double r6456198 = r6456189 * r6456197;
double r6456199 = t;
double r6456200 = r6456199 * r6456193;
double r6456201 = r6456198 - r6456200;
double r6456202 = b;
double r6456203 = r6456185 * r6456202;
double r6456204 = r6456183 * r6456179;
double r6456205 = r6456203 - r6456204;
double r6456206 = r6456201 * r6456205;
double r6456207 = k;
double r6456208 = r6456207 * r6456186;
double r6456209 = r6456193 * r6456208;
double r6456210 = r6456179 * r6456209;
double r6456211 = r6456189 * r6456186;
double r6456212 = j;
double r6456213 = r6456211 * r6456212;
double r6456214 = r6456179 * r6456213;
double r6456215 = r6456202 * r6456182;
double r6456216 = r6456215 * r6456193;
double r6456217 = r6456216 * r6456207;
double r6456218 = r6456214 + r6456217;
double r6456219 = r6456210 - r6456218;
double r6456220 = r6456206 - r6456219;
double r6456221 = r6456196 + r6456220;
double r6456222 = y4;
double r6456223 = r6456202 * r6456222;
double r6456224 = y5;
double r6456225 = r6456224 * r6456179;
double r6456226 = r6456223 - r6456225;
double r6456227 = r6456212 * r6456199;
double r6456228 = r6456197 * r6456207;
double r6456229 = r6456227 - r6456228;
double r6456230 = r6456226 * r6456229;
double r6456231 = r6456221 + r6456230;
double r6456232 = r6456183 * r6456222;
double r6456233 = r6456185 * r6456224;
double r6456234 = r6456232 - r6456233;
double r6456235 = r6456199 * r6456190;
double r6456236 = r6456192 * r6456197;
double r6456237 = r6456235 - r6456236;
double r6456238 = r6456234 * r6456237;
double r6456239 = r6456231 - r6456238;
double r6456240 = r6456190 * r6456207;
double r6456241 = r6456192 * r6456212;
double r6456242 = r6456240 - r6456241;
double r6456243 = r6456186 * r6456222;
double r6456244 = r6456224 * r6456182;
double r6456245 = r6456243 - r6456244;
double r6456246 = r6456242 * r6456245;
double r6456247 = r6456239 + r6456246;
double r6456248 = -7.571795248650015e-30;
bool r6456249 = r6456179 <= r6456248;
double r6456250 = r6456189 * r6456212;
double r6456251 = r6456207 * r6456193;
double r6456252 = r6456250 - r6456251;
double r6456253 = r6456179 * r6456186;
double r6456254 = r6456215 - r6456253;
double r6456255 = r6456252 * r6456254;
double r6456256 = r6456206 - r6456255;
double r6456257 = r6456256 + r6456196;
double r6456258 = r6456230 + r6456257;
double r6456259 = r6456197 * r6456224;
double r6456260 = r6456192 * r6456259;
double r6456261 = r6456185 * r6456260;
double r6456262 = r6456224 * r6456199;
double r6456263 = r6456190 * r6456262;
double r6456264 = r6456185 * r6456263;
double r6456265 = r6456197 * r6456183;
double r6456266 = r6456222 * r6456265;
double r6456267 = r6456266 * r6456192;
double r6456268 = r6456264 + r6456267;
double r6456269 = r6456261 - r6456268;
double r6456270 = r6456258 - r6456269;
double r6456271 = r6456246 + r6456270;
double r6456272 = -4.0262567889326796e-265;
bool r6456273 = r6456179 <= r6456272;
double r6456274 = r6456212 * r6456244;
double r6456275 = r6456274 * r6456192;
double r6456276 = r6456212 * r6456222;
double r6456277 = r6456276 * r6456192;
double r6456278 = r6456277 * r6456186;
double r6456279 = r6456244 * r6456190;
double r6456280 = r6456279 * r6456207;
double r6456281 = r6456278 + r6456280;
double r6456282 = r6456275 - r6456281;
double r6456283 = r6456258 - r6456238;
double r6456284 = r6456282 + r6456283;
double r6456285 = 7.006953314056233e-244;
bool r6456286 = r6456179 <= r6456285;
double r6456287 = r6456257 - r6456238;
double r6456288 = r6456287 + r6456246;
double r6456289 = cbrt(r6456234);
double r6456290 = r6456289 * r6456289;
double r6456291 = r6456290 * r6456237;
double r6456292 = r6456289 * r6456291;
double r6456293 = r6456258 - r6456292;
double r6456294 = r6456246 + r6456293;
double r6456295 = r6456286 ? r6456288 : r6456294;
double r6456296 = r6456273 ? r6456284 : r6456295;
double r6456297 = r6456249 ? r6456271 : r6456296;
double r6456298 = r6456181 ? r6456247 : r6456297;
return r6456298;
}



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
Results
if i < -1.885291032353528e+143Initial program 29.8
Taylor expanded around inf 27.9
if -1.885291032353528e+143 < i < -7.571795248650015e-30Initial program 26.0
Taylor expanded around inf 27.1
if -7.571795248650015e-30 < i < -4.0262567889326796e-265Initial program 24.5
Taylor expanded around inf 26.4
if -4.0262567889326796e-265 < i < 7.006953314056233e-244Initial program 26.5
Taylor expanded around 0 29.4
if 7.006953314056233e-244 < i Initial program 25.8
rmApplied add-cube-cbrt25.9
Applied associate-*r*25.9
Final simplification26.7
herbie shell --seed 2019168
(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"
(+ (- (+ (+ (- (* (- (* 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)))))