#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "math.square on complex, imaginary part";

double f_if(float re, float im) {
        float r20828 = re;
        float r20829 = im;
        float r20830 = r20828 * r20829;
        float r20831 = r20829 * r20828;
        float r20832 = r20830 + r20831;
        return r20832;
}

double f_id(double re, double im) {
        double r20833 = re;
        double r20834 = im;
        double r20835 = r20833 * r20834;
        double r20836 = r20834 * r20833;
        double r20837 = r20835 + r20836;
        return r20837;
}


double f_of(float re, float im) {
        float r20838 = re;
        float r20839 = r20838 + r20838;
        float r20840 = im;
        float r20841 = r20839 * r20840;
        return r20841;
}

double f_od(double re, double im) {
        double r20842 = re;
        double r20843 = r20842 + r20842;
        double r20844 = im;
        double r20845 = r20843 * r20844;
        return r20845;
}

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 r20846, r20847, r20848, r20849, r20850;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r20846);
        mpfr_init(r20847);
        mpfr_init(r20848);
        mpfr_init(r20849);
        mpfr_init(r20850);
}

double f_im(double re, double im) {
        mpfr_set_d(r20846, re, MPFR_RNDN);
        mpfr_set_d(r20847, im, MPFR_RNDN);
        mpfr_mul(r20848, r20846, r20847, MPFR_RNDN);
        mpfr_mul(r20849, r20847, r20846, MPFR_RNDN);
        mpfr_add(r20850, r20848, r20849, MPFR_RNDN);
        return mpfr_get_d(r20850, MPFR_RNDN);
}

static mpfr_t r20851, r20852, r20853, r20854;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20851);
        mpfr_init(r20852);
        mpfr_init(r20853);
        mpfr_init(r20854);
}

double f_fm(double re, double im) {
        mpfr_set_d(r20851, re, MPFR_RNDN);
        mpfr_add(r20852, r20851, r20851, MPFR_RNDN);
        mpfr_set_d(r20853, im, MPFR_RNDN);
        mpfr_mul(r20854, r20852, r20853, MPFR_RNDN);
        return mpfr_get_d(r20854, MPFR_RNDN);
}

static mpfr_t r20855, r20856, r20857, r20858;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20855);
        mpfr_init(r20856);
        mpfr_init(r20857);
        mpfr_init(r20858);
}

double f_dm(double re, double im) {
        mpfr_set_d(r20855, re, MPFR_RNDN);
        mpfr_add(r20856, r20855, r20855, MPFR_RNDN);
        mpfr_set_d(r20857, im, MPFR_RNDN);
        mpfr_mul(r20858, r20856, r20857, MPFR_RNDN);
        return mpfr_get_d(r20858, MPFR_RNDN);
}