#include #include #include #include #include char *name = "Complex division, imag part"; double f_if(float a, float b, float c, float d) { float r26687 = b; float r26688 = c; float r26689 = r26687 * r26688; float r26690 = a; float r26691 = d; float r26692 = r26690 * r26691; float r26693 = r26689 - r26692; float r26694 = r26688 * r26688; float r26695 = r26691 * r26691; float r26696 = r26694 + r26695; float r26697 = r26693 / r26696; return r26697; } double f_id(double a, double b, double c, double d) { double r26698 = b; double r26699 = c; double r26700 = r26698 * r26699; double r26701 = a; double r26702 = d; double r26703 = r26701 * r26702; double r26704 = r26700 - r26703; double r26705 = r26699 * r26699; double r26706 = r26702 * r26702; double r26707 = r26705 + r26706; double r26708 = r26704 / r26707; return r26708; } double f_of(float a, float b, float c, float d) { float r26709 = c; float r26710 = -1.019095860764281e+69; bool r26711 = r26709 <= r26710; float r26712 = b; float r26713 = -r26712; float r26714 = r26709 * r26709; float r26715 = d; float r26716 = r26715 * r26715; float r26717 = r26714 + r26716; float r26718 = sqrt(r26717); float r26719 = r26713 / r26718; float r26720 = 1.7442473790874188e+37; bool r26721 = r26709 <= r26720; float r26722 = r26712 * r26709; float r26723 = a; float r26724 = r26723 * r26715; float r26725 = r26722 - r26724; float r26726 = r26725 / r26718; float r26727 = r26726 / r26718; float r26728 = r26709 / r26723; float r26729 = r26715 / r26728; float r26730 = r26712 - r26729; float r26731 = r26716 + r26714; float r26732 = sqrt(r26731); float r26733 = r26730 / r26732; float r26734 = r26721 ? r26727 : r26733; float r26735 = r26711 ? r26719 : r26734; return r26735; } double f_od(double a, double b, double c, double d) { double r26736 = c; double r26737 = -1.019095860764281e+69; bool r26738 = r26736 <= r26737; double r26739 = b; double r26740 = -r26739; double r26741 = r26736 * r26736; double r26742 = d; double r26743 = r26742 * r26742; double r26744 = r26741 + r26743; double r26745 = sqrt(r26744); double r26746 = r26740 / r26745; double r26747 = 1.7442473790874188e+37; bool r26748 = r26736 <= r26747; double r26749 = r26739 * r26736; double r26750 = a; double r26751 = r26750 * r26742; double r26752 = r26749 - r26751; double r26753 = r26752 / r26745; double r26754 = r26753 / r26745; double r26755 = r26736 / r26750; double r26756 = r26742 / r26755; double r26757 = r26739 - r26756; double r26758 = r26743 + r26741; double r26759 = sqrt(r26758); double r26760 = r26757 / r26759; double r26761 = r26748 ? r26754 : r26760; double r26762 = r26738 ? r26746 : r26761; return r26762; } 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 r26763, r26764, r26765, r26766, r26767, r26768, r26769, r26770, r26771, r26772, r26773; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26763); mpfr_init(r26764); mpfr_init(r26765); mpfr_init(r26766); mpfr_init(r26767); mpfr_init(r26768); mpfr_init(r26769); mpfr_init(r26770); mpfr_init(r26771); mpfr_init(r26772); mpfr_init(r26773); } double f_im(double a, double b, double c, double d) { mpfr_set_d(r26763, b, MPFR_RNDN); mpfr_set_d(r26764, c, MPFR_RNDN); mpfr_mul(r26765, r26763, r26764, MPFR_RNDN); mpfr_set_d(r26766, a, MPFR_RNDN); mpfr_set_d(r26767, d, MPFR_RNDN); mpfr_mul(r26768, r26766, r26767, MPFR_RNDN); mpfr_sub(r26769, r26765, r26768, MPFR_RNDN); mpfr_mul(r26770, r26764, r26764, MPFR_RNDN); mpfr_mul(r26771, r26767, r26767, MPFR_RNDN); mpfr_add(r26772, r26770, r26771, MPFR_RNDN); mpfr_div(r26773, r26769, r26772, MPFR_RNDN); return mpfr_get_d(r26773, MPFR_RNDN); } static mpfr_t r26774, r26775, r26776, r26777, r26778, r26779, r26780, r26781, r26782, r26783, r26784, r26785, r26786, r26787, r26788, r26789, r26790, r26791, r26792, r26793, r26794, r26795, r26796, r26797, r26798, r26799, r26800; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26774); mpfr_init_set_str(r26775, "-1.019095860764281e+69", 10, MPFR_RNDN); mpfr_init(r26776); mpfr_init(r26777); mpfr_init(r26778); mpfr_init(r26779); mpfr_init(r26780); mpfr_init(r26781); mpfr_init(r26782); mpfr_init(r26783); mpfr_init(r26784); mpfr_init_set_str(r26785, "1.7442473790874188e+37", 10, MPFR_RNDN); mpfr_init(r26786); mpfr_init(r26787); mpfr_init(r26788); mpfr_init(r26789); mpfr_init(r26790); mpfr_init(r26791); mpfr_init(r26792); mpfr_init(r26793); mpfr_init(r26794); mpfr_init(r26795); mpfr_init(r26796); mpfr_init(r26797); mpfr_init(r26798); mpfr_init(r26799); mpfr_init(r26800); } double f_fm(double a, double b, double c, double d) { mpfr_set_d(r26774, c, MPFR_RNDN); ; mpfr_set_si(r26776, mpfr_cmp(r26774, r26775) <= 0, MPFR_RNDN); mpfr_set_d(r26777, b, MPFR_RNDN); mpfr_neg(r26778, r26777, MPFR_RNDN); mpfr_mul(r26779, r26774, r26774, MPFR_RNDN); mpfr_set_d(r26780, d, MPFR_RNDN); mpfr_mul(r26781, r26780, r26780, MPFR_RNDN); mpfr_add(r26782, r26779, r26781, MPFR_RNDN); mpfr_sqrt(r26783, r26782, MPFR_RNDN); mpfr_div(r26784, r26778, r26783, MPFR_RNDN); ; mpfr_set_si(r26786, mpfr_cmp(r26774, r26785) <= 0, MPFR_RNDN); mpfr_mul(r26787, r26777, r26774, MPFR_RNDN); mpfr_set_d(r26788, a, MPFR_RNDN); mpfr_mul(r26789, r26788, r26780, MPFR_RNDN); mpfr_sub(r26790, r26787, r26789, MPFR_RNDN); mpfr_div(r26791, r26790, r26783, MPFR_RNDN); mpfr_div(r26792, r26791, r26783, MPFR_RNDN); mpfr_div(r26793, r26774, r26788, MPFR_RNDN); mpfr_div(r26794, r26780, r26793, MPFR_RNDN); mpfr_sub(r26795, r26777, r26794, MPFR_RNDN); mpfr_add(r26796, r26781, r26779, MPFR_RNDN); mpfr_sqrt(r26797, r26796, MPFR_RNDN); mpfr_div(r26798, r26795, r26797, MPFR_RNDN); if (mpfr_get_si(r26786, MPFR_RNDN)) { mpfr_set(r26799, r26792, MPFR_RNDN); } else { mpfr_set(r26799, r26798, MPFR_RNDN); }; if (mpfr_get_si(r26776, MPFR_RNDN)) { mpfr_set(r26800, r26784, MPFR_RNDN); } else { mpfr_set(r26800, r26799, MPFR_RNDN); }; return mpfr_get_d(r26800, MPFR_RNDN); } static mpfr_t r26801, r26802, r26803, r26804, r26805, r26806, r26807, r26808, r26809, r26810, r26811, r26812, r26813, r26814, r26815, r26816, r26817, r26818, r26819, r26820, r26821, r26822, r26823, r26824, r26825, r26826, r26827; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26801); mpfr_init_set_str(r26802, "-1.019095860764281e+69", 10, MPFR_RNDN); mpfr_init(r26803); mpfr_init(r26804); mpfr_init(r26805); mpfr_init(r26806); mpfr_init(r26807); mpfr_init(r26808); mpfr_init(r26809); mpfr_init(r26810); mpfr_init(r26811); mpfr_init_set_str(r26812, "1.7442473790874188e+37", 10, MPFR_RNDN); mpfr_init(r26813); mpfr_init(r26814); mpfr_init(r26815); mpfr_init(r26816); mpfr_init(r26817); mpfr_init(r26818); mpfr_init(r26819); mpfr_init(r26820); mpfr_init(r26821); mpfr_init(r26822); mpfr_init(r26823); mpfr_init(r26824); mpfr_init(r26825); mpfr_init(r26826); mpfr_init(r26827); } double f_dm(double a, double b, double c, double d) { mpfr_set_d(r26801, c, MPFR_RNDN); ; mpfr_set_si(r26803, mpfr_cmp(r26801, r26802) <= 0, MPFR_RNDN); mpfr_set_d(r26804, b, MPFR_RNDN); mpfr_neg(r26805, r26804, MPFR_RNDN); mpfr_mul(r26806, r26801, r26801, MPFR_RNDN); mpfr_set_d(r26807, d, MPFR_RNDN); mpfr_mul(r26808, r26807, r26807, MPFR_RNDN); mpfr_add(r26809, r26806, r26808, MPFR_RNDN); mpfr_sqrt(r26810, r26809, MPFR_RNDN); mpfr_div(r26811, r26805, r26810, MPFR_RNDN); ; mpfr_set_si(r26813, mpfr_cmp(r26801, r26812) <= 0, MPFR_RNDN); mpfr_mul(r26814, r26804, r26801, MPFR_RNDN); mpfr_set_d(r26815, a, MPFR_RNDN); mpfr_mul(r26816, r26815, r26807, MPFR_RNDN); mpfr_sub(r26817, r26814, r26816, MPFR_RNDN); mpfr_div(r26818, r26817, r26810, MPFR_RNDN); mpfr_div(r26819, r26818, r26810, MPFR_RNDN); mpfr_div(r26820, r26801, r26815, MPFR_RNDN); mpfr_div(r26821, r26807, r26820, MPFR_RNDN); mpfr_sub(r26822, r26804, r26821, MPFR_RNDN); mpfr_add(r26823, r26808, r26806, MPFR_RNDN); mpfr_sqrt(r26824, r26823, MPFR_RNDN); mpfr_div(r26825, r26822, r26824, MPFR_RNDN); if (mpfr_get_si(r26813, MPFR_RNDN)) { mpfr_set(r26826, r26819, MPFR_RNDN); } else { mpfr_set(r26826, r26825, MPFR_RNDN); }; if (mpfr_get_si(r26803, MPFR_RNDN)) { mpfr_set(r26827, r26811, MPFR_RNDN); } else { mpfr_set(r26827, r26826, MPFR_RNDN); }; return mpfr_get_d(r26827, MPFR_RNDN); }