#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 r18044 = x;
        float r18045 = r18044 / r18044;
        float r18046 = 1.0f;
        float r18047 = r18046 / r18044;
        float r18048 = r18044 * r18044;
        float r18049 = sqrt(r18048);
        float r18050 = r18047 * r18049;
        float r18051 = r18045 - r18050;
        return r18051;
}

double f_id(double x) {
        double r18052 = x;
        double r18053 = r18052 / r18052;
        double r18054 = 1.0;
        double r18055 = r18054 / r18052;
        double r18056 = r18052 * r18052;
        double r18057 = sqrt(r18056);
        double r18058 = r18055 * r18057;
        double r18059 = r18053 - r18058;
        return r18059;
}


double f_of(float x) {
        float r18060 = 1.0f;
        float r18061 = x;
        float r18062 = fabs(r18061);
        float r18063 = r18062 / r18061;
        float r18064 = r18060 - r18063;
        return r18064;
}

double f_od(double x) {
        double r18065 = 1.0;
        double r18066 = x;
        double r18067 = fabs(r18066);
        double r18068 = r18067 / r18066;
        double r18069 = r18065 - r18068;
        return r18069;
}

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 r18070, r18071, r18072, r18073, r18074, r18075, r18076, r18077;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r18070);
        mpfr_init(r18071);
        mpfr_init_set_str(r18072, "1", 10, MPFR_RNDN);
        mpfr_init(r18073);
        mpfr_init(r18074);
        mpfr_init(r18075);
        mpfr_init(r18076);
        mpfr_init(r18077);
}

double f_im(double x) {
        mpfr_set_d(r18070, x, MPFR_RNDN);
        mpfr_div(r18071, r18070, r18070, MPFR_RNDN);
        ;
        mpfr_div(r18073, r18072, r18070, MPFR_RNDN);
        mpfr_mul(r18074, r18070, r18070, MPFR_RNDN);
        mpfr_sqrt(r18075, r18074, MPFR_RNDN);
        mpfr_mul(r18076, r18073, r18075, MPFR_RNDN);
        mpfr_sub(r18077, r18071, r18076, MPFR_RNDN);
        return mpfr_get_d(r18077, MPFR_RNDN);
}

static mpfr_t r18078, r18079, r18080, r18081, r18082;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18078, "1", 10, MPFR_RNDN);
        mpfr_init(r18079);
        mpfr_init(r18080);
        mpfr_init(r18081);
        mpfr_init(r18082);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r18079, x, MPFR_RNDN);
        mpfr_abs(r18080, r18079, MPFR_RNDN);
        mpfr_div(r18081, r18080, r18079, MPFR_RNDN);
        mpfr_sub(r18082, r18078, r18081, MPFR_RNDN);
        return mpfr_get_d(r18082, MPFR_RNDN);
}

static mpfr_t r18083, r18084, r18085, r18086, r18087;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18083, "1", 10, MPFR_RNDN);
        mpfr_init(r18084);
        mpfr_init(r18085);
        mpfr_init(r18086);
        mpfr_init(r18087);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r18084, x, MPFR_RNDN);
        mpfr_abs(r18085, r18084, MPFR_RNDN);
        mpfr_div(r18086, r18085, r18084, MPFR_RNDN);
        mpfr_sub(r18087, r18083, r18086, MPFR_RNDN);
        return mpfr_get_d(r18087, MPFR_RNDN);
}