#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 r25657 = x; float r25658 = y; float r25659 = z; float r25660 = r25658 + r25659; float r25661 = tan(r25660); float r25662 = a; float r25663 = tan(r25662); float r25664 = r25661 - r25663; float r25665 = r25657 + r25664; return r25665; } double f_id(double x, double y, double z, double a) { double r25666 = x; double r25667 = y; double r25668 = z; double r25669 = r25667 + r25668; double r25670 = tan(r25669); double r25671 = a; double r25672 = tan(r25671); double r25673 = r25670 - r25672; double r25674 = r25666 + r25673; return r25674; } double f_of(float x, float y, float z, float a) { float r25675 = x; float r25676 = y; float r25677 = tan(r25676); float r25678 = z; float r25679 = tan(r25678); float r25680 = r25677 + r25679; float r25681 = a; float r25682 = cos(r25681); float r25683 = r25680 * r25682; float r25684 = 1; float r25685 = r25677 * r25679; float r25686 = r25684 - r25685; float r25687 = sin(r25681); float r25688 = r25686 * r25687; float r25689 = r25683 - r25688; float r25690 = r25686 * r25682; float r25691 = r25689 / r25690; float r25692 = r25675 + r25691; return r25692; } double f_od(double x, double y, double z, double a) { double r25693 = x; double r25694 = y; double r25695 = tan(r25694); double r25696 = z; double r25697 = tan(r25696); double r25698 = r25695 + r25697; double r25699 = a; double r25700 = cos(r25699); double r25701 = r25698 * r25700; double r25702 = 1; double r25703 = r25695 * r25697; double r25704 = r25702 - r25703; double r25705 = sin(r25699); double r25706 = r25704 * r25705; double r25707 = r25701 - r25706; double r25708 = r25704 * r25700; double r25709 = r25707 / r25708; double r25710 = r25693 + r25709; return r25710; } 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 r25711, r25712, r25713, r25714, r25715, r25716, r25717, r25718, r25719; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r25711); mpfr_init(r25712); mpfr_init(r25713); mpfr_init(r25714); mpfr_init(r25715); mpfr_init(r25716); mpfr_init(r25717); mpfr_init(r25718); mpfr_init(r25719); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r25711, x, MPFR_RNDN); mpfr_set_d(r25712, y, MPFR_RNDN); mpfr_set_d(r25713, z, MPFR_RNDN); mpfr_add(r25714, r25712, r25713, MPFR_RNDN); mpfr_tan(r25715, r25714, MPFR_RNDN); mpfr_set_d(r25716, a, MPFR_RNDN); mpfr_tan(r25717, r25716, MPFR_RNDN); mpfr_sub(r25718, r25715, r25717, MPFR_RNDN); mpfr_add(r25719, r25711, r25718, MPFR_RNDN); return mpfr_get_d(r25719, MPFR_RNDN); } static mpfr_t r25720, r25721, r25722, r25723, r25724, r25725, r25726, r25727, r25728, r25729, r25730, r25731, r25732, r25733, r25734, r25735, r25736, r25737; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r25720); mpfr_init(r25721); mpfr_init(r25722); mpfr_init(r25723); mpfr_init(r25724); mpfr_init(r25725); mpfr_init(r25726); mpfr_init(r25727); mpfr_init(r25728); mpfr_init_set_str(r25729, "1", 10, MPFR_RNDN); mpfr_init(r25730); mpfr_init(r25731); mpfr_init(r25732); mpfr_init(r25733); mpfr_init(r25734); mpfr_init(r25735); mpfr_init(r25736); mpfr_init(r25737); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r25720, x, MPFR_RNDN); mpfr_set_d(r25721, y, MPFR_RNDN); mpfr_tan(r25722, r25721, MPFR_RNDN); mpfr_set_d(r25723, z, MPFR_RNDN); mpfr_tan(r25724, r25723, MPFR_RNDN); mpfr_add(r25725, r25722, r25724, MPFR_RNDN); mpfr_set_d(r25726, a, MPFR_RNDN); mpfr_cos(r25727, r25726, MPFR_RNDN); mpfr_mul(r25728, r25725, r25727, MPFR_RNDN); ; mpfr_mul(r25730, r25722, r25724, MPFR_RNDN); mpfr_sub(r25731, r25729, r25730, MPFR_RNDN); mpfr_sin(r25732, r25726, MPFR_RNDN); mpfr_mul(r25733, r25731, r25732, MPFR_RNDN); mpfr_sub(r25734, r25728, r25733, MPFR_RNDN); mpfr_mul(r25735, r25731, r25727, MPFR_RNDN); mpfr_div(r25736, r25734, r25735, MPFR_RNDN); mpfr_add(r25737, r25720, r25736, MPFR_RNDN); return mpfr_get_d(r25737, MPFR_RNDN); } static mpfr_t r25738, r25739, r25740, r25741, r25742, r25743, r25744, r25745, r25746, r25747, r25748, r25749, r25750, r25751, r25752, r25753, r25754, r25755; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r25738); mpfr_init(r25739); mpfr_init(r25740); mpfr_init(r25741); mpfr_init(r25742); mpfr_init(r25743); mpfr_init(r25744); mpfr_init(r25745); mpfr_init(r25746); mpfr_init_set_str(r25747, "1", 10, MPFR_RNDN); mpfr_init(r25748); mpfr_init(r25749); mpfr_init(r25750); mpfr_init(r25751); mpfr_init(r25752); mpfr_init(r25753); mpfr_init(r25754); mpfr_init(r25755); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r25738, x, MPFR_RNDN); mpfr_set_d(r25739, y, MPFR_RNDN); mpfr_tan(r25740, r25739, MPFR_RNDN); mpfr_set_d(r25741, z, MPFR_RNDN); mpfr_tan(r25742, r25741, MPFR_RNDN); mpfr_add(r25743, r25740, r25742, MPFR_RNDN); mpfr_set_d(r25744, a, MPFR_RNDN); mpfr_cos(r25745, r25744, MPFR_RNDN); mpfr_mul(r25746, r25743, r25745, MPFR_RNDN); ; mpfr_mul(r25748, r25740, r25742, MPFR_RNDN); mpfr_sub(r25749, r25747, r25748, MPFR_RNDN); mpfr_sin(r25750, r25744, MPFR_RNDN); mpfr_mul(r25751, r25749, r25750, MPFR_RNDN); mpfr_sub(r25752, r25746, r25751, MPFR_RNDN); mpfr_mul(r25753, r25749, r25745, MPFR_RNDN); mpfr_div(r25754, r25752, r25753, MPFR_RNDN); mpfr_add(r25755, r25738, r25754, MPFR_RNDN); return mpfr_get_d(r25755, MPFR_RNDN); }