#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r30465252 = a; float r30465253 = r30465252 * r30465252; float r30465254 = b; float r30465255 = r30465254 * r30465254; float r30465256 = r30465253 + r30465255; float r30465257 = 2; float r30465258 = pow(r30465256, r30465257); float r30465259 = 4; float r30465260 = 1; float r30465261 = r30465260 + r30465252; float r30465262 = r30465253 * r30465261; float r30465263 = 3; float r30465264 = r30465263 * r30465252; float r30465265 = r30465260 - r30465264; float r30465266 = r30465255 * r30465265; float r30465267 = r30465262 + r30465266; float r30465268 = r30465259 * r30465267; float r30465269 = r30465258 + r30465268; float r30465270 = r30465269 - r30465260; return r30465270; } double f_id(double a, double b) { double r30465271 = a; double r30465272 = r30465271 * r30465271; double r30465273 = b; double r30465274 = r30465273 * r30465273; double r30465275 = r30465272 + r30465274; double r30465276 = 2; double r30465277 = pow(r30465275, r30465276); double r30465278 = 4; double r30465279 = 1; double r30465280 = r30465279 + r30465271; double r30465281 = r30465272 * r30465280; double r30465282 = 3; double r30465283 = r30465282 * r30465271; double r30465284 = r30465279 - r30465283; double r30465285 = r30465274 * r30465284; double r30465286 = r30465281 + r30465285; double r30465287 = r30465278 * r30465286; double r30465288 = r30465277 + r30465287; double r30465289 = r30465288 - r30465279; return r30465289; } double f_of(float a, float b) { float r30465290 = a; float r30465291 = r30465290 * r30465290; float r30465292 = b; float r30465293 = r30465292 * r30465292; float r30465294 = r30465291 + r30465293; float r30465295 = 2; float r30465296 = pow(r30465294, r30465295); float r30465297 = 4; float r30465298 = 1; float r30465299 = r30465298 + r30465290; float r30465300 = r30465291 * r30465299; float r30465301 = 3; float r30465302 = r30465301 * r30465290; float r30465303 = r30465298 - r30465302; float r30465304 = r30465293 * r30465303; float r30465305 = r30465300 + r30465304; float r30465306 = r30465297 * r30465305; float r30465307 = r30465296 + r30465306; float r30465308 = sqrt(r30465307); float r30465309 = r30465308 * r30465308; float r30465310 = r30465309 - r30465298; return r30465310; } double f_od(double a, double b) { double r30465311 = a; double r30465312 = r30465311 * r30465311; double r30465313 = b; double r30465314 = r30465313 * r30465313; double r30465315 = r30465312 + r30465314; double r30465316 = 2; double r30465317 = pow(r30465315, r30465316); double r30465318 = 4; double r30465319 = 1; double r30465320 = r30465319 + r30465311; double r30465321 = r30465312 * r30465320; double r30465322 = 3; double r30465323 = r30465322 * r30465311; double r30465324 = r30465319 - r30465323; double r30465325 = r30465314 * r30465324; double r30465326 = r30465321 + r30465325; double r30465327 = r30465318 * r30465326; double r30465328 = r30465317 + r30465327; double r30465329 = sqrt(r30465328); double r30465330 = r30465329 * r30465329; double r30465331 = r30465330 - r30465319; return r30465331; } 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 r30465332, r30465333, r30465334, r30465335, r30465336, r30465337, r30465338, r30465339, r30465340, r30465341, r30465342, r30465343, r30465344, r30465345, r30465346, r30465347, r30465348, r30465349, r30465350; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r30465332); mpfr_init(r30465333); mpfr_init(r30465334); mpfr_init(r30465335); mpfr_init(r30465336); mpfr_init_set_str(r30465337, "2", 10, MPFR_RNDN); mpfr_init(r30465338); mpfr_init_set_str(r30465339, "4", 10, MPFR_RNDN); mpfr_init_set_str(r30465340, "1", 10, MPFR_RNDN); mpfr_init(r30465341); mpfr_init(r30465342); mpfr_init_set_str(r30465343, "3", 10, MPFR_RNDN); mpfr_init(r30465344); mpfr_init(r30465345); mpfr_init(r30465346); mpfr_init(r30465347); mpfr_init(r30465348); mpfr_init(r30465349); mpfr_init(r30465350); } double f_im(double a, double b) { mpfr_set_d(r30465332, a, MPFR_RNDN); mpfr_mul(r30465333, r30465332, r30465332, MPFR_RNDN); mpfr_set_d(r30465334, b, MPFR_RNDN); mpfr_mul(r30465335, r30465334, r30465334, MPFR_RNDN); mpfr_add(r30465336, r30465333, r30465335, MPFR_RNDN); ; mpfr_pow(r30465338, r30465336, r30465337, MPFR_RNDN); ; ; mpfr_add(r30465341, r30465340, r30465332, MPFR_RNDN); mpfr_mul(r30465342, r30465333, r30465341, MPFR_RNDN); ; mpfr_mul(r30465344, r30465343, r30465332, MPFR_RNDN); mpfr_sub(r30465345, r30465340, r30465344, MPFR_RNDN); mpfr_mul(r30465346, r30465335, r30465345, MPFR_RNDN); mpfr_add(r30465347, r30465342, r30465346, MPFR_RNDN); mpfr_mul(r30465348, r30465339, r30465347, MPFR_RNDN); mpfr_add(r30465349, r30465338, r30465348, MPFR_RNDN); mpfr_sub(r30465350, r30465349, r30465340, MPFR_RNDN); return mpfr_get_d(r30465350, MPFR_RNDN); } static mpfr_t r30465351, r30465352, r30465353, r30465354, r30465355, r30465356, r30465357, r30465358, r30465359, r30465360, r30465361, r30465362, r30465363, r30465364, r30465365, r30465366, r30465367, r30465368, r30465369, r30465370, r30465371; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r30465351); mpfr_init(r30465352); mpfr_init(r30465353); mpfr_init(r30465354); mpfr_init(r30465355); mpfr_init_set_str(r30465356, "2", 10, MPFR_RNDN); mpfr_init(r30465357); mpfr_init_set_str(r30465358, "4", 10, MPFR_RNDN); mpfr_init_set_str(r30465359, "1", 10, MPFR_RNDN); mpfr_init(r30465360); mpfr_init(r30465361); mpfr_init_set_str(r30465362, "3", 10, MPFR_RNDN); mpfr_init(r30465363); mpfr_init(r30465364); mpfr_init(r30465365); mpfr_init(r30465366); mpfr_init(r30465367); mpfr_init(r30465368); mpfr_init(r30465369); mpfr_init(r30465370); mpfr_init(r30465371); } double f_fm(double a, double b) { mpfr_set_d(r30465351, a, MPFR_RNDN); mpfr_mul(r30465352, r30465351, r30465351, MPFR_RNDN); mpfr_set_d(r30465353, b, MPFR_RNDN); mpfr_mul(r30465354, r30465353, r30465353, MPFR_RNDN); mpfr_add(r30465355, r30465352, r30465354, MPFR_RNDN); ; mpfr_pow(r30465357, r30465355, r30465356, MPFR_RNDN); ; ; mpfr_add(r30465360, r30465359, r30465351, MPFR_RNDN); mpfr_mul(r30465361, r30465352, r30465360, MPFR_RNDN); ; mpfr_mul(r30465363, r30465362, r30465351, MPFR_RNDN); mpfr_sub(r30465364, r30465359, r30465363, MPFR_RNDN); mpfr_mul(r30465365, r30465354, r30465364, MPFR_RNDN); mpfr_add(r30465366, r30465361, r30465365, MPFR_RNDN); mpfr_mul(r30465367, r30465358, r30465366, MPFR_RNDN); mpfr_add(r30465368, r30465357, r30465367, MPFR_RNDN); mpfr_sqrt(r30465369, r30465368, MPFR_RNDN); mpfr_mul(r30465370, r30465369, r30465369, MPFR_RNDN); mpfr_sub(r30465371, r30465370, r30465359, MPFR_RNDN); return mpfr_get_d(r30465371, MPFR_RNDN); } static mpfr_t r30465372, r30465373, r30465374, r30465375, r30465376, r30465377, r30465378, r30465379, r30465380, r30465381, r30465382, r30465383, r30465384, r30465385, r30465386, r30465387, r30465388, r30465389, r30465390, r30465391, r30465392; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r30465372); mpfr_init(r30465373); mpfr_init(r30465374); mpfr_init(r30465375); mpfr_init(r30465376); mpfr_init_set_str(r30465377, "2", 10, MPFR_RNDN); mpfr_init(r30465378); mpfr_init_set_str(r30465379, "4", 10, MPFR_RNDN); mpfr_init_set_str(r30465380, "1", 10, MPFR_RNDN); mpfr_init(r30465381); mpfr_init(r30465382); mpfr_init_set_str(r30465383, "3", 10, MPFR_RNDN); mpfr_init(r30465384); mpfr_init(r30465385); mpfr_init(r30465386); mpfr_init(r30465387); mpfr_init(r30465388); mpfr_init(r30465389); mpfr_init(r30465390); mpfr_init(r30465391); mpfr_init(r30465392); } double f_dm(double a, double b) { mpfr_set_d(r30465372, a, MPFR_RNDN); mpfr_mul(r30465373, r30465372, r30465372, MPFR_RNDN); mpfr_set_d(r30465374, b, MPFR_RNDN); mpfr_mul(r30465375, r30465374, r30465374, MPFR_RNDN); mpfr_add(r30465376, r30465373, r30465375, MPFR_RNDN); ; mpfr_pow(r30465378, r30465376, r30465377, MPFR_RNDN); ; ; mpfr_add(r30465381, r30465380, r30465372, MPFR_RNDN); mpfr_mul(r30465382, r30465373, r30465381, MPFR_RNDN); ; mpfr_mul(r30465384, r30465383, r30465372, MPFR_RNDN); mpfr_sub(r30465385, r30465380, r30465384, MPFR_RNDN); mpfr_mul(r30465386, r30465375, r30465385, MPFR_RNDN); mpfr_add(r30465387, r30465382, r30465386, MPFR_RNDN); mpfr_mul(r30465388, r30465379, r30465387, MPFR_RNDN); mpfr_add(r30465389, r30465378, r30465388, MPFR_RNDN); mpfr_sqrt(r30465390, r30465389, MPFR_RNDN); mpfr_mul(r30465391, r30465390, r30465390, MPFR_RNDN); mpfr_sub(r30465392, r30465391, r30465380, MPFR_RNDN); return mpfr_get_d(r30465392, MPFR_RNDN); }