#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 r17921 = x;
        float r17922 = r17921 / r17921;
        float r17923 = 1.0f;
        float r17924 = r17923 / r17921;
        float r17925 = r17921 * r17921;
        float r17926 = sqrt(r17925);
        float r17927 = r17924 * r17926;
        float r17928 = r17922 - r17927;
        return r17928;
}

double f_id(double x) {
        double r17929 = x;
        double r17930 = r17929 / r17929;
        double r17931 = 1.0;
        double r17932 = r17931 / r17929;
        double r17933 = r17929 * r17929;
        double r17934 = sqrt(r17933);
        double r17935 = r17932 * r17934;
        double r17936 = r17930 - r17935;
        return r17936;
}


double f_of(float x) {
        float r17937 = 1.0f;
        float r17938 = x;
        float r17939 = fabs(r17938);
        float r17940 = r17939 / r17938;
        float r17941 = r17937 - r17940;
        return r17941;
}

double f_od(double x) {
        double r17942 = 1.0;
        double r17943 = x;
        double r17944 = fabs(r17943);
        double r17945 = r17944 / r17943;
        double r17946 = r17942 - r17945;
        return r17946;
}

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 r17947, r17948, r17949, r17950, r17951, r17952, r17953, r17954;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17947);
        mpfr_init(r17948);
        mpfr_init_set_str(r17949, "1", 10, MPFR_RNDN);
        mpfr_init(r17950);
        mpfr_init(r17951);
        mpfr_init(r17952);
        mpfr_init(r17953);
        mpfr_init(r17954);
}

double f_im(double x) {
        mpfr_set_d(r17947, x, MPFR_RNDN);
        mpfr_div(r17948, r17947, r17947, MPFR_RNDN);
        ;
        mpfr_div(r17950, r17949, r17947, MPFR_RNDN);
        mpfr_mul(r17951, r17947, r17947, MPFR_RNDN);
        mpfr_sqrt(r17952, r17951, MPFR_RNDN);
        mpfr_mul(r17953, r17950, r17952, MPFR_RNDN);
        mpfr_sub(r17954, r17948, r17953, MPFR_RNDN);
        return mpfr_get_d(r17954, MPFR_RNDN);
}

static mpfr_t r17955, r17956, r17957, r17958, r17959;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r17955, "1", 10, MPFR_RNDN);
        mpfr_init(r17956);
        mpfr_init(r17957);
        mpfr_init(r17958);
        mpfr_init(r17959);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r17956, x, MPFR_RNDN);
        mpfr_abs(r17957, r17956, MPFR_RNDN);
        mpfr_div(r17958, r17957, r17956, MPFR_RNDN);
        mpfr_sub(r17959, r17955, r17958, MPFR_RNDN);
        return mpfr_get_d(r17959, MPFR_RNDN);
}

static mpfr_t r17960, r17961, r17962, r17963, r17964;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r17960, "1", 10, MPFR_RNDN);
        mpfr_init(r17961);
        mpfr_init(r17962);
        mpfr_init(r17963);
        mpfr_init(r17964);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r17961, x, MPFR_RNDN);
        mpfr_abs(r17962, r17961, MPFR_RNDN);
        mpfr_div(r17963, r17962, r17961, MPFR_RNDN);
        mpfr_sub(r17964, r17960, r17963, MPFR_RNDN);
        return mpfr_get_d(r17964, MPFR_RNDN);
}