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

char *name = "sqrt sqr";

double f_if(float x) {
        float r10681 = x;
        float r10682 = r10681 / r10681;
        float r10683 = 1.0f;
        float r10684 = r10683 / r10681;
        float r10685 = r10681 * r10681;
        float r10686 = sqrt(r10685);
        float r10687 = r10684 * r10686;
        float r10688 = r10682 - r10687;
        return r10688;
}

double f_id(double x) {
        double r10689 = x;
        double r10690 = r10689 / r10689;
        double r10691 = 1.0;
        double r10692 = r10691 / r10689;
        double r10693 = r10689 * r10689;
        double r10694 = sqrt(r10693);
        double r10695 = r10692 * r10694;
        double r10696 = r10690 - r10695;
        return r10696;
}


double f_of(float x) {
        float r10697 = 1.0f;
        float r10698 = x;
        float r10699 = fabs(r10698);
        float r10700 = r10699 / r10698;
        float r10701 = r10697 - r10700;
        return r10701;
}

double f_od(double x) {
        double r10702 = 1.0;
        double r10703 = x;
        double r10704 = fabs(r10703);
        double r10705 = r10704 / r10703;
        double r10706 = r10702 - r10705;
        return r10706;
}

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 r10707, r10708, r10709, r10710, r10711, r10712, r10713, r10714;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r10707);
        mpfr_init(r10708);
        mpfr_init_set_str(r10709, "1", 10, MPFR_RNDN);
        mpfr_init(r10710);
        mpfr_init(r10711);
        mpfr_init(r10712);
        mpfr_init(r10713);
        mpfr_init(r10714);
}

double f_im(double x) {
        mpfr_set_d(r10707, x, MPFR_RNDN);
        mpfr_div(r10708, r10707, r10707, MPFR_RNDN);
        ;
        mpfr_div(r10710, r10709, r10707, MPFR_RNDN);
        mpfr_mul(r10711, r10707, r10707, MPFR_RNDN);
        mpfr_sqrt(r10712, r10711, MPFR_RNDN);
        mpfr_mul(r10713, r10710, r10712, MPFR_RNDN);
        mpfr_sub(r10714, r10708, r10713, MPFR_RNDN);
        return mpfr_get_d(r10714, MPFR_RNDN);
}

static mpfr_t r10715, r10716, r10717, r10718, r10719;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r10715, "1", 10, MPFR_RNDN);
        mpfr_init(r10716);
        mpfr_init(r10717);
        mpfr_init(r10718);
        mpfr_init(r10719);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r10716, x, MPFR_RNDN);
        mpfr_abs(r10717, r10716, MPFR_RNDN);
        mpfr_div(r10718, r10717, r10716, MPFR_RNDN);
        mpfr_sub(r10719, r10715, r10718, MPFR_RNDN);
        return mpfr_get_d(r10719, MPFR_RNDN);
}

static mpfr_t r10720, r10721, r10722, r10723, r10724;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r10720, "1", 10, MPFR_RNDN);
        mpfr_init(r10721);
        mpfr_init(r10722);
        mpfr_init(r10723);
        mpfr_init(r10724);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r10721, x, MPFR_RNDN);
        mpfr_abs(r10722, r10721, MPFR_RNDN);
        mpfr_div(r10723, r10722, r10721, MPFR_RNDN);
        mpfr_sub(r10724, r10720, r10723, MPFR_RNDN);
        return mpfr_get_d(r10724, MPFR_RNDN);
}