#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 r27546 = x; float r27547 = y; float r27548 = z; float r27549 = r27547 + r27548; float r27550 = tan(r27549); float r27551 = a; float r27552 = tan(r27551); float r27553 = r27550 - r27552; float r27554 = r27546 + r27553; return r27554; } double f_id(double x, double y, double z, double a) { double r27555 = x; double r27556 = y; double r27557 = z; double r27558 = r27556 + r27557; double r27559 = tan(r27558); double r27560 = a; double r27561 = tan(r27560); double r27562 = r27559 - r27561; double r27563 = r27555 + r27562; return r27563; } double f_of(float x, float y, float z, float a) { float r27564 = x; float r27565 = y; float r27566 = tan(r27565); float r27567 = z; float r27568 = tan(r27567); float r27569 = r27566 + r27568; float r27570 = 1; float r27571 = 3; float r27572 = pow(r27570, r27571); float r27573 = r27566 * r27568; float r27574 = pow(r27573, r27571); float r27575 = r27572 - r27574; float r27576 = r27569 / r27575; float r27577 = r27570 * r27570; float r27578 = r27573 * r27573; float r27579 = r27570 * r27573; float r27580 = r27578 + r27579; float r27581 = r27577 + r27580; float r27582 = a; float r27583 = tan(r27582); float r27584 = -r27583; float r27585 = fma(r27576, r27581, r27584); float r27586 = r27564 + r27585; return r27586; } double f_od(double x, double y, double z, double a) { double r27587 = x; double r27588 = y; double r27589 = tan(r27588); double r27590 = z; double r27591 = tan(r27590); double r27592 = r27589 + r27591; double r27593 = 1; double r27594 = 3; double r27595 = pow(r27593, r27594); double r27596 = r27589 * r27591; double r27597 = pow(r27596, r27594); double r27598 = r27595 - r27597; double r27599 = r27592 / r27598; double r27600 = r27593 * r27593; double r27601 = r27596 * r27596; double r27602 = r27593 * r27596; double r27603 = r27601 + r27602; double r27604 = r27600 + r27603; double r27605 = a; double r27606 = tan(r27605); double r27607 = -r27606; double r27608 = fma(r27599, r27604, r27607); double r27609 = r27587 + r27608; return r27609; } 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 r27610, r27611, r27612, r27613, r27614, r27615, r27616, r27617, r27618; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r27610); mpfr_init(r27611); mpfr_init(r27612); mpfr_init(r27613); mpfr_init(r27614); mpfr_init(r27615); mpfr_init(r27616); mpfr_init(r27617); mpfr_init(r27618); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r27610, x, MPFR_RNDN); mpfr_set_d(r27611, y, MPFR_RNDN); mpfr_set_d(r27612, z, MPFR_RNDN); mpfr_add(r27613, r27611, r27612, MPFR_RNDN); mpfr_tan(r27614, r27613, MPFR_RNDN); mpfr_set_d(r27615, a, MPFR_RNDN); mpfr_tan(r27616, r27615, MPFR_RNDN); mpfr_sub(r27617, r27614, r27616, MPFR_RNDN); mpfr_add(r27618, r27610, r27617, MPFR_RNDN); return mpfr_get_d(r27618, MPFR_RNDN); } static mpfr_t r27619, r27620, r27621, r27622, r27623, r27624, r27625, r27626, r27627, r27628, r27629, r27630, r27631, r27632, r27633, r27634, r27635, r27636, r27637, r27638, r27639, r27640, r27641; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r27619); mpfr_init(r27620); mpfr_init(r27621); mpfr_init(r27622); mpfr_init(r27623); mpfr_init(r27624); mpfr_init_set_str(r27625, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27626, "3", 10, MPFR_RNDN); mpfr_init(r27627); mpfr_init(r27628); mpfr_init(r27629); mpfr_init(r27630); mpfr_init(r27631); mpfr_init(r27632); mpfr_init(r27633); mpfr_init(r27634); mpfr_init(r27635); mpfr_init(r27636); mpfr_init(r27637); mpfr_init(r27638); mpfr_init(r27639); mpfr_init(r27640); mpfr_init(r27641); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r27619, x, MPFR_RNDN); mpfr_set_d(r27620, y, MPFR_RNDN); mpfr_tan(r27621, r27620, MPFR_RNDN); mpfr_set_d(r27622, z, MPFR_RNDN); mpfr_tan(r27623, r27622, MPFR_RNDN); mpfr_add(r27624, r27621, r27623, MPFR_RNDN); ; ; mpfr_pow(r27627, r27625, r27626, MPFR_RNDN); mpfr_mul(r27628, r27621, r27623, MPFR_RNDN); mpfr_pow(r27629, r27628, r27626, MPFR_RNDN); mpfr_sub(r27630, r27627, r27629, MPFR_RNDN); mpfr_div(r27631, r27624, r27630, MPFR_RNDN); mpfr_mul(r27632, r27625, r27625, MPFR_RNDN); mpfr_mul(r27633, r27628, r27628, MPFR_RNDN); mpfr_mul(r27634, r27625, r27628, MPFR_RNDN); mpfr_add(r27635, r27633, r27634, MPFR_RNDN); mpfr_add(r27636, r27632, r27635, MPFR_RNDN); mpfr_set_d(r27637, a, MPFR_RNDN); mpfr_tan(r27638, r27637, MPFR_RNDN); mpfr_neg(r27639, r27638, MPFR_RNDN); mpfr_fma(r27640, r27631, r27636, r27639, MPFR_RNDN); mpfr_add(r27641, r27619, r27640, MPFR_RNDN); return mpfr_get_d(r27641, MPFR_RNDN); } static mpfr_t r27642, r27643, r27644, r27645, r27646, r27647, r27648, r27649, r27650, r27651, r27652, r27653, r27654, r27655, r27656, r27657, r27658, r27659, r27660, r27661, r27662, r27663, r27664; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r27642); mpfr_init(r27643); mpfr_init(r27644); mpfr_init(r27645); mpfr_init(r27646); mpfr_init(r27647); mpfr_init_set_str(r27648, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27649, "3", 10, MPFR_RNDN); mpfr_init(r27650); mpfr_init(r27651); mpfr_init(r27652); mpfr_init(r27653); mpfr_init(r27654); mpfr_init(r27655); mpfr_init(r27656); mpfr_init(r27657); mpfr_init(r27658); mpfr_init(r27659); mpfr_init(r27660); mpfr_init(r27661); mpfr_init(r27662); mpfr_init(r27663); mpfr_init(r27664); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r27642, x, MPFR_RNDN); mpfr_set_d(r27643, y, MPFR_RNDN); mpfr_tan(r27644, r27643, MPFR_RNDN); mpfr_set_d(r27645, z, MPFR_RNDN); mpfr_tan(r27646, r27645, MPFR_RNDN); mpfr_add(r27647, r27644, r27646, MPFR_RNDN); ; ; mpfr_pow(r27650, r27648, r27649, MPFR_RNDN); mpfr_mul(r27651, r27644, r27646, MPFR_RNDN); mpfr_pow(r27652, r27651, r27649, MPFR_RNDN); mpfr_sub(r27653, r27650, r27652, MPFR_RNDN); mpfr_div(r27654, r27647, r27653, MPFR_RNDN); mpfr_mul(r27655, r27648, r27648, MPFR_RNDN); mpfr_mul(r27656, r27651, r27651, MPFR_RNDN); mpfr_mul(r27657, r27648, r27651, MPFR_RNDN); mpfr_add(r27658, r27656, r27657, MPFR_RNDN); mpfr_add(r27659, r27655, r27658, MPFR_RNDN); mpfr_set_d(r27660, a, MPFR_RNDN); mpfr_tan(r27661, r27660, MPFR_RNDN); mpfr_neg(r27662, r27661, MPFR_RNDN); mpfr_fma(r27663, r27654, r27659, r27662, MPFR_RNDN); mpfr_add(r27664, r27642, r27663, MPFR_RNDN); return mpfr_get_d(r27664, MPFR_RNDN); }