#include #include #include #include #include char *name = "(+ x (- (tan (+ y z)) (tan a)))"; double f_if(float x, float y, float z, float a) { float r24939 = x; float r24940 = y; float r24941 = z; float r24942 = r24940 + r24941; float r24943 = tan(r24942); float r24944 = a; float r24945 = tan(r24944); float r24946 = r24943 - r24945; float r24947 = r24939 + r24946; return r24947; } double f_id(double x, double y, double z, double a) { double r24948 = x; double r24949 = y; double r24950 = z; double r24951 = r24949 + r24950; double r24952 = tan(r24951); double r24953 = a; double r24954 = tan(r24953); double r24955 = r24952 - r24954; double r24956 = r24948 + r24955; return r24956; } double f_of(float x, float y, float z, float a) { float r24957 = x; float r24958 = y; float r24959 = tan(r24958); float r24960 = z; float r24961 = tan(r24960); float r24962 = r24959 + r24961; float r24963 = r24959 * r24961; float r24964 = fma(r24963, r24963, r24963); float r24965 = fma(r24962, r24964, r24962); float r24966 = 1; float r24967 = 3; float r24968 = pow(r24963, r24967); float r24969 = r24966 - r24968; float r24970 = r24965 / r24969; float r24971 = a; float r24972 = tan(r24971); float r24973 = r24970 - r24972; float r24974 = 0; float r24975 = r24973 + r24974; float r24976 = r24957 + r24975; return r24976; } double f_od(double x, double y, double z, double a) { double r24977 = x; double r24978 = y; double r24979 = tan(r24978); double r24980 = z; double r24981 = tan(r24980); double r24982 = r24979 + r24981; double r24983 = r24979 * r24981; double r24984 = fma(r24983, r24983, r24983); double r24985 = fma(r24982, r24984, r24982); double r24986 = 1; double r24987 = 3; double r24988 = pow(r24983, r24987); double r24989 = r24986 - r24988; double r24990 = r24985 / r24989; double r24991 = a; double r24992 = tan(r24991); double r24993 = r24990 - r24992; double r24994 = 0; double r24995 = r24993 + r24994; double r24996 = r24977 + r24995; return r24996; } 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 r24997, r24998, r24999, r25000, r25001, r25002, r25003, r25004, r25005; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); 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); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r24997, x, MPFR_RNDN); mpfr_set_d(r24998, y, MPFR_RNDN); mpfr_set_d(r24999, z, MPFR_RNDN); mpfr_add(r25000, r24998, r24999, MPFR_RNDN); mpfr_tan(r25001, r25000, MPFR_RNDN); mpfr_set_d(r25002, a, MPFR_RNDN); mpfr_tan(r25003, r25002, MPFR_RNDN); mpfr_sub(r25004, r25001, r25003, MPFR_RNDN); mpfr_add(r25005, r24997, r25004, MPFR_RNDN); return mpfr_get_d(r25005, MPFR_RNDN); } static mpfr_t r25006, r25007, r25008, r25009, r25010, r25011, r25012, r25013, r25014, r25015, r25016, r25017, r25018, r25019, r25020, r25021, r25022, r25023, r25024, r25025; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25006); mpfr_init(r25007); mpfr_init(r25008); mpfr_init(r25009); mpfr_init(r25010); mpfr_init(r25011); mpfr_init(r25012); mpfr_init(r25013); mpfr_init(r25014); mpfr_init_set_str(r25015, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25016, "3", 10, MPFR_RNDN); mpfr_init(r25017); mpfr_init(r25018); mpfr_init(r25019); mpfr_init(r25020); mpfr_init(r25021); mpfr_init(r25022); mpfr_init_set_str(r25023, "0", 10, MPFR_RNDN); mpfr_init(r25024); mpfr_init(r25025); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r25006, x, MPFR_RNDN); mpfr_set_d(r25007, y, MPFR_RNDN); mpfr_tan(r25008, r25007, MPFR_RNDN); mpfr_set_d(r25009, z, MPFR_RNDN); mpfr_tan(r25010, r25009, MPFR_RNDN); mpfr_add(r25011, r25008, r25010, MPFR_RNDN); mpfr_mul(r25012, r25008, r25010, MPFR_RNDN); mpfr_fma(r25013, r25012, r25012, r25012, MPFR_RNDN); mpfr_fma(r25014, r25011, r25013, r25011, MPFR_RNDN); ; ; mpfr_pow(r25017, r25012, r25016, MPFR_RNDN); mpfr_sub(r25018, r25015, r25017, MPFR_RNDN); mpfr_div(r25019, r25014, r25018, MPFR_RNDN); mpfr_set_d(r25020, a, MPFR_RNDN); mpfr_tan(r25021, r25020, MPFR_RNDN); mpfr_sub(r25022, r25019, r25021, MPFR_RNDN); ; mpfr_add(r25024, r25022, r25023, MPFR_RNDN); mpfr_add(r25025, r25006, r25024, MPFR_RNDN); return mpfr_get_d(r25025, MPFR_RNDN); } static mpfr_t r25026, r25027, r25028, r25029, r25030, r25031, r25032, r25033, r25034, r25035, r25036, r25037, r25038, r25039, r25040, r25041, r25042, r25043, r25044, r25045; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); 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); mpfr_init_set_str(r25035, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25036, "3", 10, MPFR_RNDN); mpfr_init(r25037); mpfr_init(r25038); mpfr_init(r25039); mpfr_init(r25040); mpfr_init(r25041); mpfr_init(r25042); mpfr_init_set_str(r25043, "0", 10, MPFR_RNDN); mpfr_init(r25044); mpfr_init(r25045); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r25026, x, MPFR_RNDN); mpfr_set_d(r25027, y, MPFR_RNDN); mpfr_tan(r25028, r25027, MPFR_RNDN); mpfr_set_d(r25029, z, MPFR_RNDN); mpfr_tan(r25030, r25029, MPFR_RNDN); mpfr_add(r25031, r25028, r25030, MPFR_RNDN); mpfr_mul(r25032, r25028, r25030, MPFR_RNDN); mpfr_fma(r25033, r25032, r25032, r25032, MPFR_RNDN); mpfr_fma(r25034, r25031, r25033, r25031, MPFR_RNDN); ; ; mpfr_pow(r25037, r25032, r25036, MPFR_RNDN); mpfr_sub(r25038, r25035, r25037, MPFR_RNDN); mpfr_div(r25039, r25034, r25038, MPFR_RNDN); mpfr_set_d(r25040, a, MPFR_RNDN); mpfr_tan(r25041, r25040, MPFR_RNDN); mpfr_sub(r25042, r25039, r25041, MPFR_RNDN); ; mpfr_add(r25044, r25042, r25043, MPFR_RNDN); mpfr_add(r25045, r25026, r25044, MPFR_RNDN); return mpfr_get_d(r25045, MPFR_RNDN); }