#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 r7802 = re;
        float r7803 = r7802 * r7802;
        float r7804 = im;
        float r7805 = r7804 * r7804;
        float r7806 = r7803 - r7805;
        return r7806;
}

double f_id(double re, double im) {
        double r7807 = re;
        double r7808 = r7807 * r7807;
        double r7809 = im;
        double r7810 = r7809 * r7809;
        double r7811 = r7808 - r7810;
        return r7811;
}


double f_of(float re, float im) {
        float r7812 = im;
        float r7813 = re;
        float r7814 = r7812 + r7813;
        float r7815 = r7813 - r7812;
        float r7816 = r7814 * r7815;
        return r7816;
}

double f_od(double re, double im) {
        double r7817 = im;
        double r7818 = re;
        double r7819 = r7817 + r7818;
        double r7820 = r7818 - r7817;
        double r7821 = r7819 * r7820;
        return r7821;
}

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 r7822, r7823, r7824, r7825, r7826;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r7822);
        mpfr_init(r7823);
        mpfr_init(r7824);
        mpfr_init(r7825);
        mpfr_init(r7826);
}

double f_im(double re, double im) {
        mpfr_set_d(r7822, re, MPFR_RNDN);
        mpfr_mul(r7823, r7822, r7822, MPFR_RNDN);
        mpfr_set_d(r7824, im, MPFR_RNDN);
        mpfr_mul(r7825, r7824, r7824, MPFR_RNDN);
        mpfr_sub(r7826, r7823, r7825, MPFR_RNDN);
        return mpfr_get_d(r7826, MPFR_RNDN);
}

static mpfr_t r7827, r7828, r7829, r7830, r7831;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r7827);
        mpfr_init(r7828);
        mpfr_init(r7829);
        mpfr_init(r7830);
        mpfr_init(r7831);
}

double f_fm(double re, double im) {
        mpfr_set_d(r7827, im, MPFR_RNDN);
        mpfr_set_d(r7828, re, MPFR_RNDN);
        mpfr_add(r7829, r7827, r7828, MPFR_RNDN);
        mpfr_sub(r7830, r7828, r7827, MPFR_RNDN);
        mpfr_mul(r7831, r7829, r7830, MPFR_RNDN);
        return mpfr_get_d(r7831, MPFR_RNDN);
}

static mpfr_t r7832, r7833, r7834, r7835, r7836;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r7832);
        mpfr_init(r7833);
        mpfr_init(r7834);
        mpfr_init(r7835);
        mpfr_init(r7836);
}

double f_dm(double re, double im) {
        mpfr_set_d(r7832, im, MPFR_RNDN);
        mpfr_set_d(r7833, re, MPFR_RNDN);
        mpfr_add(r7834, r7832, r7833, MPFR_RNDN);
        mpfr_sub(r7835, r7833, r7832, MPFR_RNDN);
        mpfr_mul(r7836, r7834, r7835, MPFR_RNDN);
        return mpfr_get_d(r7836, MPFR_RNDN);
}