#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r27377 = a; float r27378 = r27377 * r27377; float r27379 = b; float r27380 = r27379 * r27379; float r27381 = r27378 + r27380; float r27382 = 2; float r27383 = pow(r27381, r27382); float r27384 = 4; float r27385 = 1; float r27386 = r27385 + r27377; float r27387 = r27378 * r27386; float r27388 = 3; float r27389 = r27388 * r27377; float r27390 = r27385 - r27389; float r27391 = r27380 * r27390; float r27392 = r27387 + r27391; float r27393 = r27384 * r27392; float r27394 = r27383 + r27393; float r27395 = r27394 - r27385; return r27395; } double f_id(double a, double b) { double r27396 = a; double r27397 = r27396 * r27396; double r27398 = b; double r27399 = r27398 * r27398; double r27400 = r27397 + r27399; double r27401 = 2; double r27402 = pow(r27400, r27401); double r27403 = 4; double r27404 = 1; double r27405 = r27404 + r27396; double r27406 = r27397 * r27405; double r27407 = 3; double r27408 = r27407 * r27396; double r27409 = r27404 - r27408; double r27410 = r27399 * r27409; double r27411 = r27406 + r27410; double r27412 = r27403 * r27411; double r27413 = r27402 + r27412; double r27414 = r27413 - r27404; return r27414; } double f_of(float a, float b) { float r27415 = b; float r27416 = a; float r27417 = r27416 * r27416; float r27418 = fma(r27415, r27415, r27417); float r27419 = fma(r27416, r27416, r27416); float r27420 = r27415 * r27415; float r27421 = fma(r27416, r27419, r27420); float r27422 = 4; float r27423 = 1; float r27424 = -r27423; float r27425 = fma(r27421, r27422, r27424); float r27426 = fma(r27418, r27418, r27425); return r27426; } double f_od(double a, double b) { double r27427 = b; double r27428 = a; double r27429 = r27428 * r27428; double r27430 = fma(r27427, r27427, r27429); double r27431 = fma(r27428, r27428, r27428); double r27432 = r27427 * r27427; double r27433 = fma(r27428, r27431, r27432); double r27434 = 4; double r27435 = 1; double r27436 = -r27435; double r27437 = fma(r27433, r27434, r27436); double r27438 = fma(r27430, r27430, r27437); return r27438; } 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 r27439, r27440, r27441, r27442, r27443, r27444, r27445, r27446, r27447, r27448, r27449, r27450, r27451, r27452, r27453, r27454, r27455, r27456, r27457; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r27439); mpfr_init(r27440); mpfr_init(r27441); mpfr_init(r27442); mpfr_init(r27443); mpfr_init_set_str(r27444, "2", 10, MPFR_RNDN); mpfr_init(r27445); mpfr_init_set_str(r27446, "4", 10, MPFR_RNDN); mpfr_init_set_str(r27447, "1", 10, MPFR_RNDN); mpfr_init(r27448); mpfr_init(r27449); mpfr_init_set_str(r27450, "3", 10, MPFR_RNDN); mpfr_init(r27451); mpfr_init(r27452); mpfr_init(r27453); mpfr_init(r27454); mpfr_init(r27455); mpfr_init(r27456); mpfr_init(r27457); } double f_im(double a, double b) { mpfr_set_d(r27439, a, MPFR_RNDN); mpfr_mul(r27440, r27439, r27439, MPFR_RNDN); mpfr_set_d(r27441, b, MPFR_RNDN); mpfr_mul(r27442, r27441, r27441, MPFR_RNDN); mpfr_add(r27443, r27440, r27442, MPFR_RNDN); ; mpfr_pow(r27445, r27443, r27444, MPFR_RNDN); ; ; mpfr_add(r27448, r27447, r27439, MPFR_RNDN); mpfr_mul(r27449, r27440, r27448, MPFR_RNDN); ; mpfr_mul(r27451, r27450, r27439, MPFR_RNDN); mpfr_sub(r27452, r27447, r27451, MPFR_RNDN); mpfr_mul(r27453, r27442, r27452, MPFR_RNDN); mpfr_add(r27454, r27449, r27453, MPFR_RNDN); mpfr_mul(r27455, r27446, r27454, MPFR_RNDN); mpfr_add(r27456, r27445, r27455, MPFR_RNDN); mpfr_sub(r27457, r27456, r27447, MPFR_RNDN); return mpfr_get_d(r27457, MPFR_RNDN); } static mpfr_t r27458, r27459, r27460, r27461, r27462, r27463, r27464, r27465, r27466, r27467, r27468, r27469; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r27458); mpfr_init(r27459); mpfr_init(r27460); mpfr_init(r27461); mpfr_init(r27462); mpfr_init(r27463); mpfr_init(r27464); mpfr_init_set_str(r27465, "4", 10, MPFR_RNDN); mpfr_init_set_str(r27466, "1", 10, MPFR_RNDN); mpfr_init(r27467); mpfr_init(r27468); mpfr_init(r27469); } double f_fm(double a, double b) { mpfr_set_d(r27458, b, MPFR_RNDN); mpfr_set_d(r27459, a, MPFR_RNDN); mpfr_mul(r27460, r27459, r27459, MPFR_RNDN); mpfr_fma(r27461, r27458, r27458, r27460, MPFR_RNDN); mpfr_fma(r27462, r27459, r27459, r27459, MPFR_RNDN); mpfr_mul(r27463, r27458, r27458, MPFR_RNDN); mpfr_fma(r27464, r27459, r27462, r27463, MPFR_RNDN); ; ; mpfr_neg(r27467, r27466, MPFR_RNDN); mpfr_fma(r27468, r27464, r27465, r27467, MPFR_RNDN); mpfr_fma(r27469, r27461, r27461, r27468, MPFR_RNDN); return mpfr_get_d(r27469, MPFR_RNDN); } static mpfr_t r27470, r27471, r27472, r27473, r27474, r27475, r27476, r27477, r27478, r27479, r27480, r27481; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r27470); mpfr_init(r27471); mpfr_init(r27472); mpfr_init(r27473); mpfr_init(r27474); mpfr_init(r27475); mpfr_init(r27476); mpfr_init_set_str(r27477, "4", 10, MPFR_RNDN); mpfr_init_set_str(r27478, "1", 10, MPFR_RNDN); mpfr_init(r27479); mpfr_init(r27480); mpfr_init(r27481); } double f_dm(double a, double b) { mpfr_set_d(r27470, b, MPFR_RNDN); mpfr_set_d(r27471, a, MPFR_RNDN); mpfr_mul(r27472, r27471, r27471, MPFR_RNDN); mpfr_fma(r27473, r27470, r27470, r27472, MPFR_RNDN); mpfr_fma(r27474, r27471, r27471, r27471, MPFR_RNDN); mpfr_mul(r27475, r27470, r27470, MPFR_RNDN); mpfr_fma(r27476, r27471, r27474, r27475, MPFR_RNDN); ; ; mpfr_neg(r27479, r27478, MPFR_RNDN); mpfr_fma(r27480, r27476, r27477, r27479, MPFR_RNDN); mpfr_fma(r27481, r27473, r27473, r27480, MPFR_RNDN); return mpfr_get_d(r27481, MPFR_RNDN); }