#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 r24692 = x; float r24693 = y; float r24694 = z; float r24695 = r24693 + r24694; float r24696 = tan(r24695); float r24697 = a; float r24698 = tan(r24697); float r24699 = r24696 - r24698; float r24700 = r24692 + r24699; return r24700; } double f_id(double x, double y, double z, double a) { double r24701 = x; double r24702 = y; double r24703 = z; double r24704 = r24702 + r24703; double r24705 = tan(r24704); double r24706 = a; double r24707 = tan(r24706); double r24708 = r24705 - r24707; double r24709 = r24701 + r24708; return r24709; } double f_of(float x, float y, float z, float a) { float r24710 = x; float r24711 = y; float r24712 = tan(r24711); float r24713 = z; float r24714 = tan(r24713); float r24715 = r24712 + r24714; float r24716 = 1; float r24717 = r24716 * r24716; float r24718 = r24712 * r24714; float r24719 = r24718 * r24718; float r24720 = r24717 - r24719; float r24721 = r24715 / r24720; float r24722 = r24716 + r24718; float r24723 = a; float r24724 = tan(r24723); float r24725 = -r24724; float r24726 = fma(r24721, r24722, r24725); float r24727 = r24710 + r24726; return r24727; } double f_od(double x, double y, double z, double a) { double r24728 = x; double r24729 = y; double r24730 = tan(r24729); double r24731 = z; double r24732 = tan(r24731); double r24733 = r24730 + r24732; double r24734 = 1; double r24735 = r24734 * r24734; double r24736 = r24730 * r24732; double r24737 = r24736 * r24736; double r24738 = r24735 - r24737; double r24739 = r24733 / r24738; double r24740 = r24734 + r24736; double r24741 = a; double r24742 = tan(r24741); double r24743 = -r24742; double r24744 = fma(r24739, r24740, r24743); double r24745 = r24728 + r24744; return r24745; } 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 r24746, r24747, r24748, r24749, r24750, r24751, r24752, r24753, r24754; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r24746); mpfr_init(r24747); mpfr_init(r24748); mpfr_init(r24749); mpfr_init(r24750); mpfr_init(r24751); mpfr_init(r24752); mpfr_init(r24753); mpfr_init(r24754); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r24746, x, MPFR_RNDN); mpfr_set_d(r24747, y, MPFR_RNDN); mpfr_set_d(r24748, z, MPFR_RNDN); mpfr_add(r24749, r24747, r24748, MPFR_RNDN); mpfr_tan(r24750, r24749, MPFR_RNDN); mpfr_set_d(r24751, a, MPFR_RNDN); mpfr_tan(r24752, r24751, MPFR_RNDN); mpfr_sub(r24753, r24750, r24752, MPFR_RNDN); mpfr_add(r24754, r24746, r24753, MPFR_RNDN); return mpfr_get_d(r24754, MPFR_RNDN); } static mpfr_t r24755, r24756, r24757, r24758, r24759, r24760, r24761, r24762, r24763, r24764, r24765, r24766, r24767, r24768, r24769, r24770, r24771, r24772; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r24755); mpfr_init(r24756); mpfr_init(r24757); mpfr_init(r24758); mpfr_init(r24759); mpfr_init(r24760); mpfr_init_set_str(r24761, "1", 10, MPFR_RNDN); mpfr_init(r24762); mpfr_init(r24763); mpfr_init(r24764); mpfr_init(r24765); mpfr_init(r24766); mpfr_init(r24767); mpfr_init(r24768); mpfr_init(r24769); mpfr_init(r24770); mpfr_init(r24771); mpfr_init(r24772); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r24755, x, MPFR_RNDN); mpfr_set_d(r24756, y, MPFR_RNDN); mpfr_tan(r24757, r24756, MPFR_RNDN); mpfr_set_d(r24758, z, MPFR_RNDN); mpfr_tan(r24759, r24758, MPFR_RNDN); mpfr_add(r24760, r24757, r24759, MPFR_RNDN); ; mpfr_mul(r24762, r24761, r24761, MPFR_RNDN); mpfr_mul(r24763, r24757, r24759, MPFR_RNDN); mpfr_mul(r24764, r24763, r24763, MPFR_RNDN); mpfr_sub(r24765, r24762, r24764, MPFR_RNDN); mpfr_div(r24766, r24760, r24765, MPFR_RNDN); mpfr_add(r24767, r24761, r24763, MPFR_RNDN); mpfr_set_d(r24768, a, MPFR_RNDN); mpfr_tan(r24769, r24768, MPFR_RNDN); mpfr_neg(r24770, r24769, MPFR_RNDN); mpfr_fma(r24771, r24766, r24767, r24770, MPFR_RNDN); mpfr_add(r24772, r24755, r24771, MPFR_RNDN); return mpfr_get_d(r24772, MPFR_RNDN); } static mpfr_t r24773, r24774, r24775, r24776, r24777, r24778, r24779, r24780, r24781, r24782, r24783, r24784, r24785, r24786, r24787, r24788, r24789, r24790; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r24773); mpfr_init(r24774); mpfr_init(r24775); mpfr_init(r24776); mpfr_init(r24777); mpfr_init(r24778); mpfr_init_set_str(r24779, "1", 10, MPFR_RNDN); mpfr_init(r24780); mpfr_init(r24781); mpfr_init(r24782); mpfr_init(r24783); mpfr_init(r24784); mpfr_init(r24785); mpfr_init(r24786); mpfr_init(r24787); mpfr_init(r24788); mpfr_init(r24789); mpfr_init(r24790); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r24773, x, MPFR_RNDN); mpfr_set_d(r24774, y, MPFR_RNDN); mpfr_tan(r24775, r24774, MPFR_RNDN); mpfr_set_d(r24776, z, MPFR_RNDN); mpfr_tan(r24777, r24776, MPFR_RNDN); mpfr_add(r24778, r24775, r24777, MPFR_RNDN); ; mpfr_mul(r24780, r24779, r24779, MPFR_RNDN); mpfr_mul(r24781, r24775, r24777, MPFR_RNDN); mpfr_mul(r24782, r24781, r24781, MPFR_RNDN); mpfr_sub(r24783, r24780, r24782, MPFR_RNDN); mpfr_div(r24784, r24778, r24783, MPFR_RNDN); mpfr_add(r24785, r24779, r24781, MPFR_RNDN); mpfr_set_d(r24786, a, MPFR_RNDN); mpfr_tan(r24787, r24786, MPFR_RNDN); mpfr_neg(r24788, r24787, MPFR_RNDN); mpfr_fma(r24789, r24784, r24785, r24788, MPFR_RNDN); mpfr_add(r24790, r24773, r24789, MPFR_RNDN); return mpfr_get_d(r24790, MPFR_RNDN); }