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

char *name = "Henrywood and Agarwal, Equation (3)";

double f_if(float c0, float A, float V, float l) {
        float r20759 = c0;
        float r20760 = A;
        float r20761 = V;
        float r20762 = l;
        float r20763 = r20761 * r20762;
        float r20764 = r20760 / r20763;
        float r20765 = sqrt(r20764);
        float r20766 = r20759 * r20765;
        return r20766;
}

double f_id(double c0, double A, double V, double l) {
        double r20767 = c0;
        double r20768 = A;
        double r20769 = V;
        double r20770 = l;
        double r20771 = r20769 * r20770;
        double r20772 = r20768 / r20771;
        double r20773 = sqrt(r20772);
        double r20774 = r20767 * r20773;
        return r20774;
}


double f_of(float c0, float A, float V, float l) {
        float r20775 = l;
        float r20776 = 6.538295639751805e-292f;
        bool r20777 = r20775 <= r20776;
        float r20778 = c0;
        float r20779 = A;
        float r20780 = V;
        float r20781 = r20779 / r20780;
        float r20782 = r20781 / r20775;
        float r20783 = sqrt(r20782);
        float r20784 = r20778 * r20783;
        float r20785 = sqrt(r20781);
        float r20786 = 1.0f;
        float r20787 = r20786 / r20775;
        float r20788 = sqrt(r20787);
        float r20789 = r20785 * r20788;
        float r20790 = r20778 * r20789;
        float r20791 = r20777 ? r20784 : r20790;
        return r20791;
}

double f_od(double c0, double A, double V, double l) {
        double r20792 = l;
        double r20793 = 6.538295639751805e-292;
        bool r20794 = r20792 <= r20793;
        double r20795 = c0;
        double r20796 = A;
        double r20797 = V;
        double r20798 = r20796 / r20797;
        double r20799 = r20798 / r20792;
        double r20800 = sqrt(r20799);
        double r20801 = r20795 * r20800;
        double r20802 = sqrt(r20798);
        double r20803 = 1.0;
        double r20804 = r20803 / r20792;
        double r20805 = sqrt(r20804);
        double r20806 = r20802 * r20805;
        double r20807 = r20795 * r20806;
        double r20808 = r20794 ? r20801 : r20807;
        return r20808;
}

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 r20809, r20810, r20811, r20812, r20813, r20814, r20815, r20816;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r20809);
        mpfr_init(r20810);
        mpfr_init(r20811);
        mpfr_init(r20812);
        mpfr_init(r20813);
        mpfr_init(r20814);
        mpfr_init(r20815);
        mpfr_init(r20816);
}

double f_im(double c0, double A, double V, double l) {
        mpfr_set_d(r20809, c0, MPFR_RNDN);
        mpfr_set_d(r20810, A, MPFR_RNDN);
        mpfr_set_d(r20811, V, MPFR_RNDN);
        mpfr_set_d(r20812, l, MPFR_RNDN);
        mpfr_mul(r20813, r20811, r20812, MPFR_RNDN);
        mpfr_div(r20814, r20810, r20813, MPFR_RNDN);
        mpfr_sqrt(r20815, r20814, MPFR_RNDN);
        mpfr_mul(r20816, r20809, r20815, MPFR_RNDN);
        return mpfr_get_d(r20816, MPFR_RNDN);
}

static mpfr_t r20817, r20818, r20819, r20820, r20821, r20822, r20823, r20824, r20825, r20826, r20827, r20828, r20829, r20830, r20831, r20832, r20833;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20817);
        mpfr_init_set_str(r20818, "6.538295639751805e-292", 10, MPFR_RNDN);
        mpfr_init(r20819);
        mpfr_init(r20820);
        mpfr_init(r20821);
        mpfr_init(r20822);
        mpfr_init(r20823);
        mpfr_init(r20824);
        mpfr_init(r20825);
        mpfr_init(r20826);
        mpfr_init(r20827);
        mpfr_init_set_str(r20828, "1", 10, MPFR_RNDN);
        mpfr_init(r20829);
        mpfr_init(r20830);
        mpfr_init(r20831);
        mpfr_init(r20832);
        mpfr_init(r20833);
}

double f_fm(double c0, double A, double V, double l) {
        mpfr_set_d(r20817, l, MPFR_RNDN);
        ;
        mpfr_set_si(r20819, mpfr_cmp(r20817, r20818) <= 0, MPFR_RNDN);
        mpfr_set_d(r20820, c0, MPFR_RNDN);
        mpfr_set_d(r20821, A, MPFR_RNDN);
        mpfr_set_d(r20822, V, MPFR_RNDN);
        mpfr_div(r20823, r20821, r20822, MPFR_RNDN);
        mpfr_div(r20824, r20823, r20817, MPFR_RNDN);
        mpfr_sqrt(r20825, r20824, MPFR_RNDN);
        mpfr_mul(r20826, r20820, r20825, MPFR_RNDN);
        mpfr_sqrt(r20827, r20823, MPFR_RNDN);
        ;
        mpfr_div(r20829, r20828, r20817, MPFR_RNDN);
        mpfr_sqrt(r20830, r20829, MPFR_RNDN);
        mpfr_mul(r20831, r20827, r20830, MPFR_RNDN);
        mpfr_mul(r20832, r20820, r20831, MPFR_RNDN);
        if (mpfr_get_si(r20819, MPFR_RNDN)) { mpfr_set(r20833, r20826, MPFR_RNDN); } else { mpfr_set(r20833, r20832, MPFR_RNDN); };
        return mpfr_get_d(r20833, MPFR_RNDN);
}

static mpfr_t r20834, r20835, r20836, r20837, r20838, r20839, r20840, r20841, r20842, r20843, r20844, r20845, r20846, r20847, r20848, r20849, r20850;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20834);
        mpfr_init_set_str(r20835, "6.538295639751805e-292", 10, MPFR_RNDN);
        mpfr_init(r20836);
        mpfr_init(r20837);
        mpfr_init(r20838);
        mpfr_init(r20839);
        mpfr_init(r20840);
        mpfr_init(r20841);
        mpfr_init(r20842);
        mpfr_init(r20843);
        mpfr_init(r20844);
        mpfr_init_set_str(r20845, "1", 10, MPFR_RNDN);
        mpfr_init(r20846);
        mpfr_init(r20847);
        mpfr_init(r20848);
        mpfr_init(r20849);
        mpfr_init(r20850);
}

double f_dm(double c0, double A, double V, double l) {
        mpfr_set_d(r20834, l, MPFR_RNDN);
        ;
        mpfr_set_si(r20836, mpfr_cmp(r20834, r20835) <= 0, MPFR_RNDN);
        mpfr_set_d(r20837, c0, MPFR_RNDN);
        mpfr_set_d(r20838, A, MPFR_RNDN);
        mpfr_set_d(r20839, V, MPFR_RNDN);
        mpfr_div(r20840, r20838, r20839, MPFR_RNDN);
        mpfr_div(r20841, r20840, r20834, MPFR_RNDN);
        mpfr_sqrt(r20842, r20841, MPFR_RNDN);
        mpfr_mul(r20843, r20837, r20842, MPFR_RNDN);
        mpfr_sqrt(r20844, r20840, MPFR_RNDN);
        ;
        mpfr_div(r20846, r20845, r20834, MPFR_RNDN);
        mpfr_sqrt(r20847, r20846, MPFR_RNDN);
        mpfr_mul(r20848, r20844, r20847, MPFR_RNDN);
        mpfr_mul(r20849, r20837, r20848, MPFR_RNDN);
        if (mpfr_get_si(r20836, MPFR_RNDN)) { mpfr_set(r20850, r20843, MPFR_RNDN); } else { mpfr_set(r20850, r20849, MPFR_RNDN); };
        return mpfr_get_d(r20850, MPFR_RNDN);
}