#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 r56916 = b;
        float r56917 = c;
        float r56918 = r56916 * r56917;
        float r56919 = a;
        float r56920 = d;
        float r56921 = r56919 * r56920;
        float r56922 = r56918 - r56921;
        float r56923 = r56917 * r56917;
        float r56924 = r56920 * r56920;
        float r56925 = r56923 + r56924;
        float r56926 = r56922 / r56925;
        return r56926;
}

double f_id(double a, double b, double c, double d) {
        double r56927 = b;
        double r56928 = c;
        double r56929 = r56927 * r56928;
        double r56930 = a;
        double r56931 = d;
        double r56932 = r56930 * r56931;
        double r56933 = r56929 - r56932;
        double r56934 = r56928 * r56928;
        double r56935 = r56931 * r56931;
        double r56936 = r56934 + r56935;
        double r56937 = r56933 / r56936;
        return r56937;
}


double f_of(float a, float b, float c, float d) {
        float r56938 = b;
        float r56939 = c;
        float r56940 = r56938 * r56939;
        float r56941 = r56939 * r56939;
        float r56942 = d;
        float r56943 = r56942 * r56942;
        float r56944 = r56941 + r56943;
        float r56945 = r56940 / r56944;
        float r56946 = a;
        float r56947 = sqrt(r56944);
        float r56948 = r56946 / r56947;
        float r56949 = r56942 / r56947;
        float r56950 = r56948 * r56949;
        float r56951 = r56945 - r56950;
        return r56951;
}

double f_od(double a, double b, double c, double d) {
        double r56952 = b;
        double r56953 = c;
        double r56954 = r56952 * r56953;
        double r56955 = r56953 * r56953;
        double r56956 = d;
        double r56957 = r56956 * r56956;
        double r56958 = r56955 + r56957;
        double r56959 = r56954 / r56958;
        double r56960 = a;
        double r56961 = sqrt(r56958);
        double r56962 = r56960 / r56961;
        double r56963 = r56956 / r56961;
        double r56964 = r56962 * r56963;
        double r56965 = r56959 - r56964;
        return r56965;
}

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 r56966, r56967, r56968, r56969, r56970, r56971, r56972, r56973, r56974, r56975, r56976;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r56966);
        mpfr_init(r56967);
        mpfr_init(r56968);
        mpfr_init(r56969);
        mpfr_init(r56970);
        mpfr_init(r56971);
        mpfr_init(r56972);
        mpfr_init(r56973);
        mpfr_init(r56974);
        mpfr_init(r56975);
        mpfr_init(r56976);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r56966, b, MPFR_RNDN);
        mpfr_set_d(r56967, c, MPFR_RNDN);
        mpfr_mul(r56968, r56966, r56967, MPFR_RNDN);
        mpfr_set_d(r56969, a, MPFR_RNDN);
        mpfr_set_d(r56970, d, MPFR_RNDN);
        mpfr_mul(r56971, r56969, r56970, MPFR_RNDN);
        mpfr_sub(r56972, r56968, r56971, MPFR_RNDN);
        mpfr_mul(r56973, r56967, r56967, MPFR_RNDN);
        mpfr_mul(r56974, r56970, r56970, MPFR_RNDN);
        mpfr_add(r56975, r56973, r56974, MPFR_RNDN);
        mpfr_div(r56976, r56972, r56975, MPFR_RNDN);
        return mpfr_get_d(r56976, MPFR_RNDN);
}

static mpfr_t r56977, r56978, r56979, r56980, r56981, r56982, r56983, r56984, r56985, r56986, r56987, r56988, r56989, r56990;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r56977);
        mpfr_init(r56978);
        mpfr_init(r56979);
        mpfr_init(r56980);
        mpfr_init(r56981);
        mpfr_init(r56982);
        mpfr_init(r56983);
        mpfr_init(r56984);
        mpfr_init(r56985);
        mpfr_init(r56986);
        mpfr_init(r56987);
        mpfr_init(r56988);
        mpfr_init(r56989);
        mpfr_init(r56990);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r56977, b, MPFR_RNDN);
        mpfr_set_d(r56978, c, MPFR_RNDN);
        mpfr_mul(r56979, r56977, r56978, MPFR_RNDN);
        mpfr_mul(r56980, r56978, r56978, MPFR_RNDN);
        mpfr_set_d(r56981, d, MPFR_RNDN);
        mpfr_mul(r56982, r56981, r56981, MPFR_RNDN);
        mpfr_add(r56983, r56980, r56982, MPFR_RNDN);
        mpfr_div(r56984, r56979, r56983, MPFR_RNDN);
        mpfr_set_d(r56985, a, MPFR_RNDN);
        mpfr_sqrt(r56986, r56983, MPFR_RNDN);
        mpfr_div(r56987, r56985, r56986, MPFR_RNDN);
        mpfr_div(r56988, r56981, r56986, MPFR_RNDN);
        mpfr_mul(r56989, r56987, r56988, MPFR_RNDN);
        mpfr_sub(r56990, r56984, r56989, MPFR_RNDN);
        return mpfr_get_d(r56990, MPFR_RNDN);
}

static mpfr_t r56991, r56992, r56993, r56994, r56995, r56996, r56997, r56998, r56999, r57000, r57001, r57002, r57003, r57004;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r56991);
        mpfr_init(r56992);
        mpfr_init(r56993);
        mpfr_init(r56994);
        mpfr_init(r56995);
        mpfr_init(r56996);
        mpfr_init(r56997);
        mpfr_init(r56998);
        mpfr_init(r56999);
        mpfr_init(r57000);
        mpfr_init(r57001);
        mpfr_init(r57002);
        mpfr_init(r57003);
        mpfr_init(r57004);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r56991, b, MPFR_RNDN);
        mpfr_set_d(r56992, c, MPFR_RNDN);
        mpfr_mul(r56993, r56991, r56992, MPFR_RNDN);
        mpfr_mul(r56994, r56992, r56992, MPFR_RNDN);
        mpfr_set_d(r56995, d, MPFR_RNDN);
        mpfr_mul(r56996, r56995, r56995, MPFR_RNDN);
        mpfr_add(r56997, r56994, r56996, MPFR_RNDN);
        mpfr_div(r56998, r56993, r56997, MPFR_RNDN);
        mpfr_set_d(r56999, a, MPFR_RNDN);
        mpfr_sqrt(r57000, r56997, MPFR_RNDN);
        mpfr_div(r57001, r56999, r57000, MPFR_RNDN);
        mpfr_div(r57002, r56995, r57000, MPFR_RNDN);
        mpfr_mul(r57003, r57001, r57002, MPFR_RNDN);
        mpfr_sub(r57004, r56998, r57003, MPFR_RNDN);
        return mpfr_get_d(r57004, MPFR_RNDN);
}