#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 r25561 = a;
        float r25562 = r25561 * r25561;
        float r25563 = b;
        float r25564 = r25563 * r25563;
        float r25565 = r25562 - r25564;
        float r25566 = r25565 / r25562;
        float r25567 = fabs(r25566);
        float r25568 = sqrt(r25567);
        return r25568;
}

double f_id(double a, double b) {
        double r25569 = a;
        double r25570 = r25569 * r25569;
        double r25571 = b;
        double r25572 = r25571 * r25571;
        double r25573 = r25570 - r25572;
        double r25574 = r25573 / r25570;
        double r25575 = fabs(r25574);
        double r25576 = sqrt(r25575);
        return r25576;
}


double f_of(float a, float b) {
        float r25577 = 1;
        float r25578 = b;
        float r25579 = a;
        float r25580 = r25578 / r25579;
        float r25581 = r25580 * r25580;
        float r25582 = r25577 - r25581;
        float r25583 = fabs(r25582);
        float r25584 = sqrt(r25583);
        return r25584;
}

double f_od(double a, double b) {
        double r25585 = 1;
        double r25586 = b;
        double r25587 = a;
        double r25588 = r25586 / r25587;
        double r25589 = r25588 * r25588;
        double r25590 = r25585 - r25589;
        double r25591 = fabs(r25590);
        double r25592 = sqrt(r25591);
        return r25592;
}

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 r25593, r25594, r25595, r25596, r25597, r25598, r25599, r25600;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r25593);
        mpfr_init(r25594);
        mpfr_init(r25595);
        mpfr_init(r25596);
        mpfr_init(r25597);
        mpfr_init(r25598);
        mpfr_init(r25599);
        mpfr_init(r25600);
}

double f_im(double a, double b) {
        mpfr_set_d(r25593, a, MPFR_RNDN);
        mpfr_mul(r25594, r25593, r25593, MPFR_RNDN);
        mpfr_set_d(r25595, b, MPFR_RNDN);
        mpfr_mul(r25596, r25595, r25595, MPFR_RNDN);
        mpfr_sub(r25597, r25594, r25596, MPFR_RNDN);
        mpfr_div(r25598, r25597, r25594, MPFR_RNDN);
        mpfr_abs(r25599, r25598, MPFR_RNDN);
        mpfr_sqrt(r25600, r25599, MPFR_RNDN);
        return mpfr_get_d(r25600, MPFR_RNDN);
}

static mpfr_t r25601, r25602, r25603, r25604, r25605, r25606, r25607, r25608;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r25601, "1", 10, MPFR_RNDN);
        mpfr_init(r25602);
        mpfr_init(r25603);
        mpfr_init(r25604);
        mpfr_init(r25605);
        mpfr_init(r25606);
        mpfr_init(r25607);
        mpfr_init(r25608);
}

double f_fm(double a, double b) {
        ;
        mpfr_set_d(r25602, b, MPFR_RNDN);
        mpfr_set_d(r25603, a, MPFR_RNDN);
        mpfr_div(r25604, r25602, r25603, MPFR_RNDN);
        mpfr_mul(r25605, r25604, r25604, MPFR_RNDN);
        mpfr_sub(r25606, r25601, r25605, MPFR_RNDN);
        mpfr_abs(r25607, r25606, MPFR_RNDN);
        mpfr_sqrt(r25608, r25607, MPFR_RNDN);
        return mpfr_get_d(r25608, MPFR_RNDN);
}

static mpfr_t r25609, r25610, r25611, r25612, r25613, r25614, r25615, r25616;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r25609, "1", 10, MPFR_RNDN);
        mpfr_init(r25610);
        mpfr_init(r25611);
        mpfr_init(r25612);
        mpfr_init(r25613);
        mpfr_init(r25614);
        mpfr_init(r25615);
        mpfr_init(r25616);
}

double f_dm(double a, double b) {
        ;
        mpfr_set_d(r25610, b, MPFR_RNDN);
        mpfr_set_d(r25611, a, MPFR_RNDN);
        mpfr_div(r25612, r25610, r25611, MPFR_RNDN);
        mpfr_mul(r25613, r25612, r25612, MPFR_RNDN);
        mpfr_sub(r25614, r25609, r25613, MPFR_RNDN);
        mpfr_abs(r25615, r25614, MPFR_RNDN);
        mpfr_sqrt(r25616, r25615, MPFR_RNDN);
        return mpfr_get_d(r25616, MPFR_RNDN);
}