#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 r24659 = x; float r24660 = y; float r24661 = z; float r24662 = r24660 + r24661; float r24663 = tan(r24662); float r24664 = a; float r24665 = tan(r24664); float r24666 = r24663 - r24665; float r24667 = r24659 + r24666; return r24667; } double f_id(double x, double y, double z, double a) { double r24668 = x; double r24669 = y; double r24670 = z; double r24671 = r24669 + r24670; double r24672 = tan(r24671); double r24673 = a; double r24674 = tan(r24673); double r24675 = r24672 - r24674; double r24676 = r24668 + r24675; return r24676; } double f_of(float x, float y, float z, float a) { float r24677 = x; float r24678 = y; float r24679 = tan(r24678); float r24680 = z; float r24681 = tan(r24680); float r24682 = r24679 + r24681; float r24683 = r24679 * r24681; float r24684 = fma(r24683, r24683, r24683); float r24685 = fma(r24682, r24684, r24682); float r24686 = 1; float r24687 = 3; float r24688 = pow(r24683, r24687); float r24689 = r24686 - r24688; float r24690 = r24685 / r24689; float r24691 = a; float r24692 = tan(r24691); float r24693 = r24690 - r24692; float r24694 = 0; float r24695 = r24693 + r24694; float r24696 = r24677 + r24695; return r24696; } double f_od(double x, double y, double z, double a) { double r24697 = x; double r24698 = y; double r24699 = tan(r24698); double r24700 = z; double r24701 = tan(r24700); double r24702 = r24699 + r24701; double r24703 = r24699 * r24701; double r24704 = fma(r24703, r24703, r24703); double r24705 = fma(r24702, r24704, r24702); double r24706 = 1; double r24707 = 3; double r24708 = pow(r24703, r24707); double r24709 = r24706 - r24708; double r24710 = r24705 / r24709; double r24711 = a; double r24712 = tan(r24711); double r24713 = r24710 - r24712; double r24714 = 0; double r24715 = r24713 + r24714; double r24716 = r24697 + r24715; return r24716; } 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 r24717, r24718, r24719, r24720, r24721, r24722, r24723, r24724, r24725; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r24717); mpfr_init(r24718); mpfr_init(r24719); mpfr_init(r24720); mpfr_init(r24721); mpfr_init(r24722); mpfr_init(r24723); mpfr_init(r24724); mpfr_init(r24725); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r24717, x, MPFR_RNDN); mpfr_set_d(r24718, y, MPFR_RNDN); mpfr_set_d(r24719, z, MPFR_RNDN); mpfr_add(r24720, r24718, r24719, MPFR_RNDN); mpfr_tan(r24721, r24720, MPFR_RNDN); mpfr_set_d(r24722, a, MPFR_RNDN); mpfr_tan(r24723, r24722, MPFR_RNDN); mpfr_sub(r24724, r24721, r24723, MPFR_RNDN); mpfr_add(r24725, r24717, r24724, MPFR_RNDN); return mpfr_get_d(r24725, MPFR_RNDN); } static mpfr_t r24726, r24727, r24728, r24729, r24730, r24731, r24732, r24733, r24734, r24735, r24736, r24737, r24738, r24739, r24740, r24741, r24742, r24743, r24744, r24745; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24726); mpfr_init(r24727); mpfr_init(r24728); mpfr_init(r24729); mpfr_init(r24730); mpfr_init(r24731); mpfr_init(r24732); mpfr_init(r24733); mpfr_init(r24734); mpfr_init_set_str(r24735, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24736, "3", 10, MPFR_RNDN); mpfr_init(r24737); mpfr_init(r24738); mpfr_init(r24739); mpfr_init(r24740); mpfr_init(r24741); mpfr_init(r24742); mpfr_init_set_str(r24743, "0", 10, MPFR_RNDN); mpfr_init(r24744); mpfr_init(r24745); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r24726, x, MPFR_RNDN); mpfr_set_d(r24727, y, MPFR_RNDN); mpfr_tan(r24728, r24727, MPFR_RNDN); mpfr_set_d(r24729, z, MPFR_RNDN); mpfr_tan(r24730, r24729, MPFR_RNDN); mpfr_add(r24731, r24728, r24730, MPFR_RNDN); mpfr_mul(r24732, r24728, r24730, MPFR_RNDN); mpfr_fma(r24733, r24732, r24732, r24732, MPFR_RNDN); mpfr_fma(r24734, r24731, r24733, r24731, MPFR_RNDN); ; ; mpfr_pow(r24737, r24732, r24736, MPFR_RNDN); mpfr_sub(r24738, r24735, r24737, MPFR_RNDN); mpfr_div(r24739, r24734, r24738, MPFR_RNDN); mpfr_set_d(r24740, a, MPFR_RNDN); mpfr_tan(r24741, r24740, MPFR_RNDN); mpfr_sub(r24742, r24739, r24741, MPFR_RNDN); ; mpfr_add(r24744, r24742, r24743, MPFR_RNDN); mpfr_add(r24745, r24726, r24744, MPFR_RNDN); return mpfr_get_d(r24745, MPFR_RNDN); } static mpfr_t r24746, r24747, r24748, r24749, r24750, r24751, r24752, r24753, r24754, r24755, r24756, r24757, r24758, r24759, r24760, r24761, r24762, r24763, r24764, r24765; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); 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); mpfr_init_set_str(r24755, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24756, "3", 10, MPFR_RNDN); mpfr_init(r24757); mpfr_init(r24758); mpfr_init(r24759); mpfr_init(r24760); mpfr_init(r24761); mpfr_init(r24762); mpfr_init_set_str(r24763, "0", 10, MPFR_RNDN); mpfr_init(r24764); mpfr_init(r24765); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r24746, x, MPFR_RNDN); mpfr_set_d(r24747, y, MPFR_RNDN); mpfr_tan(r24748, r24747, MPFR_RNDN); mpfr_set_d(r24749, z, MPFR_RNDN); mpfr_tan(r24750, r24749, MPFR_RNDN); mpfr_add(r24751, r24748, r24750, MPFR_RNDN); mpfr_mul(r24752, r24748, r24750, MPFR_RNDN); mpfr_fma(r24753, r24752, r24752, r24752, MPFR_RNDN); mpfr_fma(r24754, r24751, r24753, r24751, MPFR_RNDN); ; ; mpfr_pow(r24757, r24752, r24756, MPFR_RNDN); mpfr_sub(r24758, r24755, r24757, MPFR_RNDN); mpfr_div(r24759, r24754, r24758, MPFR_RNDN); mpfr_set_d(r24760, a, MPFR_RNDN); mpfr_tan(r24761, r24760, MPFR_RNDN); mpfr_sub(r24762, r24759, r24761, MPFR_RNDN); ; mpfr_add(r24764, r24762, r24763, MPFR_RNDN); mpfr_add(r24765, r24746, r24764, MPFR_RNDN); return mpfr_get_d(r24765, MPFR_RNDN); }