#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 r25677 = x; float r25678 = y; float r25679 = z; float r25680 = r25678 + r25679; float r25681 = tan(r25680); float r25682 = a; float r25683 = tan(r25682); float r25684 = r25681 - r25683; float r25685 = r25677 + r25684; return r25685; } double f_id(double x, double y, double z, double a) { double r25686 = x; double r25687 = y; double r25688 = z; double r25689 = r25687 + r25688; double r25690 = tan(r25689); double r25691 = a; double r25692 = tan(r25691); double r25693 = r25690 - r25692; double r25694 = r25686 + r25693; return r25694; } double f_of(float x, float y, float z, float a) { float r25695 = x; float r25696 = z; float r25697 = tan(r25696); float r25698 = y; float r25699 = tan(r25698); float r25700 = r25697 + r25699; float r25701 = -r25700; float r25702 = 1; float r25703 = -r25702; float r25704 = fma(r25697, r25699, r25703); float r25705 = cbrt(r25704); float r25706 = r25701 / r25705; float r25707 = r25705 * r25705; float r25708 = r25706 / r25707; float r25709 = a; float r25710 = tan(r25709); float r25711 = r25708 - r25710; float r25712 = 0; float r25713 = r25711 + r25712; float r25714 = r25695 + r25713; return r25714; } double f_od(double x, double y, double z, double a) { double r25715 = x; double r25716 = z; double r25717 = tan(r25716); double r25718 = y; double r25719 = tan(r25718); double r25720 = r25717 + r25719; double r25721 = -r25720; double r25722 = 1; double r25723 = -r25722; double r25724 = fma(r25717, r25719, r25723); double r25725 = cbrt(r25724); double r25726 = r25721 / r25725; double r25727 = r25725 * r25725; double r25728 = r25726 / r25727; double r25729 = a; double r25730 = tan(r25729); double r25731 = r25728 - r25730; double r25732 = 0; double r25733 = r25731 + r25732; double r25734 = r25715 + r25733; return r25734; } 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 r25735, r25736, r25737, r25738, r25739, r25740, r25741, r25742, r25743; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r25735); mpfr_init(r25736); mpfr_init(r25737); mpfr_init(r25738); mpfr_init(r25739); mpfr_init(r25740); mpfr_init(r25741); mpfr_init(r25742); mpfr_init(r25743); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r25735, x, MPFR_RNDN); mpfr_set_d(r25736, y, MPFR_RNDN); mpfr_set_d(r25737, z, MPFR_RNDN); mpfr_add(r25738, r25736, r25737, MPFR_RNDN); mpfr_tan(r25739, r25738, MPFR_RNDN); mpfr_set_d(r25740, a, MPFR_RNDN); mpfr_tan(r25741, r25740, MPFR_RNDN); mpfr_sub(r25742, r25739, r25741, MPFR_RNDN); mpfr_add(r25743, r25735, r25742, MPFR_RNDN); return mpfr_get_d(r25743, MPFR_RNDN); } static mpfr_t r25744, r25745, r25746, r25747, r25748, r25749, r25750, r25751, r25752, r25753, r25754, r25755, r25756, r25757, r25758, r25759, r25760, r25761, r25762, r25763; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25744); mpfr_init(r25745); mpfr_init(r25746); mpfr_init(r25747); mpfr_init(r25748); mpfr_init(r25749); mpfr_init(r25750); mpfr_init_set_str(r25751, "1", 10, MPFR_RNDN); mpfr_init(r25752); mpfr_init(r25753); mpfr_init(r25754); mpfr_init(r25755); mpfr_init(r25756); mpfr_init(r25757); mpfr_init(r25758); mpfr_init(r25759); mpfr_init(r25760); mpfr_init_set_str(r25761, "0", 10, MPFR_RNDN); mpfr_init(r25762); mpfr_init(r25763); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r25744, x, MPFR_RNDN); mpfr_set_d(r25745, z, MPFR_RNDN); mpfr_tan(r25746, r25745, MPFR_RNDN); mpfr_set_d(r25747, y, MPFR_RNDN); mpfr_tan(r25748, r25747, MPFR_RNDN); mpfr_add(r25749, r25746, r25748, MPFR_RNDN); mpfr_neg(r25750, r25749, MPFR_RNDN); ; mpfr_neg(r25752, r25751, MPFR_RNDN); mpfr_fma(r25753, r25746, r25748, r25752, MPFR_RNDN); mpfr_cbrt(r25754, r25753, MPFR_RNDN); mpfr_div(r25755, r25750, r25754, MPFR_RNDN); mpfr_mul(r25756, r25754, r25754, MPFR_RNDN); mpfr_div(r25757, r25755, r25756, MPFR_RNDN); mpfr_set_d(r25758, a, MPFR_RNDN); mpfr_tan(r25759, r25758, MPFR_RNDN); mpfr_sub(r25760, r25757, r25759, MPFR_RNDN); ; mpfr_add(r25762, r25760, r25761, MPFR_RNDN); mpfr_add(r25763, r25744, r25762, MPFR_RNDN); return mpfr_get_d(r25763, MPFR_RNDN); } static mpfr_t r25764, r25765, r25766, r25767, r25768, r25769, r25770, r25771, r25772, r25773, r25774, r25775, r25776, r25777, r25778, r25779, r25780, r25781, r25782, r25783; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r25764); mpfr_init(r25765); mpfr_init(r25766); mpfr_init(r25767); mpfr_init(r25768); mpfr_init(r25769); mpfr_init(r25770); mpfr_init_set_str(r25771, "1", 10, MPFR_RNDN); mpfr_init(r25772); mpfr_init(r25773); mpfr_init(r25774); mpfr_init(r25775); mpfr_init(r25776); mpfr_init(r25777); mpfr_init(r25778); mpfr_init(r25779); mpfr_init(r25780); mpfr_init_set_str(r25781, "0", 10, MPFR_RNDN); mpfr_init(r25782); mpfr_init(r25783); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r25764, x, MPFR_RNDN); mpfr_set_d(r25765, z, MPFR_RNDN); mpfr_tan(r25766, r25765, MPFR_RNDN); mpfr_set_d(r25767, y, MPFR_RNDN); mpfr_tan(r25768, r25767, MPFR_RNDN); mpfr_add(r25769, r25766, r25768, MPFR_RNDN); mpfr_neg(r25770, r25769, MPFR_RNDN); ; mpfr_neg(r25772, r25771, MPFR_RNDN); mpfr_fma(r25773, r25766, r25768, r25772, MPFR_RNDN); mpfr_cbrt(r25774, r25773, MPFR_RNDN); mpfr_div(r25775, r25770, r25774, MPFR_RNDN); mpfr_mul(r25776, r25774, r25774, MPFR_RNDN); mpfr_div(r25777, r25775, r25776, MPFR_RNDN); mpfr_set_d(r25778, a, MPFR_RNDN); mpfr_tan(r25779, r25778, MPFR_RNDN); mpfr_sub(r25780, r25777, r25779, MPFR_RNDN); ; mpfr_add(r25782, r25780, r25781, MPFR_RNDN); mpfr_add(r25783, r25764, r25782, MPFR_RNDN); return mpfr_get_d(r25783, MPFR_RNDN); }