#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 r18478 = x;
        float r18479 = r18478 / r18478;
        float r18480 = 1.0f;
        float r18481 = r18480 / r18478;
        float r18482 = r18478 * r18478;
        float r18483 = sqrt(r18482);
        float r18484 = r18481 * r18483;
        float r18485 = r18479 - r18484;
        return r18485;
}

double f_id(double x) {
        double r18486 = x;
        double r18487 = r18486 / r18486;
        double r18488 = 1.0;
        double r18489 = r18488 / r18486;
        double r18490 = r18486 * r18486;
        double r18491 = sqrt(r18490);
        double r18492 = r18489 * r18491;
        double r18493 = r18487 - r18492;
        return r18493;
}


double f_of(float x) {
        float r18494 = 1.0f;
        float r18495 = x;
        float r18496 = fabs(r18495);
        float r18497 = r18496 / r18495;
        float r18498 = r18494 - r18497;
        return r18498;
}

double f_od(double x) {
        double r18499 = 1.0;
        double r18500 = x;
        double r18501 = fabs(r18500);
        double r18502 = r18501 / r18500;
        double r18503 = r18499 - r18502;
        return r18503;
}

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 r18504, r18505, r18506, r18507, r18508, r18509, r18510, r18511;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r18504);
        mpfr_init(r18505);
        mpfr_init_set_str(r18506, "1", 10, MPFR_RNDN);
        mpfr_init(r18507);
        mpfr_init(r18508);
        mpfr_init(r18509);
        mpfr_init(r18510);
        mpfr_init(r18511);
}

double f_im(double x) {
        mpfr_set_d(r18504, x, MPFR_RNDN);
        mpfr_div(r18505, r18504, r18504, MPFR_RNDN);
        ;
        mpfr_div(r18507, r18506, r18504, MPFR_RNDN);
        mpfr_mul(r18508, r18504, r18504, MPFR_RNDN);
        mpfr_sqrt(r18509, r18508, MPFR_RNDN);
        mpfr_mul(r18510, r18507, r18509, MPFR_RNDN);
        mpfr_sub(r18511, r18505, r18510, MPFR_RNDN);
        return mpfr_get_d(r18511, MPFR_RNDN);
}

static mpfr_t r18512, r18513, r18514, r18515, r18516;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18512, "1", 10, MPFR_RNDN);
        mpfr_init(r18513);
        mpfr_init(r18514);
        mpfr_init(r18515);
        mpfr_init(r18516);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r18513, x, MPFR_RNDN);
        mpfr_abs(r18514, r18513, MPFR_RNDN);
        mpfr_div(r18515, r18514, r18513, MPFR_RNDN);
        mpfr_sub(r18516, r18512, r18515, MPFR_RNDN);
        return mpfr_get_d(r18516, MPFR_RNDN);
}

static mpfr_t r18517, r18518, r18519, r18520, r18521;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18517, "1", 10, MPFR_RNDN);
        mpfr_init(r18518);
        mpfr_init(r18519);
        mpfr_init(r18520);
        mpfr_init(r18521);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r18518, x, MPFR_RNDN);
        mpfr_abs(r18519, r18518, MPFR_RNDN);
        mpfr_div(r18520, r18519, r18518, MPFR_RNDN);
        mpfr_sub(r18521, r18517, r18520, MPFR_RNDN);
        return mpfr_get_d(r18521, MPFR_RNDN);
}