#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r18210 = a; float r18211 = r18210 * r18210; float r18212 = b; float r18213 = r18212 * r18212; float r18214 = r18211 + r18213; float r18215 = r18214 * r18214; float r18216 = 4.0f; float r18217 = 1.0f; float r18218 = r18217 + r18210; float r18219 = r18211 * r18218; float r18220 = 3.0f; float r18221 = r18220 * r18210; float r18222 = r18217 - r18221; float r18223 = r18213 * r18222; float r18224 = r18219 + r18223; float r18225 = r18216 * r18224; float r18226 = r18215 + r18225; float r18227 = r18226 - r18217; return r18227; } double f_id(double a, double b) { double r18228 = a; double r18229 = r18228 * r18228; double r18230 = b; double r18231 = r18230 * r18230; double r18232 = r18229 + r18231; double r18233 = r18232 * r18232; double r18234 = 4.0; double r18235 = 1.0; double r18236 = r18235 + r18228; double r18237 = r18229 * r18236; double r18238 = 3.0; double r18239 = r18238 * r18228; double r18240 = r18235 - r18239; double r18241 = r18231 * r18240; double r18242 = r18237 + r18241; double r18243 = r18234 * r18242; double r18244 = r18233 + r18243; double r18245 = r18244 - r18235; return r18245; } double f_of(float a, float b) { float r18246 = 1.0f; float r18247 = a; float r18248 = 3.0f; float r18249 = r18247 * r18248; float r18250 = r18246 - r18249; float r18251 = b; float r18252 = r18251 * r18251; float r18253 = fma(r18247, r18247, r18247); float r18254 = r18253 * r18247; float r18255 = fma(r18250, r18252, r18254); float r18256 = 4.0f; float r18257 = 2.0f; float r18258 = r18251 * r18257; float r18259 = r18258 * r18251; float r18260 = r18247 * r18247; float r18261 = pow(r18247, r18256); float r18262 = pow(r18251, r18256); float r18263 = r18261 + r18262; float r18264 = fma(r18259, r18260, r18263); float r18265 = fma(r18255, r18256, r18264); float r18266 = r18265 - r18246; return r18266; } double f_od(double a, double b) { double r18267 = 1.0; double r18268 = a; double r18269 = 3.0; double r18270 = r18268 * r18269; double r18271 = r18267 - r18270; double r18272 = b; double r18273 = r18272 * r18272; double r18274 = fma(r18268, r18268, r18268); double r18275 = r18274 * r18268; double r18276 = fma(r18271, r18273, r18275); double r18277 = 4.0; double r18278 = 2.0; double r18279 = r18272 * r18278; double r18280 = r18279 * r18272; double r18281 = r18268 * r18268; double r18282 = pow(r18268, r18277); double r18283 = pow(r18272, r18277); double r18284 = r18282 + r18283; double r18285 = fma(r18280, r18281, r18284); double r18286 = fma(r18276, r18277, r18285); double r18287 = r18286 - r18267; return r18287; } 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 r18288, r18289, r18290, r18291, r18292, r18293, r18294, r18295, r18296, r18297, r18298, r18299, r18300, r18301, r18302, r18303, r18304, r18305; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r18288); mpfr_init(r18289); mpfr_init(r18290); mpfr_init(r18291); mpfr_init(r18292); mpfr_init(r18293); mpfr_init_set_str(r18294, "4", 10, MPFR_RNDN); mpfr_init_set_str(r18295, "1", 10, MPFR_RNDN); mpfr_init(r18296); mpfr_init(r18297); mpfr_init_set_str(r18298, "3", 10, MPFR_RNDN); mpfr_init(r18299); mpfr_init(r18300); mpfr_init(r18301); mpfr_init(r18302); mpfr_init(r18303); mpfr_init(r18304); mpfr_init(r18305); } double f_im(double a, double b) { mpfr_set_d(r18288, a, MPFR_RNDN); mpfr_sqr(r18289, r18288, MPFR_RNDN); mpfr_set_d(r18290, b, MPFR_RNDN); mpfr_sqr(r18291, r18290, MPFR_RNDN); mpfr_add(r18292, r18289, r18291, MPFR_RNDN); mpfr_sqr(r18293, r18292, MPFR_RNDN); ; ; mpfr_add(r18296, r18295, r18288, MPFR_RNDN); mpfr_mul(r18297, r18289, r18296, MPFR_RNDN); ; mpfr_mul(r18299, r18298, r18288, MPFR_RNDN); mpfr_sub(r18300, r18295, r18299, MPFR_RNDN); mpfr_mul(r18301, r18291, r18300, MPFR_RNDN); mpfr_add(r18302, r18297, r18301, MPFR_RNDN); mpfr_mul(r18303, r18294, r18302, MPFR_RNDN); mpfr_add(r18304, r18293, r18303, MPFR_RNDN); mpfr_sub(r18305, r18304, r18295, MPFR_RNDN); return mpfr_get_d(r18305, MPFR_RNDN); } static mpfr_t r18306, r18307, r18308, r18309, r18310, r18311, r18312, r18313, r18314, r18315, r18316, r18317, r18318, r18319, r18320, r18321, r18322, r18323, r18324, r18325, r18326; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18306, "1", 10, MPFR_RNDN); mpfr_init(r18307); mpfr_init_set_str(r18308, "3", 10, MPFR_RNDN); mpfr_init(r18309); mpfr_init(r18310); mpfr_init(r18311); mpfr_init(r18312); mpfr_init(r18313); mpfr_init(r18314); mpfr_init(r18315); mpfr_init_set_str(r18316, "4", 10, MPFR_RNDN); mpfr_init_set_str(r18317, "2", 10, MPFR_RNDN); mpfr_init(r18318); mpfr_init(r18319); mpfr_init(r18320); mpfr_init(r18321); mpfr_init(r18322); mpfr_init(r18323); mpfr_init(r18324); mpfr_init(r18325); mpfr_init(r18326); } double f_fm(double a, double b) { ; mpfr_set_d(r18307, a, MPFR_RNDN); ; mpfr_mul(r18309, r18307, r18308, MPFR_RNDN); mpfr_sub(r18310, r18306, r18309, MPFR_RNDN); mpfr_set_d(r18311, b, MPFR_RNDN); mpfr_sqr(r18312, r18311, MPFR_RNDN); mpfr_fma(r18313, r18307, r18307, r18307, MPFR_RNDN); mpfr_mul(r18314, r18313, r18307, MPFR_RNDN); mpfr_fma(r18315, r18310, r18312, r18314, MPFR_RNDN); ; ; mpfr_mul(r18318, r18311, r18317, MPFR_RNDN); mpfr_mul(r18319, r18318, r18311, MPFR_RNDN); mpfr_mul(r18320, r18307, r18307, MPFR_RNDN); mpfr_pow(r18321, r18307, r18316, MPFR_RNDN); mpfr_pow(r18322, r18311, r18316, MPFR_RNDN); mpfr_add(r18323, r18321, r18322, MPFR_RNDN); mpfr_fma(r18324, r18319, r18320, r18323, MPFR_RNDN); mpfr_fma(r18325, r18315, r18316, r18324, MPFR_RNDN); mpfr_sub(r18326, r18325, r18306, MPFR_RNDN); return mpfr_get_d(r18326, MPFR_RNDN); } static mpfr_t r18327, r18328, r18329, r18330, r18331, r18332, r18333, r18334, r18335, r18336, r18337, r18338, r18339, r18340, r18341, r18342, r18343, r18344, r18345, r18346, r18347; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18327, "1", 10, MPFR_RNDN); mpfr_init(r18328); mpfr_init_set_str(r18329, "3", 10, MPFR_RNDN); mpfr_init(r18330); mpfr_init(r18331); mpfr_init(r18332); mpfr_init(r18333); mpfr_init(r18334); mpfr_init(r18335); mpfr_init(r18336); mpfr_init_set_str(r18337, "4", 10, MPFR_RNDN); mpfr_init_set_str(r18338, "2", 10, MPFR_RNDN); mpfr_init(r18339); mpfr_init(r18340); mpfr_init(r18341); mpfr_init(r18342); mpfr_init(r18343); mpfr_init(r18344); mpfr_init(r18345); mpfr_init(r18346); mpfr_init(r18347); } double f_dm(double a, double b) { ; mpfr_set_d(r18328, a, MPFR_RNDN); ; mpfr_mul(r18330, r18328, r18329, MPFR_RNDN); mpfr_sub(r18331, r18327, r18330, MPFR_RNDN); mpfr_set_d(r18332, b, MPFR_RNDN); mpfr_sqr(r18333, r18332, MPFR_RNDN); mpfr_fma(r18334, r18328, r18328, r18328, MPFR_RNDN); mpfr_mul(r18335, r18334, r18328, MPFR_RNDN); mpfr_fma(r18336, r18331, r18333, r18335, MPFR_RNDN); ; ; mpfr_mul(r18339, r18332, r18338, MPFR_RNDN); mpfr_mul(r18340, r18339, r18332, MPFR_RNDN); mpfr_mul(r18341, r18328, r18328, MPFR_RNDN); mpfr_pow(r18342, r18328, r18337, MPFR_RNDN); mpfr_pow(r18343, r18332, r18337, MPFR_RNDN); mpfr_add(r18344, r18342, r18343, MPFR_RNDN); mpfr_fma(r18345, r18340, r18341, r18344, MPFR_RNDN); mpfr_fma(r18346, r18336, r18337, r18345, MPFR_RNDN); mpfr_sub(r18347, r18346, r18327, MPFR_RNDN); return mpfr_get_d(r18347, MPFR_RNDN); }