#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 = -6.666946136074473e-165f; bool r15260 = r15258 <= r15259; float r15261 = r15258 * (r15258 * r15258); float r15262 = asin(r15258); float r15263 = 1.0 / tan(r15258); float r15264 = a; float r15265 = pow(r15263, r15264); float r15266 = r15262 + r15265; float r15267 = r15266 * (r15266 * r15266); float r15268 = r15261 - r15267; float r15269 = r15258 + r15262; float r15270 = r15265 + r15269; float r15271 = r15258 * r15258; float r15272 = fma(r15270, r15266, r15271); float r15273 = r15268 / r15272; float r15274 = 7.6539915406140805e-165f; bool r15275 = r15258 <= r15274; float r15276 = r15265 + r15262; float r15277 = r15276 * r15276; float r15278 = r15277 * (r15277 * r15277); float r15279 = cbrt(r15278); float r15280 = r15271 - r15279; float r15281 = r15258 + r15276; float r15282 = r15280 / r15281; float r15283 = r15275 ? r15282 : r15273; float r15284 = r15260 ? r15273 : r15283; return r15284; } double f_od(double a, double b) { double r15285 = b; double r15286 = -6.666946136074473e-165; bool r15287 = r15285 <= r15286; double r15288 = r15285 * (r15285 * r15285); double r15289 = asin(r15285); double r15290 = 1.0 / tan(r15285); double r15291 = a; double r15292 = pow(r15290, r15291); double r15293 = r15289 + r15292; double r15294 = r15293 * (r15293 * r15293); double r15295 = r15288 - r15294; double r15296 = r15285 + r15289; double r15297 = r15292 + r15296; double r15298 = r15285 * r15285; double r15299 = fma(r15297, r15293, r15298); double r15300 = r15295 / r15299; double r15301 = 7.6539915406140805e-165; bool r15302 = r15285 <= r15301; double r15303 = r15292 + r15289; double r15304 = r15303 * r15303; double r15305 = r15304 * (r15304 * r15304); double r15306 = cbrt(r15305); double r15307 = r15298 - r15306; double r15308 = r15285 + r15303; double r15309 = r15307 / r15308; double r15310 = r15302 ? r15309 : r15300; double r15311 = r15287 ? r15300 : r15310; return r15311; } 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 r15312, r15313, r15314, r15315, r15316, r15317, r15318; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r15312); mpfr_init(r15313); mpfr_init(r15314); mpfr_init(r15315); mpfr_init(r15316); mpfr_init(r15317); mpfr_init(r15318); } double f_im(double a, double b) { mpfr_set_d(r15312, b, MPFR_RNDN); mpfr_cot(r15313, r15312, MPFR_RNDN); mpfr_set_d(r15314, a, MPFR_RNDN); mpfr_pow(r15315, r15313, r15314, MPFR_RNDN); mpfr_asin(r15316, r15312, MPFR_RNDN); mpfr_add(r15317, r15315, r15316, MPFR_RNDN); mpfr_sub(r15318, r15312, 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, r15337, r15338, r15339, r15340, r15341, r15342, r15343, r15344, r15345; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r15319); mpfr_init_set_str(r15320, "-6.666946136074473e-165", 10, MPFR_RNDN); 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_set_str(r15335, "7.6539915406140805e-165", 10, MPFR_RNDN); mpfr_init(r15336); mpfr_init(r15337); mpfr_init(r15338); mpfr_init(r15339); mpfr_init(r15340); mpfr_init(r15341); mpfr_init(r15342); mpfr_init(r15343); mpfr_init(r15344); mpfr_init(r15345); } double f_fm(double a, double b) { mpfr_set_d(r15319, b, MPFR_RNDN); ; mpfr_set_si(r15321, mpfr_cmp(r15319, r15320) <= 0, MPFR_RNDN); mpfr_mul(r15322, r15319, r15319, MPFR_RNDN); mpfr_mul(r15322, r15322, r15319, MPFR_RNDN); mpfr_asin(r15323, r15319, MPFR_RNDN); mpfr_cot(r15324, r15319, MPFR_RNDN); mpfr_set_d(r15325, a, MPFR_RNDN); mpfr_pow(r15326, r15324, r15325, MPFR_RNDN); mpfr_add(r15327, r15323, r15326, MPFR_RNDN); mpfr_mul(r15328, r15327, r15327, MPFR_RNDN); mpfr_mul(r15328, r15328, r15327, MPFR_RNDN); mpfr_sub(r15329, r15322, r15328, MPFR_RNDN); mpfr_add(r15330, r15319, r15323, MPFR_RNDN); mpfr_add(r15331, r15326, r15330, MPFR_RNDN); mpfr_sqr(r15332, r15319, MPFR_RNDN); mpfr_fma(r15333, r15331, r15327, r15332, MPFR_RNDN); mpfr_div(r15334, r15329, r15333, MPFR_RNDN); ; mpfr_set_si(r15336, mpfr_cmp(r15319, r15335) <= 0, MPFR_RNDN); mpfr_add(r15337, r15326, r15323, MPFR_RNDN); mpfr_sqr(r15338, r15337, MPFR_RNDN); mpfr_mul(r15339, r15338, r15338, MPFR_RNDN); mpfr_mul(r15339, r15339, r15338, MPFR_RNDN); mpfr_cbrt(r15340, r15339, MPFR_RNDN); mpfr_sub(r15341, r15332, r15340, MPFR_RNDN); mpfr_add(r15342, r15319, r15337, MPFR_RNDN); mpfr_div(r15343, r15341, r15342, MPFR_RNDN); if (mpfr_get_si(r15336, MPFR_RNDN)) { mpfr_set(r15344, r15343, MPFR_RNDN); } else { mpfr_set(r15344, r15334, MPFR_RNDN); }; if (mpfr_get_si(r15321, MPFR_RNDN)) { mpfr_set(r15345, r15334, MPFR_RNDN); } else { mpfr_set(r15345, r15344, MPFR_RNDN); }; return mpfr_get_d(r15345, MPFR_RNDN); } static mpfr_t r15346, r15347, r15348, r15349, r15350, r15351, r15352, r15353, r15354, r15355, r15356, r15357, r15358, r15359, r15360, r15361, r15362, r15363, r15364, r15365, r15366, r15367, r15368, r15369, r15370, r15371, r15372; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r15346); mpfr_init_set_str(r15347, "-6.666946136074473e-165", 10, MPFR_RNDN); mpfr_init(r15348); mpfr_init(r15349); mpfr_init(r15350); mpfr_init(r15351); mpfr_init(r15352); mpfr_init(r15353); mpfr_init(r15354); mpfr_init(r15355); mpfr_init(r15356); mpfr_init(r15357); mpfr_init(r15358); mpfr_init(r15359); mpfr_init(r15360); mpfr_init(r15361); mpfr_init_set_str(r15362, "7.6539915406140805e-165", 10, MPFR_RNDN); mpfr_init(r15363); mpfr_init(r15364); mpfr_init(r15365); mpfr_init(r15366); mpfr_init(r15367); mpfr_init(r15368); mpfr_init(r15369); mpfr_init(r15370); mpfr_init(r15371); mpfr_init(r15372); } double f_dm(double a, double b) { mpfr_set_d(r15346, b, MPFR_RNDN); ; mpfr_set_si(r15348, mpfr_cmp(r15346, r15347) <= 0, MPFR_RNDN); mpfr_mul(r15349, r15346, r15346, MPFR_RNDN); mpfr_mul(r15349, r15349, r15346, MPFR_RNDN); mpfr_asin(r15350, r15346, MPFR_RNDN); mpfr_cot(r15351, r15346, MPFR_RNDN); mpfr_set_d(r15352, a, MPFR_RNDN); mpfr_pow(r15353, r15351, r15352, MPFR_RNDN); mpfr_add(r15354, r15350, r15353, MPFR_RNDN); mpfr_mul(r15355, r15354, r15354, MPFR_RNDN); mpfr_mul(r15355, r15355, r15354, MPFR_RNDN); mpfr_sub(r15356, r15349, r15355, MPFR_RNDN); mpfr_add(r15357, r15346, r15350, MPFR_RNDN); mpfr_add(r15358, r15353, r15357, MPFR_RNDN); mpfr_sqr(r15359, r15346, MPFR_RNDN); mpfr_fma(r15360, r15358, r15354, r15359, MPFR_RNDN); mpfr_div(r15361, r15356, r15360, MPFR_RNDN); ; mpfr_set_si(r15363, mpfr_cmp(r15346, r15362) <= 0, MPFR_RNDN); mpfr_add(r15364, r15353, r15350, MPFR_RNDN); mpfr_sqr(r15365, r15364, MPFR_RNDN); mpfr_mul(r15366, r15365, r15365, MPFR_RNDN); mpfr_mul(r15366, r15366, r15365, MPFR_RNDN); mpfr_cbrt(r15367, r15366, MPFR_RNDN); mpfr_sub(r15368, r15359, r15367, MPFR_RNDN); mpfr_add(r15369, r15346, r15364, MPFR_RNDN); mpfr_div(r15370, r15368, r15369, MPFR_RNDN); if (mpfr_get_si(r15363, MPFR_RNDN)) { mpfr_set(r15371, r15370, MPFR_RNDN); } else { mpfr_set(r15371, r15361, MPFR_RNDN); }; if (mpfr_get_si(r15348, MPFR_RNDN)) { mpfr_set(r15372, r15361, MPFR_RNDN); } else { mpfr_set(r15372, r15371, MPFR_RNDN); }; return mpfr_get_d(r15372, MPFR_RNDN); }