#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 r18300 = x;
        float r18301 = r18300 / r18300;
        float r18302 = 1.0f;
        float r18303 = r18302 / r18300;
        float r18304 = r18300 * r18300;
        float r18305 = sqrt(r18304);
        float r18306 = r18303 * r18305;
        float r18307 = r18301 - r18306;
        return r18307;
}

double f_id(double x) {
        double r18308 = x;
        double r18309 = r18308 / r18308;
        double r18310 = 1.0;
        double r18311 = r18310 / r18308;
        double r18312 = r18308 * r18308;
        double r18313 = sqrt(r18312);
        double r18314 = r18311 * r18313;
        double r18315 = r18309 - r18314;
        return r18315;
}


double f_of(float x) {
        float r18316 = 1.0f;
        float r18317 = x;
        float r18318 = fabs(r18317);
        float r18319 = r18318 / r18317;
        float r18320 = r18316 - r18319;
        return r18320;
}

double f_od(double x) {
        double r18321 = 1.0;
        double r18322 = x;
        double r18323 = fabs(r18322);
        double r18324 = r18323 / r18322;
        double r18325 = r18321 - r18324;
        return r18325;
}

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 r18326, r18327, r18328, r18329, r18330, r18331, r18332, r18333;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r18326);
        mpfr_init(r18327);
        mpfr_init_set_str(r18328, "1", 10, MPFR_RNDN);
        mpfr_init(r18329);
        mpfr_init(r18330);
        mpfr_init(r18331);
        mpfr_init(r18332);
        mpfr_init(r18333);
}

double f_im(double x) {
        mpfr_set_d(r18326, x, MPFR_RNDN);
        mpfr_div(r18327, r18326, r18326, MPFR_RNDN);
        ;
        mpfr_div(r18329, r18328, r18326, MPFR_RNDN);
        mpfr_mul(r18330, r18326, r18326, MPFR_RNDN);
        mpfr_sqrt(r18331, r18330, MPFR_RNDN);
        mpfr_mul(r18332, r18329, r18331, MPFR_RNDN);
        mpfr_sub(r18333, r18327, r18332, MPFR_RNDN);
        return mpfr_get_d(r18333, MPFR_RNDN);
}

static mpfr_t r18334, r18335, r18336, r18337, r18338;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18334, "1", 10, MPFR_RNDN);
        mpfr_init(r18335);
        mpfr_init(r18336);
        mpfr_init(r18337);
        mpfr_init(r18338);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r18335, x, MPFR_RNDN);
        mpfr_abs(r18336, r18335, MPFR_RNDN);
        mpfr_div(r18337, r18336, r18335, MPFR_RNDN);
        mpfr_sub(r18338, r18334, r18337, MPFR_RNDN);
        return mpfr_get_d(r18338, MPFR_RNDN);
}

static mpfr_t r18339, r18340, r18341, r18342, r18343;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18339, "1", 10, MPFR_RNDN);
        mpfr_init(r18340);
        mpfr_init(r18341);
        mpfr_init(r18342);
        mpfr_init(r18343);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r18340, x, MPFR_RNDN);
        mpfr_abs(r18341, r18340, MPFR_RNDN);
        mpfr_div(r18342, r18341, r18340, MPFR_RNDN);
        mpfr_sub(r18343, r18339, r18342, MPFR_RNDN);
        return mpfr_get_d(r18343, MPFR_RNDN);
}