#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 r18048 = x;
        float r18049 = r18048 / r18048;
        float r18050 = 1.0f;
        float r18051 = r18050 / r18048;
        float r18052 = r18048 * r18048;
        float r18053 = sqrt(r18052);
        float r18054 = r18051 * r18053;
        float r18055 = r18049 - r18054;
        return r18055;
}

double f_id(double x) {
        double r18056 = x;
        double r18057 = r18056 / r18056;
        double r18058 = 1.0;
        double r18059 = r18058 / r18056;
        double r18060 = r18056 * r18056;
        double r18061 = sqrt(r18060);
        double r18062 = r18059 * r18061;
        double r18063 = r18057 - r18062;
        return r18063;
}


double f_of(float x) {
        float r18064 = 1.0f;
        float r18065 = x;
        float r18066 = fabs(r18065);
        float r18067 = r18066 / r18065;
        float r18068 = r18064 - r18067;
        return r18068;
}

double f_od(double x) {
        double r18069 = 1.0;
        double r18070 = x;
        double r18071 = fabs(r18070);
        double r18072 = r18071 / r18070;
        double r18073 = r18069 - r18072;
        return r18073;
}

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 r18074, r18075, r18076, r18077, r18078, r18079, r18080, r18081;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r18074);
        mpfr_init(r18075);
        mpfr_init_set_str(r18076, "1", 10, MPFR_RNDN);
        mpfr_init(r18077);
        mpfr_init(r18078);
        mpfr_init(r18079);
        mpfr_init(r18080);
        mpfr_init(r18081);
}

double f_im(double x) {
        mpfr_set_d(r18074, x, MPFR_RNDN);
        mpfr_div(r18075, r18074, r18074, MPFR_RNDN);
        ;
        mpfr_div(r18077, r18076, r18074, MPFR_RNDN);
        mpfr_mul(r18078, r18074, r18074, MPFR_RNDN);
        mpfr_sqrt(r18079, r18078, MPFR_RNDN);
        mpfr_mul(r18080, r18077, r18079, MPFR_RNDN);
        mpfr_sub(r18081, r18075, r18080, MPFR_RNDN);
        return mpfr_get_d(r18081, MPFR_RNDN);
}

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

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

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

static mpfr_t r18087, r18088, r18089, r18090, r18091;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18087, "1", 10, MPFR_RNDN);
        mpfr_init(r18088);
        mpfr_init(r18089);
        mpfr_init(r18090);
        mpfr_init(r18091);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r18088, x, MPFR_RNDN);
        mpfr_abs(r18089, r18088, MPFR_RNDN);
        mpfr_div(r18090, r18089, r18088, MPFR_RNDN);
        mpfr_sub(r18091, r18087, r18090, MPFR_RNDN);
        return mpfr_get_d(r18091, MPFR_RNDN);
}