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

char *name = "VandenBroeck and Keller, Equation (6)";

double f_if(float F, float l) {
        float r20925 = atan2(1.0, 0.0);
        float r20926 = l;
        float r20927 = r20925 * r20926;
        float r20928 = 1;
        float r20929 = F;
        float r20930 = r20929 * r20929;
        float r20931 = r20928 / r20930;
        float r20932 = tan(r20927);
        float r20933 = r20931 * r20932;
        float r20934 = r20927 - r20933;
        return r20934;
}

double f_id(double F, double l) {
        double r20935 = atan2(1.0, 0.0);
        double r20936 = l;
        double r20937 = r20935 * r20936;
        double r20938 = 1;
        double r20939 = F;
        double r20940 = r20939 * r20939;
        double r20941 = r20938 / r20940;
        double r20942 = tan(r20937);
        double r20943 = r20941 * r20942;
        double r20944 = r20937 - r20943;
        return r20944;
}


double f_of(float F, float l) {
        float r20945 = atan2(1.0, 0.0);
        float r20946 = l;
        float r20947 = r20945 * r20946;
        float r20948 = -1.2391563621069606e+121;
        bool r20949 = r20947 <= r20948;
        float r20950 = tan(r20947);
        float r20951 = F;
        float r20952 = r20950 / r20951;
        float r20953 = cbrt(r20952);
        float r20954 = r20953 * r20953;
        float r20955 = cbrt(r20954);
        float r20956 = r20955 * r20954;
        float r20957 = r20945 / r20946;
        float r20958 = sin(r20957);
        float r20959 = cos(r20957);
        float r20960 = r20951 / r20959;
        float r20961 = r20958 * r20960;
        float r20962 = cbrt(r20961);
        float r20963 = cbrt(r20962);
        float r20964 = r20951 / r20963;
        float r20965 = r20956 / r20964;
        float r20966 = r20947 - r20965;
        float r20967 = 2.8029951490048387e+50;
        bool r20968 = r20947 <= r20967;
        float r20969 = r20946 * r20945;
        float r20970 = 1;
        float r20971 = r20970 / r20951;
        float r20972 = tan(r20969);
        float r20973 = r20972 / r20951;
        float r20974 = r20971 * r20973;
        float r20975 = r20969 - r20974;
        float r20976 = r20968 ? r20975 : r20966;
        float r20977 = r20949 ? r20966 : r20976;
        return r20977;
}

double f_od(double F, double l) {
        double r20978 = atan2(1.0, 0.0);
        double r20979 = l;
        double r20980 = r20978 * r20979;
        double r20981 = -1.2391563621069606e+121;
        bool r20982 = r20980 <= r20981;
        double r20983 = tan(r20980);
        double r20984 = F;
        double r20985 = r20983 / r20984;
        double r20986 = cbrt(r20985);
        double r20987 = r20986 * r20986;
        double r20988 = cbrt(r20987);
        double r20989 = r20988 * r20987;
        double r20990 = r20978 / r20979;
        double r20991 = sin(r20990);
        double r20992 = cos(r20990);
        double r20993 = r20984 / r20992;
        double r20994 = r20991 * r20993;
        double r20995 = cbrt(r20994);
        double r20996 = cbrt(r20995);
        double r20997 = r20984 / r20996;
        double r20998 = r20989 / r20997;
        double r20999 = r20980 - r20998;
        double r21000 = 2.8029951490048387e+50;
        bool r21001 = r20980 <= r21000;
        double r21002 = r20979 * r20978;
        double r21003 = 1;
        double r21004 = r21003 / r20984;
        double r21005 = tan(r21002);
        double r21006 = r21005 / r20984;
        double r21007 = r21004 * r21006;
        double r21008 = r21002 - r21007;
        double r21009 = r21001 ? r21008 : r20999;
        double r21010 = r20982 ? r20999 : r21009;
        return r21010;
}

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 r21011, r21012, r21013, r21014, r21015, r21016, r21017, r21018, r21019, r21020;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3984);
        mpfr_init(r21011);
        mpfr_init(r21012);
        mpfr_init(r21013);
        mpfr_init_set_str(r21014, "1", 10, MPFR_RNDN);
        mpfr_init(r21015);
        mpfr_init(r21016);
        mpfr_init(r21017);
        mpfr_init(r21018);
        mpfr_init(r21019);
        mpfr_init(r21020);
}

double f_im(double F, double l) {
        mpfr_const_pi(r21011, MPFR_RNDN);
        mpfr_set_d(r21012, l, MPFR_RNDN);
        mpfr_mul(r21013, r21011, r21012, MPFR_RNDN);
        ;
        mpfr_set_d(r21015, F, MPFR_RNDN);
        mpfr_mul(r21016, r21015, r21015, MPFR_RNDN);
        mpfr_div(r21017, r21014, r21016, MPFR_RNDN);
        mpfr_tan(r21018, r21013, MPFR_RNDN);
        mpfr_mul(r21019, r21017, r21018, MPFR_RNDN);
        mpfr_sub(r21020, r21013, r21019, MPFR_RNDN);
        return mpfr_get_d(r21020, MPFR_RNDN);
}

static mpfr_t r21021, r21022, r21023, r21024, r21025, r21026, r21027, r21028, r21029, r21030, r21031, r21032, r21033, r21034, r21035, r21036, r21037, r21038, r21039, r21040, r21041, r21042, r21043, r21044, r21045, r21046, r21047, r21048, r21049, r21050, r21051, r21052, r21053;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3984);
        mpfr_init(r21021);
        mpfr_init(r21022);
        mpfr_init(r21023);
        mpfr_init_set_str(r21024, "-1.2391563621069606e+121", 10, MPFR_RNDN);
        mpfr_init(r21025);
        mpfr_init(r21026);
        mpfr_init(r21027);
        mpfr_init(r21028);
        mpfr_init(r21029);
        mpfr_init(r21030);
        mpfr_init(r21031);
        mpfr_init(r21032);
        mpfr_init(r21033);
        mpfr_init(r21034);
        mpfr_init(r21035);
        mpfr_init(r21036);
        mpfr_init(r21037);
        mpfr_init(r21038);
        mpfr_init(r21039);
        mpfr_init(r21040);
        mpfr_init(r21041);
        mpfr_init(r21042);
        mpfr_init_set_str(r21043, "2.8029951490048387e+50", 10, MPFR_RNDN);
        mpfr_init(r21044);
        mpfr_init(r21045);
        mpfr_init_set_str(r21046, "1", 10, MPFR_RNDN);
        mpfr_init(r21047);
        mpfr_init(r21048);
        mpfr_init(r21049);
        mpfr_init(r21050);
        mpfr_init(r21051);
        mpfr_init(r21052);
        mpfr_init(r21053);
}

double f_fm(double F, double l) {
        mpfr_const_pi(r21021, MPFR_RNDN);
        mpfr_set_d(r21022, l, MPFR_RNDN);
        mpfr_mul(r21023, r21021, r21022, MPFR_RNDN);
        ;
        mpfr_set_si(r21025, mpfr_cmp(r21023, r21024) <= 0, MPFR_RNDN);
        mpfr_tan(r21026, r21023, MPFR_RNDN);
        mpfr_set_d(r21027, F, MPFR_RNDN);
        mpfr_div(r21028, r21026, r21027, MPFR_RNDN);
        mpfr_cbrt(r21029, r21028, MPFR_RNDN);
        mpfr_mul(r21030, r21029, r21029, MPFR_RNDN);
        mpfr_cbrt(r21031, r21030, MPFR_RNDN);
        mpfr_mul(r21032, r21031, r21030, MPFR_RNDN);
        mpfr_div(r21033, r21021, r21022, MPFR_RNDN);
        mpfr_sin(r21034, r21033, MPFR_RNDN);
        mpfr_cos(r21035, r21033, MPFR_RNDN);
        mpfr_div(r21036, r21027, r21035, MPFR_RNDN);
        mpfr_mul(r21037, r21034, r21036, MPFR_RNDN);
        mpfr_cbrt(r21038, r21037, MPFR_RNDN);
        mpfr_cbrt(r21039, r21038, MPFR_RNDN);
        mpfr_div(r21040, r21027, r21039, MPFR_RNDN);
        mpfr_div(r21041, r21032, r21040, MPFR_RNDN);
        mpfr_sub(r21042, r21023, r21041, MPFR_RNDN);
        ;
        mpfr_set_si(r21044, mpfr_cmp(r21023, r21043) <= 0, MPFR_RNDN);
        mpfr_mul(r21045, r21022, r21021, MPFR_RNDN);
        ;
        mpfr_div(r21047, r21046, r21027, MPFR_RNDN);
        mpfr_tan(r21048, r21045, MPFR_RNDN);
        mpfr_div(r21049, r21048, r21027, MPFR_RNDN);
        mpfr_mul(r21050, r21047, r21049, MPFR_RNDN);
        mpfr_sub(r21051, r21045, r21050, MPFR_RNDN);
        if (mpfr_get_si(r21044, MPFR_RNDN)) { mpfr_set(r21052, r21051, MPFR_RNDN); } else { mpfr_set(r21052, r21042, MPFR_RNDN); };
        if (mpfr_get_si(r21025, MPFR_RNDN)) { mpfr_set(r21053, r21042, MPFR_RNDN); } else { mpfr_set(r21053, r21052, MPFR_RNDN); };
        return mpfr_get_d(r21053, MPFR_RNDN);
}

static mpfr_t r21054, r21055, r21056, r21057, r21058, r21059, r21060, r21061, r21062, r21063, r21064, r21065, r21066, r21067, r21068, r21069, r21070, r21071, r21072, r21073, r21074, r21075, r21076, r21077, r21078, r21079, r21080, r21081, r21082, r21083, r21084, r21085, r21086;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3984);
        mpfr_init(r21054);
        mpfr_init(r21055);
        mpfr_init(r21056);
        mpfr_init_set_str(r21057, "-1.2391563621069606e+121", 10, MPFR_RNDN);
        mpfr_init(r21058);
        mpfr_init(r21059);
        mpfr_init(r21060);
        mpfr_init(r21061);
        mpfr_init(r21062);
        mpfr_init(r21063);
        mpfr_init(r21064);
        mpfr_init(r21065);
        mpfr_init(r21066);
        mpfr_init(r21067);
        mpfr_init(r21068);
        mpfr_init(r21069);
        mpfr_init(r21070);
        mpfr_init(r21071);
        mpfr_init(r21072);
        mpfr_init(r21073);
        mpfr_init(r21074);
        mpfr_init(r21075);
        mpfr_init_set_str(r21076, "2.8029951490048387e+50", 10, MPFR_RNDN);
        mpfr_init(r21077);
        mpfr_init(r21078);
        mpfr_init_set_str(r21079, "1", 10, MPFR_RNDN);
        mpfr_init(r21080);
        mpfr_init(r21081);
        mpfr_init(r21082);
        mpfr_init(r21083);
        mpfr_init(r21084);
        mpfr_init(r21085);
        mpfr_init(r21086);
}

double f_dm(double F, double l) {
        mpfr_const_pi(r21054, MPFR_RNDN);
        mpfr_set_d(r21055, l, MPFR_RNDN);
        mpfr_mul(r21056, r21054, r21055, MPFR_RNDN);
        ;
        mpfr_set_si(r21058, mpfr_cmp(r21056, r21057) <= 0, MPFR_RNDN);
        mpfr_tan(r21059, r21056, MPFR_RNDN);
        mpfr_set_d(r21060, F, MPFR_RNDN);
        mpfr_div(r21061, r21059, r21060, MPFR_RNDN);
        mpfr_cbrt(r21062, r21061, MPFR_RNDN);
        mpfr_mul(r21063, r21062, r21062, MPFR_RNDN);
        mpfr_cbrt(r21064, r21063, MPFR_RNDN);
        mpfr_mul(r21065, r21064, r21063, MPFR_RNDN);
        mpfr_div(r21066, r21054, r21055, MPFR_RNDN);
        mpfr_sin(r21067, r21066, MPFR_RNDN);
        mpfr_cos(r21068, r21066, MPFR_RNDN);
        mpfr_div(r21069, r21060, r21068, MPFR_RNDN);
        mpfr_mul(r21070, r21067, r21069, MPFR_RNDN);
        mpfr_cbrt(r21071, r21070, MPFR_RNDN);
        mpfr_cbrt(r21072, r21071, MPFR_RNDN);
        mpfr_div(r21073, r21060, r21072, MPFR_RNDN);
        mpfr_div(r21074, r21065, r21073, MPFR_RNDN);
        mpfr_sub(r21075, r21056, r21074, MPFR_RNDN);
        ;
        mpfr_set_si(r21077, mpfr_cmp(r21056, r21076) <= 0, MPFR_RNDN);
        mpfr_mul(r21078, r21055, r21054, MPFR_RNDN);
        ;
        mpfr_div(r21080, r21079, r21060, MPFR_RNDN);
        mpfr_tan(r21081, r21078, MPFR_RNDN);
        mpfr_div(r21082, r21081, r21060, MPFR_RNDN);
        mpfr_mul(r21083, r21080, r21082, MPFR_RNDN);
        mpfr_sub(r21084, r21078, r21083, MPFR_RNDN);
        if (mpfr_get_si(r21077, MPFR_RNDN)) { mpfr_set(r21085, r21084, MPFR_RNDN); } else { mpfr_set(r21085, r21075, MPFR_RNDN); };
        if (mpfr_get_si(r21058, MPFR_RNDN)) { mpfr_set(r21086, r21075, MPFR_RNDN); } else { mpfr_set(r21086, r21085, MPFR_RNDN); };
        return mpfr_get_d(r21086, MPFR_RNDN);
}