#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 r25492 = a;
        float r25493 = r25492 * r25492;
        float r25494 = b;
        float r25495 = r25494 * r25494;
        float r25496 = r25493 - r25495;
        float r25497 = r25496 / r25493;
        float r25498 = fabs(r25497);
        float r25499 = sqrt(r25498);
        return r25499;
}

double f_id(double a, double b) {
        double r25500 = a;
        double r25501 = r25500 * r25500;
        double r25502 = b;
        double r25503 = r25502 * r25502;
        double r25504 = r25501 - r25503;
        double r25505 = r25504 / r25501;
        double r25506 = fabs(r25505);
        double r25507 = sqrt(r25506);
        return r25507;
}


double f_of(float a, float b) {
        float r25508 = 1;
        float r25509 = b;
        float r25510 = a;
        float r25511 = r25509 / r25510;
        float r25512 = r25511 * r25511;
        float r25513 = cbrt(r25512);
        float r25514 = r25513 * r25513;
        float r25515 = r25514 * r25513;
        float r25516 = r25508 - r25515;
        float r25517 = fabs(r25516);
        float r25518 = sqrt(r25517);
        return r25518;
}

double f_od(double a, double b) {
        double r25519 = 1;
        double r25520 = b;
        double r25521 = a;
        double r25522 = r25520 / r25521;
        double r25523 = r25522 * r25522;
        double r25524 = cbrt(r25523);
        double r25525 = r25524 * r25524;
        double r25526 = r25525 * r25524;
        double r25527 = r25519 - r25526;
        double r25528 = fabs(r25527);
        double r25529 = sqrt(r25528);
        return r25529;
}

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 r25530, r25531, r25532, r25533, r25534, r25535, r25536, r25537;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r25530);
        mpfr_init(r25531);
        mpfr_init(r25532);
        mpfr_init(r25533);
        mpfr_init(r25534);
        mpfr_init(r25535);
        mpfr_init(r25536);
        mpfr_init(r25537);
}

double f_im(double a, double b) {
        mpfr_set_d(r25530, a, MPFR_RNDN);
        mpfr_mul(r25531, r25530, r25530, MPFR_RNDN);
        mpfr_set_d(r25532, b, MPFR_RNDN);
        mpfr_mul(r25533, r25532, r25532, MPFR_RNDN);
        mpfr_sub(r25534, r25531, r25533, MPFR_RNDN);
        mpfr_div(r25535, r25534, r25531, MPFR_RNDN);
        mpfr_abs(r25536, r25535, MPFR_RNDN);
        mpfr_sqrt(r25537, r25536, MPFR_RNDN);
        return mpfr_get_d(r25537, MPFR_RNDN);
}

static mpfr_t r25538, r25539, r25540, r25541, r25542, r25543, r25544, r25545, r25546, r25547, r25548;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r25538, "1", 10, MPFR_RNDN);
        mpfr_init(r25539);
        mpfr_init(r25540);
        mpfr_init(r25541);
        mpfr_init(r25542);
        mpfr_init(r25543);
        mpfr_init(r25544);
        mpfr_init(r25545);
        mpfr_init(r25546);
        mpfr_init(r25547);
        mpfr_init(r25548);
}

double f_fm(double a, double b) {
        ;
        mpfr_set_d(r25539, b, MPFR_RNDN);
        mpfr_set_d(r25540, a, MPFR_RNDN);
        mpfr_div(r25541, r25539, r25540, MPFR_RNDN);
        mpfr_mul(r25542, r25541, r25541, MPFR_RNDN);
        mpfr_cbrt(r25543, r25542, MPFR_RNDN);
        mpfr_mul(r25544, r25543, r25543, MPFR_RNDN);
        mpfr_mul(r25545, r25544, r25543, MPFR_RNDN);
        mpfr_sub(r25546, r25538, r25545, MPFR_RNDN);
        mpfr_abs(r25547, r25546, MPFR_RNDN);
        mpfr_sqrt(r25548, r25547, MPFR_RNDN);
        return mpfr_get_d(r25548, MPFR_RNDN);
}

static mpfr_t r25549, r25550, r25551, r25552, r25553, r25554, r25555, r25556, r25557, r25558, r25559;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r25549, "1", 10, MPFR_RNDN);
        mpfr_init(r25550);
        mpfr_init(r25551);
        mpfr_init(r25552);
        mpfr_init(r25553);
        mpfr_init(r25554);
        mpfr_init(r25555);
        mpfr_init(r25556);
        mpfr_init(r25557);
        mpfr_init(r25558);
        mpfr_init(r25559);
}

double f_dm(double a, double b) {
        ;
        mpfr_set_d(r25550, b, MPFR_RNDN);
        mpfr_set_d(r25551, a, MPFR_RNDN);
        mpfr_div(r25552, r25550, r25551, MPFR_RNDN);
        mpfr_mul(r25553, r25552, r25552, MPFR_RNDN);
        mpfr_cbrt(r25554, r25553, MPFR_RNDN);
        mpfr_mul(r25555, r25554, r25554, MPFR_RNDN);
        mpfr_mul(r25556, r25555, r25554, MPFR_RNDN);
        mpfr_sub(r25557, r25549, r25556, MPFR_RNDN);
        mpfr_abs(r25558, r25557, MPFR_RNDN);
        mpfr_sqrt(r25559, r25558, MPFR_RNDN);
        return mpfr_get_d(r25559, MPFR_RNDN);
}