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

char *name = "_multiplyComplex, real part";

double f_if(float x_re, float x_im, float y_re, float y_im) {
        float r25470 = x_re;
        float r25471 = y_re;
        float r25472 = r25470 * r25471;
        float r25473 = x_im;
        float r25474 = y_im;
        float r25475 = r25473 * r25474;
        float r25476 = r25472 - r25475;
        return r25476;
}

double f_id(double x_re, double x_im, double y_re, double y_im) {
        double r25477 = x_re;
        double r25478 = y_re;
        double r25479 = r25477 * r25478;
        double r25480 = x_im;
        double r25481 = y_im;
        double r25482 = r25480 * r25481;
        double r25483 = r25479 - r25482;
        return r25483;
}


double f_of(float x_re, float x_im, float y_re, float y_im) {
        float r25484 = x_re;
        float r25485 = y_re;
        float r25486 = r25484 * r25485;
        float r25487 = x_im;
        float r25488 = y_im;
        float r25489 = r25487 * r25488;
        float r25490 = r25486 - r25489;
        return r25490;
}

double f_od(double x_re, double x_im, double y_re, double y_im) {
        double r25491 = x_re;
        double r25492 = y_re;
        double r25493 = r25491 * r25492;
        double r25494 = x_im;
        double r25495 = y_im;
        double r25496 = r25494 * r25495;
        double r25497 = r25493 - r25496;
        return r25497;
}

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 r25498, r25499, r25500, r25501, r25502, r25503, r25504;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r25498);
        mpfr_init(r25499);
        mpfr_init(r25500);
        mpfr_init(r25501);
        mpfr_init(r25502);
        mpfr_init(r25503);
        mpfr_init(r25504);
}

double f_im(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r25498, x_re, MPFR_RNDN);
        mpfr_set_d(r25499, y_re, MPFR_RNDN);
        mpfr_mul(r25500, r25498, r25499, MPFR_RNDN);
        mpfr_set_d(r25501, x_im, MPFR_RNDN);
        mpfr_set_d(r25502, y_im, MPFR_RNDN);
        mpfr_mul(r25503, r25501, r25502, MPFR_RNDN);
        mpfr_sub(r25504, r25500, r25503, MPFR_RNDN);
        return mpfr_get_d(r25504, MPFR_RNDN);
}

static mpfr_t r25505, r25506, r25507, r25508, r25509, r25510, r25511;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25505);
        mpfr_init(r25506);
        mpfr_init(r25507);
        mpfr_init(r25508);
        mpfr_init(r25509);
        mpfr_init(r25510);
        mpfr_init(r25511);
}

double f_fm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r25505, x_re, MPFR_RNDN);
        mpfr_set_d(r25506, y_re, MPFR_RNDN);
        mpfr_mul(r25507, r25505, r25506, MPFR_RNDN);
        mpfr_set_d(r25508, x_im, MPFR_RNDN);
        mpfr_set_d(r25509, y_im, MPFR_RNDN);
        mpfr_mul(r25510, r25508, r25509, MPFR_RNDN);
        mpfr_sub(r25511, r25507, r25510, MPFR_RNDN);
        return mpfr_get_d(r25511, MPFR_RNDN);
}

static mpfr_t r25512, r25513, r25514, r25515, r25516, r25517, r25518;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25512);
        mpfr_init(r25513);
        mpfr_init(r25514);
        mpfr_init(r25515);
        mpfr_init(r25516);
        mpfr_init(r25517);
        mpfr_init(r25518);
}

double f_dm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r25512, x_re, MPFR_RNDN);
        mpfr_set_d(r25513, y_re, MPFR_RNDN);
        mpfr_mul(r25514, r25512, r25513, MPFR_RNDN);
        mpfr_set_d(r25515, x_im, MPFR_RNDN);
        mpfr_set_d(r25516, y_im, MPFR_RNDN);
        mpfr_mul(r25517, r25515, r25516, MPFR_RNDN);
        mpfr_sub(r25518, r25514, r25517, MPFR_RNDN);
        return mpfr_get_d(r25518, MPFR_RNDN);
}