#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 r9111 = x;
        float r9112 = 1;
        float r9113 = B;
        float r9114 = tan(r9113);
        float r9115 = r9112 / r9114;
        float r9116 = r9111 * r9115;
        float r9117 = -r9116;
        float r9118 = sin(r9113);
        float r9119 = r9112 / r9118;
        float r9120 = r9117 + r9119;
        return r9120;
}

double f_id(double B, double x) {
        double r9121 = x;
        double r9122 = 1;
        double r9123 = B;
        double r9124 = tan(r9123);
        double r9125 = r9122 / r9124;
        double r9126 = r9121 * r9125;
        double r9127 = -r9126;
        double r9128 = sin(r9123);
        double r9129 = r9122 / r9128;
        double r9130 = r9127 + r9129;
        return r9130;
}


double f_of(float B, float x) {
        float r9131 = 1;
        float r9132 = B;
        float r9133 = sin(r9132);
        float r9134 = r9131 / r9133;
        float r9135 = x;
        float r9136 = r9135 / r9133;
        float r9137 = cos(r9132);
        float r9138 = r9136 * r9137;
        float r9139 = r9134 - r9138;
        return r9139;
}

double f_od(double B, double x) {
        double r9140 = 1;
        double r9141 = B;
        double r9142 = sin(r9141);
        double r9143 = r9140 / r9142;
        double r9144 = x;
        double r9145 = r9144 / r9142;
        double r9146 = cos(r9141);
        double r9147 = r9145 * r9146;
        double r9148 = r9143 - r9147;
        return r9148;
}

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 r9149, r9150, r9151, r9152, r9153, r9154, r9155, r9156, r9157, r9158;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9149);
        mpfr_init_set_str(r9150, "1", 10, MPFR_RNDN);
        mpfr_init(r9151);
        mpfr_init(r9152);
        mpfr_init(r9153);
        mpfr_init(r9154);
        mpfr_init(r9155);
        mpfr_init(r9156);
        mpfr_init(r9157);
        mpfr_init(r9158);
}

double f_im(double B, double x) {
        mpfr_set_d(r9149, x, MPFR_RNDN);
        ;
        mpfr_set_d(r9151, B, MPFR_RNDN);
        mpfr_tan(r9152, r9151, MPFR_RNDN);
        mpfr_div(r9153, r9150, r9152, MPFR_RNDN);
        mpfr_mul(r9154, r9149, r9153, MPFR_RNDN);
        mpfr_neg(r9155, r9154, MPFR_RNDN);
        mpfr_sin(r9156, r9151, MPFR_RNDN);
        mpfr_div(r9157, r9150, r9156, MPFR_RNDN);
        mpfr_add(r9158, r9155, r9157, MPFR_RNDN);
        return mpfr_get_d(r9158, MPFR_RNDN);
}

static mpfr_t r9159, r9160, r9161, r9162, r9163, r9164, r9165, r9166, r9167;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9159, "1", 10, MPFR_RNDN);
        mpfr_init(r9160);
        mpfr_init(r9161);
        mpfr_init(r9162);
        mpfr_init(r9163);
        mpfr_init(r9164);
        mpfr_init(r9165);
        mpfr_init(r9166);
        mpfr_init(r9167);
}

double f_fm(double B, double x) {
        ;
        mpfr_set_d(r9160, B, MPFR_RNDN);
        mpfr_sin(r9161, r9160, MPFR_RNDN);
        mpfr_div(r9162, r9159, r9161, MPFR_RNDN);
        mpfr_set_d(r9163, x, MPFR_RNDN);
        mpfr_div(r9164, r9163, r9161, MPFR_RNDN);
        mpfr_cos(r9165, r9160, MPFR_RNDN);
        mpfr_mul(r9166, r9164, r9165, MPFR_RNDN);
        mpfr_sub(r9167, r9162, r9166, MPFR_RNDN);
        return mpfr_get_d(r9167, MPFR_RNDN);
}

static mpfr_t r9168, r9169, r9170, r9171, r9172, r9173, r9174, r9175, r9176;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9168, "1", 10, MPFR_RNDN);
        mpfr_init(r9169);
        mpfr_init(r9170);
        mpfr_init(r9171);
        mpfr_init(r9172);
        mpfr_init(r9173);
        mpfr_init(r9174);
        mpfr_init(r9175);
        mpfr_init(r9176);
}

double f_dm(double B, double x) {
        ;
        mpfr_set_d(r9169, B, MPFR_RNDN);
        mpfr_sin(r9170, r9169, MPFR_RNDN);
        mpfr_div(r9171, r9168, r9170, MPFR_RNDN);
        mpfr_set_d(r9172, x, MPFR_RNDN);
        mpfr_div(r9173, r9172, r9170, MPFR_RNDN);
        mpfr_cos(r9174, r9169, MPFR_RNDN);
        mpfr_mul(r9175, r9173, r9174, MPFR_RNDN);
        mpfr_sub(r9176, r9171, r9175, MPFR_RNDN);
        return mpfr_get_d(r9176, MPFR_RNDN);
}