#include #include #include #include #include char *name = "Linear.Matrix:det44 from linear-1.19.1.3"; double f_if(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k, float y0, float y1, float y2, float y3, float y4, float y5) { float r26174 = x; float r26175 = y; float r26176 = r26174 * r26175; float r26177 = z; float r26178 = t; float r26179 = r26177 * r26178; float r26180 = r26176 - r26179; float r26181 = a; float r26182 = b; float r26183 = r26181 * r26182; float r26184 = c; float r26185 = i; float r26186 = r26184 * r26185; float r26187 = r26183 - r26186; float r26188 = r26180 * r26187; float r26189 = j; float r26190 = r26174 * r26189; float r26191 = k; float r26192 = r26177 * r26191; float r26193 = r26190 - r26192; float r26194 = y0; float r26195 = r26194 * r26182; float r26196 = y1; float r26197 = r26196 * r26185; float r26198 = r26195 - r26197; float r26199 = r26193 * r26198; float r26200 = r26188 - r26199; float r26201 = y2; float r26202 = r26174 * r26201; float r26203 = y3; float r26204 = r26177 * r26203; float r26205 = r26202 - r26204; float r26206 = r26194 * r26184; float r26207 = r26196 * r26181; float r26208 = r26206 - r26207; float r26209 = r26205 * r26208; float r26210 = r26200 + r26209; float r26211 = r26178 * r26189; float r26212 = r26175 * r26191; float r26213 = r26211 - r26212; float r26214 = y4; float r26215 = r26214 * r26182; float r26216 = y5; float r26217 = r26216 * r26185; float r26218 = r26215 - r26217; float r26219 = r26213 * r26218; float r26220 = r26210 + r26219; float r26221 = r26178 * r26201; float r26222 = r26175 * r26203; float r26223 = r26221 - r26222; float r26224 = r26214 * r26184; float r26225 = r26216 * r26181; float r26226 = r26224 - r26225; float r26227 = r26223 * r26226; float r26228 = r26220 - r26227; float r26229 = r26191 * r26201; float r26230 = r26189 * r26203; float r26231 = r26229 - r26230; float r26232 = r26214 * r26196; float r26233 = r26216 * r26194; float r26234 = r26232 - r26233; float r26235 = r26231 * r26234; float r26236 = r26228 + r26235; return r26236; } double f_id(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 r26237 = x; double r26238 = y; double r26239 = r26237 * r26238; double r26240 = z; double r26241 = t; double r26242 = r26240 * r26241; double r26243 = r26239 - r26242; double r26244 = a; double r26245 = b; double r26246 = r26244 * r26245; double r26247 = c; double r26248 = i; double r26249 = r26247 * r26248; double r26250 = r26246 - r26249; double r26251 = r26243 * r26250; double r26252 = j; double r26253 = r26237 * r26252; double r26254 = k; double r26255 = r26240 * r26254; double r26256 = r26253 - r26255; double r26257 = y0; double r26258 = r26257 * r26245; double r26259 = y1; double r26260 = r26259 * r26248; double r26261 = r26258 - r26260; double r26262 = r26256 * r26261; double r26263 = r26251 - r26262; double r26264 = y2; double r26265 = r26237 * r26264; double r26266 = y3; double r26267 = r26240 * r26266; double r26268 = r26265 - r26267; double r26269 = r26257 * r26247; double r26270 = r26259 * r26244; double r26271 = r26269 - r26270; double r26272 = r26268 * r26271; double r26273 = r26263 + r26272; double r26274 = r26241 * r26252; double r26275 = r26238 * r26254; double r26276 = r26274 - r26275; double r26277 = y4; double r26278 = r26277 * r26245; double r26279 = y5; double r26280 = r26279 * r26248; double r26281 = r26278 - r26280; double r26282 = r26276 * r26281; double r26283 = r26273 + r26282; double r26284 = r26241 * r26264; double r26285 = r26238 * r26266; double r26286 = r26284 - r26285; double r26287 = r26277 * r26247; double r26288 = r26279 * r26244; double r26289 = r26287 - r26288; double r26290 = r26286 * r26289; double r26291 = r26283 - r26290; double r26292 = r26254 * r26264; double r26293 = r26252 * r26266; double r26294 = r26292 - r26293; double r26295 = r26277 * r26259; double r26296 = r26279 * r26257; double r26297 = r26295 - r26296; double r26298 = r26294 * r26297; double r26299 = r26291 + r26298; return r26299; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k, float y0, float y1, float y2, float y3, float y4, float y5) { float r26300 = b; float r26301 = a; float r26302 = r26300 * r26301; float r26303 = i; float r26304 = c; float r26305 = r26303 * r26304; float r26306 = r26302 - r26305; float r26307 = y; float r26308 = x; float r26309 = r26307 * r26308; float r26310 = z; float r26311 = t; float r26312 = r26310 * r26311; float r26313 = r26309 - r26312; float r26314 = r26306 * r26313; float r26315 = j; float r26316 = r26315 * r26308; float r26317 = k; float r26318 = r26310 * r26317; float r26319 = r26316 - r26318; float r26320 = y0; float r26321 = r26300 * r26320; float r26322 = y1; float r26323 = r26303 * r26322; float r26324 = r26321 - r26323; float r26325 = r26319 * r26324; float r26326 = r26314 - r26325; float r26327 = r26315 * r26311; float r26328 = r26317 * r26307; float r26329 = r26327 - r26328; float r26330 = y4; float r26331 = r26330 * r26300; float r26332 = y5; float r26333 = r26332 * r26303; float r26334 = r26331 - r26333; float r26335 = r26329 * r26334; float r26336 = y2; float r26337 = r26336 * r26308; float r26338 = y3; float r26339 = r26310 * r26338; float r26340 = r26337 - r26339; float r26341 = r26320 * r26304; float r26342 = r26322 * r26301; float r26343 = r26341 - r26342; float r26344 = r26340 * r26343; float r26345 = r26335 + r26344; float r26346 = r26326 + r26345; float r26347 = r26336 * r26311; float r26348 = r26307 * r26338; float r26349 = r26347 - r26348; float r26350 = cbrt(r26349); float r26351 = r26330 * r26304; float r26352 = r26332 * r26301; float r26353 = r26351 - r26352; float r26354 = r26350 * r26353; float r26355 = r26350 * r26350; float r26356 = r26354 * r26355; float r26357 = r26346 - r26356; float r26358 = -3.756239411623325e+307; bool r26359 = r26357 <= r26358; float r26360 = r26308 * r26307; float r26361 = r26360 - r26312; float r26362 = r26301 * r26300; float r26363 = r26304 * r26303; float r26364 = r26362 - r26363; float r26365 = r26361 * r26364; float r26366 = r26308 * r26315; float r26367 = r26366 - r26318; float r26368 = r26320 * r26300; float r26369 = r26322 * r26303; float r26370 = r26368 - r26369; float r26371 = r26367 * r26370; float r26372 = r26365 - r26371; float r26373 = r26308 * r26336; float r26374 = r26373 - r26339; float r26375 = r26374 * r26343; float r26376 = r26372 + r26375; float r26377 = r26300 * r26317; float r26378 = r26307 * r26377; float r26379 = r26330 * r26378; float r26380 = r26317 * r26303; float r26381 = r26332 * r26380; float r26382 = r26307 * r26381; float r26383 = r26379 + r26382; float r26384 = r26311 * r26303; float r26385 = r26332 * r26384; float r26386 = r26315 * r26385; float r26387 = r26383 - r26386; float r26388 = r26376 + r26387; float r26389 = r26311 * r26336; float r26390 = r26389 - r26348; float r26391 = r26390 * r26353; float r26392 = r26388 - r26391; float r26393 = r26317 * r26336; float r26394 = r26315 * r26338; float r26395 = r26393 - r26394; float r26396 = r26330 * r26322; float r26397 = r26332 * r26320; float r26398 = r26396 - r26397; float r26399 = r26395 * r26398; float r26400 = r26392 + r26399; float r26401 = 3.445107154068823e+266; bool r26402 = r26357 <= r26401; float r26403 = r26311 * r26315; float r26404 = r26307 * r26317; float r26405 = r26403 - r26404; float r26406 = r26405 * r26334; float r26407 = r26376 + r26406; float r26408 = cbrt(r26391); float r26409 = r26408 * r26408; float r26410 = r26409 * r26408; float r26411 = r26407 - r26410; float r26412 = r26411 + r26399; float r26413 = r26322 * r26330; float r26414 = r26320 * r26332; float r26415 = r26413 - r26414; float r26416 = r26338 * r26315; float r26417 = r26393 - r26416; float r26418 = r26415 * r26417; float r26419 = r26418 + r26335; float r26420 = r26338 * r26310; float r26421 = r26337 - r26420; float r26422 = r26343 * r26421; float r26423 = r26338 * r26307; float r26424 = r26389 - r26423; float r26425 = r26304 * r26330; float r26426 = r26301 * r26332; float r26427 = r26425 - r26426; float r26428 = r26424 * r26427; float r26429 = r26314 - r26428; float r26430 = r26422 + r26429; float r26431 = r26419 + r26430; float r26432 = r26402 ? r26412 : r26431; float r26433 = r26359 ? r26400 : r26432; return r26433; } double f_od(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 r26434 = b; double r26435 = a; double r26436 = r26434 * r26435; double r26437 = i; double r26438 = c; double r26439 = r26437 * r26438; double r26440 = r26436 - r26439; double r26441 = y; double r26442 = x; double r26443 = r26441 * r26442; double r26444 = z; double r26445 = t; double r26446 = r26444 * r26445; double r26447 = r26443 - r26446; double r26448 = r26440 * r26447; double r26449 = j; double r26450 = r26449 * r26442; double r26451 = k; double r26452 = r26444 * r26451; double r26453 = r26450 - r26452; double r26454 = y0; double r26455 = r26434 * r26454; double r26456 = y1; double r26457 = r26437 * r26456; double r26458 = r26455 - r26457; double r26459 = r26453 * r26458; double r26460 = r26448 - r26459; double r26461 = r26449 * r26445; double r26462 = r26451 * r26441; double r26463 = r26461 - r26462; double r26464 = y4; double r26465 = r26464 * r26434; double r26466 = y5; double r26467 = r26466 * r26437; double r26468 = r26465 - r26467; double r26469 = r26463 * r26468; double r26470 = y2; double r26471 = r26470 * r26442; double r26472 = y3; double r26473 = r26444 * r26472; double r26474 = r26471 - r26473; double r26475 = r26454 * r26438; double r26476 = r26456 * r26435; double r26477 = r26475 - r26476; double r26478 = r26474 * r26477; double r26479 = r26469 + r26478; double r26480 = r26460 + r26479; double r26481 = r26470 * r26445; double r26482 = r26441 * r26472; double r26483 = r26481 - r26482; double r26484 = cbrt(r26483); double r26485 = r26464 * r26438; double r26486 = r26466 * r26435; double r26487 = r26485 - r26486; double r26488 = r26484 * r26487; double r26489 = r26484 * r26484; double r26490 = r26488 * r26489; double r26491 = r26480 - r26490; double r26492 = -3.756239411623325e+307; bool r26493 = r26491 <= r26492; double r26494 = r26442 * r26441; double r26495 = r26494 - r26446; double r26496 = r26435 * r26434; double r26497 = r26438 * r26437; double r26498 = r26496 - r26497; double r26499 = r26495 * r26498; double r26500 = r26442 * r26449; double r26501 = r26500 - r26452; double r26502 = r26454 * r26434; double r26503 = r26456 * r26437; double r26504 = r26502 - r26503; double r26505 = r26501 * r26504; double r26506 = r26499 - r26505; double r26507 = r26442 * r26470; double r26508 = r26507 - r26473; double r26509 = r26508 * r26477; double r26510 = r26506 + r26509; double r26511 = r26434 * r26451; double r26512 = r26441 * r26511; double r26513 = r26464 * r26512; double r26514 = r26451 * r26437; double r26515 = r26466 * r26514; double r26516 = r26441 * r26515; double r26517 = r26513 + r26516; double r26518 = r26445 * r26437; double r26519 = r26466 * r26518; double r26520 = r26449 * r26519; double r26521 = r26517 - r26520; double r26522 = r26510 + r26521; double r26523 = r26445 * r26470; double r26524 = r26523 - r26482; double r26525 = r26524 * r26487; double r26526 = r26522 - r26525; double r26527 = r26451 * r26470; double r26528 = r26449 * r26472; double r26529 = r26527 - r26528; double r26530 = r26464 * r26456; double r26531 = r26466 * r26454; double r26532 = r26530 - r26531; double r26533 = r26529 * r26532; double r26534 = r26526 + r26533; double r26535 = 3.445107154068823e+266; bool r26536 = r26491 <= r26535; double r26537 = r26445 * r26449; double r26538 = r26441 * r26451; double r26539 = r26537 - r26538; double r26540 = r26539 * r26468; double r26541 = r26510 + r26540; double r26542 = cbrt(r26525); double r26543 = r26542 * r26542; double r26544 = r26543 * r26542; double r26545 = r26541 - r26544; double r26546 = r26545 + r26533; double r26547 = r26456 * r26464; double r26548 = r26454 * r26466; double r26549 = r26547 - r26548; double r26550 = r26472 * r26449; double r26551 = r26527 - r26550; double r26552 = r26549 * r26551; double r26553 = r26552 + r26469; double r26554 = r26472 * r26444; double r26555 = r26471 - r26554; double r26556 = r26477 * r26555; double r26557 = r26472 * r26441; double r26558 = r26523 - r26557; double r26559 = r26438 * r26464; double r26560 = r26435 * r26466; double r26561 = r26559 - r26560; double r26562 = r26558 * r26561; double r26563 = r26448 - r26562; double r26564 = r26556 + r26563; double r26565 = r26553 + r26564; double r26566 = r26536 ? r26546 : r26565; double r26567 = r26493 ? r26534 : r26566; return r26567; } void mpfr_fmod2(mpfr_t r, mpfr_t n, mpfr_t d, mpfr_rnd_t rmd) { mpfr_fmod(r, n, d, rmd); if (mpfr_cmp_ui(r, 0) < 0) mpfr_add(r, r, d, rmd); } static mpfr_t r26568, r26569, r26570, r26571, r26572, r26573, r26574, r26575, r26576, r26577, r26578, r26579, r26580, r26581, r26582, r26583, r26584, r26585, r26586, r26587, r26588, r26589, r26590, r26591, r26592, r26593, r26594, r26595, r26596, r26597, r26598, r26599, r26600, r26601, r26602, r26603, r26604, r26605, r26606, r26607, r26608, r26609, r26610, r26611, r26612, r26613, r26614, r26615, r26616, r26617, r26618, r26619, r26620, r26621, r26622, r26623, r26624, r26625, r26626, r26627, r26628, r26629, r26630; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26568); mpfr_init(r26569); mpfr_init(r26570); mpfr_init(r26571); mpfr_init(r26572); mpfr_init(r26573); mpfr_init(r26574); mpfr_init(r26575); mpfr_init(r26576); mpfr_init(r26577); mpfr_init(r26578); mpfr_init(r26579); mpfr_init(r26580); mpfr_init(r26581); mpfr_init(r26582); mpfr_init(r26583); mpfr_init(r26584); mpfr_init(r26585); mpfr_init(r26586); mpfr_init(r26587); mpfr_init(r26588); mpfr_init(r26589); mpfr_init(r26590); mpfr_init(r26591); mpfr_init(r26592); mpfr_init(r26593); mpfr_init(r26594); mpfr_init(r26595); mpfr_init(r26596); mpfr_init(r26597); mpfr_init(r26598); mpfr_init(r26599); mpfr_init(r26600); mpfr_init(r26601); mpfr_init(r26602); mpfr_init(r26603); mpfr_init(r26604); mpfr_init(r26605); mpfr_init(r26606); mpfr_init(r26607); mpfr_init(r26608); mpfr_init(r26609); mpfr_init(r26610); mpfr_init(r26611); mpfr_init(r26612); mpfr_init(r26613); mpfr_init(r26614); mpfr_init(r26615); mpfr_init(r26616); mpfr_init(r26617); mpfr_init(r26618); mpfr_init(r26619); mpfr_init(r26620); mpfr_init(r26621); mpfr_init(r26622); mpfr_init(r26623); mpfr_init(r26624); mpfr_init(r26625); mpfr_init(r26626); mpfr_init(r26627); mpfr_init(r26628); mpfr_init(r26629); mpfr_init(r26630); } double f_im(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) { mpfr_set_d(r26568, x, MPFR_RNDN); mpfr_set_d(r26569, y, MPFR_RNDN); mpfr_mul(r26570, r26568, r26569, MPFR_RNDN); mpfr_set_d(r26571, z, MPFR_RNDN); mpfr_set_d(r26572, t, MPFR_RNDN); mpfr_mul(r26573, r26571, r26572, MPFR_RNDN); mpfr_sub(r26574, r26570, r26573, MPFR_RNDN); mpfr_set_d(r26575, a, MPFR_RNDN); mpfr_set_d(r26576, b, MPFR_RNDN); mpfr_mul(r26577, r26575, r26576, MPFR_RNDN); mpfr_set_d(r26578, c, MPFR_RNDN); mpfr_set_d(r26579, i, MPFR_RNDN); mpfr_mul(r26580, r26578, r26579, MPFR_RNDN); mpfr_sub(r26581, r26577, r26580, MPFR_RNDN); mpfr_mul(r26582, r26574, r26581, MPFR_RNDN); mpfr_set_d(r26583, j, MPFR_RNDN); mpfr_mul(r26584, r26568, r26583, MPFR_RNDN); mpfr_set_d(r26585, k, MPFR_RNDN); mpfr_mul(r26586, r26571, r26585, MPFR_RNDN); mpfr_sub(r26587, r26584, r26586, MPFR_RNDN); mpfr_set_d(r26588, y0, MPFR_RNDN); mpfr_mul(r26589, r26588, r26576, MPFR_RNDN); mpfr_set_d(r26590, y1, MPFR_RNDN); mpfr_mul(r26591, r26590, r26579, MPFR_RNDN); mpfr_sub(r26592, r26589, r26591, MPFR_RNDN); mpfr_mul(r26593, r26587, r26592, MPFR_RNDN); mpfr_sub(r26594, r26582, r26593, MPFR_RNDN); mpfr_set_d(r26595, y2, MPFR_RNDN); mpfr_mul(r26596, r26568, r26595, MPFR_RNDN); mpfr_set_d(r26597, y3, MPFR_RNDN); mpfr_mul(r26598, r26571, r26597, MPFR_RNDN); mpfr_sub(r26599, r26596, r26598, MPFR_RNDN); mpfr_mul(r26600, r26588, r26578, MPFR_RNDN); mpfr_mul(r26601, r26590, r26575, MPFR_RNDN); mpfr_sub(r26602, r26600, r26601, MPFR_RNDN); mpfr_mul(r26603, r26599, r26602, MPFR_RNDN); mpfr_add(r26604, r26594, r26603, MPFR_RNDN); mpfr_mul(r26605, r26572, r26583, MPFR_RNDN); mpfr_mul(r26606, r26569, r26585, MPFR_RNDN); mpfr_sub(r26607, r26605, r26606, MPFR_RNDN); mpfr_set_d(r26608, y4, MPFR_RNDN); mpfr_mul(r26609, r26608, r26576, MPFR_RNDN); mpfr_set_d(r26610, y5, MPFR_RNDN); mpfr_mul(r26611, r26610, r26579, MPFR_RNDN); mpfr_sub(r26612, r26609, r26611, MPFR_RNDN); mpfr_mul(r26613, r26607, r26612, MPFR_RNDN); mpfr_add(r26614, r26604, r26613, MPFR_RNDN); mpfr_mul(r26615, r26572, r26595, MPFR_RNDN); mpfr_mul(r26616, r26569, r26597, MPFR_RNDN); mpfr_sub(r26617, r26615, r26616, MPFR_RNDN); mpfr_mul(r26618, r26608, r26578, MPFR_RNDN); mpfr_mul(r26619, r26610, r26575, MPFR_RNDN); mpfr_sub(r26620, r26618, r26619, MPFR_RNDN); mpfr_mul(r26621, r26617, r26620, MPFR_RNDN); mpfr_sub(r26622, r26614, r26621, MPFR_RNDN); mpfr_mul(r26623, r26585, r26595, MPFR_RNDN); mpfr_mul(r26624, r26583, r26597, MPFR_RNDN); mpfr_sub(r26625, r26623, r26624, MPFR_RNDN); mpfr_mul(r26626, r26608, r26590, MPFR_RNDN); mpfr_mul(r26627, r26610, r26588, MPFR_RNDN); mpfr_sub(r26628, r26626, r26627, MPFR_RNDN); mpfr_mul(r26629, r26625, r26628, MPFR_RNDN); mpfr_add(r26630, r26622, r26629, MPFR_RNDN); return mpfr_get_d(r26630, MPFR_RNDN); } static mpfr_t r26631, r26632, r26633, r26634, r26635, r26636, r26637, r26638, r26639, r26640, r26641, r26642, r26643, r26644, r26645, r26646, r26647, r26648, r26649, r26650, r26651, r26652, r26653, r26654, r26655, r26656, r26657, r26658, r26659, r26660, r26661, r26662, r26663, r26664, r26665, r26666, r26667, r26668, r26669, r26670, r26671, r26672, r26673, r26674, r26675, r26676, r26677, r26678, r26679, r26680, r26681, r26682, r26683, r26684, r26685, r26686, r26687, r26688, r26689, r26690, r26691, r26692, r26693, r26694, r26695, r26696, r26697, r26698, r26699, r26700, r26701, r26702, r26703, r26704, r26705, r26706, r26707, r26708, r26709, r26710, r26711, r26712, r26713, r26714, r26715, r26716, r26717, r26718, r26719, r26720, r26721, r26722, r26723, r26724, r26725, r26726, r26727, r26728, r26729, r26730, r26731, r26732, r26733, r26734, r26735, r26736, r26737, r26738, r26739, r26740, r26741, r26742, r26743, r26744, r26745, r26746, r26747, r26748, r26749, r26750, r26751, r26752, r26753, r26754, r26755, r26756, r26757, r26758, r26759, r26760, r26761, r26762, r26763, r26764; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26631); mpfr_init(r26632); mpfr_init(r26633); mpfr_init(r26634); mpfr_init(r26635); mpfr_init(r26636); mpfr_init(r26637); mpfr_init(r26638); mpfr_init(r26639); mpfr_init(r26640); mpfr_init(r26641); mpfr_init(r26642); mpfr_init(r26643); mpfr_init(r26644); mpfr_init(r26645); mpfr_init(r26646); mpfr_init(r26647); mpfr_init(r26648); mpfr_init(r26649); mpfr_init(r26650); mpfr_init(r26651); mpfr_init(r26652); mpfr_init(r26653); mpfr_init(r26654); mpfr_init(r26655); mpfr_init(r26656); mpfr_init(r26657); mpfr_init(r26658); mpfr_init(r26659); mpfr_init(r26660); mpfr_init(r26661); mpfr_init(r26662); mpfr_init(r26663); mpfr_init(r26664); mpfr_init(r26665); mpfr_init(r26666); mpfr_init(r26667); mpfr_init(r26668); mpfr_init(r26669); mpfr_init(r26670); mpfr_init(r26671); mpfr_init(r26672); mpfr_init(r26673); mpfr_init(r26674); mpfr_init(r26675); mpfr_init(r26676); mpfr_init(r26677); mpfr_init(r26678); mpfr_init(r26679); mpfr_init(r26680); mpfr_init(r26681); mpfr_init(r26682); mpfr_init(r26683); mpfr_init(r26684); mpfr_init(r26685); mpfr_init(r26686); mpfr_init(r26687); mpfr_init(r26688); mpfr_init_set_str(r26689, "-3.756239411623325e+307", 10, MPFR_RNDN); mpfr_init(r26690); mpfr_init(r26691); mpfr_init(r26692); mpfr_init(r26693); mpfr_init(r26694); mpfr_init(r26695); mpfr_init(r26696); mpfr_init(r26697); mpfr_init(r26698); mpfr_init(r26699); mpfr_init(r26700); mpfr_init(r26701); mpfr_init(r26702); mpfr_init(r26703); mpfr_init(r26704); mpfr_init(r26705); mpfr_init(r26706); mpfr_init(r26707); mpfr_init(r26708); mpfr_init(r26709); mpfr_init(r26710); mpfr_init(r26711); mpfr_init(r26712); mpfr_init(r26713); mpfr_init(r26714); mpfr_init(r26715); mpfr_init(r26716); mpfr_init(r26717); mpfr_init(r26718); mpfr_init(r26719); mpfr_init(r26720); mpfr_init(r26721); mpfr_init(r26722); mpfr_init(r26723); mpfr_init(r26724); mpfr_init(r26725); mpfr_init(r26726); mpfr_init(r26727); mpfr_init(r26728); mpfr_init(r26729); mpfr_init(r26730); mpfr_init(r26731); mpfr_init_set_str(r26732, "3.445107154068823e+266", 10, MPFR_RNDN); mpfr_init(r26733); mpfr_init(r26734); mpfr_init(r26735); mpfr_init(r26736); mpfr_init(r26737); mpfr_init(r26738); mpfr_init(r26739); mpfr_init(r26740); mpfr_init(r26741); mpfr_init(r26742); mpfr_init(r26743); mpfr_init(r26744); mpfr_init(r26745); mpfr_init(r26746); mpfr_init(r26747); mpfr_init(r26748); mpfr_init(r26749); mpfr_init(r26750); mpfr_init(r26751); mpfr_init(r26752); mpfr_init(r26753); mpfr_init(r26754); mpfr_init(r26755); mpfr_init(r26756); mpfr_init(r26757); mpfr_init(r26758); mpfr_init(r26759); mpfr_init(r26760); mpfr_init(r26761); mpfr_init(r26762); mpfr_init(r26763); mpfr_init(r26764); } double f_fm(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) { mpfr_set_d(r26631, b, MPFR_RNDN); mpfr_set_d(r26632, a, MPFR_RNDN); mpfr_mul(r26633, r26631, r26632, MPFR_RNDN); mpfr_set_d(r26634, i, MPFR_RNDN); mpfr_set_d(r26635, c, MPFR_RNDN); mpfr_mul(r26636, r26634, r26635, MPFR_RNDN); mpfr_sub(r26637, r26633, r26636, MPFR_RNDN); mpfr_set_d(r26638, y, MPFR_RNDN); mpfr_set_d(r26639, x, MPFR_RNDN); mpfr_mul(r26640, r26638, r26639, MPFR_RNDN); mpfr_set_d(r26641, z, MPFR_RNDN); mpfr_set_d(r26642, t, MPFR_RNDN); mpfr_mul(r26643, r26641, r26642, MPFR_RNDN); mpfr_sub(r26644, r26640, r26643, MPFR_RNDN); mpfr_mul(r26645, r26637, r26644, MPFR_RNDN); mpfr_set_d(r26646, j, MPFR_RNDN); mpfr_mul(r26647, r26646, r26639, MPFR_RNDN); mpfr_set_d(r26648, k, MPFR_RNDN); mpfr_mul(r26649, r26641, r26648, MPFR_RNDN); mpfr_sub(r26650, r26647, r26649, MPFR_RNDN); mpfr_set_d(r26651, y0, MPFR_RNDN); mpfr_mul(r26652, r26631, r26651, MPFR_RNDN); mpfr_set_d(r26653, y1, MPFR_RNDN); mpfr_mul(r26654, r26634, r26653, MPFR_RNDN); mpfr_sub(r26655, r26652, r26654, MPFR_RNDN); mpfr_mul(r26656, r26650, r26655, MPFR_RNDN); mpfr_sub(r26657, r26645, r26656, MPFR_RNDN); mpfr_mul(r26658, r26646, r26642, MPFR_RNDN); mpfr_mul(r26659, r26648, r26638, MPFR_RNDN); mpfr_sub(r26660, r26658, r26659, MPFR_RNDN); mpfr_set_d(r26661, y4, MPFR_RNDN); mpfr_mul(r26662, r26661, r26631, MPFR_RNDN); mpfr_set_d(r26663, y5, MPFR_RNDN); mpfr_mul(r26664, r26663, r26634, MPFR_RNDN); mpfr_sub(r26665, r26662, r26664, MPFR_RNDN); mpfr_mul(r26666, r26660, r26665, MPFR_RNDN); mpfr_set_d(r26667, y2, MPFR_RNDN); mpfr_mul(r26668, r26667, r26639, MPFR_RNDN); mpfr_set_d(r26669, y3, MPFR_RNDN); mpfr_mul(r26670, r26641, r26669, MPFR_RNDN); mpfr_sub(r26671, r26668, r26670, MPFR_RNDN); mpfr_mul(r26672, r26651, r26635, MPFR_RNDN); mpfr_mul(r26673, r26653, r26632, MPFR_RNDN); mpfr_sub(r26674, r26672, r26673, MPFR_RNDN); mpfr_mul(r26675, r26671, r26674, MPFR_RNDN); mpfr_add(r26676, r26666, r26675, MPFR_RNDN); mpfr_add(r26677, r26657, r26676, MPFR_RNDN); mpfr_mul(r26678, r26667, r26642, MPFR_RNDN); mpfr_mul(r26679, r26638, r26669, MPFR_RNDN); mpfr_sub(r26680, r26678, r26679, MPFR_RNDN); mpfr_cbrt(r26681, r26680, MPFR_RNDN); mpfr_mul(r26682, r26661, r26635, MPFR_RNDN); mpfr_mul(r26683, r26663, r26632, MPFR_RNDN); mpfr_sub(r26684, r26682, r26683, MPFR_RNDN); mpfr_mul(r26685, r26681, r26684, MPFR_RNDN); mpfr_mul(r26686, r26681, r26681, MPFR_RNDN); mpfr_mul(r26687, r26685, r26686, MPFR_RNDN); mpfr_sub(r26688, r26677, r26687, MPFR_RNDN); ; mpfr_set_si(r26690, mpfr_cmp(r26688, r26689) <= 0, MPFR_RNDN); mpfr_mul(r26691, r26639, r26638, MPFR_RNDN); mpfr_sub(r26692, r26691, r26643, MPFR_RNDN); mpfr_mul(r26693, r26632, r26631, MPFR_RNDN); mpfr_mul(r26694, r26635, r26634, MPFR_RNDN); mpfr_sub(r26695, r26693, r26694, MPFR_RNDN); mpfr_mul(r26696, r26692, r26695, MPFR_RNDN); mpfr_mul(r26697, r26639, r26646, MPFR_RNDN); mpfr_sub(r26698, r26697, r26649, MPFR_RNDN); mpfr_mul(r26699, r26651, r26631, MPFR_RNDN); mpfr_mul(r26700, r26653, r26634, MPFR_RNDN); mpfr_sub(r26701, r26699, r26700, MPFR_RNDN); mpfr_mul(r26702, r26698, r26701, MPFR_RNDN); mpfr_sub(r26703, r26696, r26702, MPFR_RNDN); mpfr_mul(r26704, r26639, r26667, MPFR_RNDN); mpfr_sub(r26705, r26704, r26670, MPFR_RNDN); mpfr_mul(r26706, r26705, r26674, MPFR_RNDN); mpfr_add(r26707, r26703, r26706, MPFR_RNDN); mpfr_mul(r26708, r26631, r26648, MPFR_RNDN); mpfr_mul(r26709, r26638, r26708, MPFR_RNDN); mpfr_mul(r26710, r26661, r26709, MPFR_RNDN); mpfr_mul(r26711, r26648, r26634, MPFR_RNDN); mpfr_mul(r26712, r26663, r26711, MPFR_RNDN); mpfr_mul(r26713, r26638, r26712, MPFR_RNDN); mpfr_add(r26714, r26710, r26713, MPFR_RNDN); mpfr_mul(r26715, r26642, r26634, MPFR_RNDN); mpfr_mul(r26716, r26663, r26715, MPFR_RNDN); mpfr_mul(r26717, r26646, r26716, MPFR_RNDN); mpfr_sub(r26718, r26714, r26717, MPFR_RNDN); mpfr_add(r26719, r26707, r26718, MPFR_RNDN); mpfr_mul(r26720, r26642, r26667, MPFR_RNDN); mpfr_sub(r26721, r26720, r26679, MPFR_RNDN); mpfr_mul(r26722, r26721, r26684, MPFR_RNDN); mpfr_sub(r26723, r26719, r26722, MPFR_RNDN); mpfr_mul(r26724, r26648, r26667, MPFR_RNDN); mpfr_mul(r26725, r26646, r26669, MPFR_RNDN); mpfr_sub(r26726, r26724, r26725, MPFR_RNDN); mpfr_mul(r26727, r26661, r26653, MPFR_RNDN); mpfr_mul(r26728, r26663, r26651, MPFR_RNDN); mpfr_sub(r26729, r26727, r26728, MPFR_RNDN); mpfr_mul(r26730, r26726, r26729, MPFR_RNDN); mpfr_add(r26731, r26723, r26730, MPFR_RNDN); ; mpfr_set_si(r26733, mpfr_cmp(r26688, r26732) <= 0, MPFR_RNDN); mpfr_mul(r26734, r26642, r26646, MPFR_RNDN); mpfr_mul(r26735, r26638, r26648, MPFR_RNDN); mpfr_sub(r26736, r26734, r26735, MPFR_RNDN); mpfr_mul(r26737, r26736, r26665, MPFR_RNDN); mpfr_add(r26738, r26707, r26737, MPFR_RNDN); mpfr_cbrt(r26739, r26722, MPFR_RNDN); mpfr_mul(r26740, r26739, r26739, MPFR_RNDN); mpfr_mul(r26741, r26740, r26739, MPFR_RNDN); mpfr_sub(r26742, r26738, r26741, MPFR_RNDN); mpfr_add(r26743, r26742, r26730, MPFR_RNDN); mpfr_mul(r26744, r26653, r26661, MPFR_RNDN); mpfr_mul(r26745, r26651, r26663, MPFR_RNDN); mpfr_sub(r26746, r26744, r26745, MPFR_RNDN); mpfr_mul(r26747, r26669, r26646, MPFR_RNDN); mpfr_sub(r26748, r26724, r26747, MPFR_RNDN); mpfr_mul(r26749, r26746, r26748, MPFR_RNDN); mpfr_add(r26750, r26749, r26666, MPFR_RNDN); mpfr_mul(r26751, r26669, r26641, MPFR_RNDN); mpfr_sub(r26752, r26668, r26751, MPFR_RNDN); mpfr_mul(r26753, r26674, r26752, MPFR_RNDN); mpfr_mul(r26754, r26669, r26638, MPFR_RNDN); mpfr_sub(r26755, r26720, r26754, MPFR_RNDN); mpfr_mul(r26756, r26635, r26661, MPFR_RNDN); mpfr_mul(r26757, r26632, r26663, MPFR_RNDN); mpfr_sub(r26758, r26756, r26757, MPFR_RNDN); mpfr_mul(r26759, r26755, r26758, MPFR_RNDN); mpfr_sub(r26760, r26645, r26759, MPFR_RNDN); mpfr_add(r26761, r26753, r26760, MPFR_RNDN); mpfr_add(r26762, r26750, r26761, MPFR_RNDN); if (mpfr_get_si(r26733, MPFR_RNDN)) { mpfr_set(r26763, r26743, MPFR_RNDN); } else { mpfr_set(r26763, r26762, MPFR_RNDN); }; if (mpfr_get_si(r26690, MPFR_RNDN)) { mpfr_set(r26764, r26731, MPFR_RNDN); } else { mpfr_set(r26764, r26763, MPFR_RNDN); }; return mpfr_get_d(r26764, MPFR_RNDN); } static mpfr_t r26765, r26766, r26767, r26768, r26769, r26770, r26771, r26772, r26773, r26774, r26775, r26776, r26777, r26778, r26779, r26780, r26781, r26782, r26783, r26784, r26785, r26786, r26787, r26788, r26789, r26790, r26791, r26792, r26793, r26794, r26795, r26796, r26797, r26798, r26799, r26800, r26801, r26802, r26803, r26804, r26805, r26806, r26807, r26808, r26809, r26810, r26811, r26812, r26813, r26814, r26815, r26816, r26817, r26818, r26819, r26820, r26821, r26822, r26823, r26824, r26825, r26826, r26827, r26828, r26829, r26830, r26831, r26832, r26833, r26834, r26835, r26836, r26837, r26838, r26839, r26840, r26841, r26842, r26843, r26844, r26845, r26846, r26847, r26848, r26849, r26850, r26851, r26852, r26853, r26854, r26855, r26856, r26857, r26858, r26859, r26860, r26861, r26862, r26863, r26864, r26865, r26866, r26867, r26868, r26869, r26870, r26871, r26872, r26873, r26874, r26875, r26876, r26877, r26878, r26879, r26880, r26881, r26882, r26883, r26884, r26885, r26886, r26887, r26888, r26889, r26890, r26891, r26892, r26893, r26894, r26895, r26896, r26897, r26898; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26765); mpfr_init(r26766); mpfr_init(r26767); mpfr_init(r26768); mpfr_init(r26769); mpfr_init(r26770); mpfr_init(r26771); mpfr_init(r26772); mpfr_init(r26773); mpfr_init(r26774); mpfr_init(r26775); mpfr_init(r26776); mpfr_init(r26777); mpfr_init(r26778); mpfr_init(r26779); mpfr_init(r26780); mpfr_init(r26781); mpfr_init(r26782); mpfr_init(r26783); mpfr_init(r26784); mpfr_init(r26785); mpfr_init(r26786); mpfr_init(r26787); mpfr_init(r26788); mpfr_init(r26789); mpfr_init(r26790); mpfr_init(r26791); mpfr_init(r26792); mpfr_init(r26793); mpfr_init(r26794); mpfr_init(r26795); mpfr_init(r26796); mpfr_init(r26797); mpfr_init(r26798); mpfr_init(r26799); mpfr_init(r26800); mpfr_init(r26801); mpfr_init(r26802); mpfr_init(r26803); mpfr_init(r26804); mpfr_init(r26805); mpfr_init(r26806); mpfr_init(r26807); mpfr_init(r26808); mpfr_init(r26809); mpfr_init(r26810); mpfr_init(r26811); mpfr_init(r26812); mpfr_init(r26813); mpfr_init(r26814); mpfr_init(r26815); mpfr_init(r26816); mpfr_init(r26817); mpfr_init(r26818); mpfr_init(r26819); mpfr_init(r26820); mpfr_init(r26821); mpfr_init(r26822); mpfr_init_set_str(r26823, "-3.756239411623325e+307", 10, MPFR_RNDN); mpfr_init(r26824); mpfr_init(r26825); mpfr_init(r26826); mpfr_init(r26827); mpfr_init(r26828); mpfr_init(r26829); mpfr_init(r26830); mpfr_init(r26831); mpfr_init(r26832); mpfr_init(r26833); mpfr_init(r26834); mpfr_init(r26835); mpfr_init(r26836); mpfr_init(r26837); mpfr_init(r26838); mpfr_init(r26839); mpfr_init(r26840); mpfr_init(r26841); mpfr_init(r26842); mpfr_init(r26843); mpfr_init(r26844); mpfr_init(r26845); mpfr_init(r26846); mpfr_init(r26847); mpfr_init(r26848); mpfr_init(r26849); mpfr_init(r26850); mpfr_init(r26851); mpfr_init(r26852); mpfr_init(r26853); mpfr_init(r26854); mpfr_init(r26855); mpfr_init(r26856); mpfr_init(r26857); mpfr_init(r26858); mpfr_init(r26859); mpfr_init(r26860); mpfr_init(r26861); mpfr_init(r26862); mpfr_init(r26863); mpfr_init(r26864); mpfr_init(r26865); mpfr_init_set_str(r26866, "3.445107154068823e+266", 10, MPFR_RNDN); mpfr_init(r26867); mpfr_init(r26868); mpfr_init(r26869); mpfr_init(r26870); mpfr_init(r26871); mpfr_init(r26872); mpfr_init(r26873); mpfr_init(r26874); mpfr_init(r26875); mpfr_init(r26876); mpfr_init(r26877); mpfr_init(r26878); mpfr_init(r26879); mpfr_init(r26880); mpfr_init(r26881); mpfr_init(r26882); mpfr_init(r26883); mpfr_init(r26884); mpfr_init(r26885); mpfr_init(r26886); mpfr_init(r26887); mpfr_init(r26888); mpfr_init(r26889); mpfr_init(r26890); mpfr_init(r26891); mpfr_init(r26892); mpfr_init(r26893); mpfr_init(r26894); mpfr_init(r26895); mpfr_init(r26896); mpfr_init(r26897); mpfr_init(r26898); } double f_dm(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) { mpfr_set_d(r26765, b, MPFR_RNDN); mpfr_set_d(r26766, a, MPFR_RNDN); mpfr_mul(r26767, r26765, r26766, MPFR_RNDN); mpfr_set_d(r26768, i, MPFR_RNDN); mpfr_set_d(r26769, c, MPFR_RNDN); mpfr_mul(r26770, r26768, r26769, MPFR_RNDN); mpfr_sub(r26771, r26767, r26770, MPFR_RNDN); mpfr_set_d(r26772, y, MPFR_RNDN); mpfr_set_d(r26773, x, MPFR_RNDN); mpfr_mul(r26774, r26772, r26773, MPFR_RNDN); mpfr_set_d(r26775, z, MPFR_RNDN); mpfr_set_d(r26776, t, MPFR_RNDN); mpfr_mul(r26777, r26775, r26776, MPFR_RNDN); mpfr_sub(r26778, r26774, r26777, MPFR_RNDN); mpfr_mul(r26779, r26771, r26778, MPFR_RNDN); mpfr_set_d(r26780, j, MPFR_RNDN); mpfr_mul(r26781, r26780, r26773, MPFR_RNDN); mpfr_set_d(r26782, k, MPFR_RNDN); mpfr_mul(r26783, r26775, r26782, MPFR_RNDN); mpfr_sub(r26784, r26781, r26783, MPFR_RNDN); mpfr_set_d(r26785, y0, MPFR_RNDN); mpfr_mul(r26786, r26765, r26785, MPFR_RNDN); mpfr_set_d(r26787, y1, MPFR_RNDN); mpfr_mul(r26788, r26768, r26787, MPFR_RNDN); mpfr_sub(r26789, r26786, r26788, MPFR_RNDN); mpfr_mul(r26790, r26784, r26789, MPFR_RNDN); mpfr_sub(r26791, r26779, r26790, MPFR_RNDN); mpfr_mul(r26792, r26780, r26776, MPFR_RNDN); mpfr_mul(r26793, r26782, r26772, MPFR_RNDN); mpfr_sub(r26794, r26792, r26793, MPFR_RNDN); mpfr_set_d(r26795, y4, MPFR_RNDN); mpfr_mul(r26796, r26795, r26765, MPFR_RNDN); mpfr_set_d(r26797, y5, MPFR_RNDN); mpfr_mul(r26798, r26797, r26768, MPFR_RNDN); mpfr_sub(r26799, r26796, r26798, MPFR_RNDN); mpfr_mul(r26800, r26794, r26799, MPFR_RNDN); mpfr_set_d(r26801, y2, MPFR_RNDN); mpfr_mul(r26802, r26801, r26773, MPFR_RNDN); mpfr_set_d(r26803, y3, MPFR_RNDN); mpfr_mul(r26804, r26775, r26803, MPFR_RNDN); mpfr_sub(r26805, r26802, r26804, MPFR_RNDN); mpfr_mul(r26806, r26785, r26769, MPFR_RNDN); mpfr_mul(r26807, r26787, r26766, MPFR_RNDN); mpfr_sub(r26808, r26806, r26807, MPFR_RNDN); mpfr_mul(r26809, r26805, r26808, MPFR_RNDN); mpfr_add(r26810, r26800, r26809, MPFR_RNDN); mpfr_add(r26811, r26791, r26810, MPFR_RNDN); mpfr_mul(r26812, r26801, r26776, MPFR_RNDN); mpfr_mul(r26813, r26772, r26803, MPFR_RNDN); mpfr_sub(r26814, r26812, r26813, MPFR_RNDN); mpfr_cbrt(r26815, r26814, MPFR_RNDN); mpfr_mul(r26816, r26795, r26769, MPFR_RNDN); mpfr_mul(r26817, r26797, r26766, MPFR_RNDN); mpfr_sub(r26818, r26816, r26817, MPFR_RNDN); mpfr_mul(r26819, r26815, r26818, MPFR_RNDN); mpfr_mul(r26820, r26815, r26815, MPFR_RNDN); mpfr_mul(r26821, r26819, r26820, MPFR_RNDN); mpfr_sub(r26822, r26811, r26821, MPFR_RNDN); ; mpfr_set_si(r26824, mpfr_cmp(r26822, r26823) <= 0, MPFR_RNDN); mpfr_mul(r26825, r26773, r26772, MPFR_RNDN); mpfr_sub(r26826, r26825, r26777, MPFR_RNDN); mpfr_mul(r26827, r26766, r26765, MPFR_RNDN); mpfr_mul(r26828, r26769, r26768, MPFR_RNDN); mpfr_sub(r26829, r26827, r26828, MPFR_RNDN); mpfr_mul(r26830, r26826, r26829, MPFR_RNDN); mpfr_mul(r26831, r26773, r26780, MPFR_RNDN); mpfr_sub(r26832, r26831, r26783, MPFR_RNDN); mpfr_mul(r26833, r26785, r26765, MPFR_RNDN); mpfr_mul(r26834, r26787, r26768, MPFR_RNDN); mpfr_sub(r26835, r26833, r26834, MPFR_RNDN); mpfr_mul(r26836, r26832, r26835, MPFR_RNDN); mpfr_sub(r26837, r26830, r26836, MPFR_RNDN); mpfr_mul(r26838, r26773, r26801, MPFR_RNDN); mpfr_sub(r26839, r26838, r26804, MPFR_RNDN); mpfr_mul(r26840, r26839, r26808, MPFR_RNDN); mpfr_add(r26841, r26837, r26840, MPFR_RNDN); mpfr_mul(r26842, r26765, r26782, MPFR_RNDN); mpfr_mul(r26843, r26772, r26842, MPFR_RNDN); mpfr_mul(r26844, r26795, r26843, MPFR_RNDN); mpfr_mul(r26845, r26782, r26768, MPFR_RNDN); mpfr_mul(r26846, r26797, r26845, MPFR_RNDN); mpfr_mul(r26847, r26772, r26846, MPFR_RNDN); mpfr_add(r26848, r26844, r26847, MPFR_RNDN); mpfr_mul(r26849, r26776, r26768, MPFR_RNDN); mpfr_mul(r26850, r26797, r26849, MPFR_RNDN); mpfr_mul(r26851, r26780, r26850, MPFR_RNDN); mpfr_sub(r26852, r26848, r26851, MPFR_RNDN); mpfr_add(r26853, r26841, r26852, MPFR_RNDN); mpfr_mul(r26854, r26776, r26801, MPFR_RNDN); mpfr_sub(r26855, r26854, r26813, MPFR_RNDN); mpfr_mul(r26856, r26855, r26818, MPFR_RNDN); mpfr_sub(r26857, r26853, r26856, MPFR_RNDN); mpfr_mul(r26858, r26782, r26801, MPFR_RNDN); mpfr_mul(r26859, r26780, r26803, MPFR_RNDN); mpfr_sub(r26860, r26858, r26859, MPFR_RNDN); mpfr_mul(r26861, r26795, r26787, MPFR_RNDN); mpfr_mul(r26862, r26797, r26785, MPFR_RNDN); mpfr_sub(r26863, r26861, r26862, MPFR_RNDN); mpfr_mul(r26864, r26860, r26863, MPFR_RNDN); mpfr_add(r26865, r26857, r26864, MPFR_RNDN); ; mpfr_set_si(r26867, mpfr_cmp(r26822, r26866) <= 0, MPFR_RNDN); mpfr_mul(r26868, r26776, r26780, MPFR_RNDN); mpfr_mul(r26869, r26772, r26782, MPFR_RNDN); mpfr_sub(r26870, r26868, r26869, MPFR_RNDN); mpfr_mul(r26871, r26870, r26799, MPFR_RNDN); mpfr_add(r26872, r26841, r26871, MPFR_RNDN); mpfr_cbrt(r26873, r26856, MPFR_RNDN); mpfr_mul(r26874, r26873, r26873, MPFR_RNDN); mpfr_mul(r26875, r26874, r26873, MPFR_RNDN); mpfr_sub(r26876, r26872, r26875, MPFR_RNDN); mpfr_add(r26877, r26876, r26864, MPFR_RNDN); mpfr_mul(r26878, r26787, r26795, MPFR_RNDN); mpfr_mul(r26879, r26785, r26797, MPFR_RNDN); mpfr_sub(r26880, r26878, r26879, MPFR_RNDN); mpfr_mul(r26881, r26803, r26780, MPFR_RNDN); mpfr_sub(r26882, r26858, r26881, MPFR_RNDN); mpfr_mul(r26883, r26880, r26882, MPFR_RNDN); mpfr_add(r26884, r26883, r26800, MPFR_RNDN); mpfr_mul(r26885, r26803, r26775, MPFR_RNDN); mpfr_sub(r26886, r26802, r26885, MPFR_RNDN); mpfr_mul(r26887, r26808, r26886, MPFR_RNDN); mpfr_mul(r26888, r26803, r26772, MPFR_RNDN); mpfr_sub(r26889, r26854, r26888, MPFR_RNDN); mpfr_mul(r26890, r26769, r26795, MPFR_RNDN); mpfr_mul(r26891, r26766, r26797, MPFR_RNDN); mpfr_sub(r26892, r26890, r26891, MPFR_RNDN); mpfr_mul(r26893, r26889, r26892, MPFR_RNDN); mpfr_sub(r26894, r26779, r26893, MPFR_RNDN); mpfr_add(r26895, r26887, r26894, MPFR_RNDN); mpfr_add(r26896, r26884, r26895, MPFR_RNDN); if (mpfr_get_si(r26867, MPFR_RNDN)) { mpfr_set(r26897, r26877, MPFR_RNDN); } else { mpfr_set(r26897, r26896, MPFR_RNDN); }; if (mpfr_get_si(r26824, MPFR_RNDN)) { mpfr_set(r26898, r26865, MPFR_RNDN); } else { mpfr_set(r26898, r26897, MPFR_RNDN); }; return mpfr_get_d(r26898, MPFR_RNDN); }