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

char *name = "_multiplyComplex, imaginary part";

double f_if(float x_re, float x_im, float y_re, float y_im) {
        float r9145 = x_re;
        float r9146 = y_im;
        float r9147 = r9145 * r9146;
        float r9148 = x_im;
        float r9149 = y_re;
        float r9150 = r9148 * r9149;
        float r9151 = r9147 + r9150;
        return r9151;
}

double f_id(double x_re, double x_im, double y_re, double y_im) {
        double r9152 = x_re;
        double r9153 = y_im;
        double r9154 = r9152 * r9153;
        double r9155 = x_im;
        double r9156 = y_re;
        double r9157 = r9155 * r9156;
        double r9158 = r9154 + r9157;
        return r9158;
}


double f_of(float x_re, float x_im, float y_re, float y_im) {
        float r9159 = y_re;
        float r9160 = x_im;
        float r9161 = y_im;
        float r9162 = x_re;
        float r9163 = r9161 * r9162;
        float r9164 = fma(r9159, r9160, r9163);
        return r9164;
}

double f_od(double x_re, double x_im, double y_re, double y_im) {
        double r9165 = y_re;
        double r9166 = x_im;
        double r9167 = y_im;
        double r9168 = x_re;
        double r9169 = r9167 * r9168;
        double r9170 = fma(r9165, r9166, r9169);
        return r9170;
}

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 r9171, r9172, r9173, r9174, r9175, r9176, r9177;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9171);
        mpfr_init(r9172);
        mpfr_init(r9173);
        mpfr_init(r9174);
        mpfr_init(r9175);
        mpfr_init(r9176);
        mpfr_init(r9177);
}

double f_im(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r9171, x_re, MPFR_RNDN);
        mpfr_set_d(r9172, y_im, MPFR_RNDN);
        mpfr_mul(r9173, r9171, r9172, MPFR_RNDN);
        mpfr_set_d(r9174, x_im, MPFR_RNDN);
        mpfr_set_d(r9175, y_re, MPFR_RNDN);
        mpfr_mul(r9176, r9174, r9175, MPFR_RNDN);
        mpfr_add(r9177, r9173, r9176, MPFR_RNDN);
        return mpfr_get_d(r9177, MPFR_RNDN);
}

static mpfr_t r9178, r9179, r9180, r9181, r9182, r9183;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9178);
        mpfr_init(r9179);
        mpfr_init(r9180);
        mpfr_init(r9181);
        mpfr_init(r9182);
        mpfr_init(r9183);
}

double f_fm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r9178, y_re, MPFR_RNDN);
        mpfr_set_d(r9179, x_im, MPFR_RNDN);
        mpfr_set_d(r9180, y_im, MPFR_RNDN);
        mpfr_set_d(r9181, x_re, MPFR_RNDN);
        mpfr_mul(r9182, r9180, r9181, MPFR_RNDN);
        mpfr_fma(r9183, r9178, r9179, r9182, MPFR_RNDN);
        return mpfr_get_d(r9183, MPFR_RNDN);
}

static mpfr_t r9184, r9185, r9186, r9187, r9188, r9189;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9184);
        mpfr_init(r9185);
        mpfr_init(r9186);
        mpfr_init(r9187);
        mpfr_init(r9188);
        mpfr_init(r9189);
}

double f_dm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r9184, y_re, MPFR_RNDN);
        mpfr_set_d(r9185, x_im, MPFR_RNDN);
        mpfr_set_d(r9186, y_im, MPFR_RNDN);
        mpfr_set_d(r9187, x_re, MPFR_RNDN);
        mpfr_mul(r9188, r9186, r9187, MPFR_RNDN);
        mpfr_fma(r9189, r9184, r9185, r9188, MPFR_RNDN);
        return mpfr_get_d(r9189, MPFR_RNDN);
}