#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 r24703 = x; float r24704 = y; float r24705 = z; float r24706 = r24704 + r24705; float r24707 = tan(r24706); float r24708 = a; float r24709 = tan(r24708); float r24710 = r24707 - r24709; float r24711 = r24703 + r24710; return r24711; } double f_id(double x, double y, double z, double a) { double r24712 = x; double r24713 = y; double r24714 = z; double r24715 = r24713 + r24714; double r24716 = tan(r24715); double r24717 = a; double r24718 = tan(r24717); double r24719 = r24716 - r24718; double r24720 = r24712 + r24719; return r24720; } double f_of(float x, float y, float z, float a) { float r24721 = x; float r24722 = z; float r24723 = tan(r24722); float r24724 = y; float r24725 = tan(r24724); float r24726 = 1; float r24727 = fma(r24723, r24725, r24726); float r24728 = r24723 * r24725; float r24729 = fma(r24727, r24728, r24726); float r24730 = r24723 + r24725; float r24731 = r24729 * r24730; float r24732 = 3; float r24733 = pow(r24728, r24732); float r24734 = r24726 - r24733; float r24735 = r24731 / r24734; float r24736 = a; float r24737 = tan(r24736); float r24738 = r24735 - r24737; float r24739 = r24737 - r24737; float r24740 = r24738 + r24739; float r24741 = r24721 + r24740; return r24741; } double f_od(double x, double y, double z, double a) { double r24742 = x; double r24743 = z; double r24744 = tan(r24743); double r24745 = y; double r24746 = tan(r24745); double r24747 = 1; double r24748 = fma(r24744, r24746, r24747); double r24749 = r24744 * r24746; double r24750 = fma(r24748, r24749, r24747); double r24751 = r24744 + r24746; double r24752 = r24750 * r24751; double r24753 = 3; double r24754 = pow(r24749, r24753); double r24755 = r24747 - r24754; double r24756 = r24752 / r24755; double r24757 = a; double r24758 = tan(r24757); double r24759 = r24756 - r24758; double r24760 = r24758 - r24758; double r24761 = r24759 + r24760; double r24762 = r24742 + r24761; return r24762; } 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 r24763, r24764, r24765, r24766, r24767, r24768, r24769, r24770, r24771; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); 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); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r24763, x, MPFR_RNDN); mpfr_set_d(r24764, y, MPFR_RNDN); mpfr_set_d(r24765, z, MPFR_RNDN); mpfr_add(r24766, r24764, r24765, MPFR_RNDN); mpfr_tan(r24767, r24766, MPFR_RNDN); mpfr_set_d(r24768, a, MPFR_RNDN); mpfr_tan(r24769, r24768, MPFR_RNDN); mpfr_sub(r24770, r24767, r24769, MPFR_RNDN); mpfr_add(r24771, r24763, r24770, MPFR_RNDN); return mpfr_get_d(r24771, MPFR_RNDN); } static mpfr_t r24772, r24773, r24774, r24775, r24776, r24777, r24778, r24779, r24780, r24781, r24782, r24783, r24784, r24785, r24786, r24787, r24788, r24789, r24790, r24791, r24792; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24772); mpfr_init(r24773); mpfr_init(r24774); mpfr_init(r24775); mpfr_init(r24776); mpfr_init_set_str(r24777, "1", 10, MPFR_RNDN); mpfr_init(r24778); mpfr_init(r24779); mpfr_init(r24780); mpfr_init(r24781); mpfr_init(r24782); mpfr_init_set_str(r24783, "3", 10, MPFR_RNDN); mpfr_init(r24784); mpfr_init(r24785); mpfr_init(r24786); mpfr_init(r24787); mpfr_init(r24788); mpfr_init(r24789); mpfr_init(r24790); mpfr_init(r24791); mpfr_init(r24792); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r24772, x, MPFR_RNDN); mpfr_set_d(r24773, z, MPFR_RNDN); mpfr_tan(r24774, r24773, MPFR_RNDN); mpfr_set_d(r24775, y, MPFR_RNDN); mpfr_tan(r24776, r24775, MPFR_RNDN); ; mpfr_fma(r24778, r24774, r24776, r24777, MPFR_RNDN); mpfr_mul(r24779, r24774, r24776, MPFR_RNDN); mpfr_fma(r24780, r24778, r24779, r24777, MPFR_RNDN); mpfr_add(r24781, r24774, r24776, MPFR_RNDN); mpfr_mul(r24782, r24780, r24781, MPFR_RNDN); ; mpfr_pow(r24784, r24779, r24783, MPFR_RNDN); mpfr_sub(r24785, r24777, r24784, MPFR_RNDN); mpfr_div(r24786, r24782, r24785, MPFR_RNDN); mpfr_set_d(r24787, a, MPFR_RNDN); mpfr_tan(r24788, r24787, MPFR_RNDN); mpfr_sub(r24789, r24786, r24788, MPFR_RNDN); mpfr_sub(r24790, r24788, r24788, MPFR_RNDN); mpfr_add(r24791, r24789, r24790, MPFR_RNDN); mpfr_add(r24792, r24772, r24791, MPFR_RNDN); return mpfr_get_d(r24792, MPFR_RNDN); } static mpfr_t r24793, r24794, r24795, r24796, r24797, r24798, r24799, r24800, r24801, r24802, r24803, r24804, r24805, r24806, r24807, r24808, r24809, r24810, r24811, r24812, r24813; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24793); mpfr_init(r24794); mpfr_init(r24795); mpfr_init(r24796); mpfr_init(r24797); mpfr_init_set_str(r24798, "1", 10, MPFR_RNDN); mpfr_init(r24799); mpfr_init(r24800); mpfr_init(r24801); mpfr_init(r24802); mpfr_init(r24803); mpfr_init_set_str(r24804, "3", 10, MPFR_RNDN); mpfr_init(r24805); mpfr_init(r24806); mpfr_init(r24807); mpfr_init(r24808); mpfr_init(r24809); mpfr_init(r24810); mpfr_init(r24811); mpfr_init(r24812); mpfr_init(r24813); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r24793, x, MPFR_RNDN); mpfr_set_d(r24794, z, MPFR_RNDN); mpfr_tan(r24795, r24794, MPFR_RNDN); mpfr_set_d(r24796, y, MPFR_RNDN); mpfr_tan(r24797, r24796, MPFR_RNDN); ; mpfr_fma(r24799, r24795, r24797, r24798, MPFR_RNDN); mpfr_mul(r24800, r24795, r24797, MPFR_RNDN); mpfr_fma(r24801, r24799, r24800, r24798, MPFR_RNDN); mpfr_add(r24802, r24795, r24797, MPFR_RNDN); mpfr_mul(r24803, r24801, r24802, MPFR_RNDN); ; mpfr_pow(r24805, r24800, r24804, MPFR_RNDN); mpfr_sub(r24806, r24798, r24805, MPFR_RNDN); mpfr_div(r24807, r24803, r24806, MPFR_RNDN); mpfr_set_d(r24808, a, MPFR_RNDN); mpfr_tan(r24809, r24808, MPFR_RNDN); mpfr_sub(r24810, r24807, r24809, MPFR_RNDN); mpfr_sub(r24811, r24809, r24809, MPFR_RNDN); mpfr_add(r24812, r24810, r24811, MPFR_RNDN); mpfr_add(r24813, r24793, r24812, MPFR_RNDN); return mpfr_get_d(r24813, MPFR_RNDN); }