#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r22180 = a; float r22181 = r22180 * r22180; float r22182 = b; float r22183 = r22182 * r22182; float r22184 = r22181 + r22183; float r22185 = 2; float r22186 = pow(r22184, r22185); float r22187 = 4; float r22188 = 1; float r22189 = r22188 + r22180; float r22190 = r22181 * r22189; float r22191 = 3; float r22192 = r22191 * r22180; float r22193 = r22188 - r22192; float r22194 = r22183 * r22193; float r22195 = r22190 + r22194; float r22196 = r22187 * r22195; float r22197 = r22186 + r22196; float r22198 = r22197 - r22188; return r22198; } double f_id(double a, double b) { double r22199 = a; double r22200 = r22199 * r22199; double r22201 = b; double r22202 = r22201 * r22201; double r22203 = r22200 + r22202; double r22204 = 2; double r22205 = pow(r22203, r22204); double r22206 = 4; double r22207 = 1; double r22208 = r22207 + r22199; double r22209 = r22200 * r22208; double r22210 = 3; double r22211 = r22210 * r22199; double r22212 = r22207 - r22211; double r22213 = r22202 * r22212; double r22214 = r22209 + r22213; double r22215 = r22206 * r22214; double r22216 = r22205 + r22215; double r22217 = r22216 - r22207; return r22217; } double f_of(float a, float b) { float r22218 = a; float r22219 = r22218 * r22218; float r22220 = b; float r22221 = r22220 * r22220; float r22222 = r22219 + r22221; float r22223 = 2; float r22224 = pow(r22222, r22223); float r22225 = 4; float r22226 = 1; float r22227 = r22226 + r22218; float r22228 = r22219 * r22227; float r22229 = 3; float r22230 = r22229 * r22218; float r22231 = r22226 - r22230; float r22232 = r22221 * r22231; float r22233 = r22228 + r22232; float r22234 = r22225 * r22233; float r22235 = r22224 + r22234; float r22236 = sqrt(r22235); float r22237 = r22236 * r22236; float r22238 = r22237 - r22226; return r22238; } double f_od(double a, double b) { double r22239 = a; double r22240 = r22239 * r22239; double r22241 = b; double r22242 = r22241 * r22241; double r22243 = r22240 + r22242; double r22244 = 2; double r22245 = pow(r22243, r22244); double r22246 = 4; double r22247 = 1; double r22248 = r22247 + r22239; double r22249 = r22240 * r22248; double r22250 = 3; double r22251 = r22250 * r22239; double r22252 = r22247 - r22251; double r22253 = r22242 * r22252; double r22254 = r22249 + r22253; double r22255 = r22246 * r22254; double r22256 = r22245 + r22255; double r22257 = sqrt(r22256); double r22258 = r22257 * r22257; double r22259 = r22258 - r22247; return r22259; } 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 r22260, r22261, r22262, r22263, r22264, r22265, r22266, r22267, r22268, r22269, r22270, r22271, r22272, r22273, r22274, r22275, r22276, r22277, r22278; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22260); mpfr_init(r22261); mpfr_init(r22262); mpfr_init(r22263); mpfr_init(r22264); mpfr_init_set_str(r22265, "2", 10, MPFR_RNDN); mpfr_init(r22266); mpfr_init_set_str(r22267, "4", 10, MPFR_RNDN); mpfr_init_set_str(r22268, "1", 10, MPFR_RNDN); mpfr_init(r22269); mpfr_init(r22270); mpfr_init_set_str(r22271, "3", 10, MPFR_RNDN); mpfr_init(r22272); mpfr_init(r22273); mpfr_init(r22274); mpfr_init(r22275); mpfr_init(r22276); mpfr_init(r22277); mpfr_init(r22278); } double f_im(double a, double b) { mpfr_set_d(r22260, a, MPFR_RNDN); mpfr_mul(r22261, r22260, r22260, MPFR_RNDN); mpfr_set_d(r22262, b, MPFR_RNDN); mpfr_mul(r22263, r22262, r22262, MPFR_RNDN); mpfr_add(r22264, r22261, r22263, MPFR_RNDN); ; mpfr_pow(r22266, r22264, r22265, MPFR_RNDN); ; ; mpfr_add(r22269, r22268, r22260, MPFR_RNDN); mpfr_mul(r22270, r22261, r22269, MPFR_RNDN); ; mpfr_mul(r22272, r22271, r22260, MPFR_RNDN); mpfr_sub(r22273, r22268, r22272, MPFR_RNDN); mpfr_mul(r22274, r22263, r22273, MPFR_RNDN); mpfr_add(r22275, r22270, r22274, MPFR_RNDN); mpfr_mul(r22276, r22267, r22275, MPFR_RNDN); mpfr_add(r22277, r22266, r22276, MPFR_RNDN); mpfr_sub(r22278, r22277, r22268, MPFR_RNDN); return mpfr_get_d(r22278, MPFR_RNDN); } static mpfr_t r22279, r22280, r22281, r22282, r22283, r22284, r22285, r22286, r22287, r22288, r22289, r22290, r22291, r22292, r22293, r22294, r22295, r22296, r22297, r22298, r22299; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22279); mpfr_init(r22280); mpfr_init(r22281); mpfr_init(r22282); mpfr_init(r22283); mpfr_init_set_str(r22284, "2", 10, MPFR_RNDN); mpfr_init(r22285); mpfr_init_set_str(r22286, "4", 10, MPFR_RNDN); mpfr_init_set_str(r22287, "1", 10, MPFR_RNDN); mpfr_init(r22288); mpfr_init(r22289); mpfr_init_set_str(r22290, "3", 10, MPFR_RNDN); mpfr_init(r22291); mpfr_init(r22292); mpfr_init(r22293); mpfr_init(r22294); mpfr_init(r22295); mpfr_init(r22296); mpfr_init(r22297); mpfr_init(r22298); mpfr_init(r22299); } double f_fm(double a, double b) { mpfr_set_d(r22279, a, MPFR_RNDN); mpfr_mul(r22280, r22279, r22279, MPFR_RNDN); mpfr_set_d(r22281, b, MPFR_RNDN); mpfr_mul(r22282, r22281, r22281, MPFR_RNDN); mpfr_add(r22283, r22280, r22282, MPFR_RNDN); ; mpfr_pow(r22285, r22283, r22284, MPFR_RNDN); ; ; mpfr_add(r22288, r22287, r22279, MPFR_RNDN); mpfr_mul(r22289, r22280, r22288, MPFR_RNDN); ; mpfr_mul(r22291, r22290, r22279, MPFR_RNDN); mpfr_sub(r22292, r22287, r22291, MPFR_RNDN); mpfr_mul(r22293, r22282, r22292, MPFR_RNDN); mpfr_add(r22294, r22289, r22293, MPFR_RNDN); mpfr_mul(r22295, r22286, r22294, MPFR_RNDN); mpfr_add(r22296, r22285, r22295, MPFR_RNDN); mpfr_sqrt(r22297, r22296, MPFR_RNDN); mpfr_mul(r22298, r22297, r22297, MPFR_RNDN); mpfr_sub(r22299, r22298, r22287, MPFR_RNDN); return mpfr_get_d(r22299, MPFR_RNDN); } static mpfr_t r22300, r22301, r22302, r22303, r22304, r22305, r22306, r22307, r22308, r22309, r22310, r22311, r22312, r22313, r22314, r22315, r22316, r22317, r22318, r22319, r22320; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22300); mpfr_init(r22301); mpfr_init(r22302); mpfr_init(r22303); mpfr_init(r22304); mpfr_init_set_str(r22305, "2", 10, MPFR_RNDN); mpfr_init(r22306); mpfr_init_set_str(r22307, "4", 10, MPFR_RNDN); mpfr_init_set_str(r22308, "1", 10, MPFR_RNDN); mpfr_init(r22309); mpfr_init(r22310); mpfr_init_set_str(r22311, "3", 10, MPFR_RNDN); mpfr_init(r22312); mpfr_init(r22313); mpfr_init(r22314); mpfr_init(r22315); mpfr_init(r22316); mpfr_init(r22317); mpfr_init(r22318); mpfr_init(r22319); mpfr_init(r22320); } double f_dm(double a, double b) { mpfr_set_d(r22300, a, MPFR_RNDN); mpfr_mul(r22301, r22300, r22300, MPFR_RNDN); mpfr_set_d(r22302, b, MPFR_RNDN); mpfr_mul(r22303, r22302, r22302, MPFR_RNDN); mpfr_add(r22304, r22301, r22303, MPFR_RNDN); ; mpfr_pow(r22306, r22304, r22305, MPFR_RNDN); ; ; mpfr_add(r22309, r22308, r22300, MPFR_RNDN); mpfr_mul(r22310, r22301, r22309, MPFR_RNDN); ; mpfr_mul(r22312, r22311, r22300, MPFR_RNDN); mpfr_sub(r22313, r22308, r22312, MPFR_RNDN); mpfr_mul(r22314, r22303, r22313, MPFR_RNDN); mpfr_add(r22315, r22310, r22314, MPFR_RNDN); mpfr_mul(r22316, r22307, r22315, MPFR_RNDN); mpfr_add(r22317, r22306, r22316, MPFR_RNDN); mpfr_sqrt(r22318, r22317, MPFR_RNDN); mpfr_mul(r22319, r22318, r22318, MPFR_RNDN); mpfr_sub(r22320, r22319, r22308, MPFR_RNDN); return mpfr_get_d(r22320, MPFR_RNDN); }