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

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

double f_if(float a, float b, float c, float d) {
        float r26602 = a;
        float r26603 = c;
        float r26604 = r26602 * r26603;
        float r26605 = b;
        float r26606 = d;
        float r26607 = r26605 * r26606;
        float r26608 = r26604 + r26607;
        float r26609 = r26603 * r26603;
        float r26610 = r26606 * r26606;
        float r26611 = r26609 + r26610;
        float r26612 = r26608 / r26611;
        return r26612;
}

double f_id(double a, double b, double c, double d) {
        double r26613 = a;
        double r26614 = c;
        double r26615 = r26613 * r26614;
        double r26616 = b;
        double r26617 = d;
        double r26618 = r26616 * r26617;
        double r26619 = r26615 + r26618;
        double r26620 = r26614 * r26614;
        double r26621 = r26617 * r26617;
        double r26622 = r26620 + r26621;
        double r26623 = r26619 / r26622;
        return r26623;
}


double f_of(float a, float b, float c, float d) {
        float r26624 = a;
        float r26625 = c;
        float r26626 = r26624 * r26625;
        float r26627 = b;
        float r26628 = d;
        float r26629 = r26627 * r26628;
        float r26630 = r26626 + r26629;
        float r26631 = r26625 * r26625;
        float r26632 = r26628 * r26628;
        float r26633 = r26631 + r26632;
        float r26634 = sqrt(r26633);
        float r26635 = r26630 / r26634;
        float r26636 = r26635 / r26634;
        return r26636;
}

double f_od(double a, double b, double c, double d) {
        double r26637 = a;
        double r26638 = c;
        double r26639 = r26637 * r26638;
        double r26640 = b;
        double r26641 = d;
        double r26642 = r26640 * r26641;
        double r26643 = r26639 + r26642;
        double r26644 = r26638 * r26638;
        double r26645 = r26641 * r26641;
        double r26646 = r26644 + r26645;
        double r26647 = sqrt(r26646);
        double r26648 = r26643 / r26647;
        double r26649 = r26648 / r26647;
        return r26649;
}

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 r26650, r26651, r26652, r26653, r26654, r26655, r26656, r26657, r26658, r26659, r26660;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26650);
        mpfr_init(r26651);
        mpfr_init(r26652);
        mpfr_init(r26653);
        mpfr_init(r26654);
        mpfr_init(r26655);
        mpfr_init(r26656);
        mpfr_init(r26657);
        mpfr_init(r26658);
        mpfr_init(r26659);
        mpfr_init(r26660);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r26650, a, MPFR_RNDN);
        mpfr_set_d(r26651, c, MPFR_RNDN);
        mpfr_mul(r26652, r26650, r26651, MPFR_RNDN);
        mpfr_set_d(r26653, b, MPFR_RNDN);
        mpfr_set_d(r26654, d, MPFR_RNDN);
        mpfr_mul(r26655, r26653, r26654, MPFR_RNDN);
        mpfr_add(r26656, r26652, r26655, MPFR_RNDN);
        mpfr_mul(r26657, r26651, r26651, MPFR_RNDN);
        mpfr_mul(r26658, r26654, r26654, MPFR_RNDN);
        mpfr_add(r26659, r26657, r26658, MPFR_RNDN);
        mpfr_div(r26660, r26656, r26659, MPFR_RNDN);
        return mpfr_get_d(r26660, MPFR_RNDN);
}

static mpfr_t r26661, r26662, r26663, r26664, r26665, r26666, r26667, r26668, r26669, r26670, r26671, r26672, r26673;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26661);
        mpfr_init(r26662);
        mpfr_init(r26663);
        mpfr_init(r26664);
        mpfr_init(r26665);
        mpfr_init(r26666);
        mpfr_init(r26667);
        mpfr_init(r26668);
        mpfr_init(r26669);
        mpfr_init(r26670);
        mpfr_init(r26671);
        mpfr_init(r26672);
        mpfr_init(r26673);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r26661, a, MPFR_RNDN);
        mpfr_set_d(r26662, c, MPFR_RNDN);
        mpfr_mul(r26663, r26661, r26662, MPFR_RNDN);
        mpfr_set_d(r26664, b, MPFR_RNDN);
        mpfr_set_d(r26665, d, MPFR_RNDN);
        mpfr_mul(r26666, r26664, r26665, MPFR_RNDN);
        mpfr_add(r26667, r26663, r26666, MPFR_RNDN);
        mpfr_mul(r26668, r26662, r26662, MPFR_RNDN);
        mpfr_mul(r26669, r26665, r26665, MPFR_RNDN);
        mpfr_add(r26670, r26668, r26669, MPFR_RNDN);
        mpfr_sqrt(r26671, r26670, MPFR_RNDN);
        mpfr_div(r26672, r26667, r26671, MPFR_RNDN);
        mpfr_div(r26673, r26672, r26671, MPFR_RNDN);
        return mpfr_get_d(r26673, MPFR_RNDN);
}

static mpfr_t r26674, r26675, r26676, r26677, r26678, r26679, r26680, r26681, r26682, r26683, r26684, r26685, r26686;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26674);
        mpfr_init(r26675);
        mpfr_init(r26676);
        mpfr_init(r26677);
        mpfr_init(r26678);
        mpfr_init(r26679);
        mpfr_init(r26680);
        mpfr_init(r26681);
        mpfr_init(r26682);
        mpfr_init(r26683);
        mpfr_init(r26684);
        mpfr_init(r26685);
        mpfr_init(r26686);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r26674, a, MPFR_RNDN);
        mpfr_set_d(r26675, c, MPFR_RNDN);
        mpfr_mul(r26676, r26674, r26675, MPFR_RNDN);
        mpfr_set_d(r26677, b, MPFR_RNDN);
        mpfr_set_d(r26678, d, MPFR_RNDN);
        mpfr_mul(r26679, r26677, r26678, MPFR_RNDN);
        mpfr_add(r26680, r26676, r26679, MPFR_RNDN);
        mpfr_mul(r26681, r26675, r26675, MPFR_RNDN);
        mpfr_mul(r26682, r26678, r26678, MPFR_RNDN);
        mpfr_add(r26683, r26681, r26682, MPFR_RNDN);
        mpfr_sqrt(r26684, r26683, MPFR_RNDN);
        mpfr_div(r26685, r26680, r26684, MPFR_RNDN);
        mpfr_div(r26686, r26685, r26684, MPFR_RNDN);
        return mpfr_get_d(r26686, MPFR_RNDN);
}