#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (6)"; double f_if(float F, float l) { float r23262 = atan2(1.0, 0.0); float r23263 = l; float r23264 = r23262 * r23263; float r23265 = 1; float r23266 = F; float r23267 = r23266 * r23266; float r23268 = r23265 / r23267; float r23269 = tan(r23264); float r23270 = r23268 * r23269; float r23271 = r23264 - r23270; return r23271; } double f_id(double F, double l) { double r23272 = atan2(1.0, 0.0); double r23273 = l; double r23274 = r23272 * r23273; double r23275 = 1; double r23276 = F; double r23277 = r23276 * r23276; double r23278 = r23275 / r23277; double r23279 = tan(r23274); double r23280 = r23278 * r23279; double r23281 = r23274 - r23280; return r23281; } double f_of(float F, float l) { float r23282 = l; float r23283 = -3.0570141152696557e+22; bool r23284 = r23282 <= r23283; float r23285 = atan2(1.0, 0.0); float r23286 = r23282 * r23285; float r23287 = r23285 / r23282; float r23288 = sin(r23287); float r23289 = cos(r23287); float r23290 = r23288 / r23289; float r23291 = F; float r23292 = r23290 / r23291; float r23293 = r23292 / r23291; float r23294 = r23286 - r23293; float r23295 = tan(r23286); float r23296 = r23295 / r23291; float r23297 = r23296 / r23291; float r23298 = r23286 - r23297; float r23299 = r23284 ? r23294 : r23298; return r23299; } double f_od(double F, double l) { double r23300 = l; double r23301 = -3.0570141152696557e+22; bool r23302 = r23300 <= r23301; double r23303 = atan2(1.0, 0.0); double r23304 = r23300 * r23303; double r23305 = r23303 / r23300; double r23306 = sin(r23305); double r23307 = cos(r23305); double r23308 = r23306 / r23307; double r23309 = F; double r23310 = r23308 / r23309; double r23311 = r23310 / r23309; double r23312 = r23304 - r23311; double r23313 = tan(r23304); double r23314 = r23313 / r23309; double r23315 = r23314 / r23309; double r23316 = r23304 - r23315; double r23317 = r23302 ? r23312 : r23316; return r23317; } 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 r23318, r23319, r23320, r23321, r23322, r23323, r23324, r23325, r23326, r23327; void setup_mpfr_f_im() { mpfr_set_default_prec(2448); mpfr_init(r23318); mpfr_init(r23319); mpfr_init(r23320); mpfr_init_set_str(r23321, "1", 10, MPFR_RNDN); mpfr_init(r23322); mpfr_init(r23323); mpfr_init(r23324); mpfr_init(r23325); mpfr_init(r23326); mpfr_init(r23327); } double f_im(double F, double l) { mpfr_const_pi(r23318, MPFR_RNDN); mpfr_set_d(r23319, l, MPFR_RNDN); mpfr_mul(r23320, r23318, r23319, MPFR_RNDN); ; mpfr_set_d(r23322, F, MPFR_RNDN); mpfr_mul(r23323, r23322, r23322, MPFR_RNDN); mpfr_div(r23324, r23321, r23323, MPFR_RNDN); mpfr_tan(r23325, r23320, MPFR_RNDN); mpfr_mul(r23326, r23324, r23325, MPFR_RNDN); mpfr_sub(r23327, r23320, r23326, MPFR_RNDN); return mpfr_get_d(r23327, MPFR_RNDN); } static mpfr_t r23328, r23329, r23330, r23331, r23332, r23333, r23334, r23335, r23336, r23337, r23338, r23339, r23340, r23341, r23342, r23343, r23344, r23345; void setup_mpfr_f_fm() { mpfr_set_default_prec(2448); mpfr_init(r23328); mpfr_init_set_str(r23329, "-3.0570141152696557e+22", 10, MPFR_RNDN); mpfr_init(r23330); mpfr_init(r23331); mpfr_init(r23332); mpfr_init(r23333); mpfr_init(r23334); mpfr_init(r23335); mpfr_init(r23336); mpfr_init(r23337); mpfr_init(r23338); mpfr_init(r23339); mpfr_init(r23340); mpfr_init(r23341); mpfr_init(r23342); mpfr_init(r23343); mpfr_init(r23344); mpfr_init(r23345); } double f_fm(double F, double l) { mpfr_set_d(r23328, l, MPFR_RNDN); ; mpfr_set_si(r23330, mpfr_cmp(r23328, r23329) <= 0, MPFR_RNDN); mpfr_const_pi(r23331, MPFR_RNDN); mpfr_mul(r23332, r23328, r23331, MPFR_RNDN); mpfr_div(r23333, r23331, r23328, MPFR_RNDN); mpfr_sin(r23334, r23333, MPFR_RNDN); mpfr_cos(r23335, r23333, MPFR_RNDN); mpfr_div(r23336, r23334, r23335, MPFR_RNDN); mpfr_set_d(r23337, F, MPFR_RNDN); mpfr_div(r23338, r23336, r23337, MPFR_RNDN); mpfr_div(r23339, r23338, r23337, MPFR_RNDN); mpfr_sub(r23340, r23332, r23339, MPFR_RNDN); mpfr_tan(r23341, r23332, MPFR_RNDN); mpfr_div(r23342, r23341, r23337, MPFR_RNDN); mpfr_div(r23343, r23342, r23337, MPFR_RNDN); mpfr_sub(r23344, r23332, r23343, MPFR_RNDN); if (mpfr_get_si(r23330, MPFR_RNDN)) { mpfr_set(r23345, r23340, MPFR_RNDN); } else { mpfr_set(r23345, r23344, MPFR_RNDN); }; return mpfr_get_d(r23345, MPFR_RNDN); } static mpfr_t r23346, r23347, r23348, r23349, r23350, r23351, r23352, r23353, r23354, r23355, r23356, r23357, r23358, r23359, r23360, r23361, r23362, r23363; void setup_mpfr_f_dm() { mpfr_set_default_prec(2448); mpfr_init(r23346); mpfr_init_set_str(r23347, "-3.0570141152696557e+22", 10, MPFR_RNDN); mpfr_init(r23348); mpfr_init(r23349); mpfr_init(r23350); mpfr_init(r23351); mpfr_init(r23352); mpfr_init(r23353); mpfr_init(r23354); mpfr_init(r23355); mpfr_init(r23356); mpfr_init(r23357); mpfr_init(r23358); mpfr_init(r23359); mpfr_init(r23360); mpfr_init(r23361); mpfr_init(r23362); mpfr_init(r23363); } double f_dm(double F, double l) { mpfr_set_d(r23346, l, MPFR_RNDN); ; mpfr_set_si(r23348, mpfr_cmp(r23346, r23347) <= 0, MPFR_RNDN); mpfr_const_pi(r23349, MPFR_RNDN); mpfr_mul(r23350, r23346, r23349, MPFR_RNDN); mpfr_div(r23351, r23349, r23346, MPFR_RNDN); mpfr_sin(r23352, r23351, MPFR_RNDN); mpfr_cos(r23353, r23351, MPFR_RNDN); mpfr_div(r23354, r23352, r23353, MPFR_RNDN); mpfr_set_d(r23355, F, MPFR_RNDN); mpfr_div(r23356, r23354, r23355, MPFR_RNDN); mpfr_div(r23357, r23356, r23355, MPFR_RNDN); mpfr_sub(r23358, r23350, r23357, MPFR_RNDN); mpfr_tan(r23359, r23350, MPFR_RNDN); mpfr_div(r23360, r23359, r23355, MPFR_RNDN); mpfr_div(r23361, r23360, r23355, MPFR_RNDN); mpfr_sub(r23362, r23350, r23361, MPFR_RNDN); if (mpfr_get_si(r23348, MPFR_RNDN)) { mpfr_set(r23363, r23358, MPFR_RNDN); } else { mpfr_set(r23363, r23362, MPFR_RNDN); }; return mpfr_get_d(r23363, MPFR_RNDN); }