#include #include #include #include #include char *name = "Harley's example"; double f_if(float c_p, float c_n, float t, float s) { float r24900 = 1; float r24901 = s; float r24902 = -r24901; float r24903 = exp(r24902); float r24904 = r24900 + r24903; float r24905 = r24900 / r24904; float r24906 = c_p; float r24907 = pow(r24905, r24906); float r24908 = r24900 - r24905; float r24909 = c_n; float r24910 = pow(r24908, r24909); float r24911 = r24907 * r24910; float r24912 = t; float r24913 = -r24912; float r24914 = exp(r24913); float r24915 = r24900 + r24914; float r24916 = r24900 / r24915; float r24917 = pow(r24916, r24906); float r24918 = r24900 - r24916; float r24919 = pow(r24918, r24909); float r24920 = r24917 * r24919; float r24921 = r24911 / r24920; return r24921; } double f_id(double c_p, double c_n, double t, double s) { double r24922 = 1; double r24923 = s; double r24924 = -r24923; double r24925 = exp(r24924); double r24926 = r24922 + r24925; double r24927 = r24922 / r24926; double r24928 = c_p; double r24929 = pow(r24927, r24928); double r24930 = r24922 - r24927; double r24931 = c_n; double r24932 = pow(r24930, r24931); double r24933 = r24929 * r24932; double r24934 = t; double r24935 = -r24934; double r24936 = exp(r24935); double r24937 = r24922 + r24936; double r24938 = r24922 / r24937; double r24939 = pow(r24938, r24928); double r24940 = r24922 - r24938; double r24941 = pow(r24940, r24931); double r24942 = r24939 * r24941; double r24943 = r24933 / r24942; return r24943; } double f_of(float c_p, float c_n, float t, float s) { float r24944 = 1; float r24945 = -r24944; float r24946 = s; float r24947 = -r24946; float r24948 = exp(r24947); float r24949 = r24944 + r24948; float r24950 = r24945 / r24949; float r24951 = log1p(r24950); float r24952 = t; float r24953 = -r24952; float r24954 = exp(r24953); float r24955 = r24954 + r24944; float r24956 = r24945 / r24955; float r24957 = log1p(r24956); float r24958 = r24951 - r24957; float r24959 = c_n; float r24960 = c_p; float r24961 = log1p(r24954); float r24962 = log1p(r24948); float r24963 = r24961 - r24962; float r24964 = r24960 * r24963; float r24965 = fma(r24958, r24959, r24964); float r24966 = exp(r24965); return r24966; } double f_od(double c_p, double c_n, double t, double s) { double r24967 = 1; double r24968 = -r24967; double r24969 = s; double r24970 = -r24969; double r24971 = exp(r24970); double r24972 = r24967 + r24971; double r24973 = r24968 / r24972; double r24974 = log1p(r24973); double r24975 = t; double r24976 = -r24975; double r24977 = exp(r24976); double r24978 = r24977 + r24967; double r24979 = r24968 / r24978; double r24980 = log1p(r24979); double r24981 = r24974 - r24980; double r24982 = c_n; double r24983 = c_p; double r24984 = log1p(r24977); double r24985 = log1p(r24971); double r24986 = r24984 - r24985; double r24987 = r24983 * r24986; double r24988 = fma(r24981, r24982, r24987); double r24989 = exp(r24988); return r24989; } 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 r24990, r24991, r24992, r24993, r24994, r24995, r24996, r24997, r24998, r24999, r25000, r25001, r25002, r25003, r25004, r25005, r25006, r25007, r25008, r25009, r25010, r25011; void setup_mpfr_f_im() { mpfr_set_default_prec(2192); mpfr_init_set_str(r24990, "1", 10, MPFR_RNDN); mpfr_init(r24991); mpfr_init(r24992); mpfr_init(r24993); mpfr_init(r24994); mpfr_init(r24995); mpfr_init(r24996); mpfr_init(r24997); mpfr_init(r24998); mpfr_init(r24999); mpfr_init(r25000); mpfr_init(r25001); mpfr_init(r25002); mpfr_init(r25003); mpfr_init(r25004); mpfr_init(r25005); mpfr_init(r25006); mpfr_init(r25007); mpfr_init(r25008); mpfr_init(r25009); mpfr_init(r25010); mpfr_init(r25011); } double f_im(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r24991, s, MPFR_RNDN); mpfr_neg(r24992, r24991, MPFR_RNDN); mpfr_exp(r24993, r24992, MPFR_RNDN); mpfr_add(r24994, r24990, r24993, MPFR_RNDN); mpfr_div(r24995, r24990, r24994, MPFR_RNDN); mpfr_set_d(r24996, c_p, MPFR_RNDN); mpfr_pow(r24997, r24995, r24996, MPFR_RNDN); mpfr_sub(r24998, r24990, r24995, MPFR_RNDN); mpfr_set_d(r24999, c_n, MPFR_RNDN); mpfr_pow(r25000, r24998, r24999, MPFR_RNDN); mpfr_mul(r25001, r24997, r25000, MPFR_RNDN); mpfr_set_d(r25002, t, MPFR_RNDN); mpfr_neg(r25003, r25002, MPFR_RNDN); mpfr_exp(r25004, r25003, MPFR_RNDN); mpfr_add(r25005, r24990, r25004, MPFR_RNDN); mpfr_div(r25006, r24990, r25005, MPFR_RNDN); mpfr_pow(r25007, r25006, r24996, MPFR_RNDN); mpfr_sub(r25008, r24990, r25006, MPFR_RNDN); mpfr_pow(r25009, r25008, r24999, MPFR_RNDN); mpfr_mul(r25010, r25007, r25009, MPFR_RNDN); mpfr_div(r25011, r25001, r25010, MPFR_RNDN); return mpfr_get_d(r25011, MPFR_RNDN); } static mpfr_t r25012, r25013, r25014, r25015, r25016, r25017, r25018, r25019, r25020, r25021, r25022, r25023, r25024, r25025, r25026, r25027, r25028, r25029, r25030, r25031, r25032, r25033, r25034; void setup_mpfr_f_fm() { mpfr_set_default_prec(2192); mpfr_init_set_str(r25012, "1", 10, MPFR_RNDN); mpfr_init(r25013); mpfr_init(r25014); mpfr_init(r25015); mpfr_init(r25016); mpfr_init(r25017); mpfr_init(r25018); mpfr_init(r25019); mpfr_init(r25020); mpfr_init(r25021); mpfr_init(r25022); mpfr_init(r25023); mpfr_init(r25024); mpfr_init(r25025); mpfr_init(r25026); mpfr_init(r25027); mpfr_init(r25028); mpfr_init(r25029); mpfr_init(r25030); mpfr_init(r25031); mpfr_init(r25032); mpfr_init(r25033); mpfr_init(r25034); } double f_fm(double c_p, double c_n, double t, double s) { ; mpfr_neg(r25013, r25012, MPFR_RNDN); mpfr_set_d(r25014, s, MPFR_RNDN); mpfr_neg(r25015, r25014, MPFR_RNDN); mpfr_exp(r25016, r25015, MPFR_RNDN); mpfr_add(r25017, r25012, r25016, MPFR_RNDN); mpfr_div(r25018, r25013, r25017, MPFR_RNDN); mpfr_log1p(r25019, r25018, MPFR_RNDN); mpfr_set_d(r25020, t, MPFR_RNDN); mpfr_neg(r25021, r25020, MPFR_RNDN); mpfr_exp(r25022, r25021, MPFR_RNDN); mpfr_add(r25023, r25022, r25012, MPFR_RNDN); mpfr_div(r25024, r25013, r25023, MPFR_RNDN); mpfr_log1p(r25025, r25024, MPFR_RNDN); mpfr_sub(r25026, r25019, r25025, MPFR_RNDN); mpfr_set_d(r25027, c_n, MPFR_RNDN); mpfr_set_d(r25028, c_p, MPFR_RNDN); mpfr_log1p(r25029, r25022, MPFR_RNDN); mpfr_log1p(r25030, r25016, MPFR_RNDN); mpfr_sub(r25031, r25029, r25030, MPFR_RNDN); mpfr_mul(r25032, r25028, r25031, MPFR_RNDN); mpfr_fma(r25033, r25026, r25027, r25032, MPFR_RNDN); mpfr_exp(r25034, r25033, MPFR_RNDN); return mpfr_get_d(r25034, MPFR_RNDN); } static mpfr_t r25035, r25036, r25037, r25038, r25039, r25040, r25041, r25042, r25043, r25044, r25045, r25046, r25047, r25048, r25049, r25050, r25051, r25052, r25053, r25054, r25055, r25056, r25057; void setup_mpfr_f_dm() { mpfr_set_default_prec(2192); mpfr_init_set_str(r25035, "1", 10, MPFR_RNDN); mpfr_init(r25036); mpfr_init(r25037); mpfr_init(r25038); mpfr_init(r25039); mpfr_init(r25040); mpfr_init(r25041); mpfr_init(r25042); mpfr_init(r25043); mpfr_init(r25044); mpfr_init(r25045); mpfr_init(r25046); mpfr_init(r25047); mpfr_init(r25048); mpfr_init(r25049); mpfr_init(r25050); mpfr_init(r25051); mpfr_init(r25052); mpfr_init(r25053); mpfr_init(r25054); mpfr_init(r25055); mpfr_init(r25056); mpfr_init(r25057); } double f_dm(double c_p, double c_n, double t, double s) { ; mpfr_neg(r25036, r25035, MPFR_RNDN); mpfr_set_d(r25037, s, MPFR_RNDN); mpfr_neg(r25038, r25037, MPFR_RNDN); mpfr_exp(r25039, r25038, MPFR_RNDN); mpfr_add(r25040, r25035, r25039, MPFR_RNDN); mpfr_div(r25041, r25036, r25040, MPFR_RNDN); mpfr_log1p(r25042, r25041, MPFR_RNDN); mpfr_set_d(r25043, t, MPFR_RNDN); mpfr_neg(r25044, r25043, MPFR_RNDN); mpfr_exp(r25045, r25044, MPFR_RNDN); mpfr_add(r25046, r25045, r25035, MPFR_RNDN); mpfr_div(r25047, r25036, r25046, MPFR_RNDN); mpfr_log1p(r25048, r25047, MPFR_RNDN); mpfr_sub(r25049, r25042, r25048, MPFR_RNDN); mpfr_set_d(r25050, c_n, MPFR_RNDN); mpfr_set_d(r25051, c_p, MPFR_RNDN); mpfr_log1p(r25052, r25045, MPFR_RNDN); mpfr_log1p(r25053, r25039, MPFR_RNDN); mpfr_sub(r25054, r25052, r25053, MPFR_RNDN); mpfr_mul(r25055, r25051, r25054, MPFR_RNDN); mpfr_fma(r25056, r25049, r25050, r25055, MPFR_RNDN); mpfr_exp(r25057, r25056, MPFR_RNDN); return mpfr_get_d(r25057, MPFR_RNDN); }