#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r17369 = a; float r17370 = r17369 * r17369; float r17371 = b; float r17372 = r17371 * r17371; float r17373 = r17370 + r17372; float r17374 = r17373 * r17373; float r17375 = 4.0f; float r17376 = 1.0f; float r17377 = r17376 + r17369; float r17378 = r17370 * r17377; float r17379 = 3.0f; float r17380 = r17379 * r17369; float r17381 = r17376 - r17380; float r17382 = r17372 * r17381; float r17383 = r17378 + r17382; float r17384 = r17375 * r17383; float r17385 = r17374 + r17384; float r17386 = r17385 - r17376; return r17386; } double f_id(double a, double b) { double r17387 = a; double r17388 = r17387 * r17387; double r17389 = b; double r17390 = r17389 * r17389; double r17391 = r17388 + r17390; double r17392 = r17391 * r17391; double r17393 = 4.0; double r17394 = 1.0; double r17395 = r17394 + r17387; double r17396 = r17388 * r17395; double r17397 = 3.0; double r17398 = r17397 * r17387; double r17399 = r17394 - r17398; double r17400 = r17390 * r17399; double r17401 = r17396 + r17400; double r17402 = r17393 * r17401; double r17403 = r17392 + r17402; double r17404 = r17403 - r17394; return r17404; } double f_of(float a, float b) { float r17405 = 1.0f; float r17406 = a; float r17407 = 3.0f; float r17408 = r17406 * r17407; float r17409 = r17405 - r17408; float r17410 = b; float r17411 = r17410 * r17410; float r17412 = fma(r17406, r17406, r17406); float r17413 = r17412 * r17406; float r17414 = fma(r17409, r17411, r17413); float r17415 = 4.0f; float r17416 = 2.0f; float r17417 = r17410 * r17416; float r17418 = r17417 * r17410; float r17419 = r17406 * r17406; float r17420 = pow(r17406, r17415); float r17421 = pow(r17410, r17415); float r17422 = r17420 + r17421; float r17423 = fma(r17418, r17419, r17422); float r17424 = fma(r17414, r17415, r17423); float r17425 = r17424 - r17405; return r17425; } double f_od(double a, double b) { double r17426 = 1.0; double r17427 = a; double r17428 = 3.0; double r17429 = r17427 * r17428; double r17430 = r17426 - r17429; double r17431 = b; double r17432 = r17431 * r17431; double r17433 = fma(r17427, r17427, r17427); double r17434 = r17433 * r17427; double r17435 = fma(r17430, r17432, r17434); double r17436 = 4.0; double r17437 = 2.0; double r17438 = r17431 * r17437; double r17439 = r17438 * r17431; double r17440 = r17427 * r17427; double r17441 = pow(r17427, r17436); double r17442 = pow(r17431, r17436); double r17443 = r17441 + r17442; double r17444 = fma(r17439, r17440, r17443); double r17445 = fma(r17435, r17436, r17444); double r17446 = r17445 - r17426; return r17446; } 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 r17447, r17448, r17449, r17450, r17451, r17452, r17453, r17454, r17455, r17456, r17457, r17458, r17459, r17460, r17461, r17462, r17463, r17464; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17447); mpfr_init(r17448); mpfr_init(r17449); mpfr_init(r17450); mpfr_init(r17451); mpfr_init(r17452); mpfr_init_set_str(r17453, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17454, "1", 10, MPFR_RNDN); mpfr_init(r17455); mpfr_init(r17456); mpfr_init_set_str(r17457, "3", 10, MPFR_RNDN); mpfr_init(r17458); mpfr_init(r17459); mpfr_init(r17460); mpfr_init(r17461); mpfr_init(r17462); mpfr_init(r17463); mpfr_init(r17464); } double f_im(double a, double b) { mpfr_set_d(r17447, a, MPFR_RNDN); mpfr_sqr(r17448, r17447, MPFR_RNDN); mpfr_set_d(r17449, b, MPFR_RNDN); mpfr_sqr(r17450, r17449, MPFR_RNDN); mpfr_add(r17451, r17448, r17450, MPFR_RNDN); mpfr_sqr(r17452, r17451, MPFR_RNDN); ; ; mpfr_add(r17455, r17454, r17447, MPFR_RNDN); mpfr_mul(r17456, r17448, r17455, MPFR_RNDN); ; mpfr_mul(r17458, r17457, r17447, MPFR_RNDN); mpfr_sub(r17459, r17454, r17458, MPFR_RNDN); mpfr_mul(r17460, r17450, r17459, MPFR_RNDN); mpfr_add(r17461, r17456, r17460, MPFR_RNDN); mpfr_mul(r17462, r17453, r17461, MPFR_RNDN); mpfr_add(r17463, r17452, r17462, MPFR_RNDN); mpfr_sub(r17464, r17463, r17454, MPFR_RNDN); return mpfr_get_d(r17464, MPFR_RNDN); } static mpfr_t r17465, r17466, r17467, r17468, r17469, r17470, r17471, r17472, r17473, r17474, r17475, r17476, r17477, r17478, r17479, r17480, r17481, r17482, r17483, r17484, r17485; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17465, "1", 10, MPFR_RNDN); mpfr_init(r17466); mpfr_init_set_str(r17467, "3", 10, MPFR_RNDN); mpfr_init(r17468); mpfr_init(r17469); mpfr_init(r17470); mpfr_init(r17471); mpfr_init(r17472); mpfr_init(r17473); mpfr_init(r17474); mpfr_init_set_str(r17475, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17476, "2", 10, MPFR_RNDN); mpfr_init(r17477); mpfr_init(r17478); mpfr_init(r17479); mpfr_init(r17480); mpfr_init(r17481); mpfr_init(r17482); mpfr_init(r17483); mpfr_init(r17484); mpfr_init(r17485); } double f_fm(double a, double b) { ; mpfr_set_d(r17466, a, MPFR_RNDN); ; mpfr_mul(r17468, r17466, r17467, MPFR_RNDN); mpfr_sub(r17469, r17465, r17468, MPFR_RNDN); mpfr_set_d(r17470, b, MPFR_RNDN); mpfr_sqr(r17471, r17470, MPFR_RNDN); mpfr_fma(r17472, r17466, r17466, r17466, MPFR_RNDN); mpfr_mul(r17473, r17472, r17466, MPFR_RNDN); mpfr_fma(r17474, r17469, r17471, r17473, MPFR_RNDN); ; ; mpfr_mul(r17477, r17470, r17476, MPFR_RNDN); mpfr_mul(r17478, r17477, r17470, MPFR_RNDN); mpfr_mul(r17479, r17466, r17466, MPFR_RNDN); mpfr_pow(r17480, r17466, r17475, MPFR_RNDN); mpfr_pow(r17481, r17470, r17475, MPFR_RNDN); mpfr_add(r17482, r17480, r17481, MPFR_RNDN); mpfr_fma(r17483, r17478, r17479, r17482, MPFR_RNDN); mpfr_fma(r17484, r17474, r17475, r17483, MPFR_RNDN); mpfr_sub(r17485, r17484, r17465, MPFR_RNDN); return mpfr_get_d(r17485, MPFR_RNDN); } static mpfr_t r17486, r17487, r17488, r17489, r17490, r17491, r17492, r17493, r17494, r17495, r17496, r17497, r17498, r17499, r17500, r17501, r17502, r17503, r17504, r17505, r17506; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17486, "1", 10, MPFR_RNDN); mpfr_init(r17487); mpfr_init_set_str(r17488, "3", 10, MPFR_RNDN); mpfr_init(r17489); mpfr_init(r17490); mpfr_init(r17491); mpfr_init(r17492); mpfr_init(r17493); mpfr_init(r17494); mpfr_init(r17495); mpfr_init_set_str(r17496, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17497, "2", 10, MPFR_RNDN); mpfr_init(r17498); mpfr_init(r17499); mpfr_init(r17500); mpfr_init(r17501); mpfr_init(r17502); mpfr_init(r17503); mpfr_init(r17504); mpfr_init(r17505); mpfr_init(r17506); } double f_dm(double a, double b) { ; mpfr_set_d(r17487, a, MPFR_RNDN); ; mpfr_mul(r17489, r17487, r17488, MPFR_RNDN); mpfr_sub(r17490, r17486, r17489, MPFR_RNDN); mpfr_set_d(r17491, b, MPFR_RNDN); mpfr_sqr(r17492, r17491, MPFR_RNDN); mpfr_fma(r17493, r17487, r17487, r17487, MPFR_RNDN); mpfr_mul(r17494, r17493, r17487, MPFR_RNDN); mpfr_fma(r17495, r17490, r17492, r17494, MPFR_RNDN); ; ; mpfr_mul(r17498, r17491, r17497, MPFR_RNDN); mpfr_mul(r17499, r17498, r17491, MPFR_RNDN); mpfr_mul(r17500, r17487, r17487, MPFR_RNDN); mpfr_pow(r17501, r17487, r17496, MPFR_RNDN); mpfr_pow(r17502, r17491, r17496, MPFR_RNDN); mpfr_add(r17503, r17501, r17502, MPFR_RNDN); mpfr_fma(r17504, r17499, r17500, r17503, MPFR_RNDN); mpfr_fma(r17505, r17495, r17496, r17504, MPFR_RNDN); mpfr_sub(r17506, r17505, r17486, MPFR_RNDN); return mpfr_get_d(r17506, MPFR_RNDN); }