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

char *name = "Eccentricity of an ellipse";

double f_if(float a, float b) {
        float r9114 = a;
        float r9115 = r9114 * r9114;
        float r9116 = b;
        float r9117 = r9116 * r9116;
        float r9118 = r9115 - r9117;
        float r9119 = r9118 / r9115;
        float r9120 = fabs(r9119);
        float r9121 = sqrt(r9120);
        return r9121;
}

double f_id(double a, double b) {
        double r9122 = a;
        double r9123 = r9122 * r9122;
        double r9124 = b;
        double r9125 = r9124 * r9124;
        double r9126 = r9123 - r9125;
        double r9127 = r9126 / r9123;
        double r9128 = fabs(r9127);
        double r9129 = sqrt(r9128);
        return r9129;
}


double f_of(float a, float b) {
        float r9130 = 1.0f;
        float r9131 = b;
        float r9132 = a;
        float r9133 = r9131 / r9132;
        float r9134 = r9133 * r9133;
        float r9135 = r9134 * r9134;
        float r9136 = r9130 - r9135;
        float r9137 = r9134 + r9130;
        float r9138 = r9136 / r9137;
        float r9139 = fabs(r9138);
        float r9140 = sqrt(r9139);
        return r9140;
}

double f_od(double a, double b) {
        double r9141 = 1.0;
        double r9142 = b;
        double r9143 = a;
        double r9144 = r9142 / r9143;
        double r9145 = r9144 * r9144;
        double r9146 = r9145 * r9145;
        double r9147 = r9141 - r9146;
        double r9148 = r9145 + r9141;
        double r9149 = r9147 / r9148;
        double r9150 = fabs(r9149);
        double r9151 = sqrt(r9150);
        return r9151;
}

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 r9152, r9153, r9154, r9155, r9156, r9157, r9158, r9159;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r9152);
        mpfr_init(r9153);
        mpfr_init(r9154);
        mpfr_init(r9155);
        mpfr_init(r9156);
        mpfr_init(r9157);
        mpfr_init(r9158);
        mpfr_init(r9159);
}

double f_im(double a, double b) {
        mpfr_set_d(r9152, a, MPFR_RNDN);
        mpfr_mul(r9153, r9152, r9152, MPFR_RNDN);
        mpfr_set_d(r9154, b, MPFR_RNDN);
        mpfr_mul(r9155, r9154, r9154, MPFR_RNDN);
        mpfr_sub(r9156, r9153, r9155, MPFR_RNDN);
        mpfr_div(r9157, r9156, r9153, MPFR_RNDN);
        mpfr_abs(r9158, r9157, MPFR_RNDN);
        mpfr_sqrt(r9159, r9158, MPFR_RNDN);
        return mpfr_get_d(r9159, MPFR_RNDN);
}

static mpfr_t r9160, r9161, r9162, r9163, r9164, r9165, r9166, r9167, r9168, r9169, r9170;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r9160, "1", 10, MPFR_RNDN);
        mpfr_init(r9161);
        mpfr_init(r9162);
        mpfr_init(r9163);
        mpfr_init(r9164);
        mpfr_init(r9165);
        mpfr_init(r9166);
        mpfr_init(r9167);
        mpfr_init(r9168);
        mpfr_init(r9169);
        mpfr_init(r9170);
}

double f_fm(double a, double b) {
        ;
        mpfr_set_d(r9161, b, MPFR_RNDN);
        mpfr_set_d(r9162, a, MPFR_RNDN);
        mpfr_div(r9163, r9161, r9162, MPFR_RNDN);
        mpfr_mul(r9164, r9163, r9163, MPFR_RNDN);
        mpfr_mul(r9165, r9164, r9164, MPFR_RNDN);
        mpfr_sub(r9166, r9160, r9165, MPFR_RNDN);
        mpfr_add(r9167, r9164, r9160, MPFR_RNDN);
        mpfr_div(r9168, r9166, r9167, MPFR_RNDN);
        mpfr_abs(r9169, r9168, MPFR_RNDN);
        mpfr_sqrt(r9170, r9169, MPFR_RNDN);
        return mpfr_get_d(r9170, MPFR_RNDN);
}

static mpfr_t r9171, r9172, r9173, r9174, r9175, r9176, r9177, r9178, r9179, r9180, r9181;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r9171, "1", 10, MPFR_RNDN);
        mpfr_init(r9172);
        mpfr_init(r9173);
        mpfr_init(r9174);
        mpfr_init(r9175);
        mpfr_init(r9176);
        mpfr_init(r9177);
        mpfr_init(r9178);
        mpfr_init(r9179);
        mpfr_init(r9180);
        mpfr_init(r9181);
}

double f_dm(double a, double b) {
        ;
        mpfr_set_d(r9172, b, MPFR_RNDN);
        mpfr_set_d(r9173, a, MPFR_RNDN);
        mpfr_div(r9174, r9172, r9173, MPFR_RNDN);
        mpfr_mul(r9175, r9174, r9174, MPFR_RNDN);
        mpfr_mul(r9176, r9175, r9175, MPFR_RNDN);
        mpfr_sub(r9177, r9171, r9176, MPFR_RNDN);
        mpfr_add(r9178, r9175, r9171, MPFR_RNDN);
        mpfr_div(r9179, r9177, r9178, MPFR_RNDN);
        mpfr_abs(r9180, r9179, MPFR_RNDN);
        mpfr_sqrt(r9181, r9180, MPFR_RNDN);
        return mpfr_get_d(r9181, MPFR_RNDN);
}