#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 r27060 = x; float r27061 = y; float r27062 = z; float r27063 = r27061 + r27062; float r27064 = tan(r27063); float r27065 = a; float r27066 = tan(r27065); float r27067 = r27064 - r27066; float r27068 = r27060 + r27067; return r27068; } double f_id(double x, double y, double z, double a) { double r27069 = x; double r27070 = y; double r27071 = z; double r27072 = r27070 + r27071; double r27073 = tan(r27072); double r27074 = a; double r27075 = tan(r27074); double r27076 = r27073 - r27075; double r27077 = r27069 + r27076; return r27077; } double f_of(float x, float y, float z, float a) { float r27078 = x; float r27079 = y; float r27080 = sin(r27079); float r27081 = z; float r27082 = cos(r27081); float r27083 = r27080 * r27082; float r27084 = cos(r27079); float r27085 = sin(r27081); float r27086 = r27084 * r27085; float r27087 = r27083 + r27086; float r27088 = 1; float r27089 = tan(r27079); float r27090 = tan(r27081); float r27091 = r27089 * r27090; float r27092 = r27088 - r27091; float r27093 = r27084 * r27082; float r27094 = r27092 * r27093; float r27095 = r27087 / r27094; float r27096 = a; float r27097 = tan(r27096); float r27098 = r27095 - r27097; float r27099 = r27078 + r27098; return r27099; } double f_od(double x, double y, double z, double a) { double r27100 = x; double r27101 = y; double r27102 = sin(r27101); double r27103 = z; double r27104 = cos(r27103); double r27105 = r27102 * r27104; double r27106 = cos(r27101); double r27107 = sin(r27103); double r27108 = r27106 * r27107; double r27109 = r27105 + r27108; double r27110 = 1; double r27111 = tan(r27101); double r27112 = tan(r27103); double r27113 = r27111 * r27112; double r27114 = r27110 - r27113; double r27115 = r27106 * r27104; double r27116 = r27114 * r27115; double r27117 = r27109 / r27116; double r27118 = a; double r27119 = tan(r27118); double r27120 = r27117 - r27119; double r27121 = r27100 + r27120; return r27121; } 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 r27122, r27123, r27124, r27125, r27126, r27127, r27128, r27129, r27130; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r27122); mpfr_init(r27123); mpfr_init(r27124); mpfr_init(r27125); mpfr_init(r27126); mpfr_init(r27127); mpfr_init(r27128); mpfr_init(r27129); mpfr_init(r27130); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r27122, x, MPFR_RNDN); mpfr_set_d(r27123, y, MPFR_RNDN); mpfr_set_d(r27124, z, MPFR_RNDN); mpfr_add(r27125, r27123, r27124, MPFR_RNDN); mpfr_tan(r27126, r27125, MPFR_RNDN); mpfr_set_d(r27127, a, MPFR_RNDN); mpfr_tan(r27128, r27127, MPFR_RNDN); mpfr_sub(r27129, r27126, r27128, MPFR_RNDN); mpfr_add(r27130, r27122, r27129, MPFR_RNDN); return mpfr_get_d(r27130, MPFR_RNDN); } static mpfr_t r27131, r27132, r27133, r27134, r27135, r27136, r27137, r27138, r27139, r27140, r27141, r27142, r27143, r27144, r27145, r27146, r27147, r27148, r27149, r27150, r27151, r27152; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r27131); mpfr_init(r27132); mpfr_init(r27133); mpfr_init(r27134); mpfr_init(r27135); mpfr_init(r27136); mpfr_init(r27137); mpfr_init(r27138); mpfr_init(r27139); mpfr_init(r27140); mpfr_init_set_str(r27141, "1", 10, MPFR_RNDN); mpfr_init(r27142); mpfr_init(r27143); mpfr_init(r27144); mpfr_init(r27145); mpfr_init(r27146); mpfr_init(r27147); mpfr_init(r27148); mpfr_init(r27149); mpfr_init(r27150); mpfr_init(r27151); mpfr_init(r27152); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r27131, x, MPFR_RNDN); mpfr_set_d(r27132, y, MPFR_RNDN); mpfr_sin(r27133, r27132, MPFR_RNDN); mpfr_set_d(r27134, z, MPFR_RNDN); mpfr_cos(r27135, r27134, MPFR_RNDN); mpfr_mul(r27136, r27133, r27135, MPFR_RNDN); mpfr_cos(r27137, r27132, MPFR_RNDN); mpfr_sin(r27138, r27134, MPFR_RNDN); mpfr_mul(r27139, r27137, r27138, MPFR_RNDN); mpfr_add(r27140, r27136, r27139, MPFR_RNDN); ; mpfr_tan(r27142, r27132, MPFR_RNDN); mpfr_tan(r27143, r27134, MPFR_RNDN); mpfr_mul(r27144, r27142, r27143, MPFR_RNDN); mpfr_sub(r27145, r27141, r27144, MPFR_RNDN); mpfr_mul(r27146, r27137, r27135, MPFR_RNDN); mpfr_mul(r27147, r27145, r27146, MPFR_RNDN); mpfr_div(r27148, r27140, r27147, MPFR_RNDN); mpfr_set_d(r27149, a, MPFR_RNDN); mpfr_tan(r27150, r27149, MPFR_RNDN); mpfr_sub(r27151, r27148, r27150, MPFR_RNDN); mpfr_add(r27152, r27131, r27151, MPFR_RNDN); return mpfr_get_d(r27152, MPFR_RNDN); } static mpfr_t r27153, r27154, r27155, r27156, r27157, r27158, r27159, r27160, r27161, r27162, r27163, r27164, r27165, r27166, r27167, r27168, r27169, r27170, r27171, r27172, r27173, r27174; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r27153); mpfr_init(r27154); mpfr_init(r27155); mpfr_init(r27156); mpfr_init(r27157); mpfr_init(r27158); mpfr_init(r27159); mpfr_init(r27160); mpfr_init(r27161); mpfr_init(r27162); mpfr_init_set_str(r27163, "1", 10, MPFR_RNDN); mpfr_init(r27164); mpfr_init(r27165); mpfr_init(r27166); mpfr_init(r27167); mpfr_init(r27168); mpfr_init(r27169); mpfr_init(r27170); mpfr_init(r27171); mpfr_init(r27172); mpfr_init(r27173); mpfr_init(r27174); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r27153, x, MPFR_RNDN); mpfr_set_d(r27154, y, MPFR_RNDN); mpfr_sin(r27155, r27154, MPFR_RNDN); mpfr_set_d(r27156, z, MPFR_RNDN); mpfr_cos(r27157, r27156, MPFR_RNDN); mpfr_mul(r27158, r27155, r27157, MPFR_RNDN); mpfr_cos(r27159, r27154, MPFR_RNDN); mpfr_sin(r27160, r27156, MPFR_RNDN); mpfr_mul(r27161, r27159, r27160, MPFR_RNDN); mpfr_add(r27162, r27158, r27161, MPFR_RNDN); ; mpfr_tan(r27164, r27154, MPFR_RNDN); mpfr_tan(r27165, r27156, MPFR_RNDN); mpfr_mul(r27166, r27164, r27165, MPFR_RNDN); mpfr_sub(r27167, r27163, r27166, MPFR_RNDN); mpfr_mul(r27168, r27159, r27157, MPFR_RNDN); mpfr_mul(r27169, r27167, r27168, MPFR_RNDN); mpfr_div(r27170, r27162, r27169, MPFR_RNDN); mpfr_set_d(r27171, a, MPFR_RNDN); mpfr_tan(r27172, r27171, MPFR_RNDN); mpfr_sub(r27173, r27170, r27172, MPFR_RNDN); mpfr_add(r27174, r27153, r27173, MPFR_RNDN); return mpfr_get_d(r27174, MPFR_RNDN); }