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