#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r17409 = a; float r17410 = r17409 * r17409; float r17411 = b; float r17412 = r17411 * r17411; float r17413 = r17410 + r17412; float r17414 = r17413 * r17413; float r17415 = 4.0f; float r17416 = 1.0f; float r17417 = r17416 + r17409; float r17418 = r17410 * r17417; float r17419 = 3.0f; float r17420 = r17419 * r17409; float r17421 = r17416 - r17420; float r17422 = r17412 * r17421; float r17423 = r17418 + r17422; float r17424 = r17415 * r17423; float r17425 = r17414 + r17424; float r17426 = r17425 - r17416; return r17426; } double f_id(double a, double b) { double r17427 = a; double r17428 = r17427 * r17427; double r17429 = b; double r17430 = r17429 * r17429; double r17431 = r17428 + r17430; double r17432 = r17431 * r17431; double r17433 = 4.0; double r17434 = 1.0; double r17435 = r17434 + r17427; double r17436 = r17428 * r17435; double r17437 = 3.0; double r17438 = r17437 * r17427; double r17439 = r17434 - r17438; double r17440 = r17430 * r17439; double r17441 = r17436 + r17440; double r17442 = r17433 * r17441; double r17443 = r17432 + r17442; double r17444 = r17443 - r17434; return r17444; } double f_of(float a, float b) { float r17445 = 1.0f; float r17446 = a; float r17447 = 3.0f; float r17448 = r17446 * r17447; float r17449 = r17445 - r17448; float r17450 = b; float r17451 = r17450 * r17450; float r17452 = fma(r17446, r17446, r17446); float r17453 = r17452 * r17446; float r17454 = fma(r17449, r17451, r17453); float r17455 = 4.0f; float r17456 = r17454 * r17455; float r17457 = 2.0f; float r17458 = r17450 * r17457; float r17459 = r17458 * r17450; float r17460 = r17446 * r17446; float r17461 = pow(r17446, r17455); float r17462 = pow(r17450, r17455); float r17463 = r17461 + r17462; float r17464 = fma(r17459, r17460, r17463); float r17465 = r17456 + r17464; float r17466 = r17465 - r17445; return r17466; } double f_od(double a, double b) { double r17467 = 1.0; double r17468 = a; double r17469 = 3.0; double r17470 = r17468 * r17469; double r17471 = r17467 - r17470; double r17472 = b; double r17473 = r17472 * r17472; double r17474 = fma(r17468, r17468, r17468); double r17475 = r17474 * r17468; double r17476 = fma(r17471, r17473, r17475); double r17477 = 4.0; double r17478 = r17476 * r17477; double r17479 = 2.0; double r17480 = r17472 * r17479; double r17481 = r17480 * r17472; double r17482 = r17468 * r17468; double r17483 = pow(r17468, r17477); double r17484 = pow(r17472, r17477); double r17485 = r17483 + r17484; double r17486 = fma(r17481, r17482, r17485); double r17487 = r17478 + r17486; double r17488 = r17487 - r17467; return r17488; } 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 r17489, r17490, r17491, r17492, r17493, r17494, r17495, r17496, r17497, r17498, r17499, r17500, r17501, r17502, r17503, r17504, r17505, r17506; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17489); mpfr_init(r17490); mpfr_init(r17491); mpfr_init(r17492); mpfr_init(r17493); mpfr_init(r17494); mpfr_init_set_str(r17495, "4", 10, MPFR_RNDN); mpfr_init_set_str(r17496, "1", 10, MPFR_RNDN); mpfr_init(r17497); mpfr_init(r17498); mpfr_init_set_str(r17499, "3", 10, MPFR_RNDN); mpfr_init(r17500); mpfr_init(r17501); mpfr_init(r17502); mpfr_init(r17503); mpfr_init(r17504); mpfr_init(r17505); mpfr_init(r17506); } double f_im(double a, double b) { mpfr_set_d(r17489, a, MPFR_RNDN); mpfr_sqr(r17490, r17489, MPFR_RNDN); mpfr_set_d(r17491, b, MPFR_RNDN); mpfr_sqr(r17492, r17491, MPFR_RNDN); mpfr_add(r17493, r17490, r17492, MPFR_RNDN); mpfr_sqr(r17494, r17493, MPFR_RNDN); ; ; mpfr_add(r17497, r17496, r17489, MPFR_RNDN); mpfr_mul(r17498, r17490, r17497, MPFR_RNDN); ; mpfr_mul(r17500, r17499, r17489, MPFR_RNDN); mpfr_sub(r17501, r17496, r17500, MPFR_RNDN); mpfr_mul(r17502, r17492, r17501, MPFR_RNDN); mpfr_add(r17503, r17498, r17502, MPFR_RNDN); mpfr_mul(r17504, r17495, r17503, MPFR_RNDN); mpfr_add(r17505, r17494, r17504, MPFR_RNDN); mpfr_sub(r17506, r17505, r17496, MPFR_RNDN); return mpfr_get_d(r17506, MPFR_RNDN); } static mpfr_t r17507, r17508, r17509, r17510, r17511, r17512, r17513, r17514, r17515, r17516, r17517, r17518, r17519, r17520, r17521, r17522, r17523, r17524, r17525, r17526, r17527, r17528; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17507, "1", 10, MPFR_RNDN); mpfr_init(r17508); mpfr_init_set_str(r17509, "3", 10, MPFR_RNDN); mpfr_init(r17510); mpfr_init(r17511); mpfr_init(r17512); mpfr_init(r17513); mpfr_init(r17514); mpfr_init(r17515); mpfr_init(r17516); mpfr_init_set_str(r17517, "4", 10, MPFR_RNDN); mpfr_init(r17518); mpfr_init_set_str(r17519, "2", 10, MPFR_RNDN); mpfr_init(r17520); mpfr_init(r17521); mpfr_init(r17522); mpfr_init(r17523); mpfr_init(r17524); mpfr_init(r17525); mpfr_init(r17526); mpfr_init(r17527); mpfr_init(r17528); } double f_fm(double a, double b) { ; mpfr_set_d(r17508, a, MPFR_RNDN); ; mpfr_mul(r17510, r17508, r17509, MPFR_RNDN); mpfr_sub(r17511, r17507, r17510, MPFR_RNDN); mpfr_set_d(r17512, b, MPFR_RNDN); mpfr_sqr(r17513, r17512, MPFR_RNDN); mpfr_fma(r17514, r17508, r17508, r17508, MPFR_RNDN); mpfr_mul(r17515, r17514, r17508, MPFR_RNDN); mpfr_fma(r17516, r17511, r17513, r17515, MPFR_RNDN); ; mpfr_mul(r17518, r17516, r17517, MPFR_RNDN); ; mpfr_mul(r17520, r17512, r17519, MPFR_RNDN); mpfr_mul(r17521, r17520, r17512, MPFR_RNDN); mpfr_mul(r17522, r17508, r17508, MPFR_RNDN); mpfr_pow(r17523, r17508, r17517, MPFR_RNDN); mpfr_pow(r17524, r17512, r17517, MPFR_RNDN); mpfr_add(r17525, r17523, r17524, MPFR_RNDN); mpfr_fma(r17526, r17521, r17522, r17525, MPFR_RNDN); mpfr_add(r17527, r17518, r17526, MPFR_RNDN); mpfr_sub(r17528, r17527, r17507, MPFR_RNDN); return mpfr_get_d(r17528, MPFR_RNDN); } static mpfr_t r17529, r17530, r17531, r17532, r17533, r17534, r17535, r17536, r17537, r17538, r17539, r17540, r17541, r17542, r17543, r17544, r17545, r17546, r17547, r17548, r17549, r17550; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17529, "1", 10, MPFR_RNDN); mpfr_init(r17530); mpfr_init_set_str(r17531, "3", 10, MPFR_RNDN); mpfr_init(r17532); mpfr_init(r17533); mpfr_init(r17534); mpfr_init(r17535); mpfr_init(r17536); mpfr_init(r17537); mpfr_init(r17538); mpfr_init_set_str(r17539, "4", 10, MPFR_RNDN); mpfr_init(r17540); mpfr_init_set_str(r17541, "2", 10, MPFR_RNDN); mpfr_init(r17542); mpfr_init(r17543); mpfr_init(r17544); mpfr_init(r17545); mpfr_init(r17546); mpfr_init(r17547); mpfr_init(r17548); mpfr_init(r17549); mpfr_init(r17550); } double f_dm(double a, double b) { ; mpfr_set_d(r17530, a, MPFR_RNDN); ; mpfr_mul(r17532, r17530, r17531, MPFR_RNDN); mpfr_sub(r17533, r17529, r17532, MPFR_RNDN); mpfr_set_d(r17534, b, MPFR_RNDN); mpfr_sqr(r17535, r17534, MPFR_RNDN); mpfr_fma(r17536, r17530, r17530, r17530, MPFR_RNDN); mpfr_mul(r17537, r17536, r17530, MPFR_RNDN); mpfr_fma(r17538, r17533, r17535, r17537, MPFR_RNDN); ; mpfr_mul(r17540, r17538, r17539, MPFR_RNDN); ; mpfr_mul(r17542, r17534, r17541, MPFR_RNDN); mpfr_mul(r17543, r17542, r17534, MPFR_RNDN); mpfr_mul(r17544, r17530, r17530, MPFR_RNDN); mpfr_pow(r17545, r17530, r17539, MPFR_RNDN); mpfr_pow(r17546, r17534, r17539, MPFR_RNDN); mpfr_add(r17547, r17545, r17546, MPFR_RNDN); mpfr_fma(r17548, r17543, r17544, r17547, MPFR_RNDN); mpfr_add(r17549, r17540, r17548, MPFR_RNDN); mpfr_sub(r17550, r17549, r17529, MPFR_RNDN); return mpfr_get_d(r17550, MPFR_RNDN); }