#include #include #include #include #include char *name = "Henrywood and Agarwal, Equation (12)"; double f_if(float d, float h, float l, float M, float D) { float r26511 = d; float r26512 = h; float r26513 = r26511 / r26512; float r26514 = 1; float r26515 = 2; float r26516 = r26514 / r26515; float r26517 = pow(r26513, r26516); float r26518 = l; float r26519 = r26511 / r26518; float r26520 = pow(r26519, r26516); float r26521 = r26517 * r26520; float r26522 = M; float r26523 = D; float r26524 = r26522 * r26523; float r26525 = r26515 * r26511; float r26526 = r26524 / r26525; float r26527 = pow(r26526, r26515); float r26528 = r26516 * r26527; float r26529 = r26512 / r26518; float r26530 = r26528 * r26529; float r26531 = r26514 - r26530; float r26532 = r26521 * r26531; return r26532; } double f_id(double d, double h, double l, double M, double D) { double r26533 = d; double r26534 = h; double r26535 = r26533 / r26534; double r26536 = 1; double r26537 = 2; double r26538 = r26536 / r26537; double r26539 = pow(r26535, r26538); double r26540 = l; double r26541 = r26533 / r26540; double r26542 = pow(r26541, r26538); double r26543 = r26539 * r26542; double r26544 = M; double r26545 = D; double r26546 = r26544 * r26545; double r26547 = r26537 * r26533; double r26548 = r26546 / r26547; double r26549 = pow(r26548, r26537); double r26550 = r26538 * r26549; double r26551 = r26534 / r26540; double r26552 = r26550 * r26551; double r26553 = r26536 - r26552; double r26554 = r26543 * r26553; return r26554; } double f_of(float d, float h, float l, float M, float D) { float r26555 = 1; float r26556 = h; float r26557 = D; float r26558 = M; float r26559 = r26557 * r26558; float r26560 = d; float r26561 = 2; float r26562 = r26560 * r26561; float r26563 = r26559 / r26562; float r26564 = pow(r26563, r26561); float r26565 = r26556 * r26564; float r26566 = l; float r26567 = r26566 * r26561; float r26568 = r26565 / r26567; float r26569 = r26555 - r26568; float r26570 = 1/2; float r26571 = -1; float r26572 = r26571 / r26566; float r26573 = log(r26572); float r26574 = r26571 / r26560; float r26575 = log(r26574); float r26576 = r26573 - r26575; float r26577 = r26570 * r26576; float r26578 = exp(r26577); float r26579 = r26560 / r26556; float r26580 = r26555 / r26561; float r26581 = pow(r26579, r26580); float r26582 = r26578 * r26581; float r26583 = r26569 * r26582; float r26584 = 2.5015268726055398e-259; bool r26585 = r26583 <= r26584; float r26586 = r26556 / r26566; float r26587 = r26580 * r26564; float r26588 = r26586 * r26587; float r26589 = r26555 - r26588; float r26590 = r26582 * r26589; float r26591 = 3.161821923353777e+174; bool r26592 = r26583 <= r26591; float r26593 = r26560 / r26566; float r26594 = cbrt(r26593); float r26595 = fabs(r26594); float r26596 = sqrt(r26594); float r26597 = r26595 * r26596; float r26598 = pow(r26597, r26580); float r26599 = sqrt(r26593); float r26600 = pow(r26599, r26580); float r26601 = r26598 * r26600; float r26602 = r26601 * r26581; float r26603 = r26569 * r26602; float r26604 = 8.188519636352656e+285; bool r26605 = r26583 <= r26604; float r26606 = r26555 / r26566; float r26607 = pow(r26606, r26580); float r26608 = pow(r26560, r26580); float r26609 = r26607 * r26608; float r26610 = r26581 * r26609; float r26611 = r26569 * r26610; float r26612 = r26605 ? r26583 : r26611; float r26613 = r26592 ? r26603 : r26612; float r26614 = r26585 ? r26590 : r26613; return r26614; } double f_od(double d, double h, double l, double M, double D) { double r26615 = 1; double r26616 = h; double r26617 = D; double r26618 = M; double r26619 = r26617 * r26618; double r26620 = d; double r26621 = 2; double r26622 = r26620 * r26621; double r26623 = r26619 / r26622; double r26624 = pow(r26623, r26621); double r26625 = r26616 * r26624; double r26626 = l; double r26627 = r26626 * r26621; double r26628 = r26625 / r26627; double r26629 = r26615 - r26628; double r26630 = 1/2; double r26631 = -1; double r26632 = r26631 / r26626; double r26633 = log(r26632); double r26634 = r26631 / r26620; double r26635 = log(r26634); double r26636 = r26633 - r26635; double r26637 = r26630 * r26636; double r26638 = exp(r26637); double r26639 = r26620 / r26616; double r26640 = r26615 / r26621; double r26641 = pow(r26639, r26640); double r26642 = r26638 * r26641; double r26643 = r26629 * r26642; double r26644 = 2.5015268726055398e-259; bool r26645 = r26643 <= r26644; double r26646 = r26616 / r26626; double r26647 = r26640 * r26624; double r26648 = r26646 * r26647; double r26649 = r26615 - r26648; double r26650 = r26642 * r26649; double r26651 = 3.161821923353777e+174; bool r26652 = r26643 <= r26651; double r26653 = r26620 / r26626; double r26654 = cbrt(r26653); double r26655 = fabs(r26654); double r26656 = sqrt(r26654); double r26657 = r26655 * r26656; double r26658 = pow(r26657, r26640); double r26659 = sqrt(r26653); double r26660 = pow(r26659, r26640); double r26661 = r26658 * r26660; double r26662 = r26661 * r26641; double r26663 = r26629 * r26662; double r26664 = 8.188519636352656e+285; bool r26665 = r26643 <= r26664; double r26666 = r26615 / r26626; double r26667 = pow(r26666, r26640); double r26668 = pow(r26620, r26640); double r26669 = r26667 * r26668; double r26670 = r26641 * r26669; double r26671 = r26629 * r26670; double r26672 = r26665 ? r26643 : r26671; double r26673 = r26652 ? r26663 : r26672; double r26674 = r26645 ? r26650 : r26673; return r26674; } 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 r26675, r26676, r26677, r26678, r26679, r26680, r26681, r26682, r26683, r26684, r26685, r26686, r26687, r26688, r26689, r26690, r26691, r26692, r26693, r26694, r26695, r26696; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26675); mpfr_init(r26676); mpfr_init(r26677); mpfr_init_set_str(r26678, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26679, "2", 10, MPFR_RNDN); 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(r26689); mpfr_init(r26690); mpfr_init(r26691); mpfr_init(r26692); mpfr_init(r26693); mpfr_init(r26694); mpfr_init(r26695); mpfr_init(r26696); } double f_im(double d, double h, double l, double M, double D) { mpfr_set_d(r26675, d, MPFR_RNDN); mpfr_set_d(r26676, h, MPFR_RNDN); mpfr_div(r26677, r26675, r26676, MPFR_RNDN); ; ; mpfr_div(r26680, r26678, r26679, MPFR_RNDN); mpfr_pow(r26681, r26677, r26680, MPFR_RNDN); mpfr_set_d(r26682, l, MPFR_RNDN); mpfr_div(r26683, r26675, r26682, MPFR_RNDN); mpfr_pow(r26684, r26683, r26680, MPFR_RNDN); mpfr_mul(r26685, r26681, r26684, MPFR_RNDN); mpfr_set_d(r26686, M, MPFR_RNDN); mpfr_set_d(r26687, D, MPFR_RNDN); mpfr_mul(r26688, r26686, r26687, MPFR_RNDN); mpfr_mul(r26689, r26679, r26675, MPFR_RNDN); mpfr_div(r26690, r26688, r26689, MPFR_RNDN); mpfr_pow(r26691, r26690, r26679, MPFR_RNDN); mpfr_mul(r26692, r26680, r26691, MPFR_RNDN); mpfr_div(r26693, r26676, r26682, MPFR_RNDN); mpfr_mul(r26694, r26692, r26693, MPFR_RNDN); mpfr_sub(r26695, r26678, r26694, MPFR_RNDN); mpfr_mul(r26696, r26685, r26695, MPFR_RNDN); return mpfr_get_d(r26696, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r26697, "1", 10, MPFR_RNDN); mpfr_init(r26698); mpfr_init(r26699); mpfr_init(r26700); mpfr_init(r26701); mpfr_init(r26702); mpfr_init_set_str(r26703, "2", 10, MPFR_RNDN); 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_set_str(r26712, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r26713, "-1", 10, MPFR_RNDN); 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_set_str(r26726, "2.5015268726055398e-259", 10, MPFR_RNDN); mpfr_init(r26727); mpfr_init(r26728); mpfr_init(r26729); mpfr_init(r26730); mpfr_init(r26731); mpfr_init(r26732); mpfr_init_set_str(r26733, "3.161821923353777e+174", 10, MPFR_RNDN); 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_set_str(r26746, "8.188519636352656e+285", 10, MPFR_RNDN); 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); } double f_fm(double d, double h, double l, double M, double D) { ; mpfr_set_d(r26698, h, MPFR_RNDN); mpfr_set_d(r26699, D, MPFR_RNDN); mpfr_set_d(r26700, M, MPFR_RNDN); mpfr_mul(r26701, r26699, r26700, MPFR_RNDN); mpfr_set_d(r26702, d, MPFR_RNDN); ; mpfr_mul(r26704, r26702, r26703, MPFR_RNDN); mpfr_div(r26705, r26701, r26704, MPFR_RNDN); mpfr_pow(r26706, r26705, r26703, MPFR_RNDN); mpfr_mul(r26707, r26698, r26706, MPFR_RNDN); mpfr_set_d(r26708, l, MPFR_RNDN); mpfr_mul(r26709, r26708, r26703, MPFR_RNDN); mpfr_div(r26710, r26707, r26709, MPFR_RNDN); mpfr_sub(r26711, r26697, r26710, MPFR_RNDN); ; ; mpfr_div(r26714, r26713, r26708, MPFR_RNDN); mpfr_log(r26715, r26714, MPFR_RNDN); mpfr_div(r26716, r26713, r26702, MPFR_RNDN); mpfr_log(r26717, r26716, MPFR_RNDN); mpfr_sub(r26718, r26715, r26717, MPFR_RNDN); mpfr_mul(r26719, r26712, r26718, MPFR_RNDN); mpfr_exp(r26720, r26719, MPFR_RNDN); mpfr_div(r26721, r26702, r26698, MPFR_RNDN); mpfr_div(r26722, r26697, r26703, MPFR_RNDN); mpfr_pow(r26723, r26721, r26722, MPFR_RNDN); mpfr_mul(r26724, r26720, r26723, MPFR_RNDN); mpfr_mul(r26725, r26711, r26724, MPFR_RNDN); ; mpfr_set_si(r26727, mpfr_cmp(r26725, r26726) <= 0, MPFR_RNDN); mpfr_div(r26728, r26698, r26708, MPFR_RNDN); mpfr_mul(r26729, r26722, r26706, MPFR_RNDN); mpfr_mul(r26730, r26728, r26729, MPFR_RNDN); mpfr_sub(r26731, r26697, r26730, MPFR_RNDN); mpfr_mul(r26732, r26724, r26731, MPFR_RNDN); ; mpfr_set_si(r26734, mpfr_cmp(r26725, r26733) <= 0, MPFR_RNDN); mpfr_div(r26735, r26702, r26708, MPFR_RNDN); mpfr_cbrt(r26736, r26735, MPFR_RNDN); mpfr_abs(r26737, r26736, MPFR_RNDN); mpfr_sqrt(r26738, r26736, MPFR_RNDN); mpfr_mul(r26739, r26737, r26738, MPFR_RNDN); mpfr_pow(r26740, r26739, r26722, MPFR_RNDN); mpfr_sqrt(r26741, r26735, MPFR_RNDN); mpfr_pow(r26742, r26741, r26722, MPFR_RNDN); mpfr_mul(r26743, r26740, r26742, MPFR_RNDN); mpfr_mul(r26744, r26743, r26723, MPFR_RNDN); mpfr_mul(r26745, r26711, r26744, MPFR_RNDN); ; mpfr_set_si(r26747, mpfr_cmp(r26725, r26746) <= 0, MPFR_RNDN); mpfr_div(r26748, r26697, r26708, MPFR_RNDN); mpfr_pow(r26749, r26748, r26722, MPFR_RNDN); mpfr_pow(r26750, r26702, r26722, MPFR_RNDN); mpfr_mul(r26751, r26749, r26750, MPFR_RNDN); mpfr_mul(r26752, r26723, r26751, MPFR_RNDN); mpfr_mul(r26753, r26711, r26752, MPFR_RNDN); if (mpfr_get_si(r26747, MPFR_RNDN)) { mpfr_set(r26754, r26725, MPFR_RNDN); } else { mpfr_set(r26754, r26753, MPFR_RNDN); }; if (mpfr_get_si(r26734, MPFR_RNDN)) { mpfr_set(r26755, r26745, MPFR_RNDN); } else { mpfr_set(r26755, r26754, MPFR_RNDN); }; if (mpfr_get_si(r26727, MPFR_RNDN)) { mpfr_set(r26756, r26732, MPFR_RNDN); } else { mpfr_set(r26756, r26755, MPFR_RNDN); }; return mpfr_get_d(r26756, MPFR_RNDN); } static mpfr_t r26757, r26758, r26759, r26760, r26761, r26762, r26763, r26764, 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; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r26757, "1", 10, MPFR_RNDN); mpfr_init(r26758); mpfr_init(r26759); mpfr_init(r26760); mpfr_init(r26761); mpfr_init(r26762); mpfr_init_set_str(r26763, "2", 10, MPFR_RNDN); mpfr_init(r26764); mpfr_init(r26765); mpfr_init(r26766); mpfr_init(r26767); mpfr_init(r26768); mpfr_init(r26769); mpfr_init(r26770); mpfr_init(r26771); mpfr_init_set_str(r26772, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r26773, "-1", 10, MPFR_RNDN); 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_set_str(r26786, "2.5015268726055398e-259", 10, MPFR_RNDN); mpfr_init(r26787); mpfr_init(r26788); mpfr_init(r26789); mpfr_init(r26790); mpfr_init(r26791); mpfr_init(r26792); mpfr_init_set_str(r26793, "3.161821923353777e+174", 10, MPFR_RNDN); 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_set_str(r26806, "8.188519636352656e+285", 10, MPFR_RNDN); 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); } double f_dm(double d, double h, double l, double M, double D) { ; mpfr_set_d(r26758, h, MPFR_RNDN); mpfr_set_d(r26759, D, MPFR_RNDN); mpfr_set_d(r26760, M, MPFR_RNDN); mpfr_mul(r26761, r26759, r26760, MPFR_RNDN); mpfr_set_d(r26762, d, MPFR_RNDN); ; mpfr_mul(r26764, r26762, r26763, MPFR_RNDN); mpfr_div(r26765, r26761, r26764, MPFR_RNDN); mpfr_pow(r26766, r26765, r26763, MPFR_RNDN); mpfr_mul(r26767, r26758, r26766, MPFR_RNDN); mpfr_set_d(r26768, l, MPFR_RNDN); mpfr_mul(r26769, r26768, r26763, MPFR_RNDN); mpfr_div(r26770, r26767, r26769, MPFR_RNDN); mpfr_sub(r26771, r26757, r26770, MPFR_RNDN); ; ; mpfr_div(r26774, r26773, r26768, MPFR_RNDN); mpfr_log(r26775, r26774, MPFR_RNDN); mpfr_div(r26776, r26773, r26762, MPFR_RNDN); mpfr_log(r26777, r26776, MPFR_RNDN); mpfr_sub(r26778, r26775, r26777, MPFR_RNDN); mpfr_mul(r26779, r26772, r26778, MPFR_RNDN); mpfr_exp(r26780, r26779, MPFR_RNDN); mpfr_div(r26781, r26762, r26758, MPFR_RNDN); mpfr_div(r26782, r26757, r26763, MPFR_RNDN); mpfr_pow(r26783, r26781, r26782, MPFR_RNDN); mpfr_mul(r26784, r26780, r26783, MPFR_RNDN); mpfr_mul(r26785, r26771, r26784, MPFR_RNDN); ; mpfr_set_si(r26787, mpfr_cmp(r26785, r26786) <= 0, MPFR_RNDN); mpfr_div(r26788, r26758, r26768, MPFR_RNDN); mpfr_mul(r26789, r26782, r26766, MPFR_RNDN); mpfr_mul(r26790, r26788, r26789, MPFR_RNDN); mpfr_sub(r26791, r26757, r26790, MPFR_RNDN); mpfr_mul(r26792, r26784, r26791, MPFR_RNDN); ; mpfr_set_si(r26794, mpfr_cmp(r26785, r26793) <= 0, MPFR_RNDN); mpfr_div(r26795, r26762, r26768, MPFR_RNDN); mpfr_cbrt(r26796, r26795, MPFR_RNDN); mpfr_abs(r26797, r26796, MPFR_RNDN); mpfr_sqrt(r26798, r26796, MPFR_RNDN); mpfr_mul(r26799, r26797, r26798, MPFR_RNDN); mpfr_pow(r26800, r26799, r26782, MPFR_RNDN); mpfr_sqrt(r26801, r26795, MPFR_RNDN); mpfr_pow(r26802, r26801, r26782, MPFR_RNDN); mpfr_mul(r26803, r26800, r26802, MPFR_RNDN); mpfr_mul(r26804, r26803, r26783, MPFR_RNDN); mpfr_mul(r26805, r26771, r26804, MPFR_RNDN); ; mpfr_set_si(r26807, mpfr_cmp(r26785, r26806) <= 0, MPFR_RNDN); mpfr_div(r26808, r26757, r26768, MPFR_RNDN); mpfr_pow(r26809, r26808, r26782, MPFR_RNDN); mpfr_pow(r26810, r26762, r26782, MPFR_RNDN); mpfr_mul(r26811, r26809, r26810, MPFR_RNDN); mpfr_mul(r26812, r26783, r26811, MPFR_RNDN); mpfr_mul(r26813, r26771, r26812, MPFR_RNDN); if (mpfr_get_si(r26807, MPFR_RNDN)) { mpfr_set(r26814, r26785, MPFR_RNDN); } else { mpfr_set(r26814, r26813, MPFR_RNDN); }; if (mpfr_get_si(r26794, MPFR_RNDN)) { mpfr_set(r26815, r26805, MPFR_RNDN); } else { mpfr_set(r26815, r26814, MPFR_RNDN); }; if (mpfr_get_si(r26787, MPFR_RNDN)) { mpfr_set(r26816, r26792, MPFR_RNDN); } else { mpfr_set(r26816, r26815, MPFR_RNDN); }; return mpfr_get_d(r26816, MPFR_RNDN); }