#include #include #include #include #include char *name = "expq3 (problem 3.4.2)"; double f_if(float a, float b, float eps) { float r26584 = eps; float r26585 = a; float r26586 = b; float r26587 = r26585 + r26586; float r26588 = r26587 * r26584; float r26589 = exp(r26588); float r26590 = 1; float r26591 = r26589 - r26590; float r26592 = r26584 * r26591; float r26593 = r26585 * r26584; float r26594 = exp(r26593); float r26595 = r26594 - r26590; float r26596 = r26586 * r26584; float r26597 = exp(r26596); float r26598 = r26597 - r26590; float r26599 = r26595 * r26598; float r26600 = r26592 / r26599; return r26600; } double f_id(double a, double b, double eps) { double r26601 = eps; double r26602 = a; double r26603 = b; double r26604 = r26602 + r26603; double r26605 = r26604 * r26601; double r26606 = exp(r26605); double r26607 = 1; double r26608 = r26606 - r26607; double r26609 = r26601 * r26608; double r26610 = r26602 * r26601; double r26611 = exp(r26610); double r26612 = r26611 - r26607; double r26613 = r26603 * r26601; double r26614 = exp(r26613); double r26615 = r26614 - r26607; double r26616 = r26612 * r26615; double r26617 = r26609 / r26616; return r26617; } double f_of(float a, float b, float eps) { float r26618 = eps; float r26619 = a; float r26620 = b; float r26621 = r26619 + r26620; float r26622 = r26621 * r26618; float r26623 = exp(r26622); float r26624 = 1; float r26625 = r26623 - r26624; float r26626 = r26618 * r26625; float r26627 = r26619 * r26618; float r26628 = exp(r26627); float r26629 = r26628 - r26624; float r26630 = r26620 * r26618; float r26631 = exp(r26630); float r26632 = r26631 - r26624; float r26633 = r26629 * r26632; float r26634 = r26626 / r26633; float r26635 = -inf.0; bool r26636 = r26634 <= r26635; float r26637 = r26624 / r26620; float r26638 = r26624 / r26619; float r26639 = r26637 + r26638; float r26640 = 2.0358482214599566e-35; bool r26641 = r26634 <= r26640; float r26642 = r26641 ? r26634 : r26639; float r26643 = r26636 ? r26639 : r26642; return r26643; } double f_od(double a, double b, double eps) { double r26644 = eps; double r26645 = a; double r26646 = b; double r26647 = r26645 + r26646; double r26648 = r26647 * r26644; double r26649 = exp(r26648); double r26650 = 1; double r26651 = r26649 - r26650; double r26652 = r26644 * r26651; double r26653 = r26645 * r26644; double r26654 = exp(r26653); double r26655 = r26654 - r26650; double r26656 = r26646 * r26644; double r26657 = exp(r26656); double r26658 = r26657 - r26650; double r26659 = r26655 * r26658; double r26660 = r26652 / r26659; double r26661 = -inf.0; bool r26662 = r26660 <= r26661; double r26663 = r26650 / r26646; double r26664 = r26650 / r26645; double r26665 = r26663 + r26664; double r26666 = 2.0358482214599566e-35; bool r26667 = r26660 <= r26666; double r26668 = r26667 ? r26660 : r26665; double r26669 = r26662 ? r26665 : r26668; return r26669; } 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 r26670, r26671, r26672, r26673, r26674, r26675, r26676, r26677, r26678, r26679, r26680, r26681, r26682, r26683, r26684, r26685, r26686; void setup_mpfr_f_im() { mpfr_set_default_prec(2384); mpfr_init(r26670); mpfr_init(r26671); mpfr_init(r26672); mpfr_init(r26673); mpfr_init(r26674); mpfr_init(r26675); mpfr_init_set_str(r26676, "1", 10, MPFR_RNDN); 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); } double f_im(double a, double b, double eps) { mpfr_set_d(r26670, eps, MPFR_RNDN); mpfr_set_d(r26671, a, MPFR_RNDN); mpfr_set_d(r26672, b, MPFR_RNDN); mpfr_add(r26673, r26671, r26672, MPFR_RNDN); mpfr_mul(r26674, r26673, r26670, MPFR_RNDN); mpfr_exp(r26675, r26674, MPFR_RNDN); ; mpfr_sub(r26677, r26675, r26676, MPFR_RNDN); mpfr_mul(r26678, r26670, r26677, MPFR_RNDN); mpfr_mul(r26679, r26671, r26670, MPFR_RNDN); mpfr_exp(r26680, r26679, MPFR_RNDN); mpfr_sub(r26681, r26680, r26676, MPFR_RNDN); mpfr_mul(r26682, r26672, r26670, MPFR_RNDN); mpfr_exp(r26683, r26682, MPFR_RNDN); mpfr_sub(r26684, r26683, r26676, MPFR_RNDN); mpfr_mul(r26685, r26681, r26684, MPFR_RNDN); mpfr_div(r26686, r26678, r26685, MPFR_RNDN); return mpfr_get_d(r26686, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(2384); mpfr_init(r26687); mpfr_init(r26688); mpfr_init(r26689); mpfr_init(r26690); mpfr_init(r26691); mpfr_init(r26692); mpfr_init_set_str(r26693, "1", 10, MPFR_RNDN); 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_set_str(r26704, "-inf.0", 10, MPFR_RNDN); mpfr_init(r26705); mpfr_init(r26706); mpfr_init(r26707); mpfr_init(r26708); mpfr_init_set_str(r26709, "2.0358482214599566e-35", 10, MPFR_RNDN); mpfr_init(r26710); mpfr_init(r26711); mpfr_init(r26712); } double f_fm(double a, double b, double eps) { mpfr_set_d(r26687, eps, MPFR_RNDN); mpfr_set_d(r26688, a, MPFR_RNDN); mpfr_set_d(r26689, b, MPFR_RNDN); mpfr_add(r26690, r26688, r26689, MPFR_RNDN); mpfr_mul(r26691, r26690, r26687, MPFR_RNDN); mpfr_exp(r26692, r26691, MPFR_RNDN); ; mpfr_sub(r26694, r26692, r26693, MPFR_RNDN); mpfr_mul(r26695, r26687, r26694, MPFR_RNDN); mpfr_mul(r26696, r26688, r26687, MPFR_RNDN); mpfr_exp(r26697, r26696, MPFR_RNDN); mpfr_sub(r26698, r26697, r26693, MPFR_RNDN); mpfr_mul(r26699, r26689, r26687, MPFR_RNDN); mpfr_exp(r26700, r26699, MPFR_RNDN); mpfr_sub(r26701, r26700, r26693, MPFR_RNDN); mpfr_mul(r26702, r26698, r26701, MPFR_RNDN); mpfr_div(r26703, r26695, r26702, MPFR_RNDN); ; mpfr_set_si(r26705, mpfr_cmp(r26703, r26704) <= 0, MPFR_RNDN); mpfr_div(r26706, r26693, r26689, MPFR_RNDN); mpfr_div(r26707, r26693, r26688, MPFR_RNDN); mpfr_add(r26708, r26706, r26707, MPFR_RNDN); ; mpfr_set_si(r26710, mpfr_cmp(r26703, r26709) <= 0, MPFR_RNDN); if (mpfr_get_si(r26710, MPFR_RNDN)) { mpfr_set(r26711, r26703, MPFR_RNDN); } else { mpfr_set(r26711, r26708, MPFR_RNDN); }; if (mpfr_get_si(r26705, MPFR_RNDN)) { mpfr_set(r26712, r26708, MPFR_RNDN); } else { mpfr_set(r26712, r26711, MPFR_RNDN); }; return mpfr_get_d(r26712, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_dm() { mpfr_set_default_prec(2384); mpfr_init(r26713); mpfr_init(r26714); mpfr_init(r26715); mpfr_init(r26716); mpfr_init(r26717); mpfr_init(r26718); mpfr_init_set_str(r26719, "1", 10, MPFR_RNDN); 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_set_str(r26730, "-inf.0", 10, MPFR_RNDN); mpfr_init(r26731); mpfr_init(r26732); mpfr_init(r26733); mpfr_init(r26734); mpfr_init_set_str(r26735, "2.0358482214599566e-35", 10, MPFR_RNDN); mpfr_init(r26736); mpfr_init(r26737); mpfr_init(r26738); } double f_dm(double a, double b, double eps) { mpfr_set_d(r26713, eps, MPFR_RNDN); mpfr_set_d(r26714, a, MPFR_RNDN); mpfr_set_d(r26715, b, MPFR_RNDN); mpfr_add(r26716, r26714, r26715, MPFR_RNDN); mpfr_mul(r26717, r26716, r26713, MPFR_RNDN); mpfr_exp(r26718, r26717, MPFR_RNDN); ; mpfr_sub(r26720, r26718, r26719, MPFR_RNDN); mpfr_mul(r26721, r26713, r26720, MPFR_RNDN); mpfr_mul(r26722, r26714, r26713, MPFR_RNDN); mpfr_exp(r26723, r26722, MPFR_RNDN); mpfr_sub(r26724, r26723, r26719, MPFR_RNDN); mpfr_mul(r26725, r26715, r26713, MPFR_RNDN); mpfr_exp(r26726, r26725, MPFR_RNDN); mpfr_sub(r26727, r26726, r26719, MPFR_RNDN); mpfr_mul(r26728, r26724, r26727, MPFR_RNDN); mpfr_div(r26729, r26721, r26728, MPFR_RNDN); ; mpfr_set_si(r26731, mpfr_cmp(r26729, r26730) <= 0, MPFR_RNDN); mpfr_div(r26732, r26719, r26715, MPFR_RNDN); mpfr_div(r26733, r26719, r26714, MPFR_RNDN); mpfr_add(r26734, r26732, r26733, MPFR_RNDN); ; mpfr_set_si(r26736, mpfr_cmp(r26729, r26735) <= 0, MPFR_RNDN); if (mpfr_get_si(r26736, MPFR_RNDN)) { mpfr_set(r26737, r26729, MPFR_RNDN); } else { mpfr_set(r26737, r26734, MPFR_RNDN); }; if (mpfr_get_si(r26731, MPFR_RNDN)) { mpfr_set(r26738, r26734, MPFR_RNDN); } else { mpfr_set(r26738, r26737, MPFR_RNDN); }; return mpfr_get_d(r26738, MPFR_RNDN); }