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

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


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

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

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

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

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

static mpfr_t r18229, r18230, r18231, r18232, r18233;

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

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

static mpfr_t r18234, r18235, r18236, r18237, r18238;

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

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