#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 r25193 = x; float r25194 = y; float r25195 = z; float r25196 = r25194 + r25195; float r25197 = tan(r25196); float r25198 = a; float r25199 = tan(r25198); float r25200 = r25197 - r25199; float r25201 = r25193 + r25200; return r25201; } double f_id(double x, double y, double z, double a) { double r25202 = x; double r25203 = y; double r25204 = z; double r25205 = r25203 + r25204; double r25206 = tan(r25205); double r25207 = a; double r25208 = tan(r25207); double r25209 = r25206 - r25208; double r25210 = r25202 + r25209; return r25210; } double f_of(float x, float y, float z, float a) { float r25211 = x; float r25212 = y; float r25213 = tan(r25212); float r25214 = z; float r25215 = tan(r25214); float r25216 = r25213 + r25215; float r25217 = a; float r25218 = cos(r25217); float r25219 = r25216 * r25218; float r25220 = 1; float r25221 = r25213 * r25215; float r25222 = r25220 - r25221; float r25223 = sin(r25217); float r25224 = r25222 * r25223; float r25225 = r25219 - r25224; float r25226 = r25222 * r25218; float r25227 = r25225 / r25226; float r25228 = r25211 + r25227; return r25228; } double f_od(double x, double y, double z, double a) { double r25229 = x; double r25230 = y; double r25231 = tan(r25230); double r25232 = z; double r25233 = tan(r25232); double r25234 = r25231 + r25233; double r25235 = a; double r25236 = cos(r25235); double r25237 = r25234 * r25236; double r25238 = 1; double r25239 = r25231 * r25233; double r25240 = r25238 - r25239; double r25241 = sin(r25235); double r25242 = r25240 * r25241; double r25243 = r25237 - r25242; double r25244 = r25240 * r25236; double r25245 = r25243 / r25244; double r25246 = r25229 + r25245; return r25246; } 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 r25247, r25248, r25249, r25250, r25251, r25252, r25253, r25254, r25255; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r25247); mpfr_init(r25248); mpfr_init(r25249); mpfr_init(r25250); mpfr_init(r25251); mpfr_init(r25252); mpfr_init(r25253); mpfr_init(r25254); mpfr_init(r25255); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r25247, x, MPFR_RNDN); mpfr_set_d(r25248, y, MPFR_RNDN); mpfr_set_d(r25249, z, MPFR_RNDN); mpfr_add(r25250, r25248, r25249, MPFR_RNDN); mpfr_tan(r25251, r25250, MPFR_RNDN); mpfr_set_d(r25252, a, MPFR_RNDN); mpfr_tan(r25253, r25252, MPFR_RNDN); mpfr_sub(r25254, r25251, r25253, MPFR_RNDN); mpfr_add(r25255, r25247, r25254, MPFR_RNDN); return mpfr_get_d(r25255, MPFR_RNDN); } static mpfr_t r25256, r25257, r25258, r25259, r25260, r25261, r25262, r25263, r25264, r25265, r25266, r25267, r25268, r25269, r25270, r25271, r25272, r25273; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25256); mpfr_init(r25257); mpfr_init(r25258); mpfr_init(r25259); mpfr_init(r25260); mpfr_init(r25261); mpfr_init(r25262); mpfr_init(r25263); mpfr_init(r25264); mpfr_init_set_str(r25265, "1", 10, MPFR_RNDN); mpfr_init(r25266); mpfr_init(r25267); mpfr_init(r25268); mpfr_init(r25269); mpfr_init(r25270); mpfr_init(r25271); mpfr_init(r25272); mpfr_init(r25273); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r25256, x, MPFR_RNDN); mpfr_set_d(r25257, y, MPFR_RNDN); mpfr_tan(r25258, r25257, MPFR_RNDN); mpfr_set_d(r25259, z, MPFR_RNDN); mpfr_tan(r25260, r25259, MPFR_RNDN); mpfr_add(r25261, r25258, r25260, MPFR_RNDN); mpfr_set_d(r25262, a, MPFR_RNDN); mpfr_cos(r25263, r25262, MPFR_RNDN); mpfr_mul(r25264, r25261, r25263, MPFR_RNDN); ; mpfr_mul(r25266, r25258, r25260, MPFR_RNDN); mpfr_sub(r25267, r25265, r25266, MPFR_RNDN); mpfr_sin(r25268, r25262, MPFR_RNDN); mpfr_mul(r25269, r25267, r25268, MPFR_RNDN); mpfr_sub(r25270, r25264, r25269, MPFR_RNDN); mpfr_mul(r25271, r25267, r25263, MPFR_RNDN); mpfr_div(r25272, r25270, r25271, MPFR_RNDN); mpfr_add(r25273, r25256, r25272, MPFR_RNDN); return mpfr_get_d(r25273, MPFR_RNDN); } static mpfr_t r25274, r25275, r25276, r25277, r25278, r25279, r25280, r25281, r25282, r25283, r25284, r25285, r25286, r25287, r25288, r25289, r25290, r25291; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r25274); mpfr_init(r25275); mpfr_init(r25276); mpfr_init(r25277); mpfr_init(r25278); mpfr_init(r25279); mpfr_init(r25280); mpfr_init(r25281); mpfr_init(r25282); mpfr_init_set_str(r25283, "1", 10, MPFR_RNDN); mpfr_init(r25284); mpfr_init(r25285); mpfr_init(r25286); mpfr_init(r25287); mpfr_init(r25288); mpfr_init(r25289); mpfr_init(r25290); mpfr_init(r25291); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r25274, x, MPFR_RNDN); mpfr_set_d(r25275, y, MPFR_RNDN); mpfr_tan(r25276, r25275, MPFR_RNDN); mpfr_set_d(r25277, z, MPFR_RNDN); mpfr_tan(r25278, r25277, MPFR_RNDN); mpfr_add(r25279, r25276, r25278, MPFR_RNDN); mpfr_set_d(r25280, a, MPFR_RNDN); mpfr_cos(r25281, r25280, MPFR_RNDN); mpfr_mul(r25282, r25279, r25281, MPFR_RNDN); ; mpfr_mul(r25284, r25276, r25278, MPFR_RNDN); mpfr_sub(r25285, r25283, r25284, MPFR_RNDN); mpfr_sin(r25286, r25280, MPFR_RNDN); mpfr_mul(r25287, r25285, r25286, MPFR_RNDN); mpfr_sub(r25288, r25282, r25287, MPFR_RNDN); mpfr_mul(r25289, r25285, r25281, MPFR_RNDN); mpfr_div(r25290, r25288, r25289, MPFR_RNDN); mpfr_add(r25291, r25274, r25290, MPFR_RNDN); return mpfr_get_d(r25291, MPFR_RNDN); }