#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r17581 = a; float r17582 = r17581 * r17581; float r17583 = b; float r17584 = r17583 * r17583; float r17585 = r17582 + r17584; float r17586 = r17585 * r17585; float r17587 = 4.0f; float r17588 = 1.0f; float r17589 = r17588 + r17581; float r17590 = r17582 * r17589; float r17591 = 3.0f; float r17592 = r17591 * r17581; float r17593 = r17588 - r17592; float r17594 = r17584 * r17593; float r17595 = r17590 + r17594; float r17596 = r17587 * r17595; float r17597 = r17586 + r17596; float r17598 = r17597 - r17588; return r17598; } double f_id(double a, double b) { double r17599 = a; double r17600 = r17599 * r17599; double r17601 = b; double r17602 = r17601 * r17601; double r17603 = r17600 + r17602; double r17604 = r17603 * r17603; double r17605 = 4.0; double r17606 = 1.0; double r17607 = r17606 + r17599; double r17608 = r17600 * r17607; double r17609 = 3.0; double r17610 = r17609 * r17599; double r17611 = r17606 - r17610; double r17612 = r17602 * r17611; double r17613 = r17608 + r17612; double r17614 = r17605 * r17613; double r17615 = r17604 + r17614; double r17616 = r17615 - r17606; return r17616; } double f_of(float a, float b) { float r17617 = b; float r17618 = 4.0f; float r17619 = pow(r17617, r17618); float r17620 = 2.0f; float r17621 = r17617 * r17617; float r17622 = a; float r17623 = r17622 * r17622; float r17624 = r17621 * r17623; float r17625 = r17620 * r17624; float r17626 = pow(r17622, r17618); float r17627 = r17625 + r17626; float r17628 = r17619 + r17627; float r17629 = 1.0f; float r17630 = r17629 + r17622; float r17631 = r17623 * r17630; float r17632 = 3.0f; float r17633 = r17632 * r17622; float r17634 = r17629 - r17633; float r17635 = r17621 * r17634; float r17636 = r17631 + r17635; float r17637 = r17618 * r17636; float r17638 = r17628 + r17637; float r17639 = r17638 - r17629; return r17639; } double f_od(double a, double b) { double r17640 = b; double r17641 = 4.0; double r17642 = pow(r17640, r17641); double r17643 = 2.0; double r17644 = r17640 * r17640; double r17645 = a; double r17646 = r17645 * r17645; double r17647 = r17644 * r17646; double r17648 = r17643 * r17647; double r17649 = pow(r17645, r17641); double r17650 = r17648 + r17649; double r17651 = r17642 + r17650; double r17652 = 1.0; double r17653 = r17652 + r17645; double r17654 = r17646 * r17653; double r17655 = 3.0; double r17656 = r17655 * r17645; double r17657 = r17652 - r17656; double r17658 = r17644 * r17657; double r17659 = r17654 + r17658; double r17660 = r17641 * r17659; double r17661 = r17651 + r17660; double r17662 = r17661 - r17652; return r17662; } 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 r17663, r17664, r17665, r17666, r17667, r17668, r17669, r17670, r17671, r17672, r17673, r17674, r17675, r17676, r17677, r17678, r17679, r17680; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17663); mpfr_init(r17664); mpfr_init(r17665); mpfr_init(r17666); mpfr_init(r17667); mpfr_init(r17668); mpfr_init_set_str(r17669, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17670, "1", 10, MPFR_RNDN); mpfr_init(r17671); mpfr_init(r17672); mpfr_init_set_str(r17673, "3", 10, MPFR_RNDN); mpfr_init(r17674); mpfr_init(r17675); mpfr_init(r17676); mpfr_init(r17677); mpfr_init(r17678); mpfr_init(r17679); mpfr_init(r17680); } double f_im(double a, double b) { mpfr_set_d(r17663, a, MPFR_RNDN); mpfr_sqr(r17664, r17663, MPFR_RNDN); mpfr_set_d(r17665, b, MPFR_RNDN); mpfr_sqr(r17666, r17665, MPFR_RNDN); mpfr_add(r17667, r17664, r17666, MPFR_RNDN); mpfr_sqr(r17668, r17667, MPFR_RNDN); ; ; mpfr_add(r17671, r17670, r17663, MPFR_RNDN); mpfr_mul(r17672, r17664, r17671, MPFR_RNDN); ; mpfr_mul(r17674, r17673, r17663, MPFR_RNDN); mpfr_sub(r17675, r17670, r17674, MPFR_RNDN); mpfr_mul(r17676, r17666, r17675, MPFR_RNDN); mpfr_add(r17677, r17672, r17676, MPFR_RNDN); mpfr_mul(r17678, r17669, r17677, MPFR_RNDN); mpfr_add(r17679, r17668, r17678, MPFR_RNDN); mpfr_sub(r17680, r17679, r17670, MPFR_RNDN); return mpfr_get_d(r17680, MPFR_RNDN); } static mpfr_t r17681, r17682, r17683, r17684, r17685, r17686, r17687, r17688, r17689, r17690, r17691, r17692, r17693, r17694, r17695, r17696, r17697, r17698, r17699, r17700, r17701, r17702, r17703; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17681); mpfr_init_set_str(r17682, "4", 10, MPFR_RNDN); mpfr_init(r17683); mpfr_init_set_str(r17684, "2", 10, MPFR_RNDN); mpfr_init(r17685); mpfr_init(r17686); mpfr_init(r17687); mpfr_init(r17688); mpfr_init(r17689); mpfr_init(r17690); mpfr_init(r17691); mpfr_init(r17692); mpfr_init_set_str(r17693, "1", 10, MPFR_RNDN); mpfr_init(r17694); mpfr_init(r17695); mpfr_init_set_str(r17696, "3", 10, MPFR_RNDN); mpfr_init(r17697); mpfr_init(r17698); mpfr_init(r17699); mpfr_init(r17700); mpfr_init(r17701); mpfr_init(r17702); mpfr_init(r17703); } double f_fm(double a, double b) { mpfr_set_d(r17681, b, MPFR_RNDN); ; mpfr_pow(r17683, r17681, r17682, MPFR_RNDN); ; mpfr_sqr(r17685, r17681, MPFR_RNDN); mpfr_set_d(r17686, a, MPFR_RNDN); mpfr_sqr(r17687, r17686, MPFR_RNDN); mpfr_mul(r17688, r17685, r17687, MPFR_RNDN); mpfr_mul(r17689, r17684, r17688, MPFR_RNDN); mpfr_pow(r17690, r17686, r17682, MPFR_RNDN); mpfr_add(r17691, r17689, r17690, MPFR_RNDN); mpfr_add(r17692, r17683, r17691, MPFR_RNDN); ; mpfr_add(r17694, r17693, r17686, MPFR_RNDN); mpfr_mul(r17695, r17687, r17694, MPFR_RNDN); ; mpfr_mul(r17697, r17696, r17686, MPFR_RNDN); mpfr_sub(r17698, r17693, r17697, MPFR_RNDN); mpfr_mul(r17699, r17685, r17698, MPFR_RNDN); mpfr_add(r17700, r17695, r17699, MPFR_RNDN); mpfr_mul(r17701, r17682, r17700, MPFR_RNDN); mpfr_add(r17702, r17692, r17701, MPFR_RNDN); mpfr_sub(r17703, r17702, r17693, MPFR_RNDN); return mpfr_get_d(r17703, MPFR_RNDN); } static mpfr_t r17704, r17705, r17706, r17707, r17708, r17709, r17710, r17711, r17712, r17713, r17714, r17715, r17716, r17717, r17718, r17719, r17720, r17721, r17722, r17723, r17724, r17725, r17726; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17704); mpfr_init_set_str(r17705, "4", 10, MPFR_RNDN); mpfr_init(r17706); mpfr_init_set_str(r17707, "2", 10, MPFR_RNDN); mpfr_init(r17708); mpfr_init(r17709); mpfr_init(r17710); mpfr_init(r17711); mpfr_init(r17712); mpfr_init(r17713); mpfr_init(r17714); mpfr_init(r17715); mpfr_init_set_str(r17716, "1", 10, MPFR_RNDN); mpfr_init(r17717); mpfr_init(r17718); mpfr_init_set_str(r17719, "3", 10, MPFR_RNDN); mpfr_init(r17720); mpfr_init(r17721); mpfr_init(r17722); mpfr_init(r17723); mpfr_init(r17724); mpfr_init(r17725); mpfr_init(r17726); } double f_dm(double a, double b) { mpfr_set_d(r17704, b, MPFR_RNDN); ; mpfr_pow(r17706, r17704, r17705, MPFR_RNDN); ; mpfr_sqr(r17708, r17704, MPFR_RNDN); mpfr_set_d(r17709, a, MPFR_RNDN); mpfr_sqr(r17710, r17709, MPFR_RNDN); mpfr_mul(r17711, r17708, r17710, MPFR_RNDN); mpfr_mul(r17712, r17707, r17711, MPFR_RNDN); mpfr_pow(r17713, r17709, r17705, MPFR_RNDN); mpfr_add(r17714, r17712, r17713, MPFR_RNDN); mpfr_add(r17715, r17706, r17714, MPFR_RNDN); ; mpfr_add(r17717, r17716, r17709, MPFR_RNDN); mpfr_mul(r17718, r17710, r17717, MPFR_RNDN); ; mpfr_mul(r17720, r17719, r17709, MPFR_RNDN); mpfr_sub(r17721, r17716, r17720, MPFR_RNDN); mpfr_mul(r17722, r17708, r17721, MPFR_RNDN); mpfr_add(r17723, r17718, r17722, MPFR_RNDN); mpfr_mul(r17724, r17705, r17723, MPFR_RNDN); mpfr_add(r17725, r17715, r17724, MPFR_RNDN); mpfr_sub(r17726, r17725, r17716, MPFR_RNDN); return mpfr_get_d(r17726, MPFR_RNDN); }