#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 r8513 = x;
        float r8514 = 1.0f;
        float r8515 = B;
        float r8516 = tan(r8515);
        float r8517 = r8514 / r8516;
        float r8518 = r8513 * r8517;
        float r8519 = -r8518;
        float r8520 = sin(r8515);
        float r8521 = r8514 / r8520;
        float r8522 = r8519 + r8521;
        return r8522;
}

double f_id(double B, double x) {
        double r8523 = x;
        double r8524 = 1.0;
        double r8525 = B;
        double r8526 = tan(r8525);
        double r8527 = r8524 / r8526;
        double r8528 = r8523 * r8527;
        double r8529 = -r8528;
        double r8530 = sin(r8525);
        double r8531 = r8524 / r8530;
        double r8532 = r8529 + r8531;
        return r8532;
}


double f_of(float B, float x) {
        float r8533 = 1.0f;
        float r8534 = B;
        float r8535 = sin(r8534);
        float r8536 = r8533 / r8535;
        float r8537 = tan(r8534);
        float r8538 = x;
        float r8539 = r8537 / r8538;
        float r8540 = r8533 / r8539;
        float r8541 = r8536 - r8540;
        return r8541;
}

double f_od(double B, double x) {
        double r8542 = 1.0;
        double r8543 = B;
        double r8544 = sin(r8543);
        double r8545 = r8542 / r8544;
        double r8546 = tan(r8543);
        double r8547 = x;
        double r8548 = r8546 / r8547;
        double r8549 = r8542 / r8548;
        double r8550 = r8545 - r8549;
        return r8550;
}

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 r8551, r8552, r8553, r8554, r8555, r8556, r8557, r8558, r8559, r8560;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r8551);
        mpfr_init_set_str(r8552, "1", 10, MPFR_RNDN);
        mpfr_init(r8553);
        mpfr_init(r8554);
        mpfr_init(r8555);
        mpfr_init(r8556);
        mpfr_init(r8557);
        mpfr_init(r8558);
        mpfr_init(r8559);
        mpfr_init(r8560);
}

double f_im(double B, double x) {
        mpfr_set_d(r8551, x, MPFR_RNDN);
        ;
        mpfr_set_d(r8553, B, MPFR_RNDN);
        mpfr_tan(r8554, r8553, MPFR_RNDN);
        mpfr_div(r8555, r8552, r8554, MPFR_RNDN);
        mpfr_mul(r8556, r8551, r8555, MPFR_RNDN);
        mpfr_neg(r8557, r8556, MPFR_RNDN);
        mpfr_sin(r8558, r8553, MPFR_RNDN);
        mpfr_div(r8559, r8552, r8558, MPFR_RNDN);
        mpfr_add(r8560, r8557, r8559, MPFR_RNDN);
        return mpfr_get_d(r8560, MPFR_RNDN);
}

static mpfr_t r8561, r8562, r8563, r8564, r8565, r8566, r8567, r8568, r8569;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r8561, "1", 10, MPFR_RNDN);
        mpfr_init(r8562);
        mpfr_init(r8563);
        mpfr_init(r8564);
        mpfr_init(r8565);
        mpfr_init(r8566);
        mpfr_init(r8567);
        mpfr_init(r8568);
        mpfr_init(r8569);
}

double f_fm(double B, double x) {
        ;
        mpfr_set_d(r8562, B, MPFR_RNDN);
        mpfr_sin(r8563, r8562, MPFR_RNDN);
        mpfr_div(r8564, r8561, r8563, MPFR_RNDN);
        mpfr_tan(r8565, r8562, MPFR_RNDN);
        mpfr_set_d(r8566, x, MPFR_RNDN);
        mpfr_div(r8567, r8565, r8566, MPFR_RNDN);
        mpfr_div(r8568, r8561, r8567, MPFR_RNDN);
        mpfr_sub(r8569, r8564, r8568, MPFR_RNDN);
        return mpfr_get_d(r8569, MPFR_RNDN);
}

static mpfr_t r8570, r8571, r8572, r8573, r8574, r8575, r8576, r8577, r8578;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r8570, "1", 10, MPFR_RNDN);
        mpfr_init(r8571);
        mpfr_init(r8572);
        mpfr_init(r8573);
        mpfr_init(r8574);
        mpfr_init(r8575);
        mpfr_init(r8576);
        mpfr_init(r8577);
        mpfr_init(r8578);
}

double f_dm(double B, double x) {
        ;
        mpfr_set_d(r8571, B, MPFR_RNDN);
        mpfr_sin(r8572, r8571, MPFR_RNDN);
        mpfr_div(r8573, r8570, r8572, MPFR_RNDN);
        mpfr_tan(r8574, r8571, MPFR_RNDN);
        mpfr_set_d(r8575, x, MPFR_RNDN);
        mpfr_div(r8576, r8574, r8575, MPFR_RNDN);
        mpfr_div(r8577, r8570, r8576, MPFR_RNDN);
        mpfr_sub(r8578, r8573, r8577, MPFR_RNDN);
        return mpfr_get_d(r8578, MPFR_RNDN);
}