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

char *name = "math.square on complex, real part";

double f_if(float re, float im) {
        float r20187 = re;
        float r20188 = r20187 * r20187;
        float r20189 = im;
        float r20190 = r20189 * r20189;
        float r20191 = r20188 - r20190;
        return r20191;
}

double f_id(double re, double im) {
        double r20192 = re;
        double r20193 = r20192 * r20192;
        double r20194 = im;
        double r20195 = r20194 * r20194;
        double r20196 = r20193 - r20195;
        return r20196;
}


double f_of(float re, float im) {
        float r20197 = re;
        float r20198 = im;
        float r20199 = r20197 + r20198;
        float r20200 = r20197 - r20198;
        float r20201 = r20199 * r20200;
        return r20201;
}

double f_od(double re, double im) {
        double r20202 = re;
        double r20203 = im;
        double r20204 = r20202 + r20203;
        double r20205 = r20202 - r20203;
        double r20206 = r20204 * r20205;
        return r20206;
}

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 r20207, r20208, r20209, r20210, r20211;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r20207);
        mpfr_init(r20208);
        mpfr_init(r20209);
        mpfr_init(r20210);
        mpfr_init(r20211);
}

double f_im(double re, double im) {
        mpfr_set_d(r20207, re, MPFR_RNDN);
        mpfr_mul(r20208, r20207, r20207, MPFR_RNDN);
        mpfr_set_d(r20209, im, MPFR_RNDN);
        mpfr_mul(r20210, r20209, r20209, MPFR_RNDN);
        mpfr_sub(r20211, r20208, r20210, MPFR_RNDN);
        return mpfr_get_d(r20211, MPFR_RNDN);
}

static mpfr_t r20212, r20213, r20214, r20215, r20216;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20212);
        mpfr_init(r20213);
        mpfr_init(r20214);
        mpfr_init(r20215);
        mpfr_init(r20216);
}

double f_fm(double re, double im) {
        mpfr_set_d(r20212, re, MPFR_RNDN);
        mpfr_set_d(r20213, im, MPFR_RNDN);
        mpfr_add(r20214, r20212, r20213, MPFR_RNDN);
        mpfr_sub(r20215, r20212, r20213, MPFR_RNDN);
        mpfr_mul(r20216, r20214, r20215, MPFR_RNDN);
        return mpfr_get_d(r20216, MPFR_RNDN);
}

static mpfr_t r20217, r20218, r20219, r20220, r20221;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20217);
        mpfr_init(r20218);
        mpfr_init(r20219);
        mpfr_init(r20220);
        mpfr_init(r20221);
}

double f_dm(double re, double im) {
        mpfr_set_d(r20217, re, MPFR_RNDN);
        mpfr_set_d(r20218, im, MPFR_RNDN);
        mpfr_add(r20219, r20217, r20218, MPFR_RNDN);
        mpfr_sub(r20220, r20217, r20218, MPFR_RNDN);
        mpfr_mul(r20221, r20219, r20220, MPFR_RNDN);
        return mpfr_get_d(r20221, MPFR_RNDN);
}