#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r18045 = a; float r18046 = r18045 * r18045; float r18047 = b; float r18048 = r18047 * r18047; float r18049 = r18046 + r18048; float r18050 = r18049 * r18049; float r18051 = 4.0f; float r18052 = 1.0f; float r18053 = r18052 + r18045; float r18054 = r18046 * r18053; float r18055 = 3.0f; float r18056 = r18055 * r18045; float r18057 = r18052 - r18056; float r18058 = r18048 * r18057; float r18059 = r18054 + r18058; float r18060 = r18051 * r18059; float r18061 = r18050 + r18060; float r18062 = r18061 - r18052; return r18062; } double f_id(double a, double b) { double r18063 = a; double r18064 = r18063 * r18063; double r18065 = b; double r18066 = r18065 * r18065; double r18067 = r18064 + r18066; double r18068 = r18067 * r18067; double r18069 = 4.0; double r18070 = 1.0; double r18071 = r18070 + r18063; double r18072 = r18064 * r18071; double r18073 = 3.0; double r18074 = r18073 * r18063; double r18075 = r18070 - r18074; double r18076 = r18066 * r18075; double r18077 = r18072 + r18076; double r18078 = r18069 * r18077; double r18079 = r18068 + r18078; double r18080 = r18079 - r18070; return r18080; } double f_of(float a, float b) { float r18081 = 1.0f; float r18082 = a; float r18083 = 3.0f; float r18084 = r18082 * r18083; float r18085 = r18081 - r18084; float r18086 = b; float r18087 = r18086 * r18086; float r18088 = fma(r18082, r18082, r18082); float r18089 = r18088 * r18082; float r18090 = fma(r18085, r18087, r18089); float r18091 = 4.0f; float r18092 = 2.0f; float r18093 = r18086 * r18092; float r18094 = r18093 * r18086; float r18095 = r18082 * r18082; float r18096 = pow(r18082, r18091); float r18097 = pow(r18086, r18091); float r18098 = r18096 + r18097; float r18099 = fma(r18094, r18095, r18098); float r18100 = fma(r18090, r18091, r18099); float r18101 = r18100 - r18081; return r18101; } double f_od(double a, double b) { double r18102 = 1.0; double r18103 = a; double r18104 = 3.0; double r18105 = r18103 * r18104; double r18106 = r18102 - r18105; double r18107 = b; double r18108 = r18107 * r18107; double r18109 = fma(r18103, r18103, r18103); double r18110 = r18109 * r18103; double r18111 = fma(r18106, r18108, r18110); double r18112 = 4.0; double r18113 = 2.0; double r18114 = r18107 * r18113; double r18115 = r18114 * r18107; double r18116 = r18103 * r18103; double r18117 = pow(r18103, r18112); double r18118 = pow(r18107, r18112); double r18119 = r18117 + r18118; double r18120 = fma(r18115, r18116, r18119); double r18121 = fma(r18111, r18112, r18120); double r18122 = r18121 - r18102; return r18122; } 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 r18123, r18124, r18125, r18126, r18127, r18128, r18129, r18130, r18131, r18132, r18133, r18134, r18135, r18136, r18137, r18138, r18139, r18140; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r18123); mpfr_init(r18124); mpfr_init(r18125); mpfr_init(r18126); mpfr_init(r18127); mpfr_init(r18128); mpfr_init_set_str(r18129, "4", 10, MPFR_RNDN); mpfr_init_set_str(r18130, "1", 10, MPFR_RNDN); mpfr_init(r18131); mpfr_init(r18132); mpfr_init_set_str(r18133, "3", 10, MPFR_RNDN); mpfr_init(r18134); mpfr_init(r18135); mpfr_init(r18136); mpfr_init(r18137); mpfr_init(r18138); mpfr_init(r18139); mpfr_init(r18140); } double f_im(double a, double b) { mpfr_set_d(r18123, a, MPFR_RNDN); mpfr_sqr(r18124, r18123, MPFR_RNDN); mpfr_set_d(r18125, b, MPFR_RNDN); mpfr_sqr(r18126, r18125, MPFR_RNDN); mpfr_add(r18127, r18124, r18126, MPFR_RNDN); mpfr_sqr(r18128, r18127, MPFR_RNDN); ; ; mpfr_add(r18131, r18130, r18123, MPFR_RNDN); mpfr_mul(r18132, r18124, r18131, MPFR_RNDN); ; mpfr_mul(r18134, r18133, r18123, MPFR_RNDN); mpfr_sub(r18135, r18130, r18134, MPFR_RNDN); mpfr_mul(r18136, r18126, r18135, MPFR_RNDN); mpfr_add(r18137, r18132, r18136, MPFR_RNDN); mpfr_mul(r18138, r18129, r18137, MPFR_RNDN); mpfr_add(r18139, r18128, r18138, MPFR_RNDN); mpfr_sub(r18140, r18139, r18130, MPFR_RNDN); return mpfr_get_d(r18140, MPFR_RNDN); } static mpfr_t r18141, r18142, r18143, r18144, r18145, r18146, r18147, r18148, r18149, r18150, r18151, r18152, r18153, r18154, r18155, r18156, r18157, r18158, r18159, r18160, r18161; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18141, "1", 10, MPFR_RNDN); mpfr_init(r18142); mpfr_init_set_str(r18143, "3", 10, MPFR_RNDN); mpfr_init(r18144); mpfr_init(r18145); mpfr_init(r18146); mpfr_init(r18147); mpfr_init(r18148); mpfr_init(r18149); mpfr_init(r18150); mpfr_init_set_str(r18151, "4", 10, MPFR_RNDN); mpfr_init_set_str(r18152, "2", 10, MPFR_RNDN); mpfr_init(r18153); mpfr_init(r18154); mpfr_init(r18155); mpfr_init(r18156); mpfr_init(r18157); mpfr_init(r18158); mpfr_init(r18159); mpfr_init(r18160); mpfr_init(r18161); } double f_fm(double a, double b) { ; mpfr_set_d(r18142, a, MPFR_RNDN); ; mpfr_mul(r18144, r18142, r18143, MPFR_RNDN); mpfr_sub(r18145, r18141, r18144, MPFR_RNDN); mpfr_set_d(r18146, b, MPFR_RNDN); mpfr_sqr(r18147, r18146, MPFR_RNDN); mpfr_fma(r18148, r18142, r18142, r18142, MPFR_RNDN); mpfr_mul(r18149, r18148, r18142, MPFR_RNDN); mpfr_fma(r18150, r18145, r18147, r18149, MPFR_RNDN); ; ; mpfr_mul(r18153, r18146, r18152, MPFR_RNDN); mpfr_mul(r18154, r18153, r18146, MPFR_RNDN); mpfr_mul(r18155, r18142, r18142, MPFR_RNDN); mpfr_pow(r18156, r18142, r18151, MPFR_RNDN); mpfr_pow(r18157, r18146, r18151, MPFR_RNDN); mpfr_add(r18158, r18156, r18157, MPFR_RNDN); mpfr_fma(r18159, r18154, r18155, r18158, MPFR_RNDN); mpfr_fma(r18160, r18150, r18151, r18159, MPFR_RNDN); mpfr_sub(r18161, r18160, r18141, MPFR_RNDN); return mpfr_get_d(r18161, MPFR_RNDN); } static mpfr_t r18162, r18163, r18164, r18165, r18166, r18167, r18168, r18169, r18170, r18171, r18172, r18173, r18174, r18175, r18176, r18177, r18178, r18179, r18180, r18181, r18182; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18162, "1", 10, MPFR_RNDN); mpfr_init(r18163); mpfr_init_set_str(r18164, "3", 10, MPFR_RNDN); mpfr_init(r18165); mpfr_init(r18166); mpfr_init(r18167); mpfr_init(r18168); mpfr_init(r18169); mpfr_init(r18170); mpfr_init(r18171); mpfr_init_set_str(r18172, "4", 10, MPFR_RNDN); mpfr_init_set_str(r18173, "2", 10, MPFR_RNDN); mpfr_init(r18174); mpfr_init(r18175); mpfr_init(r18176); mpfr_init(r18177); mpfr_init(r18178); mpfr_init(r18179); mpfr_init(r18180); mpfr_init(r18181); mpfr_init(r18182); } double f_dm(double a, double b) { ; mpfr_set_d(r18163, a, MPFR_RNDN); ; mpfr_mul(r18165, r18163, r18164, MPFR_RNDN); mpfr_sub(r18166, r18162, r18165, MPFR_RNDN); mpfr_set_d(r18167, b, MPFR_RNDN); mpfr_sqr(r18168, r18167, MPFR_RNDN); mpfr_fma(r18169, r18163, r18163, r18163, MPFR_RNDN); mpfr_mul(r18170, r18169, r18163, MPFR_RNDN); mpfr_fma(r18171, r18166, r18168, r18170, MPFR_RNDN); ; ; mpfr_mul(r18174, r18167, r18173, MPFR_RNDN); mpfr_mul(r18175, r18174, r18167, MPFR_RNDN); mpfr_mul(r18176, r18163, r18163, MPFR_RNDN); mpfr_pow(r18177, r18163, r18172, MPFR_RNDN); mpfr_pow(r18178, r18167, r18172, MPFR_RNDN); mpfr_add(r18179, r18177, r18178, MPFR_RNDN); mpfr_fma(r18180, r18175, r18176, r18179, MPFR_RNDN); mpfr_fma(r18181, r18171, r18172, r18180, MPFR_RNDN); mpfr_sub(r18182, r18181, r18162, MPFR_RNDN); return mpfr_get_d(r18182, MPFR_RNDN); }