#include #include #include #include #include char *name = "Jmat.Real.erfi, branch x less than or equal to 0.5"; double f_if(float x) { float r24251 = 1; float r24252 = atan2(1.0, 0.0); float r24253 = sqrt(r24252); float r24254 = r24251 / r24253; float r24255 = 2; float r24256 = x; float r24257 = fabs(r24256); float r24258 = r24255 * r24257; float r24259 = 3; float r24260 = r24255 / r24259; float r24261 = r24257 * r24257; float r24262 = r24261 * r24257; float r24263 = r24260 * r24262; float r24264 = r24258 + r24263; float r24265 = 5; float r24266 = r24251 / r24265; float r24267 = r24262 * r24257; float r24268 = r24267 * r24257; float r24269 = r24266 * r24268; float r24270 = r24264 + r24269; float r24271 = 21; float r24272 = r24251 / r24271; float r24273 = r24268 * r24257; float r24274 = r24273 * r24257; float r24275 = r24272 * r24274; float r24276 = r24270 + r24275; float r24277 = r24254 * r24276; float r24278 = fabs(r24277); return r24278; } double f_id(double x) { double r24279 = 1; double r24280 = atan2(1.0, 0.0); double r24281 = sqrt(r24280); double r24282 = r24279 / r24281; double r24283 = 2; double r24284 = x; double r24285 = fabs(r24284); double r24286 = r24283 * r24285; double r24287 = 3; double r24288 = r24283 / r24287; double r24289 = r24285 * r24285; double r24290 = r24289 * r24285; double r24291 = r24288 * r24290; double r24292 = r24286 + r24291; double r24293 = 5; double r24294 = r24279 / r24293; double r24295 = r24290 * r24285; double r24296 = r24295 * r24285; double r24297 = r24294 * r24296; double r24298 = r24292 + r24297; double r24299 = 21; double r24300 = r24279 / r24299; double r24301 = r24296 * r24285; double r24302 = r24301 * r24285; double r24303 = r24300 * r24302; double r24304 = r24298 + r24303; double r24305 = r24282 * r24304; double r24306 = fabs(r24305); return r24306; } double f_of(float x) { float r24307 = 1; float r24308 = atan2(1.0, 0.0); float r24309 = sqrt(r24308); float r24310 = r24307 / r24309; float r24311 = 2; float r24312 = x; float r24313 = fabs(r24312); float r24314 = r24311 * r24313; float r24315 = 3; float r24316 = r24311 / r24315; float r24317 = r24313 * r24313; float r24318 = r24317 * r24313; float r24319 = r24316 * r24318; float r24320 = r24314 + r24319; float r24321 = 5; float r24322 = r24307 / r24321; float r24323 = r24318 * r24313; float r24324 = r24323 * r24313; float r24325 = r24322 * r24324; float r24326 = r24320 + r24325; float r24327 = r24310 * r24326; float r24328 = 21; float r24329 = r24313 / r24328; float r24330 = r24329 / r24309; float r24331 = pow(r24313, r24315); float r24332 = r24331 * r24331; float r24333 = r24330 * r24332; float r24334 = r24327 + r24333; float r24335 = fabs(r24334); return r24335; } double f_od(double x) { double r24336 = 1; double r24337 = atan2(1.0, 0.0); double r24338 = sqrt(r24337); double r24339 = r24336 / r24338; double r24340 = 2; double r24341 = x; double r24342 = fabs(r24341); double r24343 = r24340 * r24342; double r24344 = 3; double r24345 = r24340 / r24344; double r24346 = r24342 * r24342; double r24347 = r24346 * r24342; double r24348 = r24345 * r24347; double r24349 = r24343 + r24348; double r24350 = 5; double r24351 = r24336 / r24350; double r24352 = r24347 * r24342; double r24353 = r24352 * r24342; double r24354 = r24351 * r24353; double r24355 = r24349 + r24354; double r24356 = r24339 * r24355; double r24357 = 21; double r24358 = r24342 / r24357; double r24359 = r24358 / r24338; double r24360 = pow(r24342, r24344); double r24361 = r24360 * r24360; double r24362 = r24359 * r24361; double r24363 = r24356 + r24362; double r24364 = fabs(r24363); return r24364; } 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 r24365, r24366, r24367, r24368, r24369, r24370, r24371, r24372, r24373, r24374, r24375, r24376, r24377, r24378, r24379, r24380, r24381, r24382, r24383, r24384, r24385, r24386, r24387, r24388, r24389, r24390, r24391, r24392; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r24365, "1", 10, MPFR_RNDN); mpfr_init(r24366); mpfr_init(r24367); mpfr_init(r24368); mpfr_init_set_str(r24369, "2", 10, MPFR_RNDN); mpfr_init(r24370); mpfr_init(r24371); mpfr_init(r24372); mpfr_init_set_str(r24373, "3", 10, MPFR_RNDN); mpfr_init(r24374); mpfr_init(r24375); mpfr_init(r24376); mpfr_init(r24377); mpfr_init(r24378); mpfr_init_set_str(r24379, "5", 10, MPFR_RNDN); mpfr_init(r24380); mpfr_init(r24381); mpfr_init(r24382); mpfr_init(r24383); mpfr_init(r24384); mpfr_init_set_str(r24385, "21", 10, MPFR_RNDN); mpfr_init(r24386); mpfr_init(r24387); mpfr_init(r24388); mpfr_init(r24389); mpfr_init(r24390); mpfr_init(r24391); mpfr_init(r24392); } double f_im(double x) { ; mpfr_const_pi(r24366, MPFR_RNDN); mpfr_sqrt(r24367, r24366, MPFR_RNDN); mpfr_div(r24368, r24365, r24367, MPFR_RNDN); ; mpfr_set_d(r24370, x, MPFR_RNDN); mpfr_abs(r24371, r24370, MPFR_RNDN); mpfr_mul(r24372, r24369, r24371, MPFR_RNDN); ; mpfr_div(r24374, r24369, r24373, MPFR_RNDN); mpfr_mul(r24375, r24371, r24371, MPFR_RNDN); mpfr_mul(r24376, r24375, r24371, MPFR_RNDN); mpfr_mul(r24377, r24374, r24376, MPFR_RNDN); mpfr_add(r24378, r24372, r24377, MPFR_RNDN); ; mpfr_div(r24380, r24365, r24379, MPFR_RNDN); mpfr_mul(r24381, r24376, r24371, MPFR_RNDN); mpfr_mul(r24382, r24381, r24371, MPFR_RNDN); mpfr_mul(r24383, r24380, r24382, MPFR_RNDN); mpfr_add(r24384, r24378, r24383, MPFR_RNDN); ; mpfr_div(r24386, r24365, r24385, MPFR_RNDN); mpfr_mul(r24387, r24382, r24371, MPFR_RNDN); mpfr_mul(r24388, r24387, r24371, MPFR_RNDN); mpfr_mul(r24389, r24386, r24388, MPFR_RNDN); mpfr_add(r24390, r24384, r24389, MPFR_RNDN); mpfr_mul(r24391, r24368, r24390, MPFR_RNDN); mpfr_abs(r24392, r24391, MPFR_RNDN); return mpfr_get_d(r24392, MPFR_RNDN); } static mpfr_t r24393, r24394, r24395, r24396, r24397, r24398, r24399, r24400, r24401, r24402, r24403, r24404, r24405, r24406, r24407, r24408, r24409, r24410, r24411, r24412, r24413, r24414, r24415, r24416, r24417, r24418, r24419, r24420, r24421; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24393, "1", 10, MPFR_RNDN); mpfr_init(r24394); mpfr_init(r24395); mpfr_init(r24396); mpfr_init_set_str(r24397, "2", 10, MPFR_RNDN); mpfr_init(r24398); mpfr_init(r24399); mpfr_init(r24400); mpfr_init_set_str(r24401, "3", 10, MPFR_RNDN); mpfr_init(r24402); mpfr_init(r24403); mpfr_init(r24404); mpfr_init(r24405); mpfr_init(r24406); mpfr_init_set_str(r24407, "5", 10, MPFR_RNDN); mpfr_init(r24408); mpfr_init(r24409); mpfr_init(r24410); mpfr_init(r24411); mpfr_init(r24412); mpfr_init(r24413); mpfr_init_set_str(r24414, "21", 10, MPFR_RNDN); mpfr_init(r24415); mpfr_init(r24416); mpfr_init(r24417); mpfr_init(r24418); mpfr_init(r24419); mpfr_init(r24420); mpfr_init(r24421); } double f_fm(double x) { ; mpfr_const_pi(r24394, MPFR_RNDN); mpfr_sqrt(r24395, r24394, MPFR_RNDN); mpfr_div(r24396, r24393, r24395, MPFR_RNDN); ; mpfr_set_d(r24398, x, MPFR_RNDN); mpfr_abs(r24399, r24398, MPFR_RNDN); mpfr_mul(r24400, r24397, r24399, MPFR_RNDN); ; mpfr_div(r24402, r24397, r24401, MPFR_RNDN); mpfr_mul(r24403, r24399, r24399, MPFR_RNDN); mpfr_mul(r24404, r24403, r24399, MPFR_RNDN); mpfr_mul(r24405, r24402, r24404, MPFR_RNDN); mpfr_add(r24406, r24400, r24405, MPFR_RNDN); ; mpfr_div(r24408, r24393, r24407, MPFR_RNDN); mpfr_mul(r24409, r24404, r24399, MPFR_RNDN); mpfr_mul(r24410, r24409, r24399, MPFR_RNDN); mpfr_mul(r24411, r24408, r24410, MPFR_RNDN); mpfr_add(r24412, r24406, r24411, MPFR_RNDN); mpfr_mul(r24413, r24396, r24412, MPFR_RNDN); ; mpfr_div(r24415, r24399, r24414, MPFR_RNDN); mpfr_div(r24416, r24415, r24395, MPFR_RNDN); mpfr_pow(r24417, r24399, r24401, MPFR_RNDN); mpfr_mul(r24418, r24417, r24417, MPFR_RNDN); mpfr_mul(r24419, r24416, r24418, MPFR_RNDN); mpfr_add(r24420, r24413, r24419, MPFR_RNDN); mpfr_abs(r24421, r24420, MPFR_RNDN); return mpfr_get_d(r24421, MPFR_RNDN); } static mpfr_t r24422, r24423, r24424, r24425, r24426, r24427, r24428, r24429, r24430, r24431, r24432, r24433, r24434, r24435, r24436, r24437, r24438, r24439, r24440, r24441, r24442, r24443, r24444, r24445, r24446, r24447, r24448, r24449, r24450; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24422, "1", 10, MPFR_RNDN); mpfr_init(r24423); mpfr_init(r24424); mpfr_init(r24425); mpfr_init_set_str(r24426, "2", 10, MPFR_RNDN); mpfr_init(r24427); mpfr_init(r24428); mpfr_init(r24429); mpfr_init_set_str(r24430, "3", 10, MPFR_RNDN); mpfr_init(r24431); mpfr_init(r24432); mpfr_init(r24433); mpfr_init(r24434); mpfr_init(r24435); mpfr_init_set_str(r24436, "5", 10, MPFR_RNDN); mpfr_init(r24437); mpfr_init(r24438); mpfr_init(r24439); mpfr_init(r24440); mpfr_init(r24441); mpfr_init(r24442); mpfr_init_set_str(r24443, "21", 10, MPFR_RNDN); mpfr_init(r24444); mpfr_init(r24445); mpfr_init(r24446); mpfr_init(r24447); mpfr_init(r24448); mpfr_init(r24449); mpfr_init(r24450); } double f_dm(double x) { ; mpfr_const_pi(r24423, MPFR_RNDN); mpfr_sqrt(r24424, r24423, MPFR_RNDN); mpfr_div(r24425, r24422, r24424, MPFR_RNDN); ; mpfr_set_d(r24427, x, MPFR_RNDN); mpfr_abs(r24428, r24427, MPFR_RNDN); mpfr_mul(r24429, r24426, r24428, MPFR_RNDN); ; mpfr_div(r24431, r24426, r24430, MPFR_RNDN); mpfr_mul(r24432, r24428, r24428, MPFR_RNDN); mpfr_mul(r24433, r24432, r24428, MPFR_RNDN); mpfr_mul(r24434, r24431, r24433, MPFR_RNDN); mpfr_add(r24435, r24429, r24434, MPFR_RNDN); ; mpfr_div(r24437, r24422, r24436, MPFR_RNDN); mpfr_mul(r24438, r24433, r24428, MPFR_RNDN); mpfr_mul(r24439, r24438, r24428, MPFR_RNDN); mpfr_mul(r24440, r24437, r24439, MPFR_RNDN); mpfr_add(r24441, r24435, r24440, MPFR_RNDN); mpfr_mul(r24442, r24425, r24441, MPFR_RNDN); ; mpfr_div(r24444, r24428, r24443, MPFR_RNDN); mpfr_div(r24445, r24444, r24424, MPFR_RNDN); mpfr_pow(r24446, r24428, r24430, MPFR_RNDN); mpfr_mul(r24447, r24446, r24446, MPFR_RNDN); mpfr_mul(r24448, r24445, r24447, MPFR_RNDN); mpfr_add(r24449, r24442, r24448, MPFR_RNDN); mpfr_abs(r24450, r24449, MPFR_RNDN); return mpfr_get_d(r24450, MPFR_RNDN); }