#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 r25386 = x; float r25387 = y; float r25388 = z; float r25389 = r25387 + r25388; float r25390 = tan(r25389); float r25391 = a; float r25392 = tan(r25391); float r25393 = r25390 - r25392; float r25394 = r25386 + r25393; return r25394; } double f_id(double x, double y, double z, double a) { double r25395 = x; double r25396 = y; double r25397 = z; double r25398 = r25396 + r25397; double r25399 = tan(r25398); double r25400 = a; double r25401 = tan(r25400); double r25402 = r25399 - r25401; double r25403 = r25395 + r25402; return r25403; } double f_of(float x, float y, float z, float a) { float r25404 = x; float r25405 = y; float r25406 = tan(r25405); float r25407 = z; float r25408 = tan(r25407); float r25409 = r25406 + r25408; float r25410 = a; float r25411 = cos(r25410); float r25412 = r25409 * r25411; float r25413 = 1; float r25414 = r25406 * r25408; float r25415 = r25413 - r25414; float r25416 = sin(r25410); float r25417 = r25415 * r25416; float r25418 = r25412 - r25417; float r25419 = r25415 * r25411; float r25420 = r25418 / r25419; float r25421 = r25404 + r25420; return r25421; } double f_od(double x, double y, double z, double a) { double r25422 = x; double r25423 = y; double r25424 = tan(r25423); double r25425 = z; double r25426 = tan(r25425); double r25427 = r25424 + r25426; double r25428 = a; double r25429 = cos(r25428); double r25430 = r25427 * r25429; double r25431 = 1; double r25432 = r25424 * r25426; double r25433 = r25431 - r25432; double r25434 = sin(r25428); double r25435 = r25433 * r25434; double r25436 = r25430 - r25435; double r25437 = r25433 * r25429; double r25438 = r25436 / r25437; double r25439 = r25422 + r25438; return r25439; } 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 r25440, r25441, r25442, r25443, r25444, r25445, r25446, r25447, r25448; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r25440); mpfr_init(r25441); mpfr_init(r25442); mpfr_init(r25443); mpfr_init(r25444); mpfr_init(r25445); mpfr_init(r25446); mpfr_init(r25447); mpfr_init(r25448); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r25440, x, MPFR_RNDN); mpfr_set_d(r25441, y, MPFR_RNDN); mpfr_set_d(r25442, z, MPFR_RNDN); mpfr_add(r25443, r25441, r25442, MPFR_RNDN); mpfr_tan(r25444, r25443, MPFR_RNDN); mpfr_set_d(r25445, a, MPFR_RNDN); mpfr_tan(r25446, r25445, MPFR_RNDN); mpfr_sub(r25447, r25444, r25446, MPFR_RNDN); mpfr_add(r25448, r25440, r25447, MPFR_RNDN); return mpfr_get_d(r25448, MPFR_RNDN); } static mpfr_t r25449, r25450, r25451, r25452, r25453, r25454, r25455, r25456, r25457, r25458, r25459, r25460, r25461, r25462, r25463, r25464, r25465, r25466; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25449); mpfr_init(r25450); mpfr_init(r25451); mpfr_init(r25452); mpfr_init(r25453); mpfr_init(r25454); mpfr_init(r25455); mpfr_init(r25456); mpfr_init(r25457); mpfr_init_set_str(r25458, "1", 10, MPFR_RNDN); mpfr_init(r25459); mpfr_init(r25460); mpfr_init(r25461); mpfr_init(r25462); mpfr_init(r25463); mpfr_init(r25464); mpfr_init(r25465); mpfr_init(r25466); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r25449, x, MPFR_RNDN); mpfr_set_d(r25450, y, MPFR_RNDN); mpfr_tan(r25451, r25450, MPFR_RNDN); mpfr_set_d(r25452, z, MPFR_RNDN); mpfr_tan(r25453, r25452, MPFR_RNDN); mpfr_add(r25454, r25451, r25453, MPFR_RNDN); mpfr_set_d(r25455, a, MPFR_RNDN); mpfr_cos(r25456, r25455, MPFR_RNDN); mpfr_mul(r25457, r25454, r25456, MPFR_RNDN); ; mpfr_mul(r25459, r25451, r25453, MPFR_RNDN); mpfr_sub(r25460, r25458, r25459, MPFR_RNDN); mpfr_sin(r25461, r25455, MPFR_RNDN); mpfr_mul(r25462, r25460, r25461, MPFR_RNDN); mpfr_sub(r25463, r25457, r25462, MPFR_RNDN); mpfr_mul(r25464, r25460, r25456, MPFR_RNDN); mpfr_div(r25465, r25463, r25464, MPFR_RNDN); mpfr_add(r25466, r25449, r25465, MPFR_RNDN); return mpfr_get_d(r25466, MPFR_RNDN); } static mpfr_t r25467, r25468, r25469, r25470, r25471, r25472, r25473, r25474, r25475, r25476, r25477, r25478, r25479, r25480, r25481, r25482, r25483, r25484; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r25467); mpfr_init(r25468); mpfr_init(r25469); mpfr_init(r25470); mpfr_init(r25471); mpfr_init(r25472); mpfr_init(r25473); mpfr_init(r25474); mpfr_init(r25475); mpfr_init_set_str(r25476, "1", 10, MPFR_RNDN); mpfr_init(r25477); mpfr_init(r25478); mpfr_init(r25479); mpfr_init(r25480); mpfr_init(r25481); mpfr_init(r25482); mpfr_init(r25483); mpfr_init(r25484); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r25467, x, MPFR_RNDN); mpfr_set_d(r25468, y, MPFR_RNDN); mpfr_tan(r25469, r25468, MPFR_RNDN); mpfr_set_d(r25470, z, MPFR_RNDN); mpfr_tan(r25471, r25470, MPFR_RNDN); mpfr_add(r25472, r25469, r25471, MPFR_RNDN); mpfr_set_d(r25473, a, MPFR_RNDN); mpfr_cos(r25474, r25473, MPFR_RNDN); mpfr_mul(r25475, r25472, r25474, MPFR_RNDN); ; mpfr_mul(r25477, r25469, r25471, MPFR_RNDN); mpfr_sub(r25478, r25476, r25477, MPFR_RNDN); mpfr_sin(r25479, r25473, MPFR_RNDN); mpfr_mul(r25480, r25478, r25479, MPFR_RNDN); mpfr_sub(r25481, r25475, r25480, MPFR_RNDN); mpfr_mul(r25482, r25478, r25474, MPFR_RNDN); mpfr_div(r25483, r25481, r25482, MPFR_RNDN); mpfr_add(r25484, r25467, r25483, MPFR_RNDN); return mpfr_get_d(r25484, MPFR_RNDN); }