#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r24957 = 0.5; float r24958 = re; float r24959 = sin(r24958); float r24960 = r24957 * r24959; float r24961 = im; float r24962 = -r24961; float r24963 = exp(r24962); float r24964 = exp(r24961); float r24965 = r24963 - r24964; float r24966 = r24960 * r24965; return r24966; } double f_id(double re, double im) { double r24967 = 0.5; double r24968 = re; double r24969 = sin(r24968); double r24970 = r24967 * r24969; double r24971 = im; double r24972 = -r24971; double r24973 = exp(r24972); double r24974 = exp(r24971); double r24975 = r24973 - r24974; double r24976 = r24970 * r24975; return r24976; } double f_of(float re, float im) { float r24977 = re; float r24978 = sin(r24977); float r24979 = im; float r24980 = 0.5; float r24981 = r24979 * r24980; float r24982 = r24978 * r24981; float r24983 = 1/3; float r24984 = r24983 * r24979; float r24985 = 2; float r24986 = fma(r24984, r24979, r24985); float r24987 = -r24986; float r24988 = r24982 * r24987; float r24989 = -r24980; float r24990 = r24989 * r24978; float r24991 = 5; float r24992 = pow(r24979, r24991); float r24993 = 1/60; float r24994 = r24992 * r24993; float r24995 = r24990 * r24994; float r24996 = r24988 + r24995; return r24996; } double f_od(double re, double im) { double r24997 = re; double r24998 = sin(r24997); double r24999 = im; double r25000 = 0.5; double r25001 = r24999 * r25000; double r25002 = r24998 * r25001; double r25003 = 1/3; double r25004 = r25003 * r24999; double r25005 = 2; double r25006 = fma(r25004, r24999, r25005); double r25007 = -r25006; double r25008 = r25002 * r25007; double r25009 = -r25000; double r25010 = r25009 * r24998; double r25011 = 5; double r25012 = pow(r24999, r25011); double r25013 = 1/60; double r25014 = r25012 * r25013; double r25015 = r25010 * r25014; double r25016 = r25008 + r25015; return r25016; } 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 r25017, r25018, r25019, r25020, r25021, r25022, r25023, r25024, r25025, r25026; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r25017, "0.5", 10, MPFR_RNDN); 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); } double f_im(double re, double im) { ; mpfr_set_d(r25018, re, MPFR_RNDN); mpfr_sin(r25019, r25018, MPFR_RNDN); mpfr_mul(r25020, r25017, r25019, MPFR_RNDN); mpfr_set_d(r25021, im, MPFR_RNDN); mpfr_neg(r25022, r25021, MPFR_RNDN); mpfr_exp(r25023, r25022, MPFR_RNDN); mpfr_exp(r25024, r25021, MPFR_RNDN); mpfr_sub(r25025, r25023, r25024, MPFR_RNDN); mpfr_mul(r25026, r25020, r25025, MPFR_RNDN); return mpfr_get_d(r25026, MPFR_RNDN); } static mpfr_t r25027, r25028, r25029, r25030, r25031, r25032, r25033, r25034, r25035, r25036, r25037, r25038, r25039, r25040, r25041, r25042, r25043, r25044, r25045, r25046; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r25027); mpfr_init(r25028); mpfr_init(r25029); mpfr_init_set_str(r25030, "0.5", 10, MPFR_RNDN); mpfr_init(r25031); mpfr_init(r25032); mpfr_init_set_str(r25033, "1/3", 10, MPFR_RNDN); mpfr_init(r25034); mpfr_init_set_str(r25035, "2", 10, MPFR_RNDN); mpfr_init(r25036); mpfr_init(r25037); mpfr_init(r25038); mpfr_init(r25039); mpfr_init(r25040); mpfr_init_set_str(r25041, "5", 10, MPFR_RNDN); mpfr_init(r25042); mpfr_init_set_str(r25043, "1/60", 10, MPFR_RNDN); mpfr_init(r25044); mpfr_init(r25045); mpfr_init(r25046); } double f_fm(double re, double im) { mpfr_set_d(r25027, re, MPFR_RNDN); mpfr_sin(r25028, r25027, MPFR_RNDN); mpfr_set_d(r25029, im, MPFR_RNDN); ; mpfr_mul(r25031, r25029, r25030, MPFR_RNDN); mpfr_mul(r25032, r25028, r25031, MPFR_RNDN); ; mpfr_mul(r25034, r25033, r25029, MPFR_RNDN); ; mpfr_fma(r25036, r25034, r25029, r25035, MPFR_RNDN); mpfr_neg(r25037, r25036, MPFR_RNDN); mpfr_mul(r25038, r25032, r25037, MPFR_RNDN); mpfr_neg(r25039, r25030, MPFR_RNDN); mpfr_mul(r25040, r25039, r25028, MPFR_RNDN); ; mpfr_pow(r25042, r25029, r25041, MPFR_RNDN); ; mpfr_mul(r25044, r25042, r25043, MPFR_RNDN); mpfr_mul(r25045, r25040, r25044, MPFR_RNDN); mpfr_add(r25046, r25038, r25045, MPFR_RNDN); return mpfr_get_d(r25046, MPFR_RNDN); } static mpfr_t r25047, r25048, r25049, r25050, r25051, r25052, r25053, r25054, r25055, r25056, r25057, r25058, r25059, r25060, r25061, r25062, r25063, r25064, r25065, r25066; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r25047); mpfr_init(r25048); mpfr_init(r25049); mpfr_init_set_str(r25050, "0.5", 10, MPFR_RNDN); mpfr_init(r25051); mpfr_init(r25052); mpfr_init_set_str(r25053, "1/3", 10, MPFR_RNDN); mpfr_init(r25054); mpfr_init_set_str(r25055, "2", 10, MPFR_RNDN); mpfr_init(r25056); mpfr_init(r25057); mpfr_init(r25058); mpfr_init(r25059); mpfr_init(r25060); mpfr_init_set_str(r25061, "5", 10, MPFR_RNDN); mpfr_init(r25062); mpfr_init_set_str(r25063, "1/60", 10, MPFR_RNDN); mpfr_init(r25064); mpfr_init(r25065); mpfr_init(r25066); } double f_dm(double re, double im) { mpfr_set_d(r25047, re, MPFR_RNDN); mpfr_sin(r25048, r25047, MPFR_RNDN); mpfr_set_d(r25049, im, MPFR_RNDN); ; mpfr_mul(r25051, r25049, r25050, MPFR_RNDN); mpfr_mul(r25052, r25048, r25051, MPFR_RNDN); ; mpfr_mul(r25054, r25053, r25049, MPFR_RNDN); ; mpfr_fma(r25056, r25054, r25049, r25055, MPFR_RNDN); mpfr_neg(r25057, r25056, MPFR_RNDN); mpfr_mul(r25058, r25052, r25057, MPFR_RNDN); mpfr_neg(r25059, r25050, MPFR_RNDN); mpfr_mul(r25060, r25059, r25048, MPFR_RNDN); ; mpfr_pow(r25062, r25049, r25061, MPFR_RNDN); ; mpfr_mul(r25064, r25062, r25063, MPFR_RNDN); mpfr_mul(r25065, r25060, r25064, MPFR_RNDN); mpfr_add(r25066, r25058, r25065, MPFR_RNDN); return mpfr_get_d(r25066, MPFR_RNDN); }