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