#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 r27012 = x; float r27013 = y; float r27014 = z; float r27015 = r27013 + r27014; float r27016 = tan(r27015); float r27017 = a; float r27018 = tan(r27017); float r27019 = r27016 - r27018; float r27020 = r27012 + r27019; return r27020; } double f_id(double x, double y, double z, double a) { double r27021 = x; double r27022 = y; double r27023 = z; double r27024 = r27022 + r27023; double r27025 = tan(r27024); double r27026 = a; double r27027 = tan(r27026); double r27028 = r27025 - r27027; double r27029 = r27021 + r27028; return r27029; } double f_of(float x, float y, float z, float a) { float r27030 = x; float r27031 = y; float r27032 = tan(r27031); float r27033 = z; float r27034 = tan(r27033); float r27035 = r27032 + r27034; float r27036 = a; float r27037 = cos(r27036); float r27038 = r27035 * r27037; float r27039 = 1; float r27040 = r27032 * r27034; float r27041 = r27039 - r27040; float r27042 = sin(r27036); float r27043 = r27041 * r27042; float r27044 = r27038 - r27043; float r27045 = r27041 * r27037; float r27046 = r27044 / r27045; float r27047 = r27030 + r27046; return r27047; } double f_od(double x, double y, double z, double a) { double r27048 = x; double r27049 = y; double r27050 = tan(r27049); double r27051 = z; double r27052 = tan(r27051); double r27053 = r27050 + r27052; double r27054 = a; double r27055 = cos(r27054); double r27056 = r27053 * r27055; double r27057 = 1; double r27058 = r27050 * r27052; double r27059 = r27057 - r27058; double r27060 = sin(r27054); double r27061 = r27059 * r27060; double r27062 = r27056 - r27061; double r27063 = r27059 * r27055; double r27064 = r27062 / r27063; double r27065 = r27048 + r27064; return r27065; } 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 r27066, r27067, r27068, r27069, r27070, r27071, r27072, r27073, r27074; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r27066); mpfr_init(r27067); mpfr_init(r27068); mpfr_init(r27069); mpfr_init(r27070); mpfr_init(r27071); mpfr_init(r27072); mpfr_init(r27073); mpfr_init(r27074); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r27066, x, MPFR_RNDN); mpfr_set_d(r27067, y, MPFR_RNDN); mpfr_set_d(r27068, z, MPFR_RNDN); mpfr_add(r27069, r27067, r27068, MPFR_RNDN); mpfr_tan(r27070, r27069, MPFR_RNDN); mpfr_set_d(r27071, a, MPFR_RNDN); mpfr_tan(r27072, r27071, MPFR_RNDN); mpfr_sub(r27073, r27070, r27072, MPFR_RNDN); mpfr_add(r27074, r27066, r27073, MPFR_RNDN); return mpfr_get_d(r27074, MPFR_RNDN); } static mpfr_t r27075, r27076, r27077, r27078, r27079, r27080, r27081, r27082, r27083, r27084, r27085, r27086, r27087, r27088, r27089, r27090, r27091, r27092; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r27075); mpfr_init(r27076); mpfr_init(r27077); mpfr_init(r27078); mpfr_init(r27079); mpfr_init(r27080); mpfr_init(r27081); mpfr_init(r27082); mpfr_init(r27083); mpfr_init_set_str(r27084, "1", 10, MPFR_RNDN); mpfr_init(r27085); mpfr_init(r27086); mpfr_init(r27087); mpfr_init(r27088); mpfr_init(r27089); mpfr_init(r27090); mpfr_init(r27091); mpfr_init(r27092); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r27075, x, MPFR_RNDN); mpfr_set_d(r27076, y, MPFR_RNDN); mpfr_tan(r27077, r27076, MPFR_RNDN); mpfr_set_d(r27078, z, MPFR_RNDN); mpfr_tan(r27079, r27078, MPFR_RNDN); mpfr_add(r27080, r27077, r27079, MPFR_RNDN); mpfr_set_d(r27081, a, MPFR_RNDN); mpfr_cos(r27082, r27081, MPFR_RNDN); mpfr_mul(r27083, r27080, r27082, MPFR_RNDN); ; mpfr_mul(r27085, r27077, r27079, MPFR_RNDN); mpfr_sub(r27086, r27084, r27085, MPFR_RNDN); mpfr_sin(r27087, r27081, MPFR_RNDN); mpfr_mul(r27088, r27086, r27087, MPFR_RNDN); mpfr_sub(r27089, r27083, r27088, MPFR_RNDN); mpfr_mul(r27090, r27086, r27082, MPFR_RNDN); mpfr_div(r27091, r27089, r27090, MPFR_RNDN); mpfr_add(r27092, r27075, r27091, MPFR_RNDN); return mpfr_get_d(r27092, MPFR_RNDN); } static mpfr_t r27093, r27094, r27095, r27096, r27097, r27098, r27099, r27100, r27101, r27102, r27103, r27104, r27105, r27106, r27107, r27108, r27109, r27110; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r27093); mpfr_init(r27094); mpfr_init(r27095); mpfr_init(r27096); mpfr_init(r27097); mpfr_init(r27098); mpfr_init(r27099); mpfr_init(r27100); mpfr_init(r27101); mpfr_init_set_str(r27102, "1", 10, MPFR_RNDN); mpfr_init(r27103); mpfr_init(r27104); mpfr_init(r27105); mpfr_init(r27106); mpfr_init(r27107); mpfr_init(r27108); mpfr_init(r27109); mpfr_init(r27110); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r27093, x, MPFR_RNDN); mpfr_set_d(r27094, y, MPFR_RNDN); mpfr_tan(r27095, r27094, MPFR_RNDN); mpfr_set_d(r27096, z, MPFR_RNDN); mpfr_tan(r27097, r27096, MPFR_RNDN); mpfr_add(r27098, r27095, r27097, MPFR_RNDN); mpfr_set_d(r27099, a, MPFR_RNDN); mpfr_cos(r27100, r27099, MPFR_RNDN); mpfr_mul(r27101, r27098, r27100, MPFR_RNDN); ; mpfr_mul(r27103, r27095, r27097, MPFR_RNDN); mpfr_sub(r27104, r27102, r27103, MPFR_RNDN); mpfr_sin(r27105, r27099, MPFR_RNDN); mpfr_mul(r27106, r27104, r27105, MPFR_RNDN); mpfr_sub(r27107, r27101, r27106, MPFR_RNDN); mpfr_mul(r27108, r27104, r27100, MPFR_RNDN); mpfr_div(r27109, r27107, r27108, MPFR_RNDN); mpfr_add(r27110, r27093, r27109, MPFR_RNDN); return mpfr_get_d(r27110, MPFR_RNDN); }