#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 = r26630 / r26633;
        return r26634;
}

double f_od(double a, double b, double c, double d) {
        double r26635 = a;
        double r26636 = c;
        double r26637 = r26635 * r26636;
        double r26638 = b;
        double r26639 = d;
        double r26640 = r26638 * r26639;
        double r26641 = r26637 + r26640;
        double r26642 = r26636 * r26636;
        double r26643 = r26639 * r26639;
        double r26644 = r26642 + r26643;
        double r26645 = r26641 / r26644;
        return r26645;
}

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 r26646, r26647, r26648, r26649, r26650, r26651, r26652, r26653, r26654, r26655, r26656;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r26646);
        mpfr_init(r26647);
        mpfr_init(r26648);
        mpfr_init(r26649);
        mpfr_init(r26650);
        mpfr_init(r26651);
        mpfr_init(r26652);
        mpfr_init(r26653);
        mpfr_init(r26654);
        mpfr_init(r26655);
        mpfr_init(r26656);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r26646, a, MPFR_RNDN);
        mpfr_set_d(r26647, c, MPFR_RNDN);
        mpfr_mul(r26648, r26646, r26647, MPFR_RNDN);
        mpfr_set_d(r26649, b, MPFR_RNDN);
        mpfr_set_d(r26650, d, MPFR_RNDN);
        mpfr_mul(r26651, r26649, r26650, MPFR_RNDN);
        mpfr_add(r26652, r26648, r26651, MPFR_RNDN);
        mpfr_mul(r26653, r26647, r26647, MPFR_RNDN);
        mpfr_mul(r26654, r26650, r26650, MPFR_RNDN);
        mpfr_add(r26655, r26653, r26654, MPFR_RNDN);
        mpfr_div(r26656, r26652, r26655, MPFR_RNDN);
        return mpfr_get_d(r26656, MPFR_RNDN);
}

static mpfr_t r26657, r26658, r26659, r26660, r26661, r26662, r26663, r26664, r26665, r26666, r26667;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r26657);
        mpfr_init(r26658);
        mpfr_init(r26659);
        mpfr_init(r26660);
        mpfr_init(r26661);
        mpfr_init(r26662);
        mpfr_init(r26663);
        mpfr_init(r26664);
        mpfr_init(r26665);
        mpfr_init(r26666);
        mpfr_init(r26667);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r26657, a, MPFR_RNDN);
        mpfr_set_d(r26658, c, MPFR_RNDN);
        mpfr_mul(r26659, r26657, r26658, MPFR_RNDN);
        mpfr_set_d(r26660, b, MPFR_RNDN);
        mpfr_set_d(r26661, d, MPFR_RNDN);
        mpfr_mul(r26662, r26660, r26661, MPFR_RNDN);
        mpfr_add(r26663, r26659, r26662, MPFR_RNDN);
        mpfr_mul(r26664, r26658, r26658, MPFR_RNDN);
        mpfr_mul(r26665, r26661, r26661, MPFR_RNDN);
        mpfr_add(r26666, r26664, r26665, MPFR_RNDN);
        mpfr_div(r26667, r26663, r26666, MPFR_RNDN);
        return mpfr_get_d(r26667, MPFR_RNDN);
}

static mpfr_t r26668, r26669, r26670, r26671, r26672, r26673, r26674, r26675, r26676, r26677, r26678;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r26668);
        mpfr_init(r26669);
        mpfr_init(r26670);
        mpfr_init(r26671);
        mpfr_init(r26672);
        mpfr_init(r26673);
        mpfr_init(r26674);
        mpfr_init(r26675);
        mpfr_init(r26676);
        mpfr_init(r26677);
        mpfr_init(r26678);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r26668, a, MPFR_RNDN);
        mpfr_set_d(r26669, c, MPFR_RNDN);
        mpfr_mul(r26670, r26668, r26669, MPFR_RNDN);
        mpfr_set_d(r26671, b, MPFR_RNDN);
        mpfr_set_d(r26672, d, MPFR_RNDN);
        mpfr_mul(r26673, r26671, r26672, MPFR_RNDN);
        mpfr_add(r26674, r26670, r26673, MPFR_RNDN);
        mpfr_mul(r26675, r26669, r26669, MPFR_RNDN);
        mpfr_mul(r26676, r26672, r26672, MPFR_RNDN);
        mpfr_add(r26677, r26675, r26676, MPFR_RNDN);
        mpfr_div(r26678, r26674, r26677, MPFR_RNDN);
        return mpfr_get_d(r26678, MPFR_RNDN);
}