#include #include #include #include #include char *name = "_divideComplex, real part"; double f_if(float x_re, float x_im, float y_re, float y_im) { float r9440 = x_re; float r9441 = y_re; float r9442 = r9440 * r9441; float r9443 = x_im; float r9444 = y_im; float r9445 = r9443 * r9444; float r9446 = r9442 + r9445; float r9447 = r9441 * r9441; float r9448 = r9444 * r9444; float r9449 = r9447 + r9448; float r9450 = r9446 / r9449; return r9450; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r9451 = x_re; double r9452 = y_re; double r9453 = r9451 * r9452; double r9454 = x_im; double r9455 = y_im; double r9456 = r9454 * r9455; double r9457 = r9453 + r9456; double r9458 = r9452 * r9452; double r9459 = r9455 * r9455; double r9460 = r9458 + r9459; double r9461 = r9457 / r9460; return r9461; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r9462 = 1; float r9463 = y_re; float r9464 = r9463 * r9463; float r9465 = y_im; float r9466 = r9465 * r9465; float r9467 = r9464 + r9466; float r9468 = sqrt(r9467); float r9469 = r9462 / r9468; float r9470 = x_re; float r9471 = r9470 * r9463; float r9472 = x_im; float r9473 = r9472 * r9465; float r9474 = r9471 + r9473; float r9475 = r9474 / r9468; float r9476 = r9469 * r9475; return r9476; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r9477 = 1; double r9478 = y_re; double r9479 = r9478 * r9478; double r9480 = y_im; double r9481 = r9480 * r9480; double r9482 = r9479 + r9481; double r9483 = sqrt(r9482); double r9484 = r9477 / r9483; double r9485 = x_re; double r9486 = r9485 * r9478; double r9487 = x_im; double r9488 = r9487 * r9480; double r9489 = r9486 + r9488; double r9490 = r9489 / r9483; double r9491 = r9484 * r9490; return r9491; } 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 r9492, r9493, r9494, r9495, r9496, r9497, r9498, r9499, r9500, r9501, r9502; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r9492); mpfr_init(r9493); mpfr_init(r9494); mpfr_init(r9495); mpfr_init(r9496); mpfr_init(r9497); mpfr_init(r9498); mpfr_init(r9499); mpfr_init(r9500); mpfr_init(r9501); mpfr_init(r9502); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r9492, x_re, MPFR_RNDN); mpfr_set_d(r9493, y_re, MPFR_RNDN); mpfr_mul(r9494, r9492, r9493, MPFR_RNDN); mpfr_set_d(r9495, x_im, MPFR_RNDN); mpfr_set_d(r9496, y_im, MPFR_RNDN); mpfr_mul(r9497, r9495, r9496, MPFR_RNDN); mpfr_add(r9498, r9494, r9497, MPFR_RNDN); mpfr_mul(r9499, r9493, r9493, MPFR_RNDN); mpfr_mul(r9500, r9496, r9496, MPFR_RNDN); mpfr_add(r9501, r9499, r9500, MPFR_RNDN); mpfr_div(r9502, r9498, r9501, MPFR_RNDN); return mpfr_get_d(r9502, MPFR_RNDN); } static mpfr_t r9503, r9504, r9505, r9506, r9507, r9508, r9509, r9510, r9511, r9512, r9513, r9514, r9515, r9516, r9517; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r9503, "1", 10, MPFR_RNDN); mpfr_init(r9504); mpfr_init(r9505); mpfr_init(r9506); mpfr_init(r9507); mpfr_init(r9508); mpfr_init(r9509); mpfr_init(r9510); mpfr_init(r9511); mpfr_init(r9512); mpfr_init(r9513); mpfr_init(r9514); mpfr_init(r9515); mpfr_init(r9516); mpfr_init(r9517); } double f_fm(double x_re, double x_im, double y_re, double y_im) { ; mpfr_set_d(r9504, y_re, MPFR_RNDN); mpfr_mul(r9505, r9504, r9504, MPFR_RNDN); mpfr_set_d(r9506, y_im, MPFR_RNDN); mpfr_mul(r9507, r9506, r9506, MPFR_RNDN); mpfr_add(r9508, r9505, r9507, MPFR_RNDN); mpfr_sqrt(r9509, r9508, MPFR_RNDN); mpfr_div(r9510, r9503, r9509, MPFR_RNDN); mpfr_set_d(r9511, x_re, MPFR_RNDN); mpfr_mul(r9512, r9511, r9504, MPFR_RNDN); mpfr_set_d(r9513, x_im, MPFR_RNDN); mpfr_mul(r9514, r9513, r9506, MPFR_RNDN); mpfr_add(r9515, r9512, r9514, MPFR_RNDN); mpfr_div(r9516, r9515, r9509, MPFR_RNDN); mpfr_mul(r9517, r9510, r9516, MPFR_RNDN); return mpfr_get_d(r9517, MPFR_RNDN); } static mpfr_t r9518, r9519, r9520, r9521, r9522, r9523, r9524, r9525, r9526, r9527, r9528, r9529, r9530, r9531, r9532; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r9518, "1", 10, MPFR_RNDN); mpfr_init(r9519); mpfr_init(r9520); mpfr_init(r9521); mpfr_init(r9522); mpfr_init(r9523); mpfr_init(r9524); mpfr_init(r9525); mpfr_init(r9526); mpfr_init(r9527); mpfr_init(r9528); mpfr_init(r9529); mpfr_init(r9530); mpfr_init(r9531); mpfr_init(r9532); } double f_dm(double x_re, double x_im, double y_re, double y_im) { ; mpfr_set_d(r9519, y_re, MPFR_RNDN); mpfr_mul(r9520, r9519, r9519, MPFR_RNDN); mpfr_set_d(r9521, y_im, MPFR_RNDN); mpfr_mul(r9522, r9521, r9521, MPFR_RNDN); mpfr_add(r9523, r9520, r9522, MPFR_RNDN); mpfr_sqrt(r9524, r9523, MPFR_RNDN); mpfr_div(r9525, r9518, r9524, MPFR_RNDN); mpfr_set_d(r9526, x_re, MPFR_RNDN); mpfr_mul(r9527, r9526, r9519, MPFR_RNDN); mpfr_set_d(r9528, x_im, MPFR_RNDN); mpfr_mul(r9529, r9528, r9521, MPFR_RNDN); mpfr_add(r9530, r9527, r9529, MPFR_RNDN); mpfr_div(r9531, r9530, r9524, MPFR_RNDN); mpfr_mul(r9532, r9525, r9531, MPFR_RNDN); return mpfr_get_d(r9532, MPFR_RNDN); }