#include #include #include #include #include char *name = "Complex division, imag part"; double f_if(float a, float b, float c, float d) { float r23367 = b; float r23368 = c; float r23369 = r23367 * r23368; float r23370 = a; float r23371 = d; float r23372 = r23370 * r23371; float r23373 = r23369 - r23372; float r23374 = r23368 * r23368; float r23375 = r23371 * r23371; float r23376 = r23374 + r23375; float r23377 = r23373 / r23376; return r23377; } double f_id(double a, double b, double c, double d) { double r23378 = b; double r23379 = c; double r23380 = r23378 * r23379; double r23381 = a; double r23382 = d; double r23383 = r23381 * r23382; double r23384 = r23380 - r23383; double r23385 = r23379 * r23379; double r23386 = r23382 * r23382; double r23387 = r23385 + r23386; double r23388 = r23384 / r23387; return r23388; } double f_of(float a, float b, float c, float d) { float r23389 = c; float r23390 = -2.9830793252751902e+141; bool r23391 = r23389 <= r23390; float r23392 = a; float r23393 = d; float r23394 = r23393 / r23389; float r23395 = b; float r23396 = -r23395; float r23397 = fma(r23392, r23394, r23396); float r23398 = hypot(r23393, r23389); float r23399 = r23397 / r23398; float r23400 = 7.202117007000063e+42; bool r23401 = r23389 <= r23400; float r23402 = r23389 * r23395; float r23403 = r23392 * r23393; float r23404 = r23402 - r23403; float r23405 = r23404 / r23398; float r23406 = r23405 / r23398; float r23407 = -r23392; float r23408 = fma(r23394, r23407, r23395); float r23409 = r23408 / r23398; float r23410 = r23401 ? r23406 : r23409; float r23411 = r23391 ? r23399 : r23410; return r23411; } double f_od(double a, double b, double c, double d) { double r23412 = c; double r23413 = -2.9830793252751902e+141; bool r23414 = r23412 <= r23413; double r23415 = a; double r23416 = d; double r23417 = r23416 / r23412; double r23418 = b; double r23419 = -r23418; double r23420 = fma(r23415, r23417, r23419); double r23421 = hypot(r23416, r23412); double r23422 = r23420 / r23421; double r23423 = 7.202117007000063e+42; bool r23424 = r23412 <= r23423; double r23425 = r23412 * r23418; double r23426 = r23415 * r23416; double r23427 = r23425 - r23426; double r23428 = r23427 / r23421; double r23429 = r23428 / r23421; double r23430 = -r23415; double r23431 = fma(r23417, r23430, r23418); double r23432 = r23431 / r23421; double r23433 = r23424 ? r23429 : r23432; double r23434 = r23414 ? r23422 : r23433; return r23434; } 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 r23435, r23436, r23437, r23438, r23439, r23440, r23441, r23442, r23443, r23444, r23445; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r23435); mpfr_init(r23436); mpfr_init(r23437); mpfr_init(r23438); mpfr_init(r23439); mpfr_init(r23440); mpfr_init(r23441); mpfr_init(r23442); mpfr_init(r23443); mpfr_init(r23444); mpfr_init(r23445); } double f_im(double a, double b, double c, double d) { mpfr_set_d(r23435, b, MPFR_RNDN); mpfr_set_d(r23436, c, MPFR_RNDN); mpfr_mul(r23437, r23435, r23436, MPFR_RNDN); mpfr_set_d(r23438, a, MPFR_RNDN); mpfr_set_d(r23439, d, MPFR_RNDN); mpfr_mul(r23440, r23438, r23439, MPFR_RNDN); mpfr_sub(r23441, r23437, r23440, MPFR_RNDN); mpfr_mul(r23442, r23436, r23436, MPFR_RNDN); mpfr_mul(r23443, r23439, r23439, MPFR_RNDN); mpfr_add(r23444, r23442, r23443, MPFR_RNDN); mpfr_div(r23445, r23441, r23444, MPFR_RNDN); return mpfr_get_d(r23445, MPFR_RNDN); } static mpfr_t 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_fm() { mpfr_set_default_prec(336); mpfr_init(r23446); mpfr_init_set_str(r23447, "-2.9830793252751902e+141", 10, MPFR_RNDN); mpfr_init(r23448); mpfr_init(r23449); mpfr_init(r23450); mpfr_init(r23451); mpfr_init(r23452); mpfr_init(r23453); mpfr_init(r23454); mpfr_init(r23455); mpfr_init(r23456); mpfr_init_set_str(r23457, "7.202117007000063e+42", 10, MPFR_RNDN); mpfr_init(r23458); mpfr_init(r23459); mpfr_init(r23460); mpfr_init(r23461); mpfr_init(r23462); mpfr_init(r23463); mpfr_init(r23464); mpfr_init(r23465); mpfr_init(r23466); mpfr_init(r23467); mpfr_init(r23468); } double f_fm(double a, double b, double c, double d) { mpfr_set_d(r23446, c, MPFR_RNDN); ; mpfr_set_si(r23448, mpfr_cmp(r23446, r23447) <= 0, MPFR_RNDN); mpfr_set_d(r23449, a, MPFR_RNDN); mpfr_set_d(r23450, d, MPFR_RNDN); mpfr_div(r23451, r23450, r23446, MPFR_RNDN); mpfr_set_d(r23452, b, MPFR_RNDN); mpfr_neg(r23453, r23452, MPFR_RNDN); mpfr_fma(r23454, r23449, r23451, r23453, MPFR_RNDN); mpfr_hypot(r23455, r23450, r23446, MPFR_RNDN); mpfr_div(r23456, r23454, r23455, MPFR_RNDN); ; mpfr_set_si(r23458, mpfr_cmp(r23446, r23457) <= 0, MPFR_RNDN); mpfr_mul(r23459, r23446, r23452, MPFR_RNDN); mpfr_mul(r23460, r23449, r23450, MPFR_RNDN); mpfr_sub(r23461, r23459, r23460, MPFR_RNDN); mpfr_div(r23462, r23461, r23455, MPFR_RNDN); mpfr_div(r23463, r23462, r23455, MPFR_RNDN); mpfr_neg(r23464, r23449, MPFR_RNDN); mpfr_fma(r23465, r23451, r23464, r23452, MPFR_RNDN); mpfr_div(r23466, r23465, r23455, MPFR_RNDN); if (mpfr_get_si(r23458, MPFR_RNDN)) { mpfr_set(r23467, r23463, MPFR_RNDN); } else { mpfr_set(r23467, r23466, MPFR_RNDN); }; if (mpfr_get_si(r23448, MPFR_RNDN)) { mpfr_set(r23468, r23456, MPFR_RNDN); } else { mpfr_set(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; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r23469); mpfr_init_set_str(r23470, "-2.9830793252751902e+141", 10, MPFR_RNDN); mpfr_init(r23471); mpfr_init(r23472); mpfr_init(r23473); mpfr_init(r23474); mpfr_init(r23475); mpfr_init(r23476); mpfr_init(r23477); mpfr_init(r23478); mpfr_init(r23479); mpfr_init_set_str(r23480, "7.202117007000063e+42", 10, MPFR_RNDN); mpfr_init(r23481); mpfr_init(r23482); mpfr_init(r23483); mpfr_init(r23484); mpfr_init(r23485); mpfr_init(r23486); mpfr_init(r23487); mpfr_init(r23488); mpfr_init(r23489); mpfr_init(r23490); mpfr_init(r23491); } double f_dm(double a, double b, double c, double d) { mpfr_set_d(r23469, c, MPFR_RNDN); ; mpfr_set_si(r23471, mpfr_cmp(r23469, r23470) <= 0, MPFR_RNDN); mpfr_set_d(r23472, a, MPFR_RNDN); mpfr_set_d(r23473, d, MPFR_RNDN); mpfr_div(r23474, r23473, r23469, MPFR_RNDN); mpfr_set_d(r23475, b, MPFR_RNDN); mpfr_neg(r23476, r23475, MPFR_RNDN); mpfr_fma(r23477, r23472, r23474, r23476, MPFR_RNDN); mpfr_hypot(r23478, r23473, r23469, MPFR_RNDN); mpfr_div(r23479, r23477, r23478, MPFR_RNDN); ; mpfr_set_si(r23481, mpfr_cmp(r23469, r23480) <= 0, MPFR_RNDN); mpfr_mul(r23482, r23469, r23475, MPFR_RNDN); mpfr_mul(r23483, r23472, r23473, MPFR_RNDN); mpfr_sub(r23484, r23482, r23483, MPFR_RNDN); mpfr_div(r23485, r23484, r23478, MPFR_RNDN); mpfr_div(r23486, r23485, r23478, MPFR_RNDN); mpfr_neg(r23487, r23472, MPFR_RNDN); mpfr_fma(r23488, r23474, r23487, r23475, MPFR_RNDN); mpfr_div(r23489, r23488, r23478, MPFR_RNDN); if (mpfr_get_si(r23481, MPFR_RNDN)) { mpfr_set(r23490, r23486, MPFR_RNDN); } else { mpfr_set(r23490, r23489, MPFR_RNDN); }; if (mpfr_get_si(r23471, MPFR_RNDN)) { mpfr_set(r23491, r23479, MPFR_RNDN); } else { mpfr_set(r23491, r23490, MPFR_RNDN); }; return mpfr_get_d(r23491, MPFR_RNDN); }