#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r18296 = a; float r18297 = r18296 * r18296; float r18298 = b; float r18299 = r18298 * r18298; float r18300 = r18297 + r18299; float r18301 = r18300 * r18300; float r18302 = 4.0f; float r18303 = 1.0f; float r18304 = r18303 + r18296; float r18305 = r18297 * r18304; float r18306 = 3.0f; float r18307 = r18306 * r18296; float r18308 = r18303 - r18307; float r18309 = r18299 * r18308; float r18310 = r18305 + r18309; float r18311 = r18302 * r18310; float r18312 = r18301 + r18311; float r18313 = r18312 - r18303; return r18313; } double f_id(double a, double b) { double r18314 = a; double r18315 = r18314 * r18314; double r18316 = b; double r18317 = r18316 * r18316; double r18318 = r18315 + r18317; double r18319 = r18318 * r18318; double r18320 = 4.0; double r18321 = 1.0; double r18322 = r18321 + r18314; double r18323 = r18315 * r18322; double r18324 = 3.0; double r18325 = r18324 * r18314; double r18326 = r18321 - r18325; double r18327 = r18317 * r18326; double r18328 = r18323 + r18327; double r18329 = r18320 * r18328; double r18330 = r18319 + r18329; double r18331 = r18330 - r18321; return r18331; } double f_of(float a, float b) { float r18332 = 1.0f; float r18333 = a; float r18334 = 3.0f; float r18335 = r18333 * r18334; float r18336 = r18332 - r18335; float r18337 = b; float r18338 = r18337 * r18337; float r18339 = fma(r18333, r18333, r18333); float r18340 = r18339 * r18333; float r18341 = fma(r18336, r18338, r18340); float r18342 = 4.0f; float r18343 = 2.0f; float r18344 = r18337 * r18343; float r18345 = r18344 * r18337; float r18346 = r18333 * r18333; float r18347 = pow(r18333, r18342); float r18348 = pow(r18337, r18342); float r18349 = r18347 + r18348; float r18350 = fma(r18345, r18346, r18349); float r18351 = fma(r18341, r18342, r18350); float r18352 = r18351 - r18332; return r18352; } double f_od(double a, double b) { double r18353 = 1.0; double r18354 = a; double r18355 = 3.0; double r18356 = r18354 * r18355; double r18357 = r18353 - r18356; double r18358 = b; double r18359 = r18358 * r18358; double r18360 = fma(r18354, r18354, r18354); double r18361 = r18360 * r18354; double r18362 = fma(r18357, r18359, r18361); double r18363 = 4.0; double r18364 = 2.0; double r18365 = r18358 * r18364; double r18366 = r18365 * r18358; double r18367 = r18354 * r18354; double r18368 = pow(r18354, r18363); double r18369 = pow(r18358, r18363); double r18370 = r18368 + r18369; double r18371 = fma(r18366, r18367, r18370); double r18372 = fma(r18362, r18363, r18371); double r18373 = r18372 - r18353; return r18373; } 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 r18374, r18375, r18376, r18377, r18378, r18379, r18380, r18381, r18382, r18383, r18384, r18385, r18386, r18387, r18388, r18389, r18390, r18391; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r18374); mpfr_init(r18375); mpfr_init(r18376); mpfr_init(r18377); mpfr_init(r18378); mpfr_init(r18379); mpfr_init_set_str(r18380, "4", 10, MPFR_RNDN); mpfr_init_set_str(r18381, "1", 10, MPFR_RNDN); mpfr_init(r18382); mpfr_init(r18383); mpfr_init_set_str(r18384, "3", 10, MPFR_RNDN); mpfr_init(r18385); mpfr_init(r18386); mpfr_init(r18387); mpfr_init(r18388); mpfr_init(r18389); mpfr_init(r18390); mpfr_init(r18391); } double f_im(double a, double b) { mpfr_set_d(r18374, a, MPFR_RNDN); mpfr_sqr(r18375, r18374, MPFR_RNDN); mpfr_set_d(r18376, b, MPFR_RNDN); mpfr_sqr(r18377, r18376, MPFR_RNDN); mpfr_add(r18378, r18375, r18377, MPFR_RNDN); mpfr_sqr(r18379, r18378, MPFR_RNDN); ; ; mpfr_add(r18382, r18381, r18374, MPFR_RNDN); mpfr_mul(r18383, r18375, r18382, MPFR_RNDN); ; mpfr_mul(r18385, r18384, r18374, MPFR_RNDN); mpfr_sub(r18386, r18381, r18385, MPFR_RNDN); mpfr_mul(r18387, r18377, r18386, MPFR_RNDN); mpfr_add(r18388, r18383, r18387, MPFR_RNDN); mpfr_mul(r18389, r18380, r18388, MPFR_RNDN); mpfr_add(r18390, r18379, r18389, MPFR_RNDN); mpfr_sub(r18391, r18390, r18381, MPFR_RNDN); return mpfr_get_d(r18391, MPFR_RNDN); } static mpfr_t r18392, r18393, r18394, r18395, r18396, r18397, r18398, r18399, r18400, r18401, r18402, r18403, r18404, r18405, r18406, r18407, r18408, r18409, r18410, r18411, r18412; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18392, "1", 10, MPFR_RNDN); mpfr_init(r18393); mpfr_init_set_str(r18394, "3", 10, MPFR_RNDN); mpfr_init(r18395); mpfr_init(r18396); mpfr_init(r18397); mpfr_init(r18398); mpfr_init(r18399); mpfr_init(r18400); mpfr_init(r18401); mpfr_init_set_str(r18402, "4", 10, MPFR_RNDN); mpfr_init_set_str(r18403, "2", 10, MPFR_RNDN); mpfr_init(r18404); mpfr_init(r18405); mpfr_init(r18406); mpfr_init(r18407); mpfr_init(r18408); mpfr_init(r18409); mpfr_init(r18410); mpfr_init(r18411); mpfr_init(r18412); } double f_fm(double a, double b) { ; mpfr_set_d(r18393, a, MPFR_RNDN); ; mpfr_mul(r18395, r18393, r18394, MPFR_RNDN); mpfr_sub(r18396, r18392, r18395, MPFR_RNDN); mpfr_set_d(r18397, b, MPFR_RNDN); mpfr_sqr(r18398, r18397, MPFR_RNDN); mpfr_fma(r18399, r18393, r18393, r18393, MPFR_RNDN); mpfr_mul(r18400, r18399, r18393, MPFR_RNDN); mpfr_fma(r18401, r18396, r18398, r18400, MPFR_RNDN); ; ; mpfr_mul(r18404, r18397, r18403, MPFR_RNDN); mpfr_mul(r18405, r18404, r18397, MPFR_RNDN); mpfr_mul(r18406, r18393, r18393, MPFR_RNDN); mpfr_pow(r18407, r18393, r18402, MPFR_RNDN); mpfr_pow(r18408, r18397, r18402, MPFR_RNDN); mpfr_add(r18409, r18407, r18408, MPFR_RNDN); mpfr_fma(r18410, r18405, r18406, r18409, MPFR_RNDN); mpfr_fma(r18411, r18401, r18402, r18410, MPFR_RNDN); mpfr_sub(r18412, r18411, r18392, MPFR_RNDN); return mpfr_get_d(r18412, MPFR_RNDN); } static mpfr_t r18413, r18414, r18415, r18416, r18417, r18418, r18419, r18420, r18421, r18422, r18423, r18424, r18425, r18426, r18427, r18428, r18429, r18430, r18431, r18432, r18433; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18413, "1", 10, MPFR_RNDN); mpfr_init(r18414); mpfr_init_set_str(r18415, "3", 10, MPFR_RNDN); mpfr_init(r18416); mpfr_init(r18417); mpfr_init(r18418); mpfr_init(r18419); mpfr_init(r18420); mpfr_init(r18421); mpfr_init(r18422); mpfr_init_set_str(r18423, "4", 10, MPFR_RNDN); mpfr_init_set_str(r18424, "2", 10, MPFR_RNDN); mpfr_init(r18425); mpfr_init(r18426); mpfr_init(r18427); mpfr_init(r18428); mpfr_init(r18429); mpfr_init(r18430); mpfr_init(r18431); mpfr_init(r18432); mpfr_init(r18433); } double f_dm(double a, double b) { ; mpfr_set_d(r18414, a, MPFR_RNDN); ; mpfr_mul(r18416, r18414, r18415, MPFR_RNDN); mpfr_sub(r18417, r18413, r18416, MPFR_RNDN); mpfr_set_d(r18418, b, MPFR_RNDN); mpfr_sqr(r18419, r18418, MPFR_RNDN); mpfr_fma(r18420, r18414, r18414, r18414, MPFR_RNDN); mpfr_mul(r18421, r18420, r18414, MPFR_RNDN); mpfr_fma(r18422, r18417, r18419, r18421, MPFR_RNDN); ; ; mpfr_mul(r18425, r18418, r18424, MPFR_RNDN); mpfr_mul(r18426, r18425, r18418, MPFR_RNDN); mpfr_mul(r18427, r18414, r18414, MPFR_RNDN); mpfr_pow(r18428, r18414, r18423, MPFR_RNDN); mpfr_pow(r18429, r18418, r18423, MPFR_RNDN); mpfr_add(r18430, r18428, r18429, MPFR_RNDN); mpfr_fma(r18431, r18426, r18427, r18430, MPFR_RNDN); mpfr_fma(r18432, r18422, r18423, r18431, MPFR_RNDN); mpfr_sub(r18433, r18432, r18413, MPFR_RNDN); return mpfr_get_d(r18433, MPFR_RNDN); }