#include #include #include #include #include char *name = "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1"; double f_if(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) { float r25427 = x; float r25428 = 18.0; float r25429 = r25427 * r25428; float r25430 = y; float r25431 = r25429 * r25430; float r25432 = z; float r25433 = r25431 * r25432; float r25434 = t; float r25435 = r25433 * r25434; float r25436 = a; float r25437 = 4.0; float r25438 = r25436 * r25437; float r25439 = r25438 * r25434; float r25440 = r25435 - r25439; float r25441 = b; float r25442 = c; float r25443 = r25441 * r25442; float r25444 = r25440 + r25443; float r25445 = r25427 * r25437; float r25446 = i; float r25447 = r25445 * r25446; float r25448 = r25444 - r25447; float r25449 = j; float r25450 = 27.0; float r25451 = r25449 * r25450; float r25452 = k; float r25453 = r25451 * r25452; float r25454 = r25448 - r25453; return r25454; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r25455 = x; double r25456 = 18.0; double r25457 = r25455 * r25456; double r25458 = y; double r25459 = r25457 * r25458; double r25460 = z; double r25461 = r25459 * r25460; double r25462 = t; double r25463 = r25461 * r25462; double r25464 = a; double r25465 = 4.0; double r25466 = r25464 * r25465; double r25467 = r25466 * r25462; double r25468 = r25463 - r25467; double r25469 = b; double r25470 = c; double r25471 = r25469 * r25470; double r25472 = r25468 + r25471; double r25473 = r25455 * r25465; double r25474 = i; double r25475 = r25473 * r25474; double r25476 = r25472 - r25475; double r25477 = j; double r25478 = 27.0; double r25479 = r25477 * r25478; double r25480 = k; double r25481 = r25479 * r25480; double r25482 = r25476 - r25481; return r25482; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) { float r25483 = c; float r25484 = b; float r25485 = r25483 * r25484; float r25486 = k; float r25487 = j; float r25488 = r25486 * r25487; float r25489 = 27.0; float r25490 = r25488 * r25489; float r25491 = r25485 - r25490; float r25492 = x; float r25493 = i; float r25494 = 4.0; float r25495 = r25493 * r25494; float r25496 = t; float r25497 = 18.0; float r25498 = r25496 * r25497; float r25499 = z; float r25500 = y; float r25501 = r25499 * r25500; float r25502 = r25498 * r25501; float r25503 = r25495 - r25502; float r25504 = r25492 * r25503; float r25505 = r25491 - r25504; float r25506 = a; float r25507 = r25494 * r25506; float r25508 = -r25507; float r25509 = r25496 * r25508; float r25510 = r25505 + r25509; float r25511 = -inf.0; bool r25512 = r25510 <= r25511; float r25513 = r25484 * r25483; float r25514 = r25492 * r25494; float r25515 = r25493 * r25514; float r25516 = r25513 - r25515; float r25517 = r25489 * r25486; float r25518 = r25487 * r25517; float r25519 = r25516 - r25518; float r25520 = r25500 * r25497; float r25521 = r25492 * r25520; float r25522 = r25499 * r25521; float r25523 = r25522 - r25507; float r25524 = r25496 * r25523; float r25525 = r25519 + r25524; float r25526 = 6.991261433239431e+287; bool r25527 = r25510 <= r25526; float r25528 = r25489 * r25487; float r25529 = r25486 * r25528; float r25530 = r25485 - r25529; float r25531 = r25496 * r25500; float r25532 = r25499 * r25497; float r25533 = r25531 * r25532; float r25534 = r25495 - r25533; float r25535 = r25492 * r25534; float r25536 = r25530 - r25535; float r25537 = r25508 * r25496; float r25538 = r25536 + r25537; float r25539 = r25527 ? r25510 : r25538; float r25540 = r25512 ? r25525 : r25539; return r25540; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r25541 = c; double r25542 = b; double r25543 = r25541 * r25542; double r25544 = k; double r25545 = j; double r25546 = r25544 * r25545; double r25547 = 27.0; double r25548 = r25546 * r25547; double r25549 = r25543 - r25548; double r25550 = x; double r25551 = i; double r25552 = 4.0; double r25553 = r25551 * r25552; double r25554 = t; double r25555 = 18.0; double r25556 = r25554 * r25555; double r25557 = z; double r25558 = y; double r25559 = r25557 * r25558; double r25560 = r25556 * r25559; double r25561 = r25553 - r25560; double r25562 = r25550 * r25561; double r25563 = r25549 - r25562; double r25564 = a; double r25565 = r25552 * r25564; double r25566 = -r25565; double r25567 = r25554 * r25566; double r25568 = r25563 + r25567; double r25569 = -inf.0; bool r25570 = r25568 <= r25569; double r25571 = r25542 * r25541; double r25572 = r25550 * r25552; double r25573 = r25551 * r25572; double r25574 = r25571 - r25573; double r25575 = r25547 * r25544; double r25576 = r25545 * r25575; double r25577 = r25574 - r25576; double r25578 = r25558 * r25555; double r25579 = r25550 * r25578; double r25580 = r25557 * r25579; double r25581 = r25580 - r25565; double r25582 = r25554 * r25581; double r25583 = r25577 + r25582; double r25584 = 6.991261433239431e+287; bool r25585 = r25568 <= r25584; double r25586 = r25547 * r25545; double r25587 = r25544 * r25586; double r25588 = r25543 - r25587; double r25589 = r25554 * r25558; double r25590 = r25557 * r25555; double r25591 = r25589 * r25590; double r25592 = r25553 - r25591; double r25593 = r25550 * r25592; double r25594 = r25588 - r25593; double r25595 = r25566 * r25554; double r25596 = r25594 + r25595; double r25597 = r25585 ? r25568 : r25596; double r25598 = r25570 ? r25583 : r25597; return r25598; } 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 r25599, r25600, r25601, r25602, r25603, r25604, r25605, r25606, r25607, r25608, r25609, r25610, r25611, r25612, r25613, r25614, r25615, r25616, r25617, r25618, r25619, r25620, r25621, r25622, r25623, r25624, r25625, r25626; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25599); mpfr_init_set_str(r25600, "18.0", 10, MPFR_RNDN); mpfr_init(r25601); mpfr_init(r25602); mpfr_init(r25603); mpfr_init(r25604); mpfr_init(r25605); mpfr_init(r25606); mpfr_init(r25607); mpfr_init(r25608); mpfr_init_set_str(r25609, "4.0", 10, MPFR_RNDN); mpfr_init(r25610); mpfr_init(r25611); mpfr_init(r25612); mpfr_init(r25613); mpfr_init(r25614); mpfr_init(r25615); mpfr_init(r25616); mpfr_init(r25617); mpfr_init(r25618); mpfr_init(r25619); mpfr_init(r25620); mpfr_init(r25621); mpfr_init_set_str(r25622, "27.0", 10, MPFR_RNDN); mpfr_init(r25623); mpfr_init(r25624); mpfr_init(r25625); mpfr_init(r25626); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { mpfr_set_d(r25599, x, MPFR_RNDN); ; mpfr_mul(r25601, r25599, r25600, MPFR_RNDN); mpfr_set_d(r25602, y, MPFR_RNDN); mpfr_mul(r25603, r25601, r25602, MPFR_RNDN); mpfr_set_d(r25604, z, MPFR_RNDN); mpfr_mul(r25605, r25603, r25604, MPFR_RNDN); mpfr_set_d(r25606, t, MPFR_RNDN); mpfr_mul(r25607, r25605, r25606, MPFR_RNDN); mpfr_set_d(r25608, a, MPFR_RNDN); ; mpfr_mul(r25610, r25608, r25609, MPFR_RNDN); mpfr_mul(r25611, r25610, r25606, MPFR_RNDN); mpfr_sub(r25612, r25607, r25611, MPFR_RNDN); mpfr_set_d(r25613, b, MPFR_RNDN); mpfr_set_d(r25614, c, MPFR_RNDN); mpfr_mul(r25615, r25613, r25614, MPFR_RNDN); mpfr_add(r25616, r25612, r25615, MPFR_RNDN); mpfr_mul(r25617, r25599, r25609, MPFR_RNDN); mpfr_set_d(r25618, i, MPFR_RNDN); mpfr_mul(r25619, r25617, r25618, MPFR_RNDN); mpfr_sub(r25620, r25616, r25619, MPFR_RNDN); mpfr_set_d(r25621, j, MPFR_RNDN); ; mpfr_mul(r25623, r25621, r25622, MPFR_RNDN); mpfr_set_d(r25624, k, MPFR_RNDN); mpfr_mul(r25625, r25623, r25624, MPFR_RNDN); mpfr_sub(r25626, r25620, r25625, MPFR_RNDN); return mpfr_get_d(r25626, MPFR_RNDN); } static mpfr_t r25627, r25628, r25629, r25630, r25631, r25632, r25633, r25634, r25635, r25636, r25637, r25638, r25639, r25640, r25641, r25642, r25643, r25644, r25645, r25646, r25647, r25648, r25649, r25650, r25651, r25652, r25653, r25654, r25655, r25656, r25657, r25658, r25659, r25660, r25661, r25662, r25663, r25664, r25665, r25666, r25667, r25668, r25669, r25670, r25671, r25672, r25673, r25674, r25675, r25676, r25677, r25678, r25679, r25680, r25681, r25682, r25683, r25684; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25627); mpfr_init(r25628); mpfr_init(r25629); mpfr_init(r25630); mpfr_init(r25631); mpfr_init(r25632); mpfr_init_set_str(r25633, "27.0", 10, MPFR_RNDN); mpfr_init(r25634); mpfr_init(r25635); mpfr_init(r25636); mpfr_init(r25637); mpfr_init_set_str(r25638, "4.0", 10, MPFR_RNDN); mpfr_init(r25639); mpfr_init(r25640); mpfr_init_set_str(r25641, "18.0", 10, MPFR_RNDN); mpfr_init(r25642); mpfr_init(r25643); mpfr_init(r25644); mpfr_init(r25645); mpfr_init(r25646); mpfr_init(r25647); mpfr_init(r25648); mpfr_init(r25649); mpfr_init(r25650); mpfr_init(r25651); mpfr_init(r25652); mpfr_init(r25653); mpfr_init(r25654); mpfr_init_set_str(r25655, "-inf.0", 10, MPFR_RNDN); mpfr_init(r25656); mpfr_init(r25657); mpfr_init(r25658); mpfr_init(r25659); mpfr_init(r25660); mpfr_init(r25661); mpfr_init(r25662); mpfr_init(r25663); mpfr_init(r25664); mpfr_init(r25665); mpfr_init(r25666); mpfr_init(r25667); mpfr_init(r25668); mpfr_init(r25669); mpfr_init_set_str(r25670, "6.991261433239431e+287", 10, MPFR_RNDN); mpfr_init(r25671); mpfr_init(r25672); mpfr_init(r25673); mpfr_init(r25674); mpfr_init(r25675); mpfr_init(r25676); mpfr_init(r25677); mpfr_init(r25678); mpfr_init(r25679); mpfr_init(r25680); mpfr_init(r25681); mpfr_init(r25682); mpfr_init(r25683); mpfr_init(r25684); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { mpfr_set_d(r25627, c, MPFR_RNDN); mpfr_set_d(r25628, b, MPFR_RNDN); mpfr_mul(r25629, r25627, r25628, MPFR_RNDN); mpfr_set_d(r25630, k, MPFR_RNDN); mpfr_set_d(r25631, j, MPFR_RNDN); mpfr_mul(r25632, r25630, r25631, MPFR_RNDN); ; mpfr_mul(r25634, r25632, r25633, MPFR_RNDN); mpfr_sub(r25635, r25629, r25634, MPFR_RNDN); mpfr_set_d(r25636, x, MPFR_RNDN); mpfr_set_d(r25637, i, MPFR_RNDN); ; mpfr_mul(r25639, r25637, r25638, MPFR_RNDN); mpfr_set_d(r25640, t, MPFR_RNDN); ; mpfr_mul(r25642, r25640, r25641, MPFR_RNDN); mpfr_set_d(r25643, z, MPFR_RNDN); mpfr_set_d(r25644, y, MPFR_RNDN); mpfr_mul(r25645, r25643, r25644, MPFR_RNDN); mpfr_mul(r25646, r25642, r25645, MPFR_RNDN); mpfr_sub(r25647, r25639, r25646, MPFR_RNDN); mpfr_mul(r25648, r25636, r25647, MPFR_RNDN); mpfr_sub(r25649, r25635, r25648, MPFR_RNDN); mpfr_set_d(r25650, a, MPFR_RNDN); mpfr_mul(r25651, r25638, r25650, MPFR_RNDN); mpfr_neg(r25652, r25651, MPFR_RNDN); mpfr_mul(r25653, r25640, r25652, MPFR_RNDN); mpfr_add(r25654, r25649, r25653, MPFR_RNDN); ; mpfr_set_si(r25656, mpfr_cmp(r25654, r25655) <= 0, MPFR_RNDN); mpfr_mul(r25657, r25628, r25627, MPFR_RNDN); mpfr_mul(r25658, r25636, r25638, MPFR_RNDN); mpfr_mul(r25659, r25637, r25658, MPFR_RNDN); mpfr_sub(r25660, r25657, r25659, MPFR_RNDN); mpfr_mul(r25661, r25633, r25630, MPFR_RNDN); mpfr_mul(r25662, r25631, r25661, MPFR_RNDN); mpfr_sub(r25663, r25660, r25662, MPFR_RNDN); mpfr_mul(r25664, r25644, r25641, MPFR_RNDN); mpfr_mul(r25665, r25636, r25664, MPFR_RNDN); mpfr_mul(r25666, r25643, r25665, MPFR_RNDN); mpfr_sub(r25667, r25666, r25651, MPFR_RNDN); mpfr_mul(r25668, r25640, r25667, MPFR_RNDN); mpfr_add(r25669, r25663, r25668, MPFR_RNDN); ; mpfr_set_si(r25671, mpfr_cmp(r25654, r25670) <= 0, MPFR_RNDN); mpfr_mul(r25672, r25633, r25631, MPFR_RNDN); mpfr_mul(r25673, r25630, r25672, MPFR_RNDN); mpfr_sub(r25674, r25629, r25673, MPFR_RNDN); mpfr_mul(r25675, r25640, r25644, MPFR_RNDN); mpfr_mul(r25676, r25643, r25641, MPFR_RNDN); mpfr_mul(r25677, r25675, r25676, MPFR_RNDN); mpfr_sub(r25678, r25639, r25677, MPFR_RNDN); mpfr_mul(r25679, r25636, r25678, MPFR_RNDN); mpfr_sub(r25680, r25674, r25679, MPFR_RNDN); mpfr_mul(r25681, r25652, r25640, MPFR_RNDN); mpfr_add(r25682, r25680, r25681, MPFR_RNDN); if (mpfr_get_si(r25671, MPFR_RNDN)) { mpfr_set(r25683, r25654, MPFR_RNDN); } else { mpfr_set(r25683, r25682, MPFR_RNDN); }; if (mpfr_get_si(r25656, MPFR_RNDN)) { mpfr_set(r25684, r25669, MPFR_RNDN); } else { mpfr_set(r25684, r25683, MPFR_RNDN); }; return mpfr_get_d(r25684, MPFR_RNDN); } static mpfr_t r25685, r25686, r25687, r25688, r25689, r25690, r25691, r25692, r25693, r25694, r25695, r25696, r25697, r25698, r25699, r25700, r25701, r25702, r25703, r25704, r25705, r25706, r25707, r25708, r25709, r25710, r25711, r25712, r25713, r25714, r25715, r25716, r25717, r25718, r25719, r25720, r25721, r25722, r25723, r25724, r25725, r25726, r25727, r25728, r25729, r25730, r25731, r25732, r25733, r25734, r25735, r25736, r25737, r25738, r25739, r25740, r25741, r25742; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25685); mpfr_init(r25686); mpfr_init(r25687); mpfr_init(r25688); mpfr_init(r25689); mpfr_init(r25690); mpfr_init_set_str(r25691, "27.0", 10, MPFR_RNDN); mpfr_init(r25692); mpfr_init(r25693); mpfr_init(r25694); mpfr_init(r25695); mpfr_init_set_str(r25696, "4.0", 10, MPFR_RNDN); mpfr_init(r25697); mpfr_init(r25698); mpfr_init_set_str(r25699, "18.0", 10, MPFR_RNDN); mpfr_init(r25700); mpfr_init(r25701); mpfr_init(r25702); mpfr_init(r25703); mpfr_init(r25704); mpfr_init(r25705); mpfr_init(r25706); mpfr_init(r25707); mpfr_init(r25708); mpfr_init(r25709); mpfr_init(r25710); mpfr_init(r25711); mpfr_init(r25712); mpfr_init_set_str(r25713, "-inf.0", 10, MPFR_RNDN); mpfr_init(r25714); mpfr_init(r25715); mpfr_init(r25716); mpfr_init(r25717); mpfr_init(r25718); mpfr_init(r25719); mpfr_init(r25720); mpfr_init(r25721); mpfr_init(r25722); mpfr_init(r25723); mpfr_init(r25724); mpfr_init(r25725); mpfr_init(r25726); mpfr_init(r25727); mpfr_init_set_str(r25728, "6.991261433239431e+287", 10, MPFR_RNDN); mpfr_init(r25729); mpfr_init(r25730); mpfr_init(r25731); mpfr_init(r25732); mpfr_init(r25733); mpfr_init(r25734); mpfr_init(r25735); mpfr_init(r25736); mpfr_init(r25737); mpfr_init(r25738); mpfr_init(r25739); mpfr_init(r25740); mpfr_init(r25741); mpfr_init(r25742); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { mpfr_set_d(r25685, c, MPFR_RNDN); mpfr_set_d(r25686, b, MPFR_RNDN); mpfr_mul(r25687, r25685, r25686, MPFR_RNDN); mpfr_set_d(r25688, k, MPFR_RNDN); mpfr_set_d(r25689, j, MPFR_RNDN); mpfr_mul(r25690, r25688, r25689, MPFR_RNDN); ; mpfr_mul(r25692, r25690, r25691, MPFR_RNDN); mpfr_sub(r25693, r25687, r25692, MPFR_RNDN); mpfr_set_d(r25694, x, MPFR_RNDN); mpfr_set_d(r25695, i, MPFR_RNDN); ; mpfr_mul(r25697, r25695, r25696, MPFR_RNDN); mpfr_set_d(r25698, t, MPFR_RNDN); ; mpfr_mul(r25700, r25698, r25699, MPFR_RNDN); mpfr_set_d(r25701, z, MPFR_RNDN); mpfr_set_d(r25702, y, MPFR_RNDN); mpfr_mul(r25703, r25701, r25702, MPFR_RNDN); mpfr_mul(r25704, r25700, r25703, MPFR_RNDN); mpfr_sub(r25705, r25697, r25704, MPFR_RNDN); mpfr_mul(r25706, r25694, r25705, MPFR_RNDN); mpfr_sub(r25707, r25693, r25706, MPFR_RNDN); mpfr_set_d(r25708, a, MPFR_RNDN); mpfr_mul(r25709, r25696, r25708, MPFR_RNDN); mpfr_neg(r25710, r25709, MPFR_RNDN); mpfr_mul(r25711, r25698, r25710, MPFR_RNDN); mpfr_add(r25712, r25707, r25711, MPFR_RNDN); ; mpfr_set_si(r25714, mpfr_cmp(r25712, r25713) <= 0, MPFR_RNDN); mpfr_mul(r25715, r25686, r25685, MPFR_RNDN); mpfr_mul(r25716, r25694, r25696, MPFR_RNDN); mpfr_mul(r25717, r25695, r25716, MPFR_RNDN); mpfr_sub(r25718, r25715, r25717, MPFR_RNDN); mpfr_mul(r25719, r25691, r25688, MPFR_RNDN); mpfr_mul(r25720, r25689, r25719, MPFR_RNDN); mpfr_sub(r25721, r25718, r25720, MPFR_RNDN); mpfr_mul(r25722, r25702, r25699, MPFR_RNDN); mpfr_mul(r25723, r25694, r25722, MPFR_RNDN); mpfr_mul(r25724, r25701, r25723, MPFR_RNDN); mpfr_sub(r25725, r25724, r25709, MPFR_RNDN); mpfr_mul(r25726, r25698, r25725, MPFR_RNDN); mpfr_add(r25727, r25721, r25726, MPFR_RNDN); ; mpfr_set_si(r25729, mpfr_cmp(r25712, r25728) <= 0, MPFR_RNDN); mpfr_mul(r25730, r25691, r25689, MPFR_RNDN); mpfr_mul(r25731, r25688, r25730, MPFR_RNDN); mpfr_sub(r25732, r25687, r25731, MPFR_RNDN); mpfr_mul(r25733, r25698, r25702, MPFR_RNDN); mpfr_mul(r25734, r25701, r25699, MPFR_RNDN); mpfr_mul(r25735, r25733, r25734, MPFR_RNDN); mpfr_sub(r25736, r25697, r25735, MPFR_RNDN); mpfr_mul(r25737, r25694, r25736, MPFR_RNDN); mpfr_sub(r25738, r25732, r25737, MPFR_RNDN); mpfr_mul(r25739, r25710, r25698, MPFR_RNDN); mpfr_add(r25740, r25738, r25739, MPFR_RNDN); if (mpfr_get_si(r25729, MPFR_RNDN)) { mpfr_set(r25741, r25712, MPFR_RNDN); } else { mpfr_set(r25741, r25740, MPFR_RNDN); }; if (mpfr_get_si(r25714, MPFR_RNDN)) { mpfr_set(r25742, r25727, MPFR_RNDN); } else { mpfr_set(r25742, r25741, MPFR_RNDN); }; return mpfr_get_d(r25742, MPFR_RNDN); }