#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 r20821 = re;
        float r20822 = r20821 * r20821;
        float r20823 = im;
        float r20824 = r20823 * r20823;
        float r20825 = r20822 - r20824;
        return r20825;
}

double f_id(double re, double im) {
        double r20826 = re;
        double r20827 = r20826 * r20826;
        double r20828 = im;
        double r20829 = r20828 * r20828;
        double r20830 = r20827 - r20829;
        return r20830;
}


double f_of(float re, float im) {
        float r20831 = re;
        float r20832 = im;
        float r20833 = r20831 + r20832;
        float r20834 = r20831 - r20832;
        float r20835 = r20833 * r20834;
        return r20835;
}

double f_od(double re, double im) {
        double r20836 = re;
        double r20837 = im;
        double r20838 = r20836 + r20837;
        double r20839 = r20836 - r20837;
        double r20840 = r20838 * r20839;
        return r20840;
}

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 r20841, r20842, r20843, r20844, r20845;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r20841);
        mpfr_init(r20842);
        mpfr_init(r20843);
        mpfr_init(r20844);
        mpfr_init(r20845);
}

double f_im(double re, double im) {
        mpfr_set_d(r20841, re, MPFR_RNDN);
        mpfr_mul(r20842, r20841, r20841, MPFR_RNDN);
        mpfr_set_d(r20843, im, MPFR_RNDN);
        mpfr_mul(r20844, r20843, r20843, MPFR_RNDN);
        mpfr_sub(r20845, r20842, r20844, MPFR_RNDN);
        return mpfr_get_d(r20845, MPFR_RNDN);
}

static mpfr_t r20846, r20847, r20848, r20849, r20850;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20846);
        mpfr_init(r20847);
        mpfr_init(r20848);
        mpfr_init(r20849);
        mpfr_init(r20850);
}

double f_fm(double re, double im) {
        mpfr_set_d(r20846, re, MPFR_RNDN);
        mpfr_set_d(r20847, im, MPFR_RNDN);
        mpfr_add(r20848, r20846, r20847, MPFR_RNDN);
        mpfr_sub(r20849, r20846, r20847, MPFR_RNDN);
        mpfr_mul(r20850, r20848, r20849, MPFR_RNDN);
        return mpfr_get_d(r20850, MPFR_RNDN);
}

static mpfr_t r20851, r20852, r20853, r20854, r20855;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20851);
        mpfr_init(r20852);
        mpfr_init(r20853);
        mpfr_init(r20854);
        mpfr_init(r20855);
}

double f_dm(double re, double im) {
        mpfr_set_d(r20851, re, MPFR_RNDN);
        mpfr_set_d(r20852, im, MPFR_RNDN);
        mpfr_add(r20853, r20851, r20852, MPFR_RNDN);
        mpfr_sub(r20854, r20851, r20852, MPFR_RNDN);
        mpfr_mul(r20855, r20853, r20854, MPFR_RNDN);
        return mpfr_get_d(r20855, MPFR_RNDN);
}