#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 r25268 = x; float r25269 = y; float r25270 = z; float r25271 = r25269 + r25270; float r25272 = tan(r25271); float r25273 = a; float r25274 = tan(r25273); float r25275 = r25272 - r25274; float r25276 = r25268 + r25275; return r25276; } double f_id(double x, double y, double z, double a) { double r25277 = x; double r25278 = y; double r25279 = z; double r25280 = r25278 + r25279; double r25281 = tan(r25280); double r25282 = a; double r25283 = tan(r25282); double r25284 = r25281 - r25283; double r25285 = r25277 + r25284; return r25285; } double f_of(float x, float y, float z, float a) { float r25286 = x; float r25287 = y; float r25288 = tan(r25287); float r25289 = z; float r25290 = tan(r25289); float r25291 = r25288 + r25290; float r25292 = a; float r25293 = cos(r25292); float r25294 = r25291 * r25293; float r25295 = 1; float r25296 = r25288 * r25290; float r25297 = r25295 - r25296; float r25298 = sin(r25292); float r25299 = r25297 * r25298; float r25300 = r25294 - r25299; float r25301 = r25297 * r25293; float r25302 = r25300 / r25301; float r25303 = r25286 + r25302; return r25303; } double f_od(double x, double y, double z, double a) { double r25304 = x; double r25305 = y; double r25306 = tan(r25305); double r25307 = z; double r25308 = tan(r25307); double r25309 = r25306 + r25308; double r25310 = a; double r25311 = cos(r25310); double r25312 = r25309 * r25311; double r25313 = 1; double r25314 = r25306 * r25308; double r25315 = r25313 - r25314; double r25316 = sin(r25310); double r25317 = r25315 * r25316; double r25318 = r25312 - r25317; double r25319 = r25315 * r25311; double r25320 = r25318 / r25319; double r25321 = r25304 + r25320; return r25321; } 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 r25322, r25323, r25324, r25325, r25326, r25327, r25328, r25329, r25330; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r25322); mpfr_init(r25323); mpfr_init(r25324); mpfr_init(r25325); mpfr_init(r25326); mpfr_init(r25327); mpfr_init(r25328); mpfr_init(r25329); mpfr_init(r25330); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r25322, x, MPFR_RNDN); mpfr_set_d(r25323, y, MPFR_RNDN); mpfr_set_d(r25324, z, MPFR_RNDN); mpfr_add(r25325, r25323, r25324, MPFR_RNDN); mpfr_tan(r25326, r25325, MPFR_RNDN); mpfr_set_d(r25327, a, MPFR_RNDN); mpfr_tan(r25328, r25327, MPFR_RNDN); mpfr_sub(r25329, r25326, r25328, MPFR_RNDN); mpfr_add(r25330, r25322, r25329, MPFR_RNDN); return mpfr_get_d(r25330, MPFR_RNDN); } static mpfr_t r25331, r25332, r25333, r25334, r25335, r25336, r25337, r25338, r25339, r25340, r25341, r25342, r25343, r25344, r25345, r25346, r25347, r25348; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r25331); mpfr_init(r25332); mpfr_init(r25333); mpfr_init(r25334); mpfr_init(r25335); mpfr_init(r25336); mpfr_init(r25337); mpfr_init(r25338); mpfr_init(r25339); mpfr_init_set_str(r25340, "1", 10, MPFR_RNDN); mpfr_init(r25341); mpfr_init(r25342); mpfr_init(r25343); mpfr_init(r25344); mpfr_init(r25345); mpfr_init(r25346); mpfr_init(r25347); mpfr_init(r25348); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r25331, x, MPFR_RNDN); mpfr_set_d(r25332, y, MPFR_RNDN); mpfr_tan(r25333, r25332, MPFR_RNDN); mpfr_set_d(r25334, z, MPFR_RNDN); mpfr_tan(r25335, r25334, MPFR_RNDN); mpfr_add(r25336, r25333, r25335, MPFR_RNDN); mpfr_set_d(r25337, a, MPFR_RNDN); mpfr_cos(r25338, r25337, MPFR_RNDN); mpfr_mul(r25339, r25336, r25338, MPFR_RNDN); ; mpfr_mul(r25341, r25333, r25335, MPFR_RNDN); mpfr_sub(r25342, r25340, r25341, MPFR_RNDN); mpfr_sin(r25343, r25337, MPFR_RNDN); mpfr_mul(r25344, r25342, r25343, MPFR_RNDN); mpfr_sub(r25345, r25339, r25344, MPFR_RNDN); mpfr_mul(r25346, r25342, r25338, MPFR_RNDN); mpfr_div(r25347, r25345, r25346, MPFR_RNDN); mpfr_add(r25348, r25331, r25347, MPFR_RNDN); return mpfr_get_d(r25348, MPFR_RNDN); } static mpfr_t r25349, r25350, r25351, r25352, r25353, r25354, r25355, r25356, r25357, r25358, r25359, r25360, r25361, r25362, r25363, r25364, r25365, r25366; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r25349); mpfr_init(r25350); mpfr_init(r25351); mpfr_init(r25352); mpfr_init(r25353); mpfr_init(r25354); mpfr_init(r25355); mpfr_init(r25356); mpfr_init(r25357); mpfr_init_set_str(r25358, "1", 10, MPFR_RNDN); mpfr_init(r25359); mpfr_init(r25360); mpfr_init(r25361); mpfr_init(r25362); mpfr_init(r25363); mpfr_init(r25364); mpfr_init(r25365); mpfr_init(r25366); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r25349, x, MPFR_RNDN); mpfr_set_d(r25350, y, MPFR_RNDN); mpfr_tan(r25351, r25350, MPFR_RNDN); mpfr_set_d(r25352, z, MPFR_RNDN); mpfr_tan(r25353, r25352, MPFR_RNDN); mpfr_add(r25354, r25351, r25353, MPFR_RNDN); mpfr_set_d(r25355, a, MPFR_RNDN); mpfr_cos(r25356, r25355, MPFR_RNDN); mpfr_mul(r25357, r25354, r25356, MPFR_RNDN); ; mpfr_mul(r25359, r25351, r25353, MPFR_RNDN); mpfr_sub(r25360, r25358, r25359, MPFR_RNDN); mpfr_sin(r25361, r25355, MPFR_RNDN); mpfr_mul(r25362, r25360, r25361, MPFR_RNDN); mpfr_sub(r25363, r25357, r25362, MPFR_RNDN); mpfr_mul(r25364, r25360, r25356, MPFR_RNDN); mpfr_div(r25365, r25363, r25364, MPFR_RNDN); mpfr_add(r25366, r25349, r25365, MPFR_RNDN); return mpfr_get_d(r25366, MPFR_RNDN); }