#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 r25316 = x; float r25317 = y; float r25318 = z; float r25319 = r25317 + r25318; float r25320 = tan(r25319); float r25321 = a; float r25322 = tan(r25321); float r25323 = r25320 - r25322; float r25324 = r25316 + r25323; return r25324; } double f_id(double x, double y, double z, double a) { double r25325 = x; double r25326 = y; double r25327 = z; double r25328 = r25326 + r25327; double r25329 = tan(r25328); double r25330 = a; double r25331 = tan(r25330); double r25332 = r25329 - r25331; double r25333 = r25325 + r25332; return r25333; } double f_of(float x, float y, float z, float a) { float r25334 = x; float r25335 = y; float r25336 = tan(r25335); float r25337 = z; float r25338 = tan(r25337); float r25339 = r25336 + r25338; float r25340 = a; float r25341 = cos(r25340); float r25342 = r25339 * r25341; float r25343 = 1; float r25344 = r25336 * r25338; float r25345 = r25343 - r25344; float r25346 = sin(r25340); float r25347 = r25345 * r25346; float r25348 = r25342 - r25347; float r25349 = r25345 * r25341; float r25350 = r25348 / r25349; float r25351 = r25334 + r25350; return r25351; } double f_od(double x, double y, double z, double a) { double r25352 = x; double r25353 = y; double r25354 = tan(r25353); double r25355 = z; double r25356 = tan(r25355); double r25357 = r25354 + r25356; double r25358 = a; double r25359 = cos(r25358); double r25360 = r25357 * r25359; double r25361 = 1; double r25362 = r25354 * r25356; double r25363 = r25361 - r25362; double r25364 = sin(r25358); double r25365 = r25363 * r25364; double r25366 = r25360 - r25365; double r25367 = r25363 * r25359; double r25368 = r25366 / r25367; double r25369 = r25352 + r25368; return r25369; } 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 r25370, r25371, r25372, r25373, r25374, r25375, r25376, r25377, r25378; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r25370); mpfr_init(r25371); mpfr_init(r25372); mpfr_init(r25373); mpfr_init(r25374); mpfr_init(r25375); mpfr_init(r25376); mpfr_init(r25377); mpfr_init(r25378); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r25370, x, MPFR_RNDN); mpfr_set_d(r25371, y, MPFR_RNDN); mpfr_set_d(r25372, z, MPFR_RNDN); mpfr_add(r25373, r25371, r25372, MPFR_RNDN); mpfr_tan(r25374, r25373, MPFR_RNDN); mpfr_set_d(r25375, a, MPFR_RNDN); mpfr_tan(r25376, r25375, MPFR_RNDN); mpfr_sub(r25377, r25374, r25376, MPFR_RNDN); mpfr_add(r25378, r25370, r25377, MPFR_RNDN); return mpfr_get_d(r25378, MPFR_RNDN); } static mpfr_t r25379, r25380, r25381, r25382, r25383, r25384, r25385, r25386, r25387, r25388, r25389, r25390, r25391, r25392, r25393, r25394, r25395, r25396; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25379); mpfr_init(r25380); mpfr_init(r25381); mpfr_init(r25382); mpfr_init(r25383); mpfr_init(r25384); mpfr_init(r25385); mpfr_init(r25386); mpfr_init(r25387); mpfr_init_set_str(r25388, "1", 10, MPFR_RNDN); mpfr_init(r25389); mpfr_init(r25390); mpfr_init(r25391); mpfr_init(r25392); mpfr_init(r25393); mpfr_init(r25394); mpfr_init(r25395); mpfr_init(r25396); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r25379, x, MPFR_RNDN); mpfr_set_d(r25380, y, MPFR_RNDN); mpfr_tan(r25381, r25380, MPFR_RNDN); mpfr_set_d(r25382, z, MPFR_RNDN); mpfr_tan(r25383, r25382, MPFR_RNDN); mpfr_add(r25384, r25381, r25383, MPFR_RNDN); mpfr_set_d(r25385, a, MPFR_RNDN); mpfr_cos(r25386, r25385, MPFR_RNDN); mpfr_mul(r25387, r25384, r25386, MPFR_RNDN); ; mpfr_mul(r25389, r25381, r25383, MPFR_RNDN); mpfr_sub(r25390, r25388, r25389, MPFR_RNDN); mpfr_sin(r25391, r25385, MPFR_RNDN); mpfr_mul(r25392, r25390, r25391, MPFR_RNDN); mpfr_sub(r25393, r25387, r25392, MPFR_RNDN); mpfr_mul(r25394, r25390, r25386, MPFR_RNDN); mpfr_div(r25395, r25393, r25394, MPFR_RNDN); mpfr_add(r25396, r25379, r25395, MPFR_RNDN); return mpfr_get_d(r25396, MPFR_RNDN); } static mpfr_t r25397, r25398, r25399, r25400, r25401, r25402, r25403, r25404, r25405, r25406, r25407, r25408, r25409, r25410, r25411, r25412, r25413, r25414; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r25397); mpfr_init(r25398); mpfr_init(r25399); mpfr_init(r25400); mpfr_init(r25401); mpfr_init(r25402); mpfr_init(r25403); mpfr_init(r25404); mpfr_init(r25405); mpfr_init_set_str(r25406, "1", 10, MPFR_RNDN); mpfr_init(r25407); mpfr_init(r25408); mpfr_init(r25409); mpfr_init(r25410); mpfr_init(r25411); mpfr_init(r25412); mpfr_init(r25413); mpfr_init(r25414); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r25397, x, MPFR_RNDN); mpfr_set_d(r25398, y, MPFR_RNDN); mpfr_tan(r25399, r25398, MPFR_RNDN); mpfr_set_d(r25400, z, MPFR_RNDN); mpfr_tan(r25401, r25400, MPFR_RNDN); mpfr_add(r25402, r25399, r25401, MPFR_RNDN); mpfr_set_d(r25403, a, MPFR_RNDN); mpfr_cos(r25404, r25403, MPFR_RNDN); mpfr_mul(r25405, r25402, r25404, MPFR_RNDN); ; mpfr_mul(r25407, r25399, r25401, MPFR_RNDN); mpfr_sub(r25408, r25406, r25407, MPFR_RNDN); mpfr_sin(r25409, r25403, MPFR_RNDN); mpfr_mul(r25410, r25408, r25409, MPFR_RNDN); mpfr_sub(r25411, r25405, r25410, MPFR_RNDN); mpfr_mul(r25412, r25408, r25404, MPFR_RNDN); mpfr_div(r25413, r25411, r25412, MPFR_RNDN); mpfr_add(r25414, r25397, r25413, MPFR_RNDN); return mpfr_get_d(r25414, MPFR_RNDN); }