#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 r8493 = x;
        float r8494 = 1.0f;
        float r8495 = B;
        float r8496 = tan(r8495);
        float r8497 = r8494 / r8496;
        float r8498 = r8493 * r8497;
        float r8499 = -r8498;
        float r8500 = sin(r8495);
        float r8501 = r8494 / r8500;
        float r8502 = r8499 + r8501;
        return r8502;
}

double f_id(double B, double x) {
        double r8503 = x;
        double r8504 = 1.0;
        double r8505 = B;
        double r8506 = tan(r8505);
        double r8507 = r8504 / r8506;
        double r8508 = r8503 * r8507;
        double r8509 = -r8508;
        double r8510 = sin(r8505);
        double r8511 = r8504 / r8510;
        double r8512 = r8509 + r8511;
        return r8512;
}


double f_of(float B, float x) {
        float r8513 = 1.0f;
        float r8514 = B;
        float r8515 = sin(r8514);
        float r8516 = r8513 / r8515;
        float r8517 = tan(r8514);
        float r8518 = x;
        float r8519 = r8517 / r8518;
        float r8520 = r8513 / r8519;
        float r8521 = r8516 - r8520;
        return r8521;
}

double f_od(double B, double x) {
        double r8522 = 1.0;
        double r8523 = B;
        double r8524 = sin(r8523);
        double r8525 = r8522 / r8524;
        double r8526 = tan(r8523);
        double r8527 = x;
        double r8528 = r8526 / r8527;
        double r8529 = r8522 / r8528;
        double r8530 = r8525 - r8529;
        return r8530;
}

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 r8531, r8532, r8533, r8534, r8535, r8536, r8537, r8538, r8539, r8540;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r8531);
        mpfr_init_set_str(r8532, "1", 10, MPFR_RNDN);
        mpfr_init(r8533);
        mpfr_init(r8534);
        mpfr_init(r8535);
        mpfr_init(r8536);
        mpfr_init(r8537);
        mpfr_init(r8538);
        mpfr_init(r8539);
        mpfr_init(r8540);
}

double f_im(double B, double x) {
        mpfr_set_d(r8531, x, MPFR_RNDN);
        ;
        mpfr_set_d(r8533, B, MPFR_RNDN);
        mpfr_tan(r8534, r8533, MPFR_RNDN);
        mpfr_div(r8535, r8532, r8534, MPFR_RNDN);
        mpfr_mul(r8536, r8531, r8535, MPFR_RNDN);
        mpfr_neg(r8537, r8536, MPFR_RNDN);
        mpfr_sin(r8538, r8533, MPFR_RNDN);
        mpfr_div(r8539, r8532, r8538, MPFR_RNDN);
        mpfr_add(r8540, r8537, r8539, MPFR_RNDN);
        return mpfr_get_d(r8540, MPFR_RNDN);
}

static mpfr_t r8541, r8542, r8543, r8544, r8545, r8546, r8547, r8548, r8549;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r8541, "1", 10, MPFR_RNDN);
        mpfr_init(r8542);
        mpfr_init(r8543);
        mpfr_init(r8544);
        mpfr_init(r8545);
        mpfr_init(r8546);
        mpfr_init(r8547);
        mpfr_init(r8548);
        mpfr_init(r8549);
}

double f_fm(double B, double x) {
        ;
        mpfr_set_d(r8542, B, MPFR_RNDN);
        mpfr_sin(r8543, r8542, MPFR_RNDN);
        mpfr_div(r8544, r8541, r8543, MPFR_RNDN);
        mpfr_tan(r8545, r8542, MPFR_RNDN);
        mpfr_set_d(r8546, x, MPFR_RNDN);
        mpfr_div(r8547, r8545, r8546, MPFR_RNDN);
        mpfr_div(r8548, r8541, r8547, MPFR_RNDN);
        mpfr_sub(r8549, r8544, r8548, MPFR_RNDN);
        return mpfr_get_d(r8549, MPFR_RNDN);
}

static mpfr_t r8550, r8551, r8552, r8553, r8554, r8555, r8556, r8557, r8558;

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

double f_dm(double B, double x) {
        ;
        mpfr_set_d(r8551, B, MPFR_RNDN);
        mpfr_sin(r8552, r8551, MPFR_RNDN);
        mpfr_div(r8553, r8550, r8552, MPFR_RNDN);
        mpfr_tan(r8554, r8551, MPFR_RNDN);
        mpfr_set_d(r8555, x, MPFR_RNDN);
        mpfr_div(r8556, r8554, r8555, MPFR_RNDN);
        mpfr_div(r8557, r8550, r8556, MPFR_RNDN);
        mpfr_sub(r8558, r8553, r8557, MPFR_RNDN);
        return mpfr_get_d(r8558, MPFR_RNDN);
}