#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 r7786 = atan2(1.0, 0.0);
        float r7787 = l;
        float r7788 = r7786 * r7787;
        float r7789 = 1;
        float r7790 = F;
        float r7791 = r7790 * r7790;
        float r7792 = r7789 / r7791;
        float r7793 = tan(r7788);
        float r7794 = r7792 * r7793;
        float r7795 = r7788 - r7794;
        return r7795;
}

double f_id(double F, double l) {
        double r7796 = atan2(1.0, 0.0);
        double r7797 = l;
        double r7798 = r7796 * r7797;
        double r7799 = 1;
        double r7800 = F;
        double r7801 = r7800 * r7800;
        double r7802 = r7799 / r7801;
        double r7803 = tan(r7798);
        double r7804 = r7802 * r7803;
        double r7805 = r7798 - r7804;
        return r7805;
}


double f_of(float F, float l) {
        float r7806 = atan2(1.0, 0.0);
        float r7807 = l;
        float r7808 = r7806 * r7807;
        float r7809 = 1;
        float r7810 = F;
        float r7811 = r7809 / r7810;
        float r7812 = r7810 / r7807;
        float r7813 = r7812 / r7806;
        float r7814 = 1/3;
        float r7815 = r7806 * r7814;
        float r7816 = 3;
        float r7817 = pow(r7815, r7816);
        float r7818 = cbrt(r7817);
        float r7819 = r7807 * r7810;
        float r7820 = r7818 * r7819;
        float r7821 = r7813 - r7820;
        float r7822 = r7811 / r7821;
        float r7823 = r7808 - r7822;
        return r7823;
}

double f_od(double F, double l) {
        double r7824 = atan2(1.0, 0.0);
        double r7825 = l;
        double r7826 = r7824 * r7825;
        double r7827 = 1;
        double r7828 = F;
        double r7829 = r7827 / r7828;
        double r7830 = r7828 / r7825;
        double r7831 = r7830 / r7824;
        double r7832 = 1/3;
        double r7833 = r7824 * r7832;
        double r7834 = 3;
        double r7835 = pow(r7833, r7834);
        double r7836 = cbrt(r7835);
        double r7837 = r7825 * r7828;
        double r7838 = r7836 * r7837;
        double r7839 = r7831 - r7838;
        double r7840 = r7829 / r7839;
        double r7841 = r7826 - r7840;
        return r7841;
}

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 r7842, r7843, r7844, r7845, r7846, r7847, r7848, r7849, r7850, r7851;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3152);
        mpfr_init(r7842);
        mpfr_init(r7843);
        mpfr_init(r7844);
        mpfr_init_set_str(r7845, "1", 10, MPFR_RNDN);
        mpfr_init(r7846);
        mpfr_init(r7847);
        mpfr_init(r7848);
        mpfr_init(r7849);
        mpfr_init(r7850);
        mpfr_init(r7851);
}

double f_im(double F, double l) {
        mpfr_const_pi(r7842, MPFR_RNDN);
        mpfr_set_d(r7843, l, MPFR_RNDN);
        mpfr_mul(r7844, r7842, r7843, MPFR_RNDN);
        ;
        mpfr_set_d(r7846, F, MPFR_RNDN);
        mpfr_mul(r7847, r7846, r7846, MPFR_RNDN);
        mpfr_div(r7848, r7845, r7847, MPFR_RNDN);
        mpfr_tan(r7849, r7844, MPFR_RNDN);
        mpfr_mul(r7850, r7848, r7849, MPFR_RNDN);
        mpfr_sub(r7851, r7844, r7850, MPFR_RNDN);
        return mpfr_get_d(r7851, MPFR_RNDN);
}

static mpfr_t r7852, r7853, r7854, r7855, r7856, r7857, r7858, r7859, r7860, r7861, r7862, r7863, r7864, r7865, r7866, r7867, r7868, r7869;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3152);
        mpfr_init(r7852);
        mpfr_init(r7853);
        mpfr_init(r7854);
        mpfr_init_set_str(r7855, "1", 10, MPFR_RNDN);
        mpfr_init(r7856);
        mpfr_init(r7857);
        mpfr_init(r7858);
        mpfr_init(r7859);
        mpfr_init_set_str(r7860, "1/3", 10, MPFR_RNDN);
        mpfr_init(r7861);
        mpfr_init_set_str(r7862, "3", 10, MPFR_RNDN);
        mpfr_init(r7863);
        mpfr_init(r7864);
        mpfr_init(r7865);
        mpfr_init(r7866);
        mpfr_init(r7867);
        mpfr_init(r7868);
        mpfr_init(r7869);
}

double f_fm(double F, double l) {
        mpfr_const_pi(r7852, MPFR_RNDN);
        mpfr_set_d(r7853, l, MPFR_RNDN);
        mpfr_mul(r7854, r7852, r7853, MPFR_RNDN);
        ;
        mpfr_set_d(r7856, F, MPFR_RNDN);
        mpfr_div(r7857, r7855, r7856, MPFR_RNDN);
        mpfr_div(r7858, r7856, r7853, MPFR_RNDN);
        mpfr_div(r7859, r7858, r7852, MPFR_RNDN);
        ;
        mpfr_mul(r7861, r7852, r7860, MPFR_RNDN);
        ;
        mpfr_pow(r7863, r7861, r7862, MPFR_RNDN);
        mpfr_cbrt(r7864, r7863, MPFR_RNDN);
        mpfr_mul(r7865, r7853, r7856, MPFR_RNDN);
        mpfr_mul(r7866, r7864, r7865, MPFR_RNDN);
        mpfr_sub(r7867, r7859, r7866, MPFR_RNDN);
        mpfr_div(r7868, r7857, r7867, MPFR_RNDN);
        mpfr_sub(r7869, r7854, r7868, MPFR_RNDN);
        return mpfr_get_d(r7869, MPFR_RNDN);
}

static mpfr_t r7870, r7871, r7872, r7873, r7874, r7875, r7876, r7877, r7878, r7879, r7880, r7881, r7882, r7883, r7884, r7885, r7886, r7887;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3152);
        mpfr_init(r7870);
        mpfr_init(r7871);
        mpfr_init(r7872);
        mpfr_init_set_str(r7873, "1", 10, MPFR_RNDN);
        mpfr_init(r7874);
        mpfr_init(r7875);
        mpfr_init(r7876);
        mpfr_init(r7877);
        mpfr_init_set_str(r7878, "1/3", 10, MPFR_RNDN);
        mpfr_init(r7879);
        mpfr_init_set_str(r7880, "3", 10, MPFR_RNDN);
        mpfr_init(r7881);
        mpfr_init(r7882);
        mpfr_init(r7883);
        mpfr_init(r7884);
        mpfr_init(r7885);
        mpfr_init(r7886);
        mpfr_init(r7887);
}

double f_dm(double F, double l) {
        mpfr_const_pi(r7870, MPFR_RNDN);
        mpfr_set_d(r7871, l, MPFR_RNDN);
        mpfr_mul(r7872, r7870, r7871, MPFR_RNDN);
        ;
        mpfr_set_d(r7874, F, MPFR_RNDN);
        mpfr_div(r7875, r7873, r7874, MPFR_RNDN);
        mpfr_div(r7876, r7874, r7871, MPFR_RNDN);
        mpfr_div(r7877, r7876, r7870, MPFR_RNDN);
        ;
        mpfr_mul(r7879, r7870, r7878, MPFR_RNDN);
        ;
        mpfr_pow(r7881, r7879, r7880, MPFR_RNDN);
        mpfr_cbrt(r7882, r7881, MPFR_RNDN);
        mpfr_mul(r7883, r7871, r7874, MPFR_RNDN);
        mpfr_mul(r7884, r7882, r7883, MPFR_RNDN);
        mpfr_sub(r7885, r7877, r7884, MPFR_RNDN);
        mpfr_div(r7886, r7875, r7885, MPFR_RNDN);
        mpfr_sub(r7887, r7872, r7886, MPFR_RNDN);
        return mpfr_get_d(r7887, MPFR_RNDN);
}