#include #include #include #include #include char *name = "(- b (+ (pow (cotan b) a) (asin b)))"; double f_if(float a, float b) { float r15220 = b; float r15221 = 1.0 / tan(r15220); float r15222 = a; float r15223 = pow(r15221, r15222); float r15224 = asin(r15220); float r15225 = r15223 + r15224; float r15226 = r15220 - r15225; return r15226; } double f_id(double a, double b) { double r15227 = b; double r15228 = 1.0 / tan(r15227); double r15229 = a; double r15230 = pow(r15228, r15229); double r15231 = asin(r15227); double r15232 = r15230 + r15231; double r15233 = r15227 - r15232; return r15233; } double f_of(float a, float b) { float r15234 = b; float r15235 = -6.666946136074473e-165f; bool r15236 = r15234 <= r15235; float r15237 = r15234 * (r15234 * r15234); float r15238 = asin(r15234); float r15239 = 1.0 / tan(r15234); float r15240 = a; float r15241 = pow(r15239, r15240); float r15242 = r15238 + r15241; float r15243 = r15242 * (r15242 * r15242); float r15244 = r15237 - r15243; float r15245 = r15234 + r15238; float r15246 = r15241 + r15245; float r15247 = r15234 * r15246; float r15248 = r15242 * r15242; float r15249 = r15247 + r15248; float r15250 = r15244 / r15249; float r15251 = 7.6539915406140805e-165f; bool r15252 = r15234 <= r15251; float r15253 = r15234 * r15234; float r15254 = r15241 + r15238; float r15255 = r15254 * r15254; float r15256 = r15255 * (r15255 * r15255); float r15257 = cbrt(r15256); float r15258 = r15253 - r15257; float r15259 = r15234 + r15254; float r15260 = r15258 / r15259; float r15261 = r15252 ? r15260 : r15250; float r15262 = r15236 ? r15250 : r15261; return r15262; } double f_od(double a, double b) { double r15263 = b; double r15264 = -6.666946136074473e-165; bool r15265 = r15263 <= r15264; double r15266 = r15263 * (r15263 * r15263); double r15267 = asin(r15263); double r15268 = 1.0 / tan(r15263); double r15269 = a; double r15270 = pow(r15268, r15269); double r15271 = r15267 + r15270; double r15272 = r15271 * (r15271 * r15271); double r15273 = r15266 - r15272; double r15274 = r15263 + r15267; double r15275 = r15270 + r15274; double r15276 = r15263 * r15275; double r15277 = r15271 * r15271; double r15278 = r15276 + r15277; double r15279 = r15273 / r15278; double r15280 = 7.6539915406140805e-165; bool r15281 = r15263 <= r15280; double r15282 = r15263 * r15263; double r15283 = r15270 + r15267; double r15284 = r15283 * r15283; double r15285 = r15284 * (r15284 * r15284); double r15286 = cbrt(r15285); double r15287 = r15282 - r15286; double r15288 = r15263 + r15283; double r15289 = r15287 / r15288; double r15290 = r15281 ? r15289 : r15279; double r15291 = r15265 ? r15279 : r15290; return r15291; } 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 r15292, r15293, r15294, r15295, r15296, r15297, r15298; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r15292); mpfr_init(r15293); mpfr_init(r15294); mpfr_init(r15295); mpfr_init(r15296); mpfr_init(r15297); mpfr_init(r15298); } double f_im(double a, double b) { mpfr_set_d(r15292, b, MPFR_RNDN); mpfr_cot(r15293, r15292, MPFR_RNDN); mpfr_set_d(r15294, a, MPFR_RNDN); mpfr_pow(r15295, r15293, r15294, MPFR_RNDN); mpfr_asin(r15296, r15292, MPFR_RNDN); mpfr_add(r15297, r15295, r15296, MPFR_RNDN); mpfr_sub(r15298, r15292, r15297, MPFR_RNDN); return mpfr_get_d(r15298, MPFR_RNDN); } static mpfr_t r15299, r15300, r15301, r15302, r15303, r15304, r15305, r15306, r15307, r15308, r15309, r15310, r15311, r15312, r15313, r15314, r15315, r15316, r15317, r15318, r15319, r15320, r15321, r15322, r15323, r15324, r15325, r15326, r15327; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r15299); mpfr_init_set_str(r15300, "-6.666946136074473e-165", 10, MPFR_RNDN); 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_set_str(r15316, "7.6539915406140805e-165", 10, MPFR_RNDN); mpfr_init(r15317); mpfr_init(r15318); 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); } double f_fm(double a, double b) { mpfr_set_d(r15299, b, MPFR_RNDN); ; mpfr_set_si(r15301, mpfr_cmp(r15299, r15300) <= 0, MPFR_RNDN); mpfr_mul(r15302, r15299, r15299, MPFR_RNDN); mpfr_mul(r15302, r15302, r15299, MPFR_RNDN); mpfr_asin(r15303, r15299, MPFR_RNDN); mpfr_cot(r15304, r15299, 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_add(r15310, r15299, r15303, MPFR_RNDN); mpfr_add(r15311, r15306, r15310, MPFR_RNDN); mpfr_mul(r15312, r15299, r15311, MPFR_RNDN); mpfr_sqr(r15313, r15307, MPFR_RNDN); mpfr_add(r15314, r15312, r15313, MPFR_RNDN); mpfr_div(r15315, r15309, r15314, MPFR_RNDN); ; mpfr_set_si(r15317, mpfr_cmp(r15299, r15316) <= 0, MPFR_RNDN); mpfr_sqr(r15318, r15299, MPFR_RNDN); mpfr_add(r15319, r15306, r15303, MPFR_RNDN); mpfr_sqr(r15320, r15319, MPFR_RNDN); mpfr_mul(r15321, r15320, r15320, MPFR_RNDN); mpfr_mul(r15321, r15321, r15320, MPFR_RNDN); mpfr_cbrt(r15322, r15321, MPFR_RNDN); mpfr_sub(r15323, r15318, r15322, MPFR_RNDN); mpfr_add(r15324, r15299, r15319, MPFR_RNDN); mpfr_div(r15325, r15323, r15324, MPFR_RNDN); if (mpfr_get_si(r15317, MPFR_RNDN)) { mpfr_set(r15326, r15325, MPFR_RNDN); } else { mpfr_set(r15326, r15315, MPFR_RNDN); }; if (mpfr_get_si(r15301, MPFR_RNDN)) { mpfr_set(r15327, r15315, MPFR_RNDN); } else { mpfr_set(r15327, r15326, MPFR_RNDN); }; return mpfr_get_d(r15327, MPFR_RNDN); } static mpfr_t r15328, r15329, r15330, r15331, r15332, r15333, r15334, r15335, r15336, r15337, r15338, r15339, r15340, r15341, r15342, r15343, r15344, r15345, r15346, r15347, r15348, r15349, r15350, r15351, r15352, r15353, r15354, r15355, r15356; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r15328); mpfr_init_set_str(r15329, "-6.666946136074473e-165", 10, MPFR_RNDN); mpfr_init(r15330); mpfr_init(r15331); mpfr_init(r15332); mpfr_init(r15333); mpfr_init(r15334); mpfr_init(r15335); 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_set_str(r15345, "7.6539915406140805e-165", 10, MPFR_RNDN); mpfr_init(r15346); mpfr_init(r15347); 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); } double f_dm(double a, double b) { mpfr_set_d(r15328, b, MPFR_RNDN); ; mpfr_set_si(r15330, mpfr_cmp(r15328, r15329) <= 0, MPFR_RNDN); mpfr_mul(r15331, r15328, r15328, MPFR_RNDN); mpfr_mul(r15331, r15331, r15328, MPFR_RNDN); mpfr_asin(r15332, r15328, MPFR_RNDN); mpfr_cot(r15333, r15328, MPFR_RNDN); mpfr_set_d(r15334, a, MPFR_RNDN); mpfr_pow(r15335, r15333, r15334, MPFR_RNDN); mpfr_add(r15336, r15332, r15335, MPFR_RNDN); mpfr_mul(r15337, r15336, r15336, MPFR_RNDN); mpfr_mul(r15337, r15337, r15336, MPFR_RNDN); mpfr_sub(r15338, r15331, r15337, MPFR_RNDN); mpfr_add(r15339, r15328, r15332, MPFR_RNDN); mpfr_add(r15340, r15335, r15339, MPFR_RNDN); mpfr_mul(r15341, r15328, r15340, MPFR_RNDN); mpfr_sqr(r15342, r15336, MPFR_RNDN); mpfr_add(r15343, r15341, r15342, MPFR_RNDN); mpfr_div(r15344, r15338, r15343, MPFR_RNDN); ; mpfr_set_si(r15346, mpfr_cmp(r15328, r15345) <= 0, MPFR_RNDN); mpfr_sqr(r15347, r15328, MPFR_RNDN); mpfr_add(r15348, r15335, r15332, MPFR_RNDN); mpfr_sqr(r15349, r15348, MPFR_RNDN); mpfr_mul(r15350, r15349, r15349, MPFR_RNDN); mpfr_mul(r15350, r15350, r15349, MPFR_RNDN); mpfr_cbrt(r15351, r15350, MPFR_RNDN); mpfr_sub(r15352, r15347, r15351, MPFR_RNDN); mpfr_add(r15353, r15328, r15348, MPFR_RNDN); mpfr_div(r15354, r15352, r15353, MPFR_RNDN); if (mpfr_get_si(r15346, MPFR_RNDN)) { mpfr_set(r15355, r15354, MPFR_RNDN); } else { mpfr_set(r15355, r15344, MPFR_RNDN); }; if (mpfr_get_si(r15330, MPFR_RNDN)) { mpfr_set(r15356, r15344, MPFR_RNDN); } else { mpfr_set(r15356, r15355, MPFR_RNDN); }; return mpfr_get_d(r15356, MPFR_RNDN); }