#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (20)"; double f_if(float f) { float r24304 = 1; float r24305 = atan2(1.0, 0.0); float r24306 = 4; float r24307 = r24305 / r24306; float r24308 = r24304 / r24307; float r24309 = f; float r24310 = r24307 * r24309; float r24311 = exp(r24310); float r24312 = -r24310; float r24313 = exp(r24312); float r24314 = r24311 + r24313; float r24315 = r24311 - r24313; float r24316 = r24314 / r24315; float r24317 = log(r24316); float r24318 = r24308 * r24317; float r24319 = -r24318; return r24319; } double f_id(double f) { double r24320 = 1; double r24321 = atan2(1.0, 0.0); double r24322 = 4; double r24323 = r24321 / r24322; double r24324 = r24320 / r24323; double r24325 = f; double r24326 = r24323 * r24325; double r24327 = exp(r24326); double r24328 = -r24326; double r24329 = exp(r24328); double r24330 = r24327 + r24329; double r24331 = r24327 - r24329; double r24332 = r24330 / r24331; double r24333 = log(r24332); double r24334 = r24324 * r24333; double r24335 = -r24334; return r24335; } double f_of(float f) { float r24336 = 4; float r24337 = atan2(1.0, 0.0); float r24338 = r24336 / r24337; float r24339 = -r24338; float r24340 = cbrt(r24338); float r24341 = r24340 * r24340; float r24342 = log(r24341); float r24343 = r24339 * r24342; float r24344 = log(r24340); float r24345 = f; float r24346 = -r24345; float r24347 = r24337 * r24345; float r24348 = 7/96; float r24349 = 3/2; float r24350 = fma(r24347, r24348, r24349); float r24351 = log(r24345); float r24352 = r24337 / r24336; float r24353 = r24351 / r24352; float r24354 = fma(r24346, r24350, r24353); float r24355 = fma(r24339, r24344, r24354); float r24356 = r24343 + r24355; return r24356; } double f_od(double f) { double r24357 = 4; double r24358 = atan2(1.0, 0.0); double r24359 = r24357 / r24358; double r24360 = -r24359; double r24361 = cbrt(r24359); double r24362 = r24361 * r24361; double r24363 = log(r24362); double r24364 = r24360 * r24363; double r24365 = log(r24361); double r24366 = f; double r24367 = -r24366; double r24368 = r24358 * r24366; double r24369 = 7/96; double r24370 = 3/2; double r24371 = fma(r24368, r24369, r24370); double r24372 = log(r24366); double r24373 = r24358 / r24357; double r24374 = r24372 / r24373; double r24375 = fma(r24367, r24371, r24374); double r24376 = fma(r24360, r24365, r24375); double r24377 = r24364 + r24376; return r24377; } 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 r24378, r24379, r24380, r24381, r24382, r24383, r24384, r24385, r24386, r24387, r24388, r24389, r24390, r24391, r24392, r24393; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init_set_str(r24378, "1", 10, MPFR_RNDN); mpfr_init(r24379); mpfr_init_set_str(r24380, "4", 10, MPFR_RNDN); mpfr_init(r24381); mpfr_init(r24382); mpfr_init(r24383); mpfr_init(r24384); mpfr_init(r24385); mpfr_init(r24386); mpfr_init(r24387); mpfr_init(r24388); mpfr_init(r24389); mpfr_init(r24390); mpfr_init(r24391); mpfr_init(r24392); mpfr_init(r24393); } double f_im(double f) { ; mpfr_const_pi(r24379, MPFR_RNDN); ; mpfr_div(r24381, r24379, r24380, MPFR_RNDN); mpfr_div(r24382, r24378, r24381, MPFR_RNDN); mpfr_set_d(r24383, f, MPFR_RNDN); mpfr_mul(r24384, r24381, r24383, MPFR_RNDN); mpfr_exp(r24385, r24384, MPFR_RNDN); mpfr_neg(r24386, r24384, MPFR_RNDN); mpfr_exp(r24387, r24386, MPFR_RNDN); mpfr_add(r24388, r24385, r24387, MPFR_RNDN); mpfr_sub(r24389, r24385, r24387, MPFR_RNDN); mpfr_div(r24390, r24388, r24389, MPFR_RNDN); mpfr_log(r24391, r24390, MPFR_RNDN); mpfr_mul(r24392, r24382, r24391, MPFR_RNDN); mpfr_neg(r24393, r24392, MPFR_RNDN); return mpfr_get_d(r24393, MPFR_RNDN); } static mpfr_t r24394, r24395, r24396, r24397, r24398, r24399, r24400, r24401, r24402, r24403, r24404, r24405, r24406, r24407, r24408, r24409, r24410, r24411, r24412, r24413, r24414; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r24394, "4", 10, MPFR_RNDN); mpfr_init(r24395); mpfr_init(r24396); mpfr_init(r24397); mpfr_init(r24398); mpfr_init(r24399); mpfr_init(r24400); mpfr_init(r24401); mpfr_init(r24402); mpfr_init(r24403); mpfr_init(r24404); mpfr_init(r24405); mpfr_init_set_str(r24406, "7/96", 10, MPFR_RNDN); mpfr_init_set_str(r24407, "3/2", 10, MPFR_RNDN); mpfr_init(r24408); mpfr_init(r24409); mpfr_init(r24410); mpfr_init(r24411); mpfr_init(r24412); mpfr_init(r24413); mpfr_init(r24414); } double f_fm(double f) { ; mpfr_const_pi(r24395, MPFR_RNDN); mpfr_div(r24396, r24394, r24395, MPFR_RNDN); mpfr_neg(r24397, r24396, MPFR_RNDN); mpfr_cbrt(r24398, r24396, MPFR_RNDN); mpfr_mul(r24399, r24398, r24398, MPFR_RNDN); mpfr_log(r24400, r24399, MPFR_RNDN); mpfr_mul(r24401, r24397, r24400, MPFR_RNDN); mpfr_log(r24402, r24398, MPFR_RNDN); mpfr_set_d(r24403, f, MPFR_RNDN); mpfr_neg(r24404, r24403, MPFR_RNDN); mpfr_mul(r24405, r24395, r24403, MPFR_RNDN); ; ; mpfr_fma(r24408, r24405, r24406, r24407, MPFR_RNDN); mpfr_log(r24409, r24403, MPFR_RNDN); mpfr_div(r24410, r24395, r24394, MPFR_RNDN); mpfr_div(r24411, r24409, r24410, MPFR_RNDN); mpfr_fma(r24412, r24404, r24408, r24411, MPFR_RNDN); mpfr_fma(r24413, r24397, r24402, r24412, MPFR_RNDN); mpfr_add(r24414, r24401, r24413, MPFR_RNDN); return mpfr_get_d(r24414, MPFR_RNDN); } static mpfr_t r24415, r24416, r24417, r24418, r24419, r24420, r24421, r24422, r24423, r24424, r24425, r24426, r24427, r24428, r24429, r24430, r24431, r24432, r24433, r24434, r24435; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r24415, "4", 10, MPFR_RNDN); mpfr_init(r24416); mpfr_init(r24417); mpfr_init(r24418); mpfr_init(r24419); mpfr_init(r24420); mpfr_init(r24421); mpfr_init(r24422); mpfr_init(r24423); mpfr_init(r24424); mpfr_init(r24425); mpfr_init(r24426); mpfr_init_set_str(r24427, "7/96", 10, MPFR_RNDN); mpfr_init_set_str(r24428, "3/2", 10, MPFR_RNDN); mpfr_init(r24429); mpfr_init(r24430); mpfr_init(r24431); mpfr_init(r24432); mpfr_init(r24433); mpfr_init(r24434); mpfr_init(r24435); } double f_dm(double f) { ; mpfr_const_pi(r24416, MPFR_RNDN); mpfr_div(r24417, r24415, r24416, MPFR_RNDN); mpfr_neg(r24418, r24417, MPFR_RNDN); mpfr_cbrt(r24419, r24417, MPFR_RNDN); mpfr_mul(r24420, r24419, r24419, MPFR_RNDN); mpfr_log(r24421, r24420, MPFR_RNDN); mpfr_mul(r24422, r24418, r24421, MPFR_RNDN); mpfr_log(r24423, r24419, MPFR_RNDN); mpfr_set_d(r24424, f, MPFR_RNDN); mpfr_neg(r24425, r24424, MPFR_RNDN); mpfr_mul(r24426, r24416, r24424, MPFR_RNDN); ; ; mpfr_fma(r24429, r24426, r24427, r24428, MPFR_RNDN); mpfr_log(r24430, r24424, MPFR_RNDN); mpfr_div(r24431, r24416, r24415, MPFR_RNDN); mpfr_div(r24432, r24430, r24431, MPFR_RNDN); mpfr_fma(r24433, r24425, r24429, r24432, MPFR_RNDN); mpfr_fma(r24434, r24418, r24423, r24433, MPFR_RNDN); mpfr_add(r24435, r24422, r24434, MPFR_RNDN); return mpfr_get_d(r24435, MPFR_RNDN); }