#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (26)"; double f_if(float a, float b) { float r24354 = a; float r24355 = r24354 * r24354; float r24356 = b; float r24357 = r24356 * r24356; float r24358 = r24355 + r24357; float r24359 = 2; float r24360 = pow(r24358, r24359); float r24361 = 4; float r24362 = r24361 * r24357; float r24363 = r24360 + r24362; float r24364 = 1; float r24365 = r24363 - r24364; return r24365; } double f_id(double a, double b) { double r24366 = a; double r24367 = r24366 * r24366; double r24368 = b; double r24369 = r24368 * r24368; double r24370 = r24367 + r24369; double r24371 = 2; double r24372 = pow(r24370, r24371); double r24373 = 4; double r24374 = r24373 * r24369; double r24375 = r24372 + r24374; double r24376 = 1; double r24377 = r24375 - r24376; return r24377; } double f_of(float a, float b) { float r24378 = b; float r24379 = r24378 * r24378; float r24380 = a; float r24381 = r24380 * r24380; float r24382 = r24379 + r24381; float r24383 = sqrt(r24382); float r24384 = cbrt(r24383); float r24385 = r24382 * r24384; float r24386 = r24385 * r24384; float r24387 = r24381 + r24379; float r24388 = cbrt(r24387); float r24389 = 2; float r24390 = pow(r24388, r24389); float r24391 = r24386 * r24390; float r24392 = 4; float r24393 = r24392 * r24379; float r24394 = r24391 + r24393; float r24395 = 1; float r24396 = r24394 - r24395; return r24396; } double f_od(double a, double b) { double r24397 = b; double r24398 = r24397 * r24397; double r24399 = a; double r24400 = r24399 * r24399; double r24401 = r24398 + r24400; double r24402 = sqrt(r24401); double r24403 = cbrt(r24402); double r24404 = r24401 * r24403; double r24405 = r24404 * r24403; double r24406 = r24400 + r24398; double r24407 = cbrt(r24406); double r24408 = 2; double r24409 = pow(r24407, r24408); double r24410 = r24405 * r24409; double r24411 = 4; double r24412 = r24411 * r24398; double r24413 = r24410 + r24412; double r24414 = 1; double r24415 = r24413 - r24414; return r24415; } 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 r24416, r24417, r24418, r24419, r24420, r24421, r24422, r24423, r24424, r24425, r24426, r24427; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r24416); mpfr_init(r24417); mpfr_init(r24418); mpfr_init(r24419); mpfr_init(r24420); mpfr_init_set_str(r24421, "2", 10, MPFR_RNDN); mpfr_init(r24422); mpfr_init_set_str(r24423, "4", 10, MPFR_RNDN); mpfr_init(r24424); mpfr_init(r24425); mpfr_init_set_str(r24426, "1", 10, MPFR_RNDN); mpfr_init(r24427); } double f_im(double a, double b) { mpfr_set_d(r24416, a, MPFR_RNDN); mpfr_mul(r24417, r24416, r24416, MPFR_RNDN); mpfr_set_d(r24418, b, MPFR_RNDN); mpfr_mul(r24419, r24418, r24418, MPFR_RNDN); mpfr_add(r24420, r24417, r24419, MPFR_RNDN); ; mpfr_pow(r24422, r24420, r24421, MPFR_RNDN); ; mpfr_mul(r24424, r24423, r24419, MPFR_RNDN); mpfr_add(r24425, r24422, r24424, MPFR_RNDN); ; mpfr_sub(r24427, r24425, r24426, MPFR_RNDN); return mpfr_get_d(r24427, MPFR_RNDN); } static mpfr_t r24428, r24429, r24430, r24431, r24432, r24433, r24434, r24435, r24436, r24437, r24438, r24439, r24440, r24441, r24442, r24443, r24444, r24445, r24446; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r24428); mpfr_init(r24429); mpfr_init(r24430); mpfr_init(r24431); mpfr_init(r24432); mpfr_init(r24433); mpfr_init(r24434); mpfr_init(r24435); mpfr_init(r24436); mpfr_init(r24437); mpfr_init(r24438); mpfr_init_set_str(r24439, "2", 10, MPFR_RNDN); mpfr_init(r24440); mpfr_init(r24441); mpfr_init_set_str(r24442, "4", 10, MPFR_RNDN); mpfr_init(r24443); mpfr_init(r24444); mpfr_init_set_str(r24445, "1", 10, MPFR_RNDN); mpfr_init(r24446); } double f_fm(double a, double b) { mpfr_set_d(r24428, b, MPFR_RNDN); mpfr_mul(r24429, r24428, r24428, MPFR_RNDN); mpfr_set_d(r24430, a, MPFR_RNDN); mpfr_mul(r24431, r24430, r24430, MPFR_RNDN); mpfr_add(r24432, r24429, r24431, MPFR_RNDN); mpfr_sqrt(r24433, r24432, MPFR_RNDN); mpfr_cbrt(r24434, r24433, MPFR_RNDN); mpfr_mul(r24435, r24432, r24434, MPFR_RNDN); mpfr_mul(r24436, r24435, r24434, MPFR_RNDN); mpfr_add(r24437, r24431, r24429, MPFR_RNDN); mpfr_cbrt(r24438, r24437, MPFR_RNDN); ; mpfr_pow(r24440, r24438, r24439, MPFR_RNDN); mpfr_mul(r24441, r24436, r24440, MPFR_RNDN); ; mpfr_mul(r24443, r24442, r24429, MPFR_RNDN); mpfr_add(r24444, r24441, r24443, MPFR_RNDN); ; mpfr_sub(r24446, r24444, r24445, MPFR_RNDN); return mpfr_get_d(r24446, MPFR_RNDN); } static mpfr_t r24447, r24448, r24449, r24450, r24451, r24452, r24453, r24454, r24455, r24456, r24457, r24458, r24459, r24460, r24461, r24462, r24463, r24464, r24465; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r24447); mpfr_init(r24448); mpfr_init(r24449); mpfr_init(r24450); mpfr_init(r24451); mpfr_init(r24452); mpfr_init(r24453); mpfr_init(r24454); mpfr_init(r24455); mpfr_init(r24456); mpfr_init(r24457); mpfr_init_set_str(r24458, "2", 10, MPFR_RNDN); mpfr_init(r24459); mpfr_init(r24460); mpfr_init_set_str(r24461, "4", 10, MPFR_RNDN); mpfr_init(r24462); mpfr_init(r24463); mpfr_init_set_str(r24464, "1", 10, MPFR_RNDN); mpfr_init(r24465); } double f_dm(double a, double b) { mpfr_set_d(r24447, b, MPFR_RNDN); mpfr_mul(r24448, r24447, r24447, MPFR_RNDN); mpfr_set_d(r24449, a, MPFR_RNDN); mpfr_mul(r24450, r24449, r24449, MPFR_RNDN); mpfr_add(r24451, r24448, r24450, MPFR_RNDN); mpfr_sqrt(r24452, r24451, MPFR_RNDN); mpfr_cbrt(r24453, r24452, MPFR_RNDN); mpfr_mul(r24454, r24451, r24453, MPFR_RNDN); mpfr_mul(r24455, r24454, r24453, MPFR_RNDN); mpfr_add(r24456, r24450, r24448, MPFR_RNDN); mpfr_cbrt(r24457, r24456, MPFR_RNDN); ; mpfr_pow(r24459, r24457, r24458, MPFR_RNDN); mpfr_mul(r24460, r24455, r24459, MPFR_RNDN); ; mpfr_mul(r24462, r24461, r24448, MPFR_RNDN); mpfr_add(r24463, r24460, r24462, MPFR_RNDN); ; mpfr_sub(r24465, r24463, r24464, MPFR_RNDN); return mpfr_get_d(r24465, MPFR_RNDN); }