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