#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 r20793 = re;
        float r20794 = r20793 * r20793;
        float r20795 = im;
        float r20796 = r20795 * r20795;
        float r20797 = r20794 - r20796;
        return r20797;
}

double f_id(double re, double im) {
        double r20798 = re;
        double r20799 = r20798 * r20798;
        double r20800 = im;
        double r20801 = r20800 * r20800;
        double r20802 = r20799 - r20801;
        return r20802;
}


double f_of(float re, float im) {
        float r20803 = re;
        float r20804 = im;
        float r20805 = r20803 + r20804;
        float r20806 = r20803 - r20804;
        float r20807 = r20805 * r20806;
        return r20807;
}

double f_od(double re, double im) {
        double r20808 = re;
        double r20809 = im;
        double r20810 = r20808 + r20809;
        double r20811 = r20808 - r20809;
        double r20812 = r20810 * r20811;
        return r20812;
}

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 r20813, r20814, r20815, r20816, r20817;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r20813);
        mpfr_init(r20814);
        mpfr_init(r20815);
        mpfr_init(r20816);
        mpfr_init(r20817);
}

double f_im(double re, double im) {
        mpfr_set_d(r20813, re, MPFR_RNDN);
        mpfr_mul(r20814, r20813, r20813, MPFR_RNDN);
        mpfr_set_d(r20815, im, MPFR_RNDN);
        mpfr_mul(r20816, r20815, r20815, MPFR_RNDN);
        mpfr_sub(r20817, r20814, r20816, MPFR_RNDN);
        return mpfr_get_d(r20817, MPFR_RNDN);
}

static mpfr_t r20818, r20819, r20820, r20821, r20822;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20818);
        mpfr_init(r20819);
        mpfr_init(r20820);
        mpfr_init(r20821);
        mpfr_init(r20822);
}

double f_fm(double re, double im) {
        mpfr_set_d(r20818, re, MPFR_RNDN);
        mpfr_set_d(r20819, im, MPFR_RNDN);
        mpfr_add(r20820, r20818, r20819, MPFR_RNDN);
        mpfr_sub(r20821, r20818, r20819, MPFR_RNDN);
        mpfr_mul(r20822, r20820, r20821, MPFR_RNDN);
        return mpfr_get_d(r20822, MPFR_RNDN);
}

static mpfr_t r20823, r20824, r20825, r20826, r20827;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20823);
        mpfr_init(r20824);
        mpfr_init(r20825);
        mpfr_init(r20826);
        mpfr_init(r20827);
}

double f_dm(double re, double im) {
        mpfr_set_d(r20823, re, MPFR_RNDN);
        mpfr_set_d(r20824, im, MPFR_RNDN);
        mpfr_add(r20825, r20823, r20824, MPFR_RNDN);
        mpfr_sub(r20826, r20823, r20824, MPFR_RNDN);
        mpfr_mul(r20827, r20825, r20826, MPFR_RNDN);
        return mpfr_get_d(r20827, MPFR_RNDN);
}