#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 r8733 = re;
        float r8734 = im;
        float r8735 = r8733 * r8734;
        float r8736 = r8734 * r8733;
        float r8737 = r8735 + r8736;
        return r8737;
}

double f_id(double re, double im) {
        double r8738 = re;
        double r8739 = im;
        double r8740 = r8738 * r8739;
        double r8741 = r8739 * r8738;
        double r8742 = r8740 + r8741;
        return r8742;
}


double f_of(float re, float im) {
        float r8743 = re;
        float r8744 = im;
        float r8745 = r8743 * r8744;
        float r8746 = r8744 * r8743;
        float r8747 = r8745 + r8746;
        return r8747;
}

double f_od(double re, double im) {
        double r8748 = re;
        double r8749 = im;
        double r8750 = r8748 * r8749;
        double r8751 = r8749 * r8748;
        double r8752 = r8750 + r8751;
        return r8752;
}

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 r8753, r8754, r8755, r8756, r8757;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r8753);
        mpfr_init(r8754);
        mpfr_init(r8755);
        mpfr_init(r8756);
        mpfr_init(r8757);
}

double f_im(double re, double im) {
        mpfr_set_d(r8753, re, MPFR_RNDN);
        mpfr_set_d(r8754, im, MPFR_RNDN);
        mpfr_mul(r8755, r8753, r8754, MPFR_RNDN);
        mpfr_mul(r8756, r8754, r8753, MPFR_RNDN);
        mpfr_add(r8757, r8755, r8756, MPFR_RNDN);
        return mpfr_get_d(r8757, MPFR_RNDN);
}

static mpfr_t r8758, r8759, r8760, r8761, r8762;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r8758);
        mpfr_init(r8759);
        mpfr_init(r8760);
        mpfr_init(r8761);
        mpfr_init(r8762);
}

double f_fm(double re, double im) {
        mpfr_set_d(r8758, re, MPFR_RNDN);
        mpfr_set_d(r8759, im, MPFR_RNDN);
        mpfr_mul(r8760, r8758, r8759, MPFR_RNDN);
        mpfr_mul(r8761, r8759, r8758, MPFR_RNDN);
        mpfr_add(r8762, r8760, r8761, MPFR_RNDN);
        return mpfr_get_d(r8762, MPFR_RNDN);
}

static mpfr_t r8763, r8764, r8765, r8766, r8767;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r8763);
        mpfr_init(r8764);
        mpfr_init(r8765);
        mpfr_init(r8766);
        mpfr_init(r8767);
}

double f_dm(double re, double im) {
        mpfr_set_d(r8763, re, MPFR_RNDN);
        mpfr_set_d(r8764, im, MPFR_RNDN);
        mpfr_mul(r8765, r8763, r8764, MPFR_RNDN);
        mpfr_mul(r8766, r8764, r8763, MPFR_RNDN);
        mpfr_add(r8767, r8765, r8766, MPFR_RNDN);
        return mpfr_get_d(r8767, MPFR_RNDN);
}