#include #include #include #include #include char *name = "Harley's example"; double f_if(float c_p, float c_n, float t, float s) { float r21949 = 1; float r21950 = s; float r21951 = -r21950; float r21952 = exp(r21951); float r21953 = r21949 + r21952; float r21954 = r21949 / r21953; float r21955 = c_p; float r21956 = pow(r21954, r21955); float r21957 = r21949 - r21954; float r21958 = c_n; float r21959 = pow(r21957, r21958); float r21960 = r21956 * r21959; float r21961 = t; float r21962 = -r21961; float r21963 = exp(r21962); float r21964 = r21949 + r21963; float r21965 = r21949 / r21964; float r21966 = pow(r21965, r21955); float r21967 = r21949 - r21965; float r21968 = pow(r21967, r21958); float r21969 = r21966 * r21968; float r21970 = r21960 / r21969; return r21970; } double f_id(double c_p, double c_n, double t, double s) { double r21971 = 1; double r21972 = s; double r21973 = -r21972; double r21974 = exp(r21973); double r21975 = r21971 + r21974; double r21976 = r21971 / r21975; double r21977 = c_p; double r21978 = pow(r21976, r21977); double r21979 = r21971 - r21976; double r21980 = c_n; double r21981 = pow(r21979, r21980); double r21982 = r21978 * r21981; double r21983 = t; double r21984 = -r21983; double r21985 = exp(r21984); double r21986 = r21971 + r21985; double r21987 = r21971 / r21986; double r21988 = pow(r21987, r21977); double r21989 = r21971 - r21987; double r21990 = pow(r21989, r21980); double r21991 = r21988 * r21990; double r21992 = r21982 / r21991; return r21992; } double f_of(float c_p, float c_n, float t, float s) { float r21993 = 1; float r21994 = s; float r21995 = -r21994; float r21996 = exp(r21995); float r21997 = r21993 + r21996; float r21998 = log(r21997); float r21999 = t; float r22000 = -r21999; float r22001 = exp(r22000); float r22002 = r22001 + r21993; float r22003 = log(r22002); float r22004 = r21998 - r22003; float r22005 = c_p; float r22006 = -r22005; float r22007 = r22004 * r22006; float r22008 = r21993 / r21997; float r22009 = r21993 - r22008; float r22010 = r21993 / r22002; float r22011 = r21993 - r22010; float r22012 = r22009 / r22011; float r22013 = log(r22012); float r22014 = 3; float r22015 = pow(r22013, r22014); float r22016 = cbrt(r22015); float r22017 = c_n; float r22018 = r22016 * r22017; float r22019 = r22007 + r22018; float r22020 = exp(r22019); return r22020; } double f_od(double c_p, double c_n, double t, double s) { double r22021 = 1; double r22022 = s; double r22023 = -r22022; double r22024 = exp(r22023); double r22025 = r22021 + r22024; double r22026 = log(r22025); double r22027 = t; double r22028 = -r22027; double r22029 = exp(r22028); double r22030 = r22029 + r22021; double r22031 = log(r22030); double r22032 = r22026 - r22031; double r22033 = c_p; double r22034 = -r22033; double r22035 = r22032 * r22034; double r22036 = r22021 / r22025; double r22037 = r22021 - r22036; double r22038 = r22021 / r22030; double r22039 = r22021 - r22038; double r22040 = r22037 / r22039; double r22041 = log(r22040); double r22042 = 3; double r22043 = pow(r22041, r22042); double r22044 = cbrt(r22043); double r22045 = c_n; double r22046 = r22044 * r22045; double r22047 = r22035 + r22046; double r22048 = exp(r22047); return r22048; } 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 r22049, r22050, r22051, r22052, r22053, r22054, r22055, r22056, r22057, r22058, r22059, r22060, r22061, r22062, r22063, r22064, r22065, r22066, r22067, r22068, r22069, r22070; void setup_mpfr_f_im() { mpfr_set_default_prec(1936); mpfr_init_set_str(r22049, "1", 10, MPFR_RNDN); mpfr_init(r22050); mpfr_init(r22051); mpfr_init(r22052); mpfr_init(r22053); mpfr_init(r22054); mpfr_init(r22055); mpfr_init(r22056); mpfr_init(r22057); mpfr_init(r22058); mpfr_init(r22059); mpfr_init(r22060); mpfr_init(r22061); mpfr_init(r22062); mpfr_init(r22063); mpfr_init(r22064); mpfr_init(r22065); mpfr_init(r22066); mpfr_init(r22067); mpfr_init(r22068); mpfr_init(r22069); mpfr_init(r22070); } double f_im(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r22050, s, MPFR_RNDN); mpfr_neg(r22051, r22050, MPFR_RNDN); mpfr_exp(r22052, r22051, MPFR_RNDN); mpfr_add(r22053, r22049, r22052, MPFR_RNDN); mpfr_div(r22054, r22049, r22053, MPFR_RNDN); mpfr_set_d(r22055, c_p, MPFR_RNDN); mpfr_pow(r22056, r22054, r22055, MPFR_RNDN); mpfr_sub(r22057, r22049, r22054, MPFR_RNDN); mpfr_set_d(r22058, c_n, MPFR_RNDN); mpfr_pow(r22059, r22057, r22058, MPFR_RNDN); mpfr_mul(r22060, r22056, r22059, MPFR_RNDN); mpfr_set_d(r22061, t, MPFR_RNDN); mpfr_neg(r22062, r22061, MPFR_RNDN); mpfr_exp(r22063, r22062, MPFR_RNDN); mpfr_add(r22064, r22049, r22063, MPFR_RNDN); mpfr_div(r22065, r22049, r22064, MPFR_RNDN); mpfr_pow(r22066, r22065, r22055, MPFR_RNDN); mpfr_sub(r22067, r22049, r22065, MPFR_RNDN); mpfr_pow(r22068, r22067, r22058, MPFR_RNDN); mpfr_mul(r22069, r22066, r22068, MPFR_RNDN); mpfr_div(r22070, r22060, r22069, MPFR_RNDN); return mpfr_get_d(r22070, MPFR_RNDN); } static mpfr_t r22071, r22072, r22073, r22074, r22075, r22076, r22077, r22078, r22079, r22080, r22081, r22082, r22083, r22084, r22085, r22086, r22087, r22088, r22089, r22090, r22091, r22092, r22093, r22094, r22095, r22096, r22097, r22098; void setup_mpfr_f_fm() { mpfr_set_default_prec(1936); mpfr_init_set_str(r22071, "1", 10, MPFR_RNDN); mpfr_init(r22072); mpfr_init(r22073); mpfr_init(r22074); mpfr_init(r22075); mpfr_init(r22076); mpfr_init(r22077); mpfr_init(r22078); mpfr_init(r22079); mpfr_init(r22080); mpfr_init(r22081); mpfr_init(r22082); mpfr_init(r22083); mpfr_init(r22084); mpfr_init(r22085); mpfr_init(r22086); mpfr_init(r22087); mpfr_init(r22088); mpfr_init(r22089); mpfr_init(r22090); mpfr_init(r22091); mpfr_init_set_str(r22092, "3", 10, MPFR_RNDN); mpfr_init(r22093); mpfr_init(r22094); mpfr_init(r22095); mpfr_init(r22096); mpfr_init(r22097); mpfr_init(r22098); } double f_fm(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r22072, s, MPFR_RNDN); mpfr_neg(r22073, r22072, MPFR_RNDN); mpfr_exp(r22074, r22073, MPFR_RNDN); mpfr_add(r22075, r22071, r22074, MPFR_RNDN); mpfr_log(r22076, r22075, MPFR_RNDN); mpfr_set_d(r22077, t, MPFR_RNDN); mpfr_neg(r22078, r22077, MPFR_RNDN); mpfr_exp(r22079, r22078, MPFR_RNDN); mpfr_add(r22080, r22079, r22071, MPFR_RNDN); mpfr_log(r22081, r22080, MPFR_RNDN); mpfr_sub(r22082, r22076, r22081, MPFR_RNDN); mpfr_set_d(r22083, c_p, MPFR_RNDN); mpfr_neg(r22084, r22083, MPFR_RNDN); mpfr_mul(r22085, r22082, r22084, MPFR_RNDN); mpfr_div(r22086, r22071, r22075, MPFR_RNDN); mpfr_sub(r22087, r22071, r22086, MPFR_RNDN); mpfr_div(r22088, r22071, r22080, MPFR_RNDN); mpfr_sub(r22089, r22071, r22088, MPFR_RNDN); mpfr_div(r22090, r22087, r22089, MPFR_RNDN); mpfr_log(r22091, r22090, MPFR_RNDN); ; mpfr_pow(r22093, r22091, r22092, MPFR_RNDN); mpfr_cbrt(r22094, r22093, MPFR_RNDN); mpfr_set_d(r22095, c_n, MPFR_RNDN); mpfr_mul(r22096, r22094, r22095, MPFR_RNDN); mpfr_add(r22097, r22085, r22096, MPFR_RNDN); mpfr_exp(r22098, r22097, MPFR_RNDN); return mpfr_get_d(r22098, MPFR_RNDN); } static mpfr_t r22099, r22100, r22101, r22102, r22103, r22104, r22105, r22106, r22107, r22108, r22109, r22110, r22111, r22112, r22113, r22114, r22115, r22116, r22117, r22118, r22119, r22120, r22121, r22122, r22123, r22124, r22125, r22126; void setup_mpfr_f_dm() { mpfr_set_default_prec(1936); mpfr_init_set_str(r22099, "1", 10, MPFR_RNDN); mpfr_init(r22100); mpfr_init(r22101); mpfr_init(r22102); mpfr_init(r22103); mpfr_init(r22104); mpfr_init(r22105); mpfr_init(r22106); mpfr_init(r22107); mpfr_init(r22108); mpfr_init(r22109); mpfr_init(r22110); mpfr_init(r22111); mpfr_init(r22112); mpfr_init(r22113); mpfr_init(r22114); mpfr_init(r22115); mpfr_init(r22116); mpfr_init(r22117); mpfr_init(r22118); mpfr_init(r22119); mpfr_init_set_str(r22120, "3", 10, MPFR_RNDN); mpfr_init(r22121); mpfr_init(r22122); mpfr_init(r22123); mpfr_init(r22124); mpfr_init(r22125); mpfr_init(r22126); } double f_dm(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r22100, s, MPFR_RNDN); mpfr_neg(r22101, r22100, MPFR_RNDN); mpfr_exp(r22102, r22101, MPFR_RNDN); mpfr_add(r22103, r22099, r22102, MPFR_RNDN); mpfr_log(r22104, r22103, MPFR_RNDN); mpfr_set_d(r22105, t, MPFR_RNDN); mpfr_neg(r22106, r22105, MPFR_RNDN); mpfr_exp(r22107, r22106, MPFR_RNDN); mpfr_add(r22108, r22107, r22099, MPFR_RNDN); mpfr_log(r22109, r22108, MPFR_RNDN); mpfr_sub(r22110, r22104, r22109, MPFR_RNDN); mpfr_set_d(r22111, c_p, MPFR_RNDN); mpfr_neg(r22112, r22111, MPFR_RNDN); mpfr_mul(r22113, r22110, r22112, MPFR_RNDN); mpfr_div(r22114, r22099, r22103, MPFR_RNDN); mpfr_sub(r22115, r22099, r22114, MPFR_RNDN); mpfr_div(r22116, r22099, r22108, MPFR_RNDN); mpfr_sub(r22117, r22099, r22116, MPFR_RNDN); mpfr_div(r22118, r22115, r22117, MPFR_RNDN); mpfr_log(r22119, r22118, MPFR_RNDN); ; mpfr_pow(r22121, r22119, r22120, MPFR_RNDN); mpfr_cbrt(r22122, r22121, MPFR_RNDN); mpfr_set_d(r22123, c_n, MPFR_RNDN); mpfr_mul(r22124, r22122, r22123, MPFR_RNDN); mpfr_add(r22125, r22113, r22124, MPFR_RNDN); mpfr_exp(r22126, r22125, MPFR_RNDN); return mpfr_get_d(r22126, MPFR_RNDN); }