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

char *name = "Complex division, imag part";

double f_if(float a, float b, float c, float d) {
        float r26862 = b;
        float r26863 = c;
        float r26864 = r26862 * r26863;
        float r26865 = a;
        float r26866 = d;
        float r26867 = r26865 * r26866;
        float r26868 = r26864 - r26867;
        float r26869 = r26863 * r26863;
        float r26870 = r26866 * r26866;
        float r26871 = r26869 + r26870;
        float r26872 = r26868 / r26871;
        return r26872;
}

double f_id(double a, double b, double c, double d) {
        double r26873 = b;
        double r26874 = c;
        double r26875 = r26873 * r26874;
        double r26876 = a;
        double r26877 = d;
        double r26878 = r26876 * r26877;
        double r26879 = r26875 - r26878;
        double r26880 = r26874 * r26874;
        double r26881 = r26877 * r26877;
        double r26882 = r26880 + r26881;
        double r26883 = r26879 / r26882;
        return r26883;
}


double f_of(float a, float b, float c, float d) {
        float r26884 = c;
        float r26885 = 6.5435474222459e+39;
        bool r26886 = r26884 <= r26885;
        float r26887 = b;
        float r26888 = r26887 * r26884;
        float r26889 = a;
        float r26890 = d;
        float r26891 = r26889 * r26890;
        float r26892 = r26888 - r26891;
        float r26893 = r26884 * r26884;
        float r26894 = r26890 * r26890;
        float r26895 = r26893 + r26894;
        float r26896 = sqrt(r26895);
        float r26897 = r26892 / r26896;
        float r26898 = r26897 / r26896;
        float r26899 = r26884 / r26889;
        float r26900 = r26890 / r26899;
        float r26901 = r26887 - r26900;
        float r26902 = r26894 + r26893;
        float r26903 = sqrt(r26902);
        float r26904 = r26901 / r26903;
        float r26905 = r26886 ? r26898 : r26904;
        return r26905;
}

double f_od(double a, double b, double c, double d) {
        double r26906 = c;
        double r26907 = 6.5435474222459e+39;
        bool r26908 = r26906 <= r26907;
        double r26909 = b;
        double r26910 = r26909 * r26906;
        double r26911 = a;
        double r26912 = d;
        double r26913 = r26911 * r26912;
        double r26914 = r26910 - r26913;
        double r26915 = r26906 * r26906;
        double r26916 = r26912 * r26912;
        double r26917 = r26915 + r26916;
        double r26918 = sqrt(r26917);
        double r26919 = r26914 / r26918;
        double r26920 = r26919 / r26918;
        double r26921 = r26906 / r26911;
        double r26922 = r26912 / r26921;
        double r26923 = r26909 - r26922;
        double r26924 = r26916 + r26915;
        double r26925 = sqrt(r26924);
        double r26926 = r26923 / r26925;
        double r26927 = r26908 ? r26920 : r26926;
        return r26927;
}

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 r26928, r26929, r26930, r26931, r26932, r26933, r26934, r26935, r26936, r26937, r26938;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26928);
        mpfr_init(r26929);
        mpfr_init(r26930);
        mpfr_init(r26931);
        mpfr_init(r26932);
        mpfr_init(r26933);
        mpfr_init(r26934);
        mpfr_init(r26935);
        mpfr_init(r26936);
        mpfr_init(r26937);
        mpfr_init(r26938);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r26928, b, MPFR_RNDN);
        mpfr_set_d(r26929, c, MPFR_RNDN);
        mpfr_mul(r26930, r26928, r26929, MPFR_RNDN);
        mpfr_set_d(r26931, a, MPFR_RNDN);
        mpfr_set_d(r26932, d, MPFR_RNDN);
        mpfr_mul(r26933, r26931, r26932, MPFR_RNDN);
        mpfr_sub(r26934, r26930, r26933, MPFR_RNDN);
        mpfr_mul(r26935, r26929, r26929, MPFR_RNDN);
        mpfr_mul(r26936, r26932, r26932, MPFR_RNDN);
        mpfr_add(r26937, r26935, r26936, MPFR_RNDN);
        mpfr_div(r26938, r26934, r26937, MPFR_RNDN);
        return mpfr_get_d(r26938, MPFR_RNDN);
}

static mpfr_t r26939, r26940, r26941, r26942, r26943, r26944, r26945, r26946, r26947, r26948, r26949, r26950, r26951, r26952, r26953, r26954, r26955, r26956, r26957, r26958, r26959, r26960;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26939);
        mpfr_init_set_str(r26940, "6.5435474222459e+39", 10, MPFR_RNDN);
        mpfr_init(r26941);
        mpfr_init(r26942);
        mpfr_init(r26943);
        mpfr_init(r26944);
        mpfr_init(r26945);
        mpfr_init(r26946);
        mpfr_init(r26947);
        mpfr_init(r26948);
        mpfr_init(r26949);
        mpfr_init(r26950);
        mpfr_init(r26951);
        mpfr_init(r26952);
        mpfr_init(r26953);
        mpfr_init(r26954);
        mpfr_init(r26955);
        mpfr_init(r26956);
        mpfr_init(r26957);
        mpfr_init(r26958);
        mpfr_init(r26959);
        mpfr_init(r26960);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r26939, c, MPFR_RNDN);
        ;
        mpfr_set_si(r26941, mpfr_cmp(r26939, r26940) <= 0, MPFR_RNDN);
        mpfr_set_d(r26942, b, MPFR_RNDN);
        mpfr_mul(r26943, r26942, r26939, MPFR_RNDN);
        mpfr_set_d(r26944, a, MPFR_RNDN);
        mpfr_set_d(r26945, d, MPFR_RNDN);
        mpfr_mul(r26946, r26944, r26945, MPFR_RNDN);
        mpfr_sub(r26947, r26943, r26946, MPFR_RNDN);
        mpfr_mul(r26948, r26939, r26939, MPFR_RNDN);
        mpfr_mul(r26949, r26945, r26945, MPFR_RNDN);
        mpfr_add(r26950, r26948, r26949, MPFR_RNDN);
        mpfr_sqrt(r26951, r26950, MPFR_RNDN);
        mpfr_div(r26952, r26947, r26951, MPFR_RNDN);
        mpfr_div(r26953, r26952, r26951, MPFR_RNDN);
        mpfr_div(r26954, r26939, r26944, MPFR_RNDN);
        mpfr_div(r26955, r26945, r26954, MPFR_RNDN);
        mpfr_sub(r26956, r26942, r26955, MPFR_RNDN);
        mpfr_add(r26957, r26949, r26948, MPFR_RNDN);
        mpfr_sqrt(r26958, r26957, MPFR_RNDN);
        mpfr_div(r26959, r26956, r26958, MPFR_RNDN);
        if (mpfr_get_si(r26941, MPFR_RNDN)) { mpfr_set(r26960, r26953, MPFR_RNDN); } else { mpfr_set(r26960, r26959, MPFR_RNDN); };
        return mpfr_get_d(r26960, MPFR_RNDN);
}

static mpfr_t r26961, r26962, r26963, r26964, r26965, r26966, r26967, r26968, r26969, r26970, r26971, r26972, r26973, r26974, r26975, r26976, r26977, r26978, r26979, r26980, r26981, r26982;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26961);
        mpfr_init_set_str(r26962, "6.5435474222459e+39", 10, MPFR_RNDN);
        mpfr_init(r26963);
        mpfr_init(r26964);
        mpfr_init(r26965);
        mpfr_init(r26966);
        mpfr_init(r26967);
        mpfr_init(r26968);
        mpfr_init(r26969);
        mpfr_init(r26970);
        mpfr_init(r26971);
        mpfr_init(r26972);
        mpfr_init(r26973);
        mpfr_init(r26974);
        mpfr_init(r26975);
        mpfr_init(r26976);
        mpfr_init(r26977);
        mpfr_init(r26978);
        mpfr_init(r26979);
        mpfr_init(r26980);
        mpfr_init(r26981);
        mpfr_init(r26982);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r26961, c, MPFR_RNDN);
        ;
        mpfr_set_si(r26963, mpfr_cmp(r26961, r26962) <= 0, MPFR_RNDN);
        mpfr_set_d(r26964, b, MPFR_RNDN);
        mpfr_mul(r26965, r26964, r26961, MPFR_RNDN);
        mpfr_set_d(r26966, a, MPFR_RNDN);
        mpfr_set_d(r26967, d, MPFR_RNDN);
        mpfr_mul(r26968, r26966, r26967, MPFR_RNDN);
        mpfr_sub(r26969, r26965, r26968, MPFR_RNDN);
        mpfr_mul(r26970, r26961, r26961, MPFR_RNDN);
        mpfr_mul(r26971, r26967, r26967, MPFR_RNDN);
        mpfr_add(r26972, r26970, r26971, MPFR_RNDN);
        mpfr_sqrt(r26973, r26972, MPFR_RNDN);
        mpfr_div(r26974, r26969, r26973, MPFR_RNDN);
        mpfr_div(r26975, r26974, r26973, MPFR_RNDN);
        mpfr_div(r26976, r26961, r26966, MPFR_RNDN);
        mpfr_div(r26977, r26967, r26976, MPFR_RNDN);
        mpfr_sub(r26978, r26964, r26977, MPFR_RNDN);
        mpfr_add(r26979, r26971, r26970, MPFR_RNDN);
        mpfr_sqrt(r26980, r26979, MPFR_RNDN);
        mpfr_div(r26981, r26978, r26980, MPFR_RNDN);
        if (mpfr_get_si(r26963, MPFR_RNDN)) { mpfr_set(r26982, r26975, MPFR_RNDN); } else { mpfr_set(r26982, r26981, MPFR_RNDN); };
        return mpfr_get_d(r26982, MPFR_RNDN);
}