#include #include #include #include #include char *name = "(- b (+ (pow (cotan b) a) (asin b)))"; double f_if(float a, float b) { float r15244 = b; float r15245 = 1.0 / tan(r15244); float r15246 = a; float r15247 = pow(r15245, r15246); float r15248 = asin(r15244); float r15249 = r15247 + r15248; float r15250 = r15244 - r15249; return r15250; } double f_id(double a, double b) { double r15251 = b; double r15252 = 1.0 / tan(r15251); double r15253 = a; double r15254 = pow(r15252, r15253); double r15255 = asin(r15251); double r15256 = r15254 + r15255; double r15257 = r15251 - r15256; return r15257; } double f_of(float a, float b) { float r15258 = b; float r15259 = r15258 * (r15258 * r15258); float r15260 = asin(r15258); float r15261 = 1.0 / tan(r15258); float r15262 = a; float r15263 = pow(r15261, r15262); float r15264 = r15260 + r15263; float r15265 = r15264 * (r15264 * r15264); float r15266 = r15259 - r15265; float r15267 = r15266 * (r15266 * r15266); float r15268 = r15258 * r15258; float r15269 = r15258 + r15260; float r15270 = r15269 + r15263; float r15271 = r15270 * r15264; float r15272 = r15268 + r15271; float r15273 = r15272 * (r15272 * r15272); float r15274 = r15267 / r15273; float r15275 = cbrt(r15274); return r15275; } double f_od(double a, double b) { double r15276 = b; double r15277 = r15276 * (r15276 * r15276); double r15278 = asin(r15276); double r15279 = 1.0 / tan(r15276); double r15280 = a; double r15281 = pow(r15279, r15280); double r15282 = r15278 + r15281; double r15283 = r15282 * (r15282 * r15282); double r15284 = r15277 - r15283; double r15285 = r15284 * (r15284 * r15284); double r15286 = r15276 * r15276; double r15287 = r15276 + r15278; double r15288 = r15287 + r15281; double r15289 = r15288 * r15282; double r15290 = r15286 + r15289; double r15291 = r15290 * (r15290 * r15290); double r15292 = r15285 / r15291; double r15293 = cbrt(r15292); return r15293; } 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 r15294, r15295, r15296, r15297, r15298, r15299, r15300; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r15294); mpfr_init(r15295); mpfr_init(r15296); mpfr_init(r15297); mpfr_init(r15298); mpfr_init(r15299); mpfr_init(r15300); } double f_im(double a, double b) { mpfr_set_d(r15294, b, MPFR_RNDN); mpfr_cot(r15295, r15294, MPFR_RNDN); mpfr_set_d(r15296, a, MPFR_RNDN); mpfr_pow(r15297, r15295, r15296, MPFR_RNDN); mpfr_asin(r15298, r15294, MPFR_RNDN); mpfr_add(r15299, r15297, r15298, MPFR_RNDN); mpfr_sub(r15300, r15294, r15299, MPFR_RNDN); return mpfr_get_d(r15300, MPFR_RNDN); } static mpfr_t r15301, r15302, r15303, r15304, r15305, r15306, r15307, r15308, r15309, r15310, r15311, r15312, r15313, r15314, r15315, r15316, r15317, r15318; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r15301); mpfr_init(r15302); mpfr_init(r15303); mpfr_init(r15304); mpfr_init(r15305); mpfr_init(r15306); mpfr_init(r15307); mpfr_init(r15308); mpfr_init(r15309); mpfr_init(r15310); mpfr_init(r15311); mpfr_init(r15312); mpfr_init(r15313); mpfr_init(r15314); mpfr_init(r15315); mpfr_init(r15316); mpfr_init(r15317); mpfr_init(r15318); } double f_fm(double a, double b) { mpfr_set_d(r15301, b, MPFR_RNDN); mpfr_mul(r15302, r15301, r15301, MPFR_RNDN); mpfr_mul(r15302, r15302, r15301, MPFR_RNDN); mpfr_asin(r15303, r15301, MPFR_RNDN); mpfr_cot(r15304, r15301, MPFR_RNDN); mpfr_set_d(r15305, a, MPFR_RNDN); mpfr_pow(r15306, r15304, r15305, MPFR_RNDN); mpfr_add(r15307, r15303, r15306, MPFR_RNDN); mpfr_mul(r15308, r15307, r15307, MPFR_RNDN); mpfr_mul(r15308, r15308, r15307, MPFR_RNDN); mpfr_sub(r15309, r15302, r15308, MPFR_RNDN); mpfr_mul(r15310, r15309, r15309, MPFR_RNDN); mpfr_mul(r15310, r15310, r15309, MPFR_RNDN); mpfr_sqr(r15311, r15301, MPFR_RNDN); mpfr_add(r15312, r15301, r15303, MPFR_RNDN); mpfr_add(r15313, r15312, r15306, MPFR_RNDN); mpfr_mul(r15314, r15313, r15307, MPFR_RNDN); mpfr_add(r15315, r15311, r15314, MPFR_RNDN); mpfr_mul(r15316, r15315, r15315, MPFR_RNDN); mpfr_mul(r15316, r15316, r15315, MPFR_RNDN); mpfr_div(r15317, r15310, r15316, MPFR_RNDN); mpfr_cbrt(r15318, r15317, MPFR_RNDN); return mpfr_get_d(r15318, MPFR_RNDN); } static mpfr_t r15319, r15320, r15321, r15322, r15323, r15324, r15325, r15326, r15327, r15328, r15329, r15330, r15331, r15332, r15333, r15334, r15335, r15336; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r15319); mpfr_init(r15320); mpfr_init(r15321); mpfr_init(r15322); mpfr_init(r15323); mpfr_init(r15324); mpfr_init(r15325); mpfr_init(r15326); mpfr_init(r15327); mpfr_init(r15328); mpfr_init(r15329); mpfr_init(r15330); mpfr_init(r15331); mpfr_init(r15332); mpfr_init(r15333); mpfr_init(r15334); mpfr_init(r15335); mpfr_init(r15336); } double f_dm(double a, double b) { mpfr_set_d(r15319, b, MPFR_RNDN); mpfr_mul(r15320, r15319, r15319, MPFR_RNDN); mpfr_mul(r15320, r15320, r15319, MPFR_RNDN); mpfr_asin(r15321, r15319, MPFR_RNDN); mpfr_cot(r15322, r15319, MPFR_RNDN); mpfr_set_d(r15323, a, MPFR_RNDN); mpfr_pow(r15324, r15322, r15323, MPFR_RNDN); mpfr_add(r15325, r15321, r15324, MPFR_RNDN); mpfr_mul(r15326, r15325, r15325, MPFR_RNDN); mpfr_mul(r15326, r15326, r15325, MPFR_RNDN); mpfr_sub(r15327, r15320, r15326, MPFR_RNDN); mpfr_mul(r15328, r15327, r15327, MPFR_RNDN); mpfr_mul(r15328, r15328, r15327, MPFR_RNDN); mpfr_sqr(r15329, r15319, MPFR_RNDN); mpfr_add(r15330, r15319, r15321, MPFR_RNDN); mpfr_add(r15331, r15330, r15324, MPFR_RNDN); mpfr_mul(r15332, r15331, r15325, MPFR_RNDN); mpfr_add(r15333, r15329, r15332, MPFR_RNDN); mpfr_mul(r15334, r15333, r15333, MPFR_RNDN); mpfr_mul(r15334, r15334, r15333, MPFR_RNDN); mpfr_div(r15335, r15328, r15334, MPFR_RNDN); mpfr_cbrt(r15336, r15335, MPFR_RNDN); return mpfr_get_d(r15336, MPFR_RNDN); }