#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "(sqrt (/ -2.839573235346269e-37 (* (* a b) (fabs a))))";

double f_if(float a, float b) {
        float r14695 = -2.839573235346269e-37f;
        float r14696 = a;
        float r14697 = b;
        float r14698 = r14696 * r14697;
        float r14699 = fabs(r14696);
        float r14700 = r14698 * r14699;
        float r14701 = r14695 / r14700;
        float r14702 = sqrt(r14701);
        return r14702;
}

double f_id(double a, double b) {
        double r14703 = -2.839573235346269e-37;
        double r14704 = a;
        double r14705 = b;
        double r14706 = r14704 * r14705;
        double r14707 = fabs(r14704);
        double r14708 = r14706 * r14707;
        double r14709 = r14703 / r14708;
        double r14710 = sqrt(r14709);
        return r14710;
}


double f_of(float a, float b) {
        float r14711 = a;
        float r14712 = b;
        float r14713 = r14711 * r14712;
        float r14714 = fabs(r14711);
        float r14715 = r14713 * r14714;
        float r14716 = -3.2336614114229e-311f;
        bool r14717 = r14715 <= r14716;
        float r14718 = -2.839573235346269e-37f;
        float r14719 = r14718 / r14713;
        float r14720 = sqrt(r14719);
        float r14721 = sqrt(r14714);
        float r14722 = r14720 / r14721;
        float r14723 = r14718 / r14712;
        float r14724 = r14723 / r14711;
        float r14725 = sqrt(r14724);
        float r14726 = r14725 / r14721;
        float r14727 = r14717 ? r14722 : r14726;
        return r14727;
}

double f_od(double a, double b) {
        double r14728 = a;
        double r14729 = b;
        double r14730 = r14728 * r14729;
        double r14731 = fabs(r14728);
        double r14732 = r14730 * r14731;
        double r14733 = -3.2336614114229e-311;
        bool r14734 = r14732 <= r14733;
        double r14735 = -2.839573235346269e-37;
        double r14736 = r14735 / r14730;
        double r14737 = sqrt(r14736);
        double r14738 = sqrt(r14731);
        double r14739 = r14737 / r14738;
        double r14740 = r14735 / r14729;
        double r14741 = r14740 / r14728;
        double r14742 = sqrt(r14741);
        double r14743 = r14742 / r14738;
        double r14744 = r14734 ? r14739 : r14743;
        return r14744;
}

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 r14745, r14746, r14747, r14748, r14749, r14750, r14751, r14752;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r14745, "-2.839573235346269e-37", 10, MPFR_RNDN);
        mpfr_init(r14746);
        mpfr_init(r14747);
        mpfr_init(r14748);
        mpfr_init(r14749);
        mpfr_init(r14750);
        mpfr_init(r14751);
        mpfr_init(r14752);
}

double f_im(double a, double b) {
        ;
        mpfr_set_d(r14746, a, MPFR_RNDN);
        mpfr_set_d(r14747, b, MPFR_RNDN);
        mpfr_mul(r14748, r14746, r14747, MPFR_RNDN);
        mpfr_abs(r14749, r14746, MPFR_RNDN);
        mpfr_mul(r14750, r14748, r14749, MPFR_RNDN);
        mpfr_div(r14751, r14745, r14750, MPFR_RNDN);
        mpfr_sqrt(r14752, r14751, MPFR_RNDN);
        return mpfr_get_d(r14752, MPFR_RNDN);
}

static mpfr_t r14753, r14754, r14755, r14756, r14757, r14758, r14759, r14760, r14761, r14762, r14763, r14764, r14765, r14766, r14767, r14768, r14769;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r14753);
        mpfr_init(r14754);
        mpfr_init(r14755);
        mpfr_init(r14756);
        mpfr_init(r14757);
        mpfr_init_set_str(r14758, "-3.2336614114229e-311", 10, MPFR_RNDN);
        mpfr_init(r14759);
        mpfr_init_set_str(r14760, "-2.839573235346269e-37", 10, MPFR_RNDN);
        mpfr_init(r14761);
        mpfr_init(r14762);
        mpfr_init(r14763);
        mpfr_init(r14764);
        mpfr_init(r14765);
        mpfr_init(r14766);
        mpfr_init(r14767);
        mpfr_init(r14768);
        mpfr_init(r14769);
}

double f_fm(double a, double b) {
        mpfr_set_d(r14753, a, MPFR_RNDN);
        mpfr_set_d(r14754, b, MPFR_RNDN);
        mpfr_mul(r14755, r14753, r14754, MPFR_RNDN);
        mpfr_abs(r14756, r14753, MPFR_RNDN);
        mpfr_mul(r14757, r14755, r14756, MPFR_RNDN);
        ;
        mpfr_set_si(r14759, mpfr_cmp(r14757, r14758) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r14761, r14760, r14755, MPFR_RNDN);
        mpfr_sqrt(r14762, r14761, MPFR_RNDN);
        mpfr_sqrt(r14763, r14756, MPFR_RNDN);
        mpfr_div(r14764, r14762, r14763, MPFR_RNDN);
        mpfr_div(r14765, r14760, r14754, MPFR_RNDN);
        mpfr_div(r14766, r14765, r14753, MPFR_RNDN);
        mpfr_sqrt(r14767, r14766, MPFR_RNDN);
        mpfr_div(r14768, r14767, r14763, MPFR_RNDN);
        if (mpfr_get_si(r14759, MPFR_RNDN)) { mpfr_set(r14769, r14764, MPFR_RNDN); } else { mpfr_set(r14769, r14768, MPFR_RNDN); };
        return mpfr_get_d(r14769, MPFR_RNDN);
}

static mpfr_t r14770, r14771, r14772, r14773, r14774, r14775, r14776, r14777, r14778, r14779, r14780, r14781, r14782, r14783, r14784, r14785, r14786;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r14770);
        mpfr_init(r14771);
        mpfr_init(r14772);
        mpfr_init(r14773);
        mpfr_init(r14774);
        mpfr_init_set_str(r14775, "-3.2336614114229e-311", 10, MPFR_RNDN);
        mpfr_init(r14776);
        mpfr_init_set_str(r14777, "-2.839573235346269e-37", 10, MPFR_RNDN);
        mpfr_init(r14778);
        mpfr_init(r14779);
        mpfr_init(r14780);
        mpfr_init(r14781);
        mpfr_init(r14782);
        mpfr_init(r14783);
        mpfr_init(r14784);
        mpfr_init(r14785);
        mpfr_init(r14786);
}

double f_dm(double a, double b) {
        mpfr_set_d(r14770, a, MPFR_RNDN);
        mpfr_set_d(r14771, b, MPFR_RNDN);
        mpfr_mul(r14772, r14770, r14771, MPFR_RNDN);
        mpfr_abs(r14773, r14770, MPFR_RNDN);
        mpfr_mul(r14774, r14772, r14773, MPFR_RNDN);
        ;
        mpfr_set_si(r14776, mpfr_cmp(r14774, r14775) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r14778, r14777, r14772, MPFR_RNDN);
        mpfr_sqrt(r14779, r14778, MPFR_RNDN);
        mpfr_sqrt(r14780, r14773, MPFR_RNDN);
        mpfr_div(r14781, r14779, r14780, MPFR_RNDN);
        mpfr_div(r14782, r14777, r14771, MPFR_RNDN);
        mpfr_div(r14783, r14782, r14770, MPFR_RNDN);
        mpfr_sqrt(r14784, r14783, MPFR_RNDN);
        mpfr_div(r14785, r14784, r14780, MPFR_RNDN);
        if (mpfr_get_si(r14776, MPFR_RNDN)) { mpfr_set(r14786, r14781, MPFR_RNDN); } else { mpfr_set(r14786, r14785, MPFR_RNDN); };
        return mpfr_get_d(r14786, MPFR_RNDN);
}

