#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 r8466 = re;
        float r8467 = r8466 * r8466;
        float r8468 = im;
        float r8469 = r8468 * r8468;
        float r8470 = r8467 - r8469;
        return r8470;
}

double f_id(double re, double im) {
        double r8471 = re;
        double r8472 = r8471 * r8471;
        double r8473 = im;
        double r8474 = r8473 * r8473;
        double r8475 = r8472 - r8474;
        return r8475;
}


double f_of(float re, float im) {
        float r8476 = re;
        float r8477 = im;
        float r8478 = r8476 + r8477;
        float r8479 = r8476 - r8477;
        float r8480 = r8478 * r8479;
        return r8480;
}

double f_od(double re, double im) {
        double r8481 = re;
        double r8482 = im;
        double r8483 = r8481 + r8482;
        double r8484 = r8481 - r8482;
        double r8485 = r8483 * r8484;
        return r8485;
}

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 r8486, r8487, r8488, r8489, r8490;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r8486);
        mpfr_init(r8487);
        mpfr_init(r8488);
        mpfr_init(r8489);
        mpfr_init(r8490);
}

double f_im(double re, double im) {
        mpfr_set_d(r8486, re, MPFR_RNDN);
        mpfr_mul(r8487, r8486, r8486, MPFR_RNDN);
        mpfr_set_d(r8488, im, MPFR_RNDN);
        mpfr_mul(r8489, r8488, r8488, MPFR_RNDN);
        mpfr_sub(r8490, r8487, r8489, MPFR_RNDN);
        return mpfr_get_d(r8490, MPFR_RNDN);
}

static mpfr_t r8491, r8492, r8493, r8494, r8495;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r8491);
        mpfr_init(r8492);
        mpfr_init(r8493);
        mpfr_init(r8494);
        mpfr_init(r8495);
}

double f_fm(double re, double im) {
        mpfr_set_d(r8491, re, MPFR_RNDN);
        mpfr_set_d(r8492, im, MPFR_RNDN);
        mpfr_add(r8493, r8491, r8492, MPFR_RNDN);
        mpfr_sub(r8494, r8491, r8492, MPFR_RNDN);
        mpfr_mul(r8495, r8493, r8494, MPFR_RNDN);
        return mpfr_get_d(r8495, MPFR_RNDN);
}

static mpfr_t r8496, r8497, r8498, r8499, r8500;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r8496);
        mpfr_init(r8497);
        mpfr_init(r8498);
        mpfr_init(r8499);
        mpfr_init(r8500);
}

double f_dm(double re, double im) {
        mpfr_set_d(r8496, re, MPFR_RNDN);
        mpfr_set_d(r8497, im, MPFR_RNDN);
        mpfr_add(r8498, r8496, r8497, MPFR_RNDN);
        mpfr_sub(r8499, r8496, r8497, MPFR_RNDN);
        mpfr_mul(r8500, r8498, r8499, MPFR_RNDN);
        return mpfr_get_d(r8500, MPFR_RNDN);
}