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