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

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

double f_if(float B, float x) {
        float r8965 = x;
        float r8966 = 1;
        float r8967 = B;
        float r8968 = tan(r8967);
        float r8969 = r8966 / r8968;
        float r8970 = r8965 * r8969;
        float r8971 = -r8970;
        float r8972 = sin(r8967);
        float r8973 = r8966 / r8972;
        float r8974 = r8971 + r8973;
        return r8974;
}

double f_id(double B, double x) {
        double r8975 = x;
        double r8976 = 1;
        double r8977 = B;
        double r8978 = tan(r8977);
        double r8979 = r8976 / r8978;
        double r8980 = r8975 * r8979;
        double r8981 = -r8980;
        double r8982 = sin(r8977);
        double r8983 = r8976 / r8982;
        double r8984 = r8981 + r8983;
        return r8984;
}


double f_of(float B, float x) {
        float r8985 = 1;
        float r8986 = B;
        float r8987 = sin(r8986);
        float r8988 = r8985 / r8987;
        float r8989 = x;
        float r8990 = cos(r8986);
        float r8991 = r8989 * r8990;
        float r8992 = r8991 / r8987;
        float r8993 = r8988 - r8992;
        return r8993;
}

double f_od(double B, double x) {
        double r8994 = 1;
        double r8995 = B;
        double r8996 = sin(r8995);
        double r8997 = r8994 / r8996;
        double r8998 = x;
        double r8999 = cos(r8995);
        double r9000 = r8998 * r8999;
        double r9001 = r9000 / r8996;
        double r9002 = r8997 - r9001;
        return r9002;
}

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 r9003, r9004, r9005, r9006, r9007, r9008, r9009, r9010, r9011, r9012;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9003);
        mpfr_init_set_str(r9004, "1", 10, MPFR_RNDN);
        mpfr_init(r9005);
        mpfr_init(r9006);
        mpfr_init(r9007);
        mpfr_init(r9008);
        mpfr_init(r9009);
        mpfr_init(r9010);
        mpfr_init(r9011);
        mpfr_init(r9012);
}

double f_im(double B, double x) {
        mpfr_set_d(r9003, x, MPFR_RNDN);
        ;
        mpfr_set_d(r9005, B, MPFR_RNDN);
        mpfr_tan(r9006, r9005, MPFR_RNDN);
        mpfr_div(r9007, r9004, r9006, MPFR_RNDN);
        mpfr_mul(r9008, r9003, r9007, MPFR_RNDN);
        mpfr_neg(r9009, r9008, MPFR_RNDN);
        mpfr_sin(r9010, r9005, MPFR_RNDN);
        mpfr_div(r9011, r9004, r9010, MPFR_RNDN);
        mpfr_add(r9012, r9009, r9011, MPFR_RNDN);
        return mpfr_get_d(r9012, MPFR_RNDN);
}

static mpfr_t r9013, r9014, r9015, r9016, r9017, r9018, r9019, r9020, r9021;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9013, "1", 10, MPFR_RNDN);
        mpfr_init(r9014);
        mpfr_init(r9015);
        mpfr_init(r9016);
        mpfr_init(r9017);
        mpfr_init(r9018);
        mpfr_init(r9019);
        mpfr_init(r9020);
        mpfr_init(r9021);
}

double f_fm(double B, double x) {
        ;
        mpfr_set_d(r9014, B, MPFR_RNDN);
        mpfr_sin(r9015, r9014, MPFR_RNDN);
        mpfr_div(r9016, r9013, r9015, MPFR_RNDN);
        mpfr_set_d(r9017, x, MPFR_RNDN);
        mpfr_cos(r9018, r9014, MPFR_RNDN);
        mpfr_mul(r9019, r9017, r9018, MPFR_RNDN);
        mpfr_div(r9020, r9019, r9015, MPFR_RNDN);
        mpfr_sub(r9021, r9016, r9020, MPFR_RNDN);
        return mpfr_get_d(r9021, MPFR_RNDN);
}

static mpfr_t r9022, r9023, r9024, r9025, r9026, r9027, r9028, r9029, r9030;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9022, "1", 10, MPFR_RNDN);
        mpfr_init(r9023);
        mpfr_init(r9024);
        mpfr_init(r9025);
        mpfr_init(r9026);
        mpfr_init(r9027);
        mpfr_init(r9028);
        mpfr_init(r9029);
        mpfr_init(r9030);
}

double f_dm(double B, double x) {
        ;
        mpfr_set_d(r9023, B, MPFR_RNDN);
        mpfr_sin(r9024, r9023, MPFR_RNDN);
        mpfr_div(r9025, r9022, r9024, MPFR_RNDN);
        mpfr_set_d(r9026, x, MPFR_RNDN);
        mpfr_cos(r9027, r9023, MPFR_RNDN);
        mpfr_mul(r9028, r9026, r9027, MPFR_RNDN);
        mpfr_div(r9029, r9028, r9024, MPFR_RNDN);
        mpfr_sub(r9030, r9025, r9029, MPFR_RNDN);
        return mpfr_get_d(r9030, MPFR_RNDN);
}