#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 r21051 = re;
        float r21052 = r21051 * r21051;
        float r21053 = im;
        float r21054 = r21053 * r21053;
        float r21055 = r21052 - r21054;
        return r21055;
}

double f_id(double re, double im) {
        double r21056 = re;
        double r21057 = r21056 * r21056;
        double r21058 = im;
        double r21059 = r21058 * r21058;
        double r21060 = r21057 - r21059;
        return r21060;
}


double f_of(float re, float im) {
        float r21061 = re;
        float r21062 = im;
        float r21063 = r21061 + r21062;
        float r21064 = r21061 - r21062;
        float r21065 = r21063 * r21064;
        return r21065;
}

double f_od(double re, double im) {
        double r21066 = re;
        double r21067 = im;
        double r21068 = r21066 + r21067;
        double r21069 = r21066 - r21067;
        double r21070 = r21068 * r21069;
        return r21070;
}

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 r21071, r21072, r21073, r21074, r21075;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r21071);
        mpfr_init(r21072);
        mpfr_init(r21073);
        mpfr_init(r21074);
        mpfr_init(r21075);
}

double f_im(double re, double im) {
        mpfr_set_d(r21071, re, MPFR_RNDN);
        mpfr_mul(r21072, r21071, r21071, MPFR_RNDN);
        mpfr_set_d(r21073, im, MPFR_RNDN);
        mpfr_mul(r21074, r21073, r21073, MPFR_RNDN);
        mpfr_sub(r21075, r21072, r21074, MPFR_RNDN);
        return mpfr_get_d(r21075, MPFR_RNDN);
}

static mpfr_t r21076, r21077, r21078, r21079, r21080;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r21076);
        mpfr_init(r21077);
        mpfr_init(r21078);
        mpfr_init(r21079);
        mpfr_init(r21080);
}

double f_fm(double re, double im) {
        mpfr_set_d(r21076, re, MPFR_RNDN);
        mpfr_set_d(r21077, im, MPFR_RNDN);
        mpfr_add(r21078, r21076, r21077, MPFR_RNDN);
        mpfr_sub(r21079, r21076, r21077, MPFR_RNDN);
        mpfr_mul(r21080, r21078, r21079, MPFR_RNDN);
        return mpfr_get_d(r21080, MPFR_RNDN);
}

static mpfr_t r21081, r21082, r21083, r21084, r21085;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r21081);
        mpfr_init(r21082);
        mpfr_init(r21083);
        mpfr_init(r21084);
        mpfr_init(r21085);
}

double f_dm(double re, double im) {
        mpfr_set_d(r21081, re, MPFR_RNDN);
        mpfr_set_d(r21082, im, MPFR_RNDN);
        mpfr_add(r21083, r21081, r21082, MPFR_RNDN);
        mpfr_sub(r21084, r21081, r21082, MPFR_RNDN);
        mpfr_mul(r21085, r21083, r21084, MPFR_RNDN);
        return mpfr_get_d(r21085, MPFR_RNDN);
}