#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 r9348 = a;
        float r9349 = r9348 * r9348;
        float r9350 = b;
        float r9351 = r9350 * r9350;
        float r9352 = r9349 - r9351;
        float r9353 = r9352 / r9349;
        float r9354 = fabs(r9353);
        float r9355 = sqrt(r9354);
        return r9355;
}

double f_id(double a, double b) {
        double r9356 = a;
        double r9357 = r9356 * r9356;
        double r9358 = b;
        double r9359 = r9358 * r9358;
        double r9360 = r9357 - r9359;
        double r9361 = r9360 / r9357;
        double r9362 = fabs(r9361);
        double r9363 = sqrt(r9362);
        return r9363;
}


double f_of(float a, float b) {
        float r9364 = 1;
        float r9365 = b;
        float r9366 = a;
        float r9367 = r9365 / r9366;
        float r9368 = r9367 * r9367;
        float r9369 = exp(r9368);
        float r9370 = log(r9369);
        float r9371 = r9364 - r9370;
        float r9372 = fabs(r9371);
        float r9373 = sqrt(r9372);
        return r9373;
}

double f_od(double a, double b) {
        double r9374 = 1;
        double r9375 = b;
        double r9376 = a;
        double r9377 = r9375 / r9376;
        double r9378 = r9377 * r9377;
        double r9379 = exp(r9378);
        double r9380 = log(r9379);
        double r9381 = r9374 - r9380;
        double r9382 = fabs(r9381);
        double r9383 = sqrt(r9382);
        return r9383;
}

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 r9384, r9385, r9386, r9387, r9388, r9389, r9390, r9391;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r9384);
        mpfr_init(r9385);
        mpfr_init(r9386);
        mpfr_init(r9387);
        mpfr_init(r9388);
        mpfr_init(r9389);
        mpfr_init(r9390);
        mpfr_init(r9391);
}

double f_im(double a, double b) {
        mpfr_set_d(r9384, a, MPFR_RNDN);
        mpfr_mul(r9385, r9384, r9384, MPFR_RNDN);
        mpfr_set_d(r9386, b, MPFR_RNDN);
        mpfr_mul(r9387, r9386, r9386, MPFR_RNDN);
        mpfr_sub(r9388, r9385, r9387, MPFR_RNDN);
        mpfr_div(r9389, r9388, r9385, MPFR_RNDN);
        mpfr_abs(r9390, r9389, MPFR_RNDN);
        mpfr_sqrt(r9391, r9390, MPFR_RNDN);
        return mpfr_get_d(r9391, MPFR_RNDN);
}

static mpfr_t r9392, r9393, r9394, r9395, r9396, r9397, r9398, r9399, r9400, r9401;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r9392, "1", 10, MPFR_RNDN);
        mpfr_init(r9393);
        mpfr_init(r9394);
        mpfr_init(r9395);
        mpfr_init(r9396);
        mpfr_init(r9397);
        mpfr_init(r9398);
        mpfr_init(r9399);
        mpfr_init(r9400);
        mpfr_init(r9401);
}

double f_fm(double a, double b) {
        ;
        mpfr_set_d(r9393, b, MPFR_RNDN);
        mpfr_set_d(r9394, a, MPFR_RNDN);
        mpfr_div(r9395, r9393, r9394, MPFR_RNDN);
        mpfr_mul(r9396, r9395, r9395, MPFR_RNDN);
        mpfr_exp(r9397, r9396, MPFR_RNDN);
        mpfr_log(r9398, r9397, MPFR_RNDN);
        mpfr_sub(r9399, r9392, r9398, MPFR_RNDN);
        mpfr_abs(r9400, r9399, MPFR_RNDN);
        mpfr_sqrt(r9401, r9400, MPFR_RNDN);
        return mpfr_get_d(r9401, MPFR_RNDN);
}

static mpfr_t r9402, r9403, r9404, r9405, r9406, r9407, r9408, r9409, r9410, r9411;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r9402, "1", 10, MPFR_RNDN);
        mpfr_init(r9403);
        mpfr_init(r9404);
        mpfr_init(r9405);
        mpfr_init(r9406);
        mpfr_init(r9407);
        mpfr_init(r9408);
        mpfr_init(r9409);
        mpfr_init(r9410);
        mpfr_init(r9411);
}

double f_dm(double a, double b) {
        ;
        mpfr_set_d(r9403, b, MPFR_RNDN);
        mpfr_set_d(r9404, a, MPFR_RNDN);
        mpfr_div(r9405, r9403, r9404, MPFR_RNDN);
        mpfr_mul(r9406, r9405, r9405, MPFR_RNDN);
        mpfr_exp(r9407, r9406, MPFR_RNDN);
        mpfr_log(r9408, r9407, MPFR_RNDN);
        mpfr_sub(r9409, r9402, r9408, MPFR_RNDN);
        mpfr_abs(r9410, r9409, MPFR_RNDN);
        mpfr_sqrt(r9411, r9410, MPFR_RNDN);
        return mpfr_get_d(r9411, MPFR_RNDN);
}