#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 r18193 = x;
        float r18194 = r18193 / r18193;
        float r18195 = 1.0f;
        float r18196 = r18195 / r18193;
        float r18197 = r18193 * r18193;
        float r18198 = sqrt(r18197);
        float r18199 = r18196 * r18198;
        float r18200 = r18194 - r18199;
        return r18200;
}

double f_id(double x) {
        double r18201 = x;
        double r18202 = r18201 / r18201;
        double r18203 = 1.0;
        double r18204 = r18203 / r18201;
        double r18205 = r18201 * r18201;
        double r18206 = sqrt(r18205);
        double r18207 = r18204 * r18206;
        double r18208 = r18202 - r18207;
        return r18208;
}


double f_of(float x) {
        float r18209 = 1.0f;
        float r18210 = x;
        float r18211 = fabs(r18210);
        float r18212 = r18211 / r18210;
        float r18213 = r18209 - r18212;
        return r18213;
}

double f_od(double x) {
        double r18214 = 1.0;
        double r18215 = x;
        double r18216 = fabs(r18215);
        double r18217 = r18216 / r18215;
        double r18218 = r18214 - r18217;
        return r18218;
}

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 r18219, r18220, r18221, r18222, r18223, r18224, r18225, r18226;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r18219);
        mpfr_init(r18220);
        mpfr_init_set_str(r18221, "1", 10, MPFR_RNDN);
        mpfr_init(r18222);
        mpfr_init(r18223);
        mpfr_init(r18224);
        mpfr_init(r18225);
        mpfr_init(r18226);
}

double f_im(double x) {
        mpfr_set_d(r18219, x, MPFR_RNDN);
        mpfr_div(r18220, r18219, r18219, MPFR_RNDN);
        ;
        mpfr_div(r18222, r18221, r18219, MPFR_RNDN);
        mpfr_mul(r18223, r18219, r18219, MPFR_RNDN);
        mpfr_sqrt(r18224, r18223, MPFR_RNDN);
        mpfr_mul(r18225, r18222, r18224, MPFR_RNDN);
        mpfr_sub(r18226, r18220, r18225, MPFR_RNDN);
        return mpfr_get_d(r18226, MPFR_RNDN);
}

static mpfr_t r18227, r18228, r18229, r18230, r18231;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18227, "1", 10, MPFR_RNDN);
        mpfr_init(r18228);
        mpfr_init(r18229);
        mpfr_init(r18230);
        mpfr_init(r18231);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r18228, x, MPFR_RNDN);
        mpfr_abs(r18229, r18228, MPFR_RNDN);
        mpfr_div(r18230, r18229, r18228, MPFR_RNDN);
        mpfr_sub(r18231, r18227, r18230, MPFR_RNDN);
        return mpfr_get_d(r18231, MPFR_RNDN);
}

static mpfr_t r18232, r18233, r18234, r18235, r18236;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18232, "1", 10, MPFR_RNDN);
        mpfr_init(r18233);
        mpfr_init(r18234);
        mpfr_init(r18235);
        mpfr_init(r18236);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r18233, x, MPFR_RNDN);
        mpfr_abs(r18234, r18233, MPFR_RNDN);
        mpfr_div(r18235, r18234, r18233, MPFR_RNDN);
        mpfr_sub(r18236, r18232, r18235, MPFR_RNDN);
        return mpfr_get_d(r18236, MPFR_RNDN);
}