#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "_divideComplex, real part";

double f_if(float x_re, float x_im, float y_re, float y_im) {
        float r9869 = x_re;
        float r9870 = y_re;
        float r9871 = r9869 * r9870;
        float r9872 = x_im;
        float r9873 = y_im;
        float r9874 = r9872 * r9873;
        float r9875 = r9871 + r9874;
        float r9876 = r9870 * r9870;
        float r9877 = r9873 * r9873;
        float r9878 = r9876 + r9877;
        float r9879 = r9875 / r9878;
        return r9879;
}

double f_id(double x_re, double x_im, double y_re, double y_im) {
        double r9880 = x_re;
        double r9881 = y_re;
        double r9882 = r9880 * r9881;
        double r9883 = x_im;
        double r9884 = y_im;
        double r9885 = r9883 * r9884;
        double r9886 = r9882 + r9885;
        double r9887 = r9881 * r9881;
        double r9888 = r9884 * r9884;
        double r9889 = r9887 + r9888;
        double r9890 = r9886 / r9889;
        return r9890;
}


double f_of(float x_re, float x_im, float y_re, float y_im) {
        float r9891 = y_im;
        float r9892 = 2.287685333116204e+71;
        bool r9893 = r9891 <= r9892;
        float r9894 = x_re;
        float r9895 = y_re;
        float r9896 = r9894 * r9895;
        float r9897 = x_im;
        float r9898 = r9897 * r9891;
        float r9899 = r9896 + r9898;
        float r9900 = r9895 * r9895;
        float r9901 = r9891 * r9891;
        float r9902 = r9900 + r9901;
        float r9903 = sqrt(r9902);
        float r9904 = r9899 / r9903;
        float r9905 = r9904 / r9903;
        float r9906 = 2.6646868450529295e+148;
        bool r9907 = r9891 <= r9906;
        float r9908 = r9897 / r9903;
        float r9909 = cbrt(r9899);
        float r9910 = r9909 * r9909;
        float r9911 = cbrt(r9902);
        float r9912 = r9911 * r9911;
        float r9913 = r9910 / r9912;
        float r9914 = r9909 / r9911;
        float r9915 = r9913 * r9914;
        float r9916 = r9907 ? r9908 : r9915;
        float r9917 = r9893 ? r9905 : r9916;
        return r9917;
}

double f_od(double x_re, double x_im, double y_re, double y_im) {
        double r9918 = y_im;
        double r9919 = 2.287685333116204e+71;
        bool r9920 = r9918 <= r9919;
        double r9921 = x_re;
        double r9922 = y_re;
        double r9923 = r9921 * r9922;
        double r9924 = x_im;
        double r9925 = r9924 * r9918;
        double r9926 = r9923 + r9925;
        double r9927 = r9922 * r9922;
        double r9928 = r9918 * r9918;
        double r9929 = r9927 + r9928;
        double r9930 = sqrt(r9929);
        double r9931 = r9926 / r9930;
        double r9932 = r9931 / r9930;
        double r9933 = 2.6646868450529295e+148;
        bool r9934 = r9918 <= r9933;
        double r9935 = r9924 / r9930;
        double r9936 = cbrt(r9926);
        double r9937 = r9936 * r9936;
        double r9938 = cbrt(r9929);
        double r9939 = r9938 * r9938;
        double r9940 = r9937 / r9939;
        double r9941 = r9936 / r9938;
        double r9942 = r9940 * r9941;
        double r9943 = r9934 ? r9935 : r9942;
        double r9944 = r9920 ? r9932 : r9943;
        return r9944;
}

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 r9945, r9946, r9947, r9948, r9949, r9950, r9951, r9952, r9953, r9954, r9955;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r9945);
        mpfr_init(r9946);
        mpfr_init(r9947);
        mpfr_init(r9948);
        mpfr_init(r9949);
        mpfr_init(r9950);
        mpfr_init(r9951);
        mpfr_init(r9952);
        mpfr_init(r9953);
        mpfr_init(r9954);
        mpfr_init(r9955);
}

double f_im(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r9945, x_re, MPFR_RNDN);
        mpfr_set_d(r9946, y_re, MPFR_RNDN);
        mpfr_mul(r9947, r9945, r9946, MPFR_RNDN);
        mpfr_set_d(r9948, x_im, MPFR_RNDN);
        mpfr_set_d(r9949, y_im, MPFR_RNDN);
        mpfr_mul(r9950, r9948, r9949, MPFR_RNDN);
        mpfr_add(r9951, r9947, r9950, MPFR_RNDN);
        mpfr_mul(r9952, r9946, r9946, MPFR_RNDN);
        mpfr_mul(r9953, r9949, r9949, MPFR_RNDN);
        mpfr_add(r9954, r9952, r9953, MPFR_RNDN);
        mpfr_div(r9955, r9951, r9954, MPFR_RNDN);
        return mpfr_get_d(r9955, MPFR_RNDN);
}

static mpfr_t r9956, r9957, r9958, r9959, r9960, r9961, r9962, r9963, r9964, r9965, r9966, r9967, r9968, r9969, r9970, r9971, r9972, r9973, r9974, r9975, r9976, r9977, r9978, r9979, r9980, r9981, r9982;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r9956);
        mpfr_init_set_str(r9957, "2.287685333116204e+71", 10, MPFR_RNDN);
        mpfr_init(r9958);
        mpfr_init(r9959);
        mpfr_init(r9960);
        mpfr_init(r9961);
        mpfr_init(r9962);
        mpfr_init(r9963);
        mpfr_init(r9964);
        mpfr_init(r9965);
        mpfr_init(r9966);
        mpfr_init(r9967);
        mpfr_init(r9968);
        mpfr_init(r9969);
        mpfr_init(r9970);
        mpfr_init_set_str(r9971, "2.6646868450529295e+148", 10, MPFR_RNDN);
        mpfr_init(r9972);
        mpfr_init(r9973);
        mpfr_init(r9974);
        mpfr_init(r9975);
        mpfr_init(r9976);
        mpfr_init(r9977);
        mpfr_init(r9978);
        mpfr_init(r9979);
        mpfr_init(r9980);
        mpfr_init(r9981);
        mpfr_init(r9982);
}

double f_fm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r9956, y_im, MPFR_RNDN);
        ;
        mpfr_set_si(r9958, mpfr_cmp(r9956, r9957) <= 0, MPFR_RNDN);
        mpfr_set_d(r9959, x_re, MPFR_RNDN);
        mpfr_set_d(r9960, y_re, MPFR_RNDN);
        mpfr_mul(r9961, r9959, r9960, MPFR_RNDN);
        mpfr_set_d(r9962, x_im, MPFR_RNDN);
        mpfr_mul(r9963, r9962, r9956, MPFR_RNDN);
        mpfr_add(r9964, r9961, r9963, MPFR_RNDN);
        mpfr_mul(r9965, r9960, r9960, MPFR_RNDN);
        mpfr_mul(r9966, r9956, r9956, MPFR_RNDN);
        mpfr_add(r9967, r9965, r9966, MPFR_RNDN);
        mpfr_sqrt(r9968, r9967, MPFR_RNDN);
        mpfr_div(r9969, r9964, r9968, MPFR_RNDN);
        mpfr_div(r9970, r9969, r9968, MPFR_RNDN);
        ;
        mpfr_set_si(r9972, mpfr_cmp(r9956, r9971) <= 0, MPFR_RNDN);
        mpfr_div(r9973, r9962, r9968, MPFR_RNDN);
        mpfr_cbrt(r9974, r9964, MPFR_RNDN);
        mpfr_mul(r9975, r9974, r9974, MPFR_RNDN);
        mpfr_cbrt(r9976, r9967, MPFR_RNDN);
        mpfr_mul(r9977, r9976, r9976, MPFR_RNDN);
        mpfr_div(r9978, r9975, r9977, MPFR_RNDN);
        mpfr_div(r9979, r9974, r9976, MPFR_RNDN);
        mpfr_mul(r9980, r9978, r9979, MPFR_RNDN);
        if (mpfr_get_si(r9972, MPFR_RNDN)) { mpfr_set(r9981, r9973, MPFR_RNDN); } else { mpfr_set(r9981, r9980, MPFR_RNDN); };
        if (mpfr_get_si(r9958, MPFR_RNDN)) { mpfr_set(r9982, r9970, MPFR_RNDN); } else { mpfr_set(r9982, r9981, MPFR_RNDN); };
        return mpfr_get_d(r9982, MPFR_RNDN);
}

static mpfr_t r9983, r9984, r9985, r9986, r9987, r9988, r9989, r9990, r9991, r9992, r9993, r9994, r9995, r9996, r9997, r9998, r9999, r10000, r10001, r10002, r10003, r10004, r10005, r10006, r10007, r10008, r10009;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r9983);
        mpfr_init_set_str(r9984, "2.287685333116204e+71", 10, MPFR_RNDN);
        mpfr_init(r9985);
        mpfr_init(r9986);
        mpfr_init(r9987);
        mpfr_init(r9988);
        mpfr_init(r9989);
        mpfr_init(r9990);
        mpfr_init(r9991);
        mpfr_init(r9992);
        mpfr_init(r9993);
        mpfr_init(r9994);
        mpfr_init(r9995);
        mpfr_init(r9996);
        mpfr_init(r9997);
        mpfr_init_set_str(r9998, "2.6646868450529295e+148", 10, MPFR_RNDN);
        mpfr_init(r9999);
        mpfr_init(r10000);
        mpfr_init(r10001);
        mpfr_init(r10002);
        mpfr_init(r10003);
        mpfr_init(r10004);
        mpfr_init(r10005);
        mpfr_init(r10006);
        mpfr_init(r10007);
        mpfr_init(r10008);
        mpfr_init(r10009);
}

double f_dm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r9983, y_im, MPFR_RNDN);
        ;
        mpfr_set_si(r9985, mpfr_cmp(r9983, r9984) <= 0, MPFR_RNDN);
        mpfr_set_d(r9986, x_re, MPFR_RNDN);
        mpfr_set_d(r9987, y_re, MPFR_RNDN);
        mpfr_mul(r9988, r9986, r9987, MPFR_RNDN);
        mpfr_set_d(r9989, x_im, MPFR_RNDN);
        mpfr_mul(r9990, r9989, r9983, MPFR_RNDN);
        mpfr_add(r9991, r9988, r9990, MPFR_RNDN);
        mpfr_mul(r9992, r9987, r9987, MPFR_RNDN);
        mpfr_mul(r9993, r9983, r9983, MPFR_RNDN);
        mpfr_add(r9994, r9992, r9993, MPFR_RNDN);
        mpfr_sqrt(r9995, r9994, MPFR_RNDN);
        mpfr_div(r9996, r9991, r9995, MPFR_RNDN);
        mpfr_div(r9997, r9996, r9995, MPFR_RNDN);
        ;
        mpfr_set_si(r9999, mpfr_cmp(r9983, r9998) <= 0, MPFR_RNDN);
        mpfr_div(r10000, r9989, r9995, MPFR_RNDN);
        mpfr_cbrt(r10001, r9991, MPFR_RNDN);
        mpfr_mul(r10002, r10001, r10001, MPFR_RNDN);
        mpfr_cbrt(r10003, r9994, MPFR_RNDN);
        mpfr_mul(r10004, r10003, r10003, MPFR_RNDN);
        mpfr_div(r10005, r10002, r10004, MPFR_RNDN);
        mpfr_div(r10006, r10001, r10003, MPFR_RNDN);
        mpfr_mul(r10007, r10005, r10006, MPFR_RNDN);
        if (mpfr_get_si(r9999, MPFR_RNDN)) { mpfr_set(r10008, r10000, MPFR_RNDN); } else { mpfr_set(r10008, r10007, MPFR_RNDN); };
        if (mpfr_get_si(r9985, MPFR_RNDN)) { mpfr_set(r10009, r9997, MPFR_RNDN); } else { mpfr_set(r10009, r10008, MPFR_RNDN); };
        return mpfr_get_d(r10009, MPFR_RNDN);
}