#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 r17917 = x;
        float r17918 = r17917 / r17917;
        float r17919 = 1.0f;
        float r17920 = r17919 / r17917;
        float r17921 = r17917 * r17917;
        float r17922 = sqrt(r17921);
        float r17923 = r17920 * r17922;
        float r17924 = r17918 - r17923;
        return r17924;
}

double f_id(double x) {
        double r17925 = x;
        double r17926 = r17925 / r17925;
        double r17927 = 1.0;
        double r17928 = r17927 / r17925;
        double r17929 = r17925 * r17925;
        double r17930 = sqrt(r17929);
        double r17931 = r17928 * r17930;
        double r17932 = r17926 - r17931;
        return r17932;
}


double f_of(float x) {
        float r17933 = 1.0f;
        float r17934 = x;
        float r17935 = fabs(r17934);
        float r17936 = r17935 / r17934;
        float r17937 = r17933 - r17936;
        return r17937;
}

double f_od(double x) {
        double r17938 = 1.0;
        double r17939 = x;
        double r17940 = fabs(r17939);
        double r17941 = r17940 / r17939;
        double r17942 = r17938 - r17941;
        return r17942;
}

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 r17943, r17944, r17945, r17946, r17947, r17948, r17949, r17950;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17943);
        mpfr_init(r17944);
        mpfr_init_set_str(r17945, "1", 10, MPFR_RNDN);
        mpfr_init(r17946);
        mpfr_init(r17947);
        mpfr_init(r17948);
        mpfr_init(r17949);
        mpfr_init(r17950);
}

double f_im(double x) {
        mpfr_set_d(r17943, x, MPFR_RNDN);
        mpfr_div(r17944, r17943, r17943, MPFR_RNDN);
        ;
        mpfr_div(r17946, r17945, r17943, MPFR_RNDN);
        mpfr_mul(r17947, r17943, r17943, MPFR_RNDN);
        mpfr_sqrt(r17948, r17947, MPFR_RNDN);
        mpfr_mul(r17949, r17946, r17948, MPFR_RNDN);
        mpfr_sub(r17950, r17944, r17949, MPFR_RNDN);
        return mpfr_get_d(r17950, MPFR_RNDN);
}

static mpfr_t r17951, r17952, r17953, r17954, r17955;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r17951, "1", 10, MPFR_RNDN);
        mpfr_init(r17952);
        mpfr_init(r17953);
        mpfr_init(r17954);
        mpfr_init(r17955);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r17952, x, MPFR_RNDN);
        mpfr_abs(r17953, r17952, MPFR_RNDN);
        mpfr_div(r17954, r17953, r17952, MPFR_RNDN);
        mpfr_sub(r17955, r17951, r17954, MPFR_RNDN);
        return mpfr_get_d(r17955, MPFR_RNDN);
}

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

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

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