#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 r20121 = re;
        float r20122 = im;
        float r20123 = r20121 * r20122;
        float r20124 = r20122 * r20121;
        float r20125 = r20123 + r20124;
        return r20125;
}

double f_id(double re, double im) {
        double r20126 = re;
        double r20127 = im;
        double r20128 = r20126 * r20127;
        double r20129 = r20127 * r20126;
        double r20130 = r20128 + r20129;
        return r20130;
}


double f_of(float re, float im) {
        float r20131 = re;
        float r20132 = im;
        float r20133 = r20132 + r20132;
        float r20134 = r20131 * r20133;
        return r20134;
}

double f_od(double re, double im) {
        double r20135 = re;
        double r20136 = im;
        double r20137 = r20136 + r20136;
        double r20138 = r20135 * r20137;
        return r20138;
}

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 r20139, r20140, r20141, r20142, r20143;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r20139);
        mpfr_init(r20140);
        mpfr_init(r20141);
        mpfr_init(r20142);
        mpfr_init(r20143);
}

double f_im(double re, double im) {
        mpfr_set_d(r20139, re, MPFR_RNDN);
        mpfr_set_d(r20140, im, MPFR_RNDN);
        mpfr_mul(r20141, r20139, r20140, MPFR_RNDN);
        mpfr_mul(r20142, r20140, r20139, MPFR_RNDN);
        mpfr_add(r20143, r20141, r20142, MPFR_RNDN);
        return mpfr_get_d(r20143, MPFR_RNDN);
}

static mpfr_t r20144, r20145, r20146, r20147;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20144);
        mpfr_init(r20145);
        mpfr_init(r20146);
        mpfr_init(r20147);
}

double f_fm(double re, double im) {
        mpfr_set_d(r20144, re, MPFR_RNDN);
        mpfr_set_d(r20145, im, MPFR_RNDN);
        mpfr_add(r20146, r20145, r20145, MPFR_RNDN);
        mpfr_mul(r20147, r20144, r20146, MPFR_RNDN);
        return mpfr_get_d(r20147, MPFR_RNDN);
}

static mpfr_t r20148, r20149, r20150, r20151;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20148);
        mpfr_init(r20149);
        mpfr_init(r20150);
        mpfr_init(r20151);
}

double f_dm(double re, double im) {
        mpfr_set_d(r20148, re, MPFR_RNDN);
        mpfr_set_d(r20149, im, MPFR_RNDN);
        mpfr_add(r20150, r20149, r20149, MPFR_RNDN);
        mpfr_mul(r20151, r20148, r20150, MPFR_RNDN);
        return mpfr_get_d(r20151, MPFR_RNDN);
}