#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 r27676 = x; float r27677 = y; float r27678 = z; float r27679 = r27677 + r27678; float r27680 = tan(r27679); float r27681 = a; float r27682 = tan(r27681); float r27683 = r27680 - r27682; float r27684 = r27676 + r27683; return r27684; } double f_id(double x, double y, double z, double a) { double r27685 = x; double r27686 = y; double r27687 = z; double r27688 = r27686 + r27687; double r27689 = tan(r27688); double r27690 = a; double r27691 = tan(r27690); double r27692 = r27689 - r27691; double r27693 = r27685 + r27692; return r27693; } double f_of(float x, float y, float z, float a) { float r27694 = x; float r27695 = y; float r27696 = tan(r27695); float r27697 = z; float r27698 = tan(r27697); float r27699 = r27696 + r27698; float r27700 = a; float r27701 = cos(r27700); float r27702 = r27699 * r27701; float r27703 = 1; float r27704 = r27696 * r27698; float r27705 = r27703 - r27704; float r27706 = sin(r27700); float r27707 = r27705 * r27706; float r27708 = r27702 - r27707; float r27709 = r27705 * r27701; float r27710 = r27708 / r27709; float r27711 = r27694 + r27710; return r27711; } double f_od(double x, double y, double z, double a) { double r27712 = x; double r27713 = y; double r27714 = tan(r27713); double r27715 = z; double r27716 = tan(r27715); double r27717 = r27714 + r27716; double r27718 = a; double r27719 = cos(r27718); double r27720 = r27717 * r27719; double r27721 = 1; double r27722 = r27714 * r27716; double r27723 = r27721 - r27722; double r27724 = sin(r27718); double r27725 = r27723 * r27724; double r27726 = r27720 - r27725; double r27727 = r27723 * r27719; double r27728 = r27726 / r27727; double r27729 = r27712 + r27728; return r27729; } 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 r27730, r27731, r27732, r27733, r27734, r27735, r27736, r27737, r27738; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r27730); mpfr_init(r27731); mpfr_init(r27732); mpfr_init(r27733); mpfr_init(r27734); mpfr_init(r27735); mpfr_init(r27736); mpfr_init(r27737); mpfr_init(r27738); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r27730, x, MPFR_RNDN); mpfr_set_d(r27731, y, MPFR_RNDN); mpfr_set_d(r27732, z, MPFR_RNDN); mpfr_add(r27733, r27731, r27732, MPFR_RNDN); mpfr_tan(r27734, r27733, MPFR_RNDN); mpfr_set_d(r27735, a, MPFR_RNDN); mpfr_tan(r27736, r27735, MPFR_RNDN); mpfr_sub(r27737, r27734, r27736, MPFR_RNDN); mpfr_add(r27738, r27730, r27737, MPFR_RNDN); return mpfr_get_d(r27738, MPFR_RNDN); } static mpfr_t r27739, r27740, r27741, r27742, r27743, r27744, r27745, r27746, r27747, r27748, r27749, r27750, r27751, r27752, r27753, r27754, r27755, r27756; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r27739); mpfr_init(r27740); mpfr_init(r27741); mpfr_init(r27742); mpfr_init(r27743); mpfr_init(r27744); mpfr_init(r27745); mpfr_init(r27746); mpfr_init(r27747); mpfr_init_set_str(r27748, "1", 10, MPFR_RNDN); mpfr_init(r27749); mpfr_init(r27750); mpfr_init(r27751); mpfr_init(r27752); mpfr_init(r27753); mpfr_init(r27754); mpfr_init(r27755); mpfr_init(r27756); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r27739, x, MPFR_RNDN); mpfr_set_d(r27740, y, MPFR_RNDN); mpfr_tan(r27741, r27740, MPFR_RNDN); mpfr_set_d(r27742, z, MPFR_RNDN); mpfr_tan(r27743, r27742, MPFR_RNDN); mpfr_add(r27744, r27741, r27743, MPFR_RNDN); mpfr_set_d(r27745, a, MPFR_RNDN); mpfr_cos(r27746, r27745, MPFR_RNDN); mpfr_mul(r27747, r27744, r27746, MPFR_RNDN); ; mpfr_mul(r27749, r27741, r27743, MPFR_RNDN); mpfr_sub(r27750, r27748, r27749, MPFR_RNDN); mpfr_sin(r27751, r27745, MPFR_RNDN); mpfr_mul(r27752, r27750, r27751, MPFR_RNDN); mpfr_sub(r27753, r27747, r27752, MPFR_RNDN); mpfr_mul(r27754, r27750, r27746, MPFR_RNDN); mpfr_div(r27755, r27753, r27754, MPFR_RNDN); mpfr_add(r27756, r27739, r27755, MPFR_RNDN); return mpfr_get_d(r27756, MPFR_RNDN); } static mpfr_t r27757, r27758, r27759, r27760, r27761, r27762, r27763, r27764, r27765, r27766, r27767, r27768, r27769, r27770, r27771, r27772, r27773, r27774; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r27757); mpfr_init(r27758); mpfr_init(r27759); mpfr_init(r27760); mpfr_init(r27761); mpfr_init(r27762); mpfr_init(r27763); mpfr_init(r27764); mpfr_init(r27765); mpfr_init_set_str(r27766, "1", 10, MPFR_RNDN); mpfr_init(r27767); mpfr_init(r27768); mpfr_init(r27769); mpfr_init(r27770); mpfr_init(r27771); mpfr_init(r27772); mpfr_init(r27773); mpfr_init(r27774); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r27757, x, MPFR_RNDN); mpfr_set_d(r27758, y, MPFR_RNDN); mpfr_tan(r27759, r27758, MPFR_RNDN); mpfr_set_d(r27760, z, MPFR_RNDN); mpfr_tan(r27761, r27760, MPFR_RNDN); mpfr_add(r27762, r27759, r27761, MPFR_RNDN); mpfr_set_d(r27763, a, MPFR_RNDN); mpfr_cos(r27764, r27763, MPFR_RNDN); mpfr_mul(r27765, r27762, r27764, MPFR_RNDN); ; mpfr_mul(r27767, r27759, r27761, MPFR_RNDN); mpfr_sub(r27768, r27766, r27767, MPFR_RNDN); mpfr_sin(r27769, r27763, MPFR_RNDN); mpfr_mul(r27770, r27768, r27769, MPFR_RNDN); mpfr_sub(r27771, r27765, r27770, MPFR_RNDN); mpfr_mul(r27772, r27768, r27764, MPFR_RNDN); mpfr_div(r27773, r27771, r27772, MPFR_RNDN); mpfr_add(r27774, r27757, r27773, MPFR_RNDN); return mpfr_get_d(r27774, MPFR_RNDN); }