#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 r24377 = x; float r24378 = y; float r24379 = z; float r24380 = r24378 + r24379; float r24381 = tan(r24380); float r24382 = a; float r24383 = tan(r24382); float r24384 = r24381 - r24383; float r24385 = r24377 + r24384; return r24385; } double f_id(double x, double y, double z, double a) { double r24386 = x; double r24387 = y; double r24388 = z; double r24389 = r24387 + r24388; double r24390 = tan(r24389); double r24391 = a; double r24392 = tan(r24391); double r24393 = r24390 - r24392; double r24394 = r24386 + r24393; return r24394; } double f_of(float x, float y, float z, float a) { float r24395 = x; float r24396 = y; float r24397 = tan(r24396); float r24398 = z; float r24399 = tan(r24398); float r24400 = r24397 + r24399; float r24401 = 1; float r24402 = sin(r24396); float r24403 = sin(r24398); float r24404 = r24402 * r24403; float r24405 = cos(r24396); float r24406 = cos(r24398); float r24407 = r24405 * r24406; float r24408 = r24404 / r24407; float r24409 = r24401 - r24408; float r24410 = r24400 / r24409; float r24411 = a; float r24412 = tan(r24411); float r24413 = r24410 - r24412; float r24414 = r24395 + r24413; return r24414; } double f_od(double x, double y, double z, double a) { double r24415 = x; double r24416 = y; double r24417 = tan(r24416); double r24418 = z; double r24419 = tan(r24418); double r24420 = r24417 + r24419; double r24421 = 1; double r24422 = sin(r24416); double r24423 = sin(r24418); double r24424 = r24422 * r24423; double r24425 = cos(r24416); double r24426 = cos(r24418); double r24427 = r24425 * r24426; double r24428 = r24424 / r24427; double r24429 = r24421 - r24428; double r24430 = r24420 / r24429; double r24431 = a; double r24432 = tan(r24431); double r24433 = r24430 - r24432; double r24434 = r24415 + r24433; return r24434; } 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 r24435, r24436, r24437, r24438, r24439, r24440, r24441, r24442, r24443; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r24435); mpfr_init(r24436); mpfr_init(r24437); mpfr_init(r24438); mpfr_init(r24439); mpfr_init(r24440); mpfr_init(r24441); mpfr_init(r24442); mpfr_init(r24443); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r24435, x, MPFR_RNDN); mpfr_set_d(r24436, y, MPFR_RNDN); mpfr_set_d(r24437, z, MPFR_RNDN); mpfr_add(r24438, r24436, r24437, MPFR_RNDN); mpfr_tan(r24439, r24438, MPFR_RNDN); mpfr_set_d(r24440, a, MPFR_RNDN); mpfr_tan(r24441, r24440, MPFR_RNDN); mpfr_sub(r24442, r24439, r24441, MPFR_RNDN); mpfr_add(r24443, r24435, r24442, MPFR_RNDN); return mpfr_get_d(r24443, MPFR_RNDN); } static mpfr_t r24444, r24445, r24446, r24447, r24448, r24449, r24450, r24451, r24452, r24453, r24454, r24455, r24456, r24457, r24458, r24459, r24460, r24461, r24462, r24463; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r24444); mpfr_init(r24445); mpfr_init(r24446); mpfr_init(r24447); mpfr_init(r24448); mpfr_init(r24449); mpfr_init_set_str(r24450, "1", 10, MPFR_RNDN); mpfr_init(r24451); mpfr_init(r24452); mpfr_init(r24453); mpfr_init(r24454); mpfr_init(r24455); mpfr_init(r24456); mpfr_init(r24457); mpfr_init(r24458); mpfr_init(r24459); mpfr_init(r24460); mpfr_init(r24461); mpfr_init(r24462); mpfr_init(r24463); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r24444, x, MPFR_RNDN); mpfr_set_d(r24445, y, MPFR_RNDN); mpfr_tan(r24446, r24445, MPFR_RNDN); mpfr_set_d(r24447, z, MPFR_RNDN); mpfr_tan(r24448, r24447, MPFR_RNDN); mpfr_add(r24449, r24446, r24448, MPFR_RNDN); ; mpfr_sin(r24451, r24445, MPFR_RNDN); mpfr_sin(r24452, r24447, MPFR_RNDN); mpfr_mul(r24453, r24451, r24452, MPFR_RNDN); mpfr_cos(r24454, r24445, MPFR_RNDN); mpfr_cos(r24455, r24447, MPFR_RNDN); mpfr_mul(r24456, r24454, r24455, MPFR_RNDN); mpfr_div(r24457, r24453, r24456, MPFR_RNDN); mpfr_sub(r24458, r24450, r24457, MPFR_RNDN); mpfr_div(r24459, r24449, r24458, MPFR_RNDN); mpfr_set_d(r24460, a, MPFR_RNDN); mpfr_tan(r24461, r24460, MPFR_RNDN); mpfr_sub(r24462, r24459, r24461, MPFR_RNDN); mpfr_add(r24463, r24444, r24462, MPFR_RNDN); return mpfr_get_d(r24463, MPFR_RNDN); } static mpfr_t r24464, r24465, r24466, r24467, r24468, r24469, r24470, r24471, r24472, r24473, r24474, r24475, r24476, r24477, r24478, r24479, r24480, r24481, r24482, r24483; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r24464); mpfr_init(r24465); mpfr_init(r24466); mpfr_init(r24467); mpfr_init(r24468); mpfr_init(r24469); mpfr_init_set_str(r24470, "1", 10, MPFR_RNDN); mpfr_init(r24471); mpfr_init(r24472); mpfr_init(r24473); mpfr_init(r24474); mpfr_init(r24475); mpfr_init(r24476); mpfr_init(r24477); mpfr_init(r24478); mpfr_init(r24479); mpfr_init(r24480); mpfr_init(r24481); mpfr_init(r24482); mpfr_init(r24483); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r24464, x, MPFR_RNDN); mpfr_set_d(r24465, y, MPFR_RNDN); mpfr_tan(r24466, r24465, MPFR_RNDN); mpfr_set_d(r24467, z, MPFR_RNDN); mpfr_tan(r24468, r24467, MPFR_RNDN); mpfr_add(r24469, r24466, r24468, MPFR_RNDN); ; mpfr_sin(r24471, r24465, MPFR_RNDN); mpfr_sin(r24472, r24467, MPFR_RNDN); mpfr_mul(r24473, r24471, r24472, MPFR_RNDN); mpfr_cos(r24474, r24465, MPFR_RNDN); mpfr_cos(r24475, r24467, MPFR_RNDN); mpfr_mul(r24476, r24474, r24475, MPFR_RNDN); mpfr_div(r24477, r24473, r24476, MPFR_RNDN); mpfr_sub(r24478, r24470, r24477, MPFR_RNDN); mpfr_div(r24479, r24469, r24478, MPFR_RNDN); mpfr_set_d(r24480, a, MPFR_RNDN); mpfr_tan(r24481, r24480, MPFR_RNDN); mpfr_sub(r24482, r24479, r24481, MPFR_RNDN); mpfr_add(r24483, r24464, r24482, MPFR_RNDN); return mpfr_get_d(r24483, MPFR_RNDN); }