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