#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 r20833 = re;
        float r20834 = r20833 * r20833;
        float r20835 = im;
        float r20836 = r20835 * r20835;
        float r20837 = r20834 - r20836;
        return r20837;
}

double f_id(double re, double im) {
        double r20838 = re;
        double r20839 = r20838 * r20838;
        double r20840 = im;
        double r20841 = r20840 * r20840;
        double r20842 = r20839 - r20841;
        return r20842;
}


double f_of(float re, float im) {
        float r20843 = re;
        float r20844 = r20843 * r20843;
        float r20845 = im;
        float r20846 = r20845 * r20845;
        float r20847 = r20844 - r20846;
        return r20847;
}

double f_od(double re, double im) {
        double r20848 = re;
        double r20849 = r20848 * r20848;
        double r20850 = im;
        double r20851 = r20850 * r20850;
        double r20852 = r20849 - r20851;
        return r20852;
}

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 r20853, r20854, r20855, r20856, r20857;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r20853);
        mpfr_init(r20854);
        mpfr_init(r20855);
        mpfr_init(r20856);
        mpfr_init(r20857);
}

double f_im(double re, double im) {
        mpfr_set_d(r20853, re, MPFR_RNDN);
        mpfr_mul(r20854, r20853, r20853, MPFR_RNDN);
        mpfr_set_d(r20855, im, MPFR_RNDN);
        mpfr_mul(r20856, r20855, r20855, MPFR_RNDN);
        mpfr_sub(r20857, r20854, r20856, MPFR_RNDN);
        return mpfr_get_d(r20857, MPFR_RNDN);
}

static mpfr_t r20858, r20859, r20860, r20861, r20862;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20858);
        mpfr_init(r20859);
        mpfr_init(r20860);
        mpfr_init(r20861);
        mpfr_init(r20862);
}

double f_fm(double re, double im) {
        mpfr_set_d(r20858, re, MPFR_RNDN);
        mpfr_mul(r20859, r20858, r20858, MPFR_RNDN);
        mpfr_set_d(r20860, im, MPFR_RNDN);
        mpfr_mul(r20861, r20860, r20860, MPFR_RNDN);
        mpfr_sub(r20862, r20859, r20861, MPFR_RNDN);
        return mpfr_get_d(r20862, MPFR_RNDN);
}

static mpfr_t r20863, r20864, r20865, r20866, r20867;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20863);
        mpfr_init(r20864);
        mpfr_init(r20865);
        mpfr_init(r20866);
        mpfr_init(r20867);
}

double f_dm(double re, double im) {
        mpfr_set_d(r20863, re, MPFR_RNDN);
        mpfr_mul(r20864, r20863, r20863, MPFR_RNDN);
        mpfr_set_d(r20865, im, MPFR_RNDN);
        mpfr_mul(r20866, r20865, r20865, MPFR_RNDN);
        mpfr_sub(r20867, r20864, r20866, MPFR_RNDN);
        return mpfr_get_d(r20867, MPFR_RNDN);
}