#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 r23329 = 1; float r23330 = atan2(1.0, 0.0); float r23331 = sqrt(r23330); float r23332 = r23329 / r23331; float r23333 = 2; float r23334 = x; float r23335 = fabs(r23334); float r23336 = r23333 * r23335; float r23337 = 3; float r23338 = r23333 / r23337; float r23339 = r23335 * r23335; float r23340 = r23339 * r23335; float r23341 = r23338 * r23340; float r23342 = r23336 + r23341; float r23343 = 5; float r23344 = r23329 / r23343; float r23345 = r23340 * r23335; float r23346 = r23345 * r23335; float r23347 = r23344 * r23346; float r23348 = r23342 + r23347; float r23349 = 21; float r23350 = r23329 / r23349; float r23351 = r23346 * r23335; float r23352 = r23351 * r23335; float r23353 = r23350 * r23352; float r23354 = r23348 + r23353; float r23355 = r23332 * r23354; float r23356 = fabs(r23355); return r23356; } double f_id(double x) { double r23357 = 1; double r23358 = atan2(1.0, 0.0); double r23359 = sqrt(r23358); double r23360 = r23357 / r23359; double r23361 = 2; double r23362 = x; double r23363 = fabs(r23362); double r23364 = r23361 * r23363; double r23365 = 3; double r23366 = r23361 / r23365; double r23367 = r23363 * r23363; double r23368 = r23367 * r23363; double r23369 = r23366 * r23368; double r23370 = r23364 + r23369; double r23371 = 5; double r23372 = r23357 / r23371; double r23373 = r23368 * r23363; double r23374 = r23373 * r23363; double r23375 = r23372 * r23374; double r23376 = r23370 + r23375; double r23377 = 21; double r23378 = r23357 / r23377; double r23379 = r23374 * r23363; double r23380 = r23379 * r23363; double r23381 = r23378 * r23380; double r23382 = r23376 + r23381; double r23383 = r23360 * r23382; double r23384 = fabs(r23383); return r23384; } double f_of(float x) { float r23385 = 1; float r23386 = atan2(1.0, 0.0); float r23387 = sqrt(r23386); float r23388 = r23385 / r23387; float r23389 = 2; float r23390 = x; float r23391 = fabs(r23390); float r23392 = r23389 * r23391; float r23393 = 3; float r23394 = r23389 / r23393; float r23395 = r23391 * r23391; float r23396 = r23395 * r23391; float r23397 = r23394 * r23396; float r23398 = r23392 + r23397; float r23399 = 5; float r23400 = r23385 / r23399; float r23401 = r23396 * r23391; float r23402 = r23401 * r23391; float r23403 = r23400 * r23402; float r23404 = r23398 + r23403; float r23405 = 21; float r23406 = r23385 / r23405; float r23407 = r23402 * r23391; float r23408 = r23407 * r23391; float r23409 = r23406 * r23408; float r23410 = r23404 + r23409; float r23411 = r23388 * r23410; float r23412 = fabs(r23411); return r23412; } double f_od(double x) { double r23413 = 1; double r23414 = atan2(1.0, 0.0); double r23415 = sqrt(r23414); double r23416 = r23413 / r23415; double r23417 = 2; double r23418 = x; double r23419 = fabs(r23418); double r23420 = r23417 * r23419; double r23421 = 3; double r23422 = r23417 / r23421; double r23423 = r23419 * r23419; double r23424 = r23423 * r23419; double r23425 = r23422 * r23424; double r23426 = r23420 + r23425; double r23427 = 5; double r23428 = r23413 / r23427; double r23429 = r23424 * r23419; double r23430 = r23429 * r23419; double r23431 = r23428 * r23430; double r23432 = r23426 + r23431; double r23433 = 21; double r23434 = r23413 / r23433; double r23435 = r23430 * r23419; double r23436 = r23435 * r23419; double r23437 = r23434 * r23436; double r23438 = r23432 + r23437; double r23439 = r23416 * r23438; double r23440 = fabs(r23439); return r23440; } 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 r23441, r23442, r23443, r23444, r23445, r23446, r23447, r23448, r23449, r23450, r23451, r23452, r23453, r23454, r23455, r23456, r23457, r23458, r23459, r23460, r23461, r23462, r23463, r23464, r23465, r23466, r23467, r23468; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r23441, "1", 10, MPFR_RNDN); mpfr_init(r23442); mpfr_init(r23443); mpfr_init(r23444); mpfr_init_set_str(r23445, "2", 10, MPFR_RNDN); mpfr_init(r23446); mpfr_init(r23447); mpfr_init(r23448); mpfr_init_set_str(r23449, "3", 10, MPFR_RNDN); mpfr_init(r23450); mpfr_init(r23451); mpfr_init(r23452); mpfr_init(r23453); mpfr_init(r23454); mpfr_init_set_str(r23455, "5", 10, MPFR_RNDN); mpfr_init(r23456); mpfr_init(r23457); mpfr_init(r23458); mpfr_init(r23459); mpfr_init(r23460); mpfr_init_set_str(r23461, "21", 10, MPFR_RNDN); mpfr_init(r23462); mpfr_init(r23463); mpfr_init(r23464); mpfr_init(r23465); mpfr_init(r23466); mpfr_init(r23467); mpfr_init(r23468); } double f_im(double x) { ; mpfr_const_pi(r23442, MPFR_RNDN); mpfr_sqrt(r23443, r23442, MPFR_RNDN); mpfr_div(r23444, r23441, r23443, MPFR_RNDN); ; mpfr_set_d(r23446, x, MPFR_RNDN); mpfr_abs(r23447, r23446, MPFR_RNDN); mpfr_mul(r23448, r23445, r23447, MPFR_RNDN); ; mpfr_div(r23450, r23445, r23449, MPFR_RNDN); mpfr_mul(r23451, r23447, r23447, MPFR_RNDN); mpfr_mul(r23452, r23451, r23447, MPFR_RNDN); mpfr_mul(r23453, r23450, r23452, MPFR_RNDN); mpfr_add(r23454, r23448, r23453, MPFR_RNDN); ; mpfr_div(r23456, r23441, r23455, MPFR_RNDN); mpfr_mul(r23457, r23452, r23447, MPFR_RNDN); mpfr_mul(r23458, r23457, r23447, MPFR_RNDN); mpfr_mul(r23459, r23456, r23458, MPFR_RNDN); mpfr_add(r23460, r23454, r23459, MPFR_RNDN); ; mpfr_div(r23462, r23441, r23461, MPFR_RNDN); mpfr_mul(r23463, r23458, r23447, MPFR_RNDN); mpfr_mul(r23464, r23463, r23447, MPFR_RNDN); mpfr_mul(r23465, r23462, r23464, MPFR_RNDN); mpfr_add(r23466, r23460, r23465, MPFR_RNDN); mpfr_mul(r23467, r23444, r23466, MPFR_RNDN); mpfr_abs(r23468, r23467, MPFR_RNDN); return mpfr_get_d(r23468, MPFR_RNDN); } static mpfr_t r23469, r23470, r23471, r23472, r23473, r23474, r23475, r23476, r23477, r23478, r23479, r23480, r23481, r23482, r23483, r23484, r23485, r23486, r23487, r23488, r23489, r23490, r23491, r23492, r23493, r23494, r23495, r23496; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23469, "1", 10, MPFR_RNDN); mpfr_init(r23470); mpfr_init(r23471); mpfr_init(r23472); mpfr_init_set_str(r23473, "2", 10, MPFR_RNDN); mpfr_init(r23474); mpfr_init(r23475); mpfr_init(r23476); mpfr_init_set_str(r23477, "3", 10, MPFR_RNDN); mpfr_init(r23478); mpfr_init(r23479); mpfr_init(r23480); mpfr_init(r23481); mpfr_init(r23482); mpfr_init_set_str(r23483, "5", 10, MPFR_RNDN); mpfr_init(r23484); mpfr_init(r23485); mpfr_init(r23486); mpfr_init(r23487); mpfr_init(r23488); mpfr_init_set_str(r23489, "21", 10, MPFR_RNDN); mpfr_init(r23490); mpfr_init(r23491); mpfr_init(r23492); mpfr_init(r23493); mpfr_init(r23494); mpfr_init(r23495); mpfr_init(r23496); } double f_fm(double x) { ; mpfr_const_pi(r23470, MPFR_RNDN); mpfr_sqrt(r23471, r23470, MPFR_RNDN); mpfr_div(r23472, r23469, r23471, MPFR_RNDN); ; mpfr_set_d(r23474, x, MPFR_RNDN); mpfr_abs(r23475, r23474, MPFR_RNDN); mpfr_mul(r23476, r23473, r23475, MPFR_RNDN); ; mpfr_div(r23478, r23473, r23477, MPFR_RNDN); mpfr_mul(r23479, r23475, r23475, MPFR_RNDN); mpfr_mul(r23480, r23479, r23475, MPFR_RNDN); mpfr_mul(r23481, r23478, r23480, MPFR_RNDN); mpfr_add(r23482, r23476, r23481, MPFR_RNDN); ; mpfr_div(r23484, r23469, r23483, MPFR_RNDN); mpfr_mul(r23485, r23480, r23475, MPFR_RNDN); mpfr_mul(r23486, r23485, r23475, MPFR_RNDN); mpfr_mul(r23487, r23484, r23486, MPFR_RNDN); mpfr_add(r23488, r23482, r23487, MPFR_RNDN); ; mpfr_div(r23490, r23469, r23489, MPFR_RNDN); mpfr_mul(r23491, r23486, r23475, MPFR_RNDN); mpfr_mul(r23492, r23491, r23475, MPFR_RNDN); mpfr_mul(r23493, r23490, r23492, MPFR_RNDN); mpfr_add(r23494, r23488, r23493, MPFR_RNDN); mpfr_mul(r23495, r23472, r23494, MPFR_RNDN); mpfr_abs(r23496, r23495, MPFR_RNDN); return mpfr_get_d(r23496, MPFR_RNDN); } static mpfr_t r23497, r23498, r23499, r23500, r23501, r23502, r23503, r23504, r23505, r23506, r23507, r23508, r23509, r23510, r23511, r23512, r23513, r23514, r23515, r23516, r23517, r23518, r23519, r23520, r23521, r23522, r23523, r23524; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23497, "1", 10, MPFR_RNDN); mpfr_init(r23498); mpfr_init(r23499); mpfr_init(r23500); mpfr_init_set_str(r23501, "2", 10, MPFR_RNDN); mpfr_init(r23502); mpfr_init(r23503); mpfr_init(r23504); mpfr_init_set_str(r23505, "3", 10, MPFR_RNDN); mpfr_init(r23506); mpfr_init(r23507); mpfr_init(r23508); mpfr_init(r23509); mpfr_init(r23510); mpfr_init_set_str(r23511, "5", 10, MPFR_RNDN); mpfr_init(r23512); mpfr_init(r23513); mpfr_init(r23514); mpfr_init(r23515); mpfr_init(r23516); mpfr_init_set_str(r23517, "21", 10, MPFR_RNDN); mpfr_init(r23518); mpfr_init(r23519); mpfr_init(r23520); mpfr_init(r23521); mpfr_init(r23522); mpfr_init(r23523); mpfr_init(r23524); } double f_dm(double x) { ; mpfr_const_pi(r23498, MPFR_RNDN); mpfr_sqrt(r23499, r23498, MPFR_RNDN); mpfr_div(r23500, r23497, r23499, MPFR_RNDN); ; mpfr_set_d(r23502, x, MPFR_RNDN); mpfr_abs(r23503, r23502, MPFR_RNDN); mpfr_mul(r23504, r23501, r23503, MPFR_RNDN); ; mpfr_div(r23506, r23501, r23505, MPFR_RNDN); mpfr_mul(r23507, r23503, r23503, MPFR_RNDN); mpfr_mul(r23508, r23507, r23503, MPFR_RNDN); mpfr_mul(r23509, r23506, r23508, MPFR_RNDN); mpfr_add(r23510, r23504, r23509, MPFR_RNDN); ; mpfr_div(r23512, r23497, r23511, MPFR_RNDN); mpfr_mul(r23513, r23508, r23503, MPFR_RNDN); mpfr_mul(r23514, r23513, r23503, MPFR_RNDN); mpfr_mul(r23515, r23512, r23514, MPFR_RNDN); mpfr_add(r23516, r23510, r23515, MPFR_RNDN); ; mpfr_div(r23518, r23497, r23517, MPFR_RNDN); mpfr_mul(r23519, r23514, r23503, MPFR_RNDN); mpfr_mul(r23520, r23519, r23503, MPFR_RNDN); mpfr_mul(r23521, r23518, r23520, MPFR_RNDN); mpfr_add(r23522, r23516, r23521, MPFR_RNDN); mpfr_mul(r23523, r23500, r23522, MPFR_RNDN); mpfr_abs(r23524, r23523, MPFR_RNDN); return mpfr_get_d(r23524, MPFR_RNDN); }