#include #include #include #include #include char *name = "Kahan's exp quotient"; double f_if(float x) { float r23719 = x; float r23720 = exp(r23719); float r23721 = 1; float r23722 = r23720 - r23721; float r23723 = r23722 / r23719; return r23723; } double f_id(double x) { double r23724 = x; double r23725 = exp(r23724); double r23726 = 1; double r23727 = r23725 - r23726; double r23728 = r23727 / r23724; return r23728; } double f_of(float x) { float r23729 = x; float r23730 = -0.00013592676717935897; bool r23731 = r23729 <= r23730; float r23732 = exp(r23729); float r23733 = 1; float r23734 = r23732 - r23733; float r23735 = cbrt(r23734); float r23736 = r23735 * r23735; float r23737 = cbrt(r23736); float r23738 = cbrt(r23735); float r23739 = r23737 * r23738; float r23740 = r23736 * r23739; float r23741 = r23740 / r23729; float r23742 = 1/6; float r23743 = 2; float r23744 = pow(r23729, r23743); float r23745 = r23742 * r23744; float r23746 = 1/2; float r23747 = r23746 * r23729; float r23748 = r23733 + r23747; float r23749 = r23745 + r23748; float r23750 = r23731 ? r23741 : r23749; return r23750; } double f_od(double x) { double r23751 = x; double r23752 = -0.00013592676717935897; bool r23753 = r23751 <= r23752; double r23754 = exp(r23751); double r23755 = 1; double r23756 = r23754 - r23755; double r23757 = cbrt(r23756); double r23758 = r23757 * r23757; double r23759 = cbrt(r23758); double r23760 = cbrt(r23757); double r23761 = r23759 * r23760; double r23762 = r23758 * r23761; double r23763 = r23762 / r23751; double r23764 = 1/6; double r23765 = 2; double r23766 = pow(r23751, r23765); double r23767 = r23764 * r23766; double r23768 = 1/2; double r23769 = r23768 * r23751; double r23770 = r23755 + r23769; double r23771 = r23767 + r23770; double r23772 = r23753 ? r23763 : r23771; return r23772; } 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 r23773, r23774, r23775, r23776, r23777; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r23773); mpfr_init(r23774); mpfr_init_set_str(r23775, "1", 10, MPFR_RNDN); mpfr_init(r23776); mpfr_init(r23777); } double f_im(double x) { mpfr_set_d(r23773, x, MPFR_RNDN); mpfr_exp(r23774, r23773, MPFR_RNDN); ; mpfr_sub(r23776, r23774, r23775, MPFR_RNDN); mpfr_div(r23777, r23776, r23773, MPFR_RNDN); return mpfr_get_d(r23777, MPFR_RNDN); } static mpfr_t r23778, r23779, r23780, r23781, r23782, r23783, r23784, r23785, r23786, r23787, r23788, r23789, r23790, r23791, r23792, r23793, r23794, r23795, r23796, r23797, r23798, r23799; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r23778); mpfr_init_set_str(r23779, "-0.00013592676717935897", 10, MPFR_RNDN); mpfr_init(r23780); mpfr_init(r23781); mpfr_init_set_str(r23782, "1", 10, MPFR_RNDN); mpfr_init(r23783); mpfr_init(r23784); mpfr_init(r23785); mpfr_init(r23786); mpfr_init(r23787); mpfr_init(r23788); mpfr_init(r23789); mpfr_init(r23790); mpfr_init_set_str(r23791, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r23792, "2", 10, MPFR_RNDN); mpfr_init(r23793); mpfr_init(r23794); mpfr_init_set_str(r23795, "1/2", 10, MPFR_RNDN); mpfr_init(r23796); mpfr_init(r23797); mpfr_init(r23798); mpfr_init(r23799); } double f_fm(double x) { mpfr_set_d(r23778, x, MPFR_RNDN); ; mpfr_set_si(r23780, mpfr_cmp(r23778, r23779) <= 0, MPFR_RNDN); mpfr_exp(r23781, r23778, MPFR_RNDN); ; mpfr_sub(r23783, r23781, r23782, MPFR_RNDN); mpfr_cbrt(r23784, r23783, MPFR_RNDN); mpfr_mul(r23785, r23784, r23784, MPFR_RNDN); mpfr_cbrt(r23786, r23785, MPFR_RNDN); mpfr_cbrt(r23787, r23784, MPFR_RNDN); mpfr_mul(r23788, r23786, r23787, MPFR_RNDN); mpfr_mul(r23789, r23785, r23788, MPFR_RNDN); mpfr_div(r23790, r23789, r23778, MPFR_RNDN); ; ; mpfr_pow(r23793, r23778, r23792, MPFR_RNDN); mpfr_mul(r23794, r23791, r23793, MPFR_RNDN); ; mpfr_mul(r23796, r23795, r23778, MPFR_RNDN); mpfr_add(r23797, r23782, r23796, MPFR_RNDN); mpfr_add(r23798, r23794, r23797, MPFR_RNDN); if (mpfr_get_si(r23780, MPFR_RNDN)) { mpfr_set(r23799, r23790, MPFR_RNDN); } else { mpfr_set(r23799, r23798, MPFR_RNDN); }; return mpfr_get_d(r23799, MPFR_RNDN); } static mpfr_t r23800, r23801, r23802, r23803, r23804, r23805, r23806, r23807, r23808, r23809, r23810, r23811, r23812, r23813, r23814, r23815, r23816, r23817, r23818, r23819, r23820, r23821; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r23800); mpfr_init_set_str(r23801, "-0.00013592676717935897", 10, MPFR_RNDN); mpfr_init(r23802); mpfr_init(r23803); mpfr_init_set_str(r23804, "1", 10, MPFR_RNDN); mpfr_init(r23805); mpfr_init(r23806); mpfr_init(r23807); mpfr_init(r23808); mpfr_init(r23809); mpfr_init(r23810); mpfr_init(r23811); mpfr_init(r23812); mpfr_init_set_str(r23813, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r23814, "2", 10, MPFR_RNDN); mpfr_init(r23815); mpfr_init(r23816); mpfr_init_set_str(r23817, "1/2", 10, MPFR_RNDN); mpfr_init(r23818); mpfr_init(r23819); mpfr_init(r23820); mpfr_init(r23821); } double f_dm(double x) { mpfr_set_d(r23800, x, MPFR_RNDN); ; mpfr_set_si(r23802, mpfr_cmp(r23800, r23801) <= 0, MPFR_RNDN); mpfr_exp(r23803, r23800, MPFR_RNDN); ; mpfr_sub(r23805, r23803, r23804, MPFR_RNDN); mpfr_cbrt(r23806, r23805, MPFR_RNDN); mpfr_mul(r23807, r23806, r23806, MPFR_RNDN); mpfr_cbrt(r23808, r23807, MPFR_RNDN); mpfr_cbrt(r23809, r23806, MPFR_RNDN); mpfr_mul(r23810, r23808, r23809, MPFR_RNDN); mpfr_mul(r23811, r23807, r23810, MPFR_RNDN); mpfr_div(r23812, r23811, r23800, MPFR_RNDN); ; ; mpfr_pow(r23815, r23800, r23814, MPFR_RNDN); mpfr_mul(r23816, r23813, r23815, MPFR_RNDN); ; mpfr_mul(r23818, r23817, r23800, MPFR_RNDN); mpfr_add(r23819, r23804, r23818, MPFR_RNDN); mpfr_add(r23820, r23816, r23819, MPFR_RNDN); if (mpfr_get_si(r23802, MPFR_RNDN)) { mpfr_set(r23821, r23812, MPFR_RNDN); } else { mpfr_set(r23821, r23820, MPFR_RNDN); }; return mpfr_get_d(r23821, MPFR_RNDN); }