#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 r25501 = x; float r25502 = y; float r25503 = z; float r25504 = r25502 + r25503; float r25505 = tan(r25504); float r25506 = a; float r25507 = tan(r25506); float r25508 = r25505 - r25507; float r25509 = r25501 + r25508; return r25509; } double f_id(double x, double y, double z, double a) { double r25510 = x; double r25511 = y; double r25512 = z; double r25513 = r25511 + r25512; double r25514 = tan(r25513); double r25515 = a; double r25516 = tan(r25515); double r25517 = r25514 - r25516; double r25518 = r25510 + r25517; return r25518; } double f_of(float x, float y, float z, float a) { float r25519 = x; float r25520 = 1; float r25521 = y; float r25522 = tan(r25521); float r25523 = z; float r25524 = tan(r25523); float r25525 = r25522 * r25524; float r25526 = r25520 - r25525; float r25527 = r25520 / r25526; float r25528 = a; float r25529 = sin(r25528); float r25530 = cos(r25528); float r25531 = r25529 / r25530; float r25532 = -r25520; float r25533 = fma(r25524, r25522, r25532); float r25534 = r25524 + r25522; float r25535 = fma(r25531, r25533, r25534); float r25536 = r25527 * r25535; float r25537 = r25519 + r25536; return r25537; } double f_od(double x, double y, double z, double a) { double r25538 = x; double r25539 = 1; double r25540 = y; double r25541 = tan(r25540); double r25542 = z; double r25543 = tan(r25542); double r25544 = r25541 * r25543; double r25545 = r25539 - r25544; double r25546 = r25539 / r25545; double r25547 = a; double r25548 = sin(r25547); double r25549 = cos(r25547); double r25550 = r25548 / r25549; double r25551 = -r25539; double r25552 = fma(r25543, r25541, r25551); double r25553 = r25543 + r25541; double r25554 = fma(r25550, r25552, r25553); double r25555 = r25546 * r25554; double r25556 = r25538 + r25555; return r25556; } 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 r25557, r25558, r25559, r25560, r25561, r25562, r25563, r25564, r25565; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r25557); mpfr_init(r25558); mpfr_init(r25559); mpfr_init(r25560); mpfr_init(r25561); mpfr_init(r25562); mpfr_init(r25563); mpfr_init(r25564); mpfr_init(r25565); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r25557, x, MPFR_RNDN); mpfr_set_d(r25558, y, MPFR_RNDN); mpfr_set_d(r25559, z, MPFR_RNDN); mpfr_add(r25560, r25558, r25559, MPFR_RNDN); mpfr_tan(r25561, r25560, MPFR_RNDN); mpfr_set_d(r25562, a, MPFR_RNDN); mpfr_tan(r25563, r25562, MPFR_RNDN); mpfr_sub(r25564, r25561, r25563, MPFR_RNDN); mpfr_add(r25565, r25557, r25564, MPFR_RNDN); return mpfr_get_d(r25565, MPFR_RNDN); } static mpfr_t r25566, r25567, r25568, r25569, r25570, r25571, r25572, r25573, r25574, r25575, r25576, r25577, r25578, r25579, r25580, r25581, r25582, r25583, r25584; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r25566); mpfr_init_set_str(r25567, "1", 10, MPFR_RNDN); mpfr_init(r25568); mpfr_init(r25569); mpfr_init(r25570); mpfr_init(r25571); mpfr_init(r25572); mpfr_init(r25573); mpfr_init(r25574); mpfr_init(r25575); mpfr_init(r25576); mpfr_init(r25577); mpfr_init(r25578); mpfr_init(r25579); mpfr_init(r25580); mpfr_init(r25581); mpfr_init(r25582); mpfr_init(r25583); mpfr_init(r25584); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r25566, x, MPFR_RNDN); ; mpfr_set_d(r25568, y, MPFR_RNDN); mpfr_tan(r25569, r25568, MPFR_RNDN); mpfr_set_d(r25570, z, MPFR_RNDN); mpfr_tan(r25571, r25570, MPFR_RNDN); mpfr_mul(r25572, r25569, r25571, MPFR_RNDN); mpfr_sub(r25573, r25567, r25572, MPFR_RNDN); mpfr_div(r25574, r25567, r25573, MPFR_RNDN); mpfr_set_d(r25575, a, MPFR_RNDN); mpfr_sin(r25576, r25575, MPFR_RNDN); mpfr_cos(r25577, r25575, MPFR_RNDN); mpfr_div(r25578, r25576, r25577, MPFR_RNDN); mpfr_neg(r25579, r25567, MPFR_RNDN); mpfr_fma(r25580, r25571, r25569, r25579, MPFR_RNDN); mpfr_add(r25581, r25571, r25569, MPFR_RNDN); mpfr_fma(r25582, r25578, r25580, r25581, MPFR_RNDN); mpfr_mul(r25583, r25574, r25582, MPFR_RNDN); mpfr_add(r25584, r25566, r25583, MPFR_RNDN); return mpfr_get_d(r25584, MPFR_RNDN); } static mpfr_t r25585, r25586, r25587, r25588, r25589, r25590, r25591, r25592, r25593, r25594, r25595, r25596, r25597, r25598, r25599, r25600, r25601, r25602, r25603; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r25585); mpfr_init_set_str(r25586, "1", 10, MPFR_RNDN); mpfr_init(r25587); mpfr_init(r25588); mpfr_init(r25589); mpfr_init(r25590); mpfr_init(r25591); mpfr_init(r25592); mpfr_init(r25593); mpfr_init(r25594); mpfr_init(r25595); mpfr_init(r25596); mpfr_init(r25597); mpfr_init(r25598); mpfr_init(r25599); mpfr_init(r25600); mpfr_init(r25601); mpfr_init(r25602); mpfr_init(r25603); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r25585, x, MPFR_RNDN); ; mpfr_set_d(r25587, y, MPFR_RNDN); mpfr_tan(r25588, r25587, MPFR_RNDN); mpfr_set_d(r25589, z, MPFR_RNDN); mpfr_tan(r25590, r25589, MPFR_RNDN); mpfr_mul(r25591, r25588, r25590, MPFR_RNDN); mpfr_sub(r25592, r25586, r25591, MPFR_RNDN); mpfr_div(r25593, r25586, r25592, MPFR_RNDN); mpfr_set_d(r25594, a, MPFR_RNDN); mpfr_sin(r25595, r25594, MPFR_RNDN); mpfr_cos(r25596, r25594, MPFR_RNDN); mpfr_div(r25597, r25595, r25596, MPFR_RNDN); mpfr_neg(r25598, r25586, MPFR_RNDN); mpfr_fma(r25599, r25590, r25588, r25598, MPFR_RNDN); mpfr_add(r25600, r25590, r25588, MPFR_RNDN); mpfr_fma(r25601, r25597, r25599, r25600, MPFR_RNDN); mpfr_mul(r25602, r25593, r25601, MPFR_RNDN); mpfr_add(r25603, r25585, r25602, MPFR_RNDN); return mpfr_get_d(r25603, MPFR_RNDN); }