#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 r27199 = x; float r27200 = y; float r27201 = z; float r27202 = r27200 + r27201; float r27203 = tan(r27202); float r27204 = a; float r27205 = tan(r27204); float r27206 = r27203 - r27205; float r27207 = r27199 + r27206; return r27207; } double f_id(double x, double y, double z, double a) { double r27208 = x; double r27209 = y; double r27210 = z; double r27211 = r27209 + r27210; double r27212 = tan(r27211); double r27213 = a; double r27214 = tan(r27213); double r27215 = r27212 - r27214; double r27216 = r27208 + r27215; return r27216; } double f_of(float x, float y, float z, float a) { float r27217 = x; float r27218 = y; float r27219 = tan(r27218); float r27220 = z; float r27221 = tan(r27220); float r27222 = r27219 + r27221; float r27223 = a; float r27224 = cos(r27223); float r27225 = r27222 * r27224; float r27226 = 1; float r27227 = r27219 * r27221; float r27228 = r27226 - r27227; float r27229 = sin(r27223); float r27230 = r27228 * r27229; float r27231 = r27225 - r27230; float r27232 = r27228 * r27224; float r27233 = r27231 / r27232; float r27234 = r27217 + r27233; return r27234; } double f_od(double x, double y, double z, double a) { double r27235 = x; double r27236 = y; double r27237 = tan(r27236); double r27238 = z; double r27239 = tan(r27238); double r27240 = r27237 + r27239; double r27241 = a; double r27242 = cos(r27241); double r27243 = r27240 * r27242; double r27244 = 1; double r27245 = r27237 * r27239; double r27246 = r27244 - r27245; double r27247 = sin(r27241); double r27248 = r27246 * r27247; double r27249 = r27243 - r27248; double r27250 = r27246 * r27242; double r27251 = r27249 / r27250; double r27252 = r27235 + r27251; return r27252; } 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 r27253, r27254, r27255, r27256, r27257, r27258, r27259, r27260, r27261; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r27253); mpfr_init(r27254); mpfr_init(r27255); mpfr_init(r27256); mpfr_init(r27257); mpfr_init(r27258); mpfr_init(r27259); mpfr_init(r27260); mpfr_init(r27261); } double f_im(double x, double y, double z, double a) { mpfr_set_d(r27253, x, MPFR_RNDN); mpfr_set_d(r27254, y, MPFR_RNDN); mpfr_set_d(r27255, z, MPFR_RNDN); mpfr_add(r27256, r27254, r27255, MPFR_RNDN); mpfr_tan(r27257, r27256, MPFR_RNDN); mpfr_set_d(r27258, a, MPFR_RNDN); mpfr_tan(r27259, r27258, MPFR_RNDN); mpfr_sub(r27260, r27257, r27259, MPFR_RNDN); mpfr_add(r27261, r27253, r27260, MPFR_RNDN); return mpfr_get_d(r27261, MPFR_RNDN); } static mpfr_t r27262, r27263, r27264, r27265, r27266, r27267, r27268, r27269, r27270, r27271, r27272, r27273, r27274, r27275, r27276, r27277, r27278, r27279; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r27262); mpfr_init(r27263); mpfr_init(r27264); mpfr_init(r27265); mpfr_init(r27266); mpfr_init(r27267); mpfr_init(r27268); mpfr_init(r27269); mpfr_init(r27270); mpfr_init_set_str(r27271, "1", 10, MPFR_RNDN); mpfr_init(r27272); mpfr_init(r27273); mpfr_init(r27274); mpfr_init(r27275); mpfr_init(r27276); mpfr_init(r27277); mpfr_init(r27278); mpfr_init(r27279); } double f_fm(double x, double y, double z, double a) { mpfr_set_d(r27262, x, MPFR_RNDN); mpfr_set_d(r27263, y, MPFR_RNDN); mpfr_tan(r27264, r27263, MPFR_RNDN); mpfr_set_d(r27265, z, MPFR_RNDN); mpfr_tan(r27266, r27265, MPFR_RNDN); mpfr_add(r27267, r27264, r27266, MPFR_RNDN); mpfr_set_d(r27268, a, MPFR_RNDN); mpfr_cos(r27269, r27268, MPFR_RNDN); mpfr_mul(r27270, r27267, r27269, MPFR_RNDN); ; mpfr_mul(r27272, r27264, r27266, MPFR_RNDN); mpfr_sub(r27273, r27271, r27272, MPFR_RNDN); mpfr_sin(r27274, r27268, MPFR_RNDN); mpfr_mul(r27275, r27273, r27274, MPFR_RNDN); mpfr_sub(r27276, r27270, r27275, MPFR_RNDN); mpfr_mul(r27277, r27273, r27269, MPFR_RNDN); mpfr_div(r27278, r27276, r27277, MPFR_RNDN); mpfr_add(r27279, r27262, r27278, MPFR_RNDN); return mpfr_get_d(r27279, MPFR_RNDN); } static mpfr_t r27280, r27281, r27282, r27283, r27284, r27285, r27286, r27287, r27288, r27289, r27290, r27291, r27292, r27293, r27294, r27295, r27296, r27297; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r27280); mpfr_init(r27281); mpfr_init(r27282); mpfr_init(r27283); mpfr_init(r27284); mpfr_init(r27285); mpfr_init(r27286); mpfr_init(r27287); mpfr_init(r27288); mpfr_init_set_str(r27289, "1", 10, MPFR_RNDN); mpfr_init(r27290); mpfr_init(r27291); mpfr_init(r27292); mpfr_init(r27293); mpfr_init(r27294); mpfr_init(r27295); mpfr_init(r27296); mpfr_init(r27297); } double f_dm(double x, double y, double z, double a) { mpfr_set_d(r27280, x, MPFR_RNDN); mpfr_set_d(r27281, y, MPFR_RNDN); mpfr_tan(r27282, r27281, MPFR_RNDN); mpfr_set_d(r27283, z, MPFR_RNDN); mpfr_tan(r27284, r27283, MPFR_RNDN); mpfr_add(r27285, r27282, r27284, MPFR_RNDN); mpfr_set_d(r27286, a, MPFR_RNDN); mpfr_cos(r27287, r27286, MPFR_RNDN); mpfr_mul(r27288, r27285, r27287, MPFR_RNDN); ; mpfr_mul(r27290, r27282, r27284, MPFR_RNDN); mpfr_sub(r27291, r27289, r27290, MPFR_RNDN); mpfr_sin(r27292, r27286, MPFR_RNDN); mpfr_mul(r27293, r27291, r27292, MPFR_RNDN); mpfr_sub(r27294, r27288, r27293, MPFR_RNDN); mpfr_mul(r27295, r27291, r27287, MPFR_RNDN); mpfr_div(r27296, r27294, r27295, MPFR_RNDN); mpfr_add(r27297, r27280, r27296, MPFR_RNDN); return mpfr_get_d(r27297, MPFR_RNDN); }