#include #include #include #include #include char *name = "expax (section 3.5)"; double f_if(float a, float x) { float r26894 = a; float r26895 = x; float r26896 = r26894 * r26895; float r26897 = exp(r26896); float r26898 = 1; float r26899 = r26897 - r26898; return r26899; } double f_id(double a, double x) { double r26900 = a; double r26901 = x; double r26902 = r26900 * r26901; double r26903 = exp(r26902); double r26904 = 1; double r26905 = r26903 - r26904; return r26905; } double f_of(float a, float x) { float r26906 = a; float r26907 = x; float r26908 = r26906 * r26907; float r26909 = exp(r26908); float r26910 = 1; float r26911 = r26909 - r26910; float r26912 = -1.8436612778273367e-08; bool r26913 = r26911 <= r26912; float r26914 = sqrt(r26909); float r26915 = r26914 + r26910; float r26916 = r26914 - r26910; float r26917 = r26915 * r26916; float r26918 = r26908 * r26908; float r26919 = 1/6; float r26920 = r26908 * r26919; float r26921 = 1/2; float r26922 = r26920 + r26921; float r26923 = r26918 * r26922; float r26924 = r26923 + r26908; float r26925 = r26913 ? r26917 : r26924; return r26925; } double f_od(double a, double x) { double r26926 = a; double r26927 = x; double r26928 = r26926 * r26927; double r26929 = exp(r26928); double r26930 = 1; double r26931 = r26929 - r26930; double r26932 = -1.8436612778273367e-08; bool r26933 = r26931 <= r26932; double r26934 = sqrt(r26929); double r26935 = r26934 + r26930; double r26936 = r26934 - r26930; double r26937 = r26935 * r26936; double r26938 = r26928 * r26928; double r26939 = 1/6; double r26940 = r26928 * r26939; double r26941 = 1/2; double r26942 = r26940 + r26941; double r26943 = r26938 * r26942; double r26944 = r26943 + r26928; double r26945 = r26933 ? r26937 : r26944; return r26945; } 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 r26946, r26947, r26948, r26949, r26950, r26951; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r26946); mpfr_init(r26947); mpfr_init(r26948); mpfr_init(r26949); mpfr_init_set_str(r26950, "1", 10, MPFR_RNDN); mpfr_init(r26951); } double f_im(double a, double x) { mpfr_set_d(r26946, a, MPFR_RNDN); mpfr_set_d(r26947, x, MPFR_RNDN); mpfr_mul(r26948, r26946, r26947, MPFR_RNDN); mpfr_exp(r26949, r26948, MPFR_RNDN); ; mpfr_sub(r26951, r26949, r26950, MPFR_RNDN); return mpfr_get_d(r26951, MPFR_RNDN); } static mpfr_t r26952, r26953, r26954, r26955, r26956, r26957, r26958, r26959, r26960, r26961, r26962, r26963, r26964, r26965, r26966, r26967, r26968, r26969, r26970, r26971; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r26952); mpfr_init(r26953); mpfr_init(r26954); mpfr_init(r26955); mpfr_init_set_str(r26956, "1", 10, MPFR_RNDN); mpfr_init(r26957); mpfr_init_set_str(r26958, "-1.8436612778273367e-08", 10, MPFR_RNDN); mpfr_init(r26959); mpfr_init(r26960); mpfr_init(r26961); mpfr_init(r26962); mpfr_init(r26963); mpfr_init(r26964); mpfr_init_set_str(r26965, "1/6", 10, MPFR_RNDN); mpfr_init(r26966); mpfr_init_set_str(r26967, "1/2", 10, MPFR_RNDN); mpfr_init(r26968); mpfr_init(r26969); mpfr_init(r26970); mpfr_init(r26971); } double f_fm(double a, double x) { mpfr_set_d(r26952, a, MPFR_RNDN); mpfr_set_d(r26953, x, MPFR_RNDN); mpfr_mul(r26954, r26952, r26953, MPFR_RNDN); mpfr_exp(r26955, r26954, MPFR_RNDN); ; mpfr_sub(r26957, r26955, r26956, MPFR_RNDN); ; mpfr_set_si(r26959, mpfr_cmp(r26957, r26958) <= 0, MPFR_RNDN); mpfr_sqrt(r26960, r26955, MPFR_RNDN); mpfr_add(r26961, r26960, r26956, MPFR_RNDN); mpfr_sub(r26962, r26960, r26956, MPFR_RNDN); mpfr_mul(r26963, r26961, r26962, MPFR_RNDN); mpfr_mul(r26964, r26954, r26954, MPFR_RNDN); ; mpfr_mul(r26966, r26954, r26965, MPFR_RNDN); ; mpfr_add(r26968, r26966, r26967, MPFR_RNDN); mpfr_mul(r26969, r26964, r26968, MPFR_RNDN); mpfr_add(r26970, r26969, r26954, MPFR_RNDN); if (mpfr_get_si(r26959, MPFR_RNDN)) { mpfr_set(r26971, r26963, MPFR_RNDN); } else { mpfr_set(r26971, r26970, MPFR_RNDN); }; return mpfr_get_d(r26971, MPFR_RNDN); } static mpfr_t r26972, r26973, r26974, r26975, r26976, r26977, r26978, r26979, r26980, r26981, r26982, r26983, r26984, r26985, r26986, r26987, r26988, r26989, r26990, r26991; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r26972); mpfr_init(r26973); mpfr_init(r26974); mpfr_init(r26975); mpfr_init_set_str(r26976, "1", 10, MPFR_RNDN); mpfr_init(r26977); mpfr_init_set_str(r26978, "-1.8436612778273367e-08", 10, MPFR_RNDN); mpfr_init(r26979); mpfr_init(r26980); mpfr_init(r26981); mpfr_init(r26982); mpfr_init(r26983); mpfr_init(r26984); mpfr_init_set_str(r26985, "1/6", 10, MPFR_RNDN); mpfr_init(r26986); mpfr_init_set_str(r26987, "1/2", 10, MPFR_RNDN); mpfr_init(r26988); mpfr_init(r26989); mpfr_init(r26990); mpfr_init(r26991); } double f_dm(double a, double x) { mpfr_set_d(r26972, a, MPFR_RNDN); mpfr_set_d(r26973, x, MPFR_RNDN); mpfr_mul(r26974, r26972, r26973, MPFR_RNDN); mpfr_exp(r26975, r26974, MPFR_RNDN); ; mpfr_sub(r26977, r26975, r26976, MPFR_RNDN); ; mpfr_set_si(r26979, mpfr_cmp(r26977, r26978) <= 0, MPFR_RNDN); mpfr_sqrt(r26980, r26975, MPFR_RNDN); mpfr_add(r26981, r26980, r26976, MPFR_RNDN); mpfr_sub(r26982, r26980, r26976, MPFR_RNDN); mpfr_mul(r26983, r26981, r26982, MPFR_RNDN); mpfr_mul(r26984, r26974, r26974, MPFR_RNDN); ; mpfr_mul(r26986, r26974, r26985, MPFR_RNDN); ; mpfr_add(r26988, r26986, r26987, MPFR_RNDN); mpfr_mul(r26989, r26984, r26988, MPFR_RNDN); mpfr_add(r26990, r26989, r26974, MPFR_RNDN); if (mpfr_get_si(r26979, MPFR_RNDN)) { mpfr_set(r26991, r26983, MPFR_RNDN); } else { mpfr_set(r26991, r26990, MPFR_RNDN); }; return mpfr_get_d(r26991, MPFR_RNDN); }