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

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

double f_if(float F, float B, float x) {
        float r24277 = x;
        float r24278 = 1;
        float r24279 = B;
        float r24280 = tan(r24279);
        float r24281 = r24278 / r24280;
        float r24282 = r24277 * r24281;
        float r24283 = -r24282;
        float r24284 = F;
        float r24285 = sin(r24279);
        float r24286 = r24284 / r24285;
        float r24287 = r24284 * r24284;
        float r24288 = 2;
        float r24289 = r24287 + r24288;
        float r24290 = r24288 * r24277;
        float r24291 = r24289 + r24290;
        float r24292 = r24278 / r24288;
        float r24293 = -r24292;
        float r24294 = pow(r24291, r24293);
        float r24295 = r24286 * r24294;
        float r24296 = r24283 + r24295;
        return r24296;
}

double f_id(double F, double B, double x) {
        double r24297 = x;
        double r24298 = 1;
        double r24299 = B;
        double r24300 = tan(r24299);
        double r24301 = r24298 / r24300;
        double r24302 = r24297 * r24301;
        double r24303 = -r24302;
        double r24304 = F;
        double r24305 = sin(r24299);
        double r24306 = r24304 / r24305;
        double r24307 = r24304 * r24304;
        double r24308 = 2;
        double r24309 = r24307 + r24308;
        double r24310 = r24308 * r24297;
        double r24311 = r24309 + r24310;
        double r24312 = r24298 / r24308;
        double r24313 = -r24312;
        double r24314 = pow(r24311, r24313);
        double r24315 = r24306 * r24314;
        double r24316 = r24303 + r24315;
        return r24316;
}


double f_of(float F, float B, float x) {
        float r24317 = x;
        float r24318 = 2;
        float r24319 = F;
        float r24320 = fma(r24319, r24319, r24318);
        float r24321 = fma(r24317, r24318, r24320);
        float r24322 = 1;
        float r24323 = r24322 / r24318;
        float r24324 = -r24323;
        float r24325 = pow(r24321, r24324);
        float r24326 = B;
        float r24327 = sin(r24326);
        float r24328 = r24322 / r24327;
        float r24329 = r24319 * r24328;
        float r24330 = -r24317;
        float r24331 = tan(r24326);
        float r24332 = r24330 / r24331;
        float r24333 = fma(r24325, r24329, r24332);
        return r24333;
}

double f_od(double F, double B, double x) {
        double r24334 = x;
        double r24335 = 2;
        double r24336 = F;
        double r24337 = fma(r24336, r24336, r24335);
        double r24338 = fma(r24334, r24335, r24337);
        double r24339 = 1;
        double r24340 = r24339 / r24335;
        double r24341 = -r24340;
        double r24342 = pow(r24338, r24341);
        double r24343 = B;
        double r24344 = sin(r24343);
        double r24345 = r24339 / r24344;
        double r24346 = r24336 * r24345;
        double r24347 = -r24334;
        double r24348 = tan(r24343);
        double r24349 = r24347 / r24348;
        double r24350 = fma(r24342, r24346, r24349);
        return r24350;
}

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 r24351, r24352, r24353, r24354, r24355, r24356, r24357, r24358, r24359, r24360, r24361, r24362, r24363, r24364, r24365, r24366, r24367, r24368, r24369, r24370;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r24351);
        mpfr_init_set_str(r24352, "1", 10, MPFR_RNDN);
        mpfr_init(r24353);
        mpfr_init(r24354);
        mpfr_init(r24355);
        mpfr_init(r24356);
        mpfr_init(r24357);
        mpfr_init(r24358);
        mpfr_init(r24359);
        mpfr_init(r24360);
        mpfr_init(r24361);
        mpfr_init_set_str(r24362, "2", 10, MPFR_RNDN);
        mpfr_init(r24363);
        mpfr_init(r24364);
        mpfr_init(r24365);
        mpfr_init(r24366);
        mpfr_init(r24367);
        mpfr_init(r24368);
        mpfr_init(r24369);
        mpfr_init(r24370);
}

double f_im(double F, double B, double x) {
        mpfr_set_d(r24351, x, MPFR_RNDN);
        ;
        mpfr_set_d(r24353, B, MPFR_RNDN);
        mpfr_tan(r24354, r24353, MPFR_RNDN);
        mpfr_div(r24355, r24352, r24354, MPFR_RNDN);
        mpfr_mul(r24356, r24351, r24355, MPFR_RNDN);
        mpfr_neg(r24357, r24356, MPFR_RNDN);
        mpfr_set_d(r24358, F, MPFR_RNDN);
        mpfr_sin(r24359, r24353, MPFR_RNDN);
        mpfr_div(r24360, r24358, r24359, MPFR_RNDN);
        mpfr_mul(r24361, r24358, r24358, MPFR_RNDN);
        ;
        mpfr_add(r24363, r24361, r24362, MPFR_RNDN);
        mpfr_mul(r24364, r24362, r24351, MPFR_RNDN);
        mpfr_add(r24365, r24363, r24364, MPFR_RNDN);
        mpfr_div(r24366, r24352, r24362, MPFR_RNDN);
        mpfr_neg(r24367, r24366, MPFR_RNDN);
        mpfr_pow(r24368, r24365, r24367, MPFR_RNDN);
        mpfr_mul(r24369, r24360, r24368, MPFR_RNDN);
        mpfr_add(r24370, r24357, r24369, MPFR_RNDN);
        return mpfr_get_d(r24370, MPFR_RNDN);
}

static mpfr_t r24371, r24372, r24373, r24374, r24375, r24376, r24377, r24378, r24379, r24380, r24381, r24382, r24383, r24384, r24385, r24386, r24387;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r24371);
        mpfr_init_set_str(r24372, "2", 10, MPFR_RNDN);
        mpfr_init(r24373);
        mpfr_init(r24374);
        mpfr_init(r24375);
        mpfr_init_set_str(r24376, "1", 10, MPFR_RNDN);
        mpfr_init(r24377);
        mpfr_init(r24378);
        mpfr_init(r24379);
        mpfr_init(r24380);
        mpfr_init(r24381);
        mpfr_init(r24382);
        mpfr_init(r24383);
        mpfr_init(r24384);
        mpfr_init(r24385);
        mpfr_init(r24386);
        mpfr_init(r24387);
}

double f_fm(double F, double B, double x) {
        mpfr_set_d(r24371, x, MPFR_RNDN);
        ;
        mpfr_set_d(r24373, F, MPFR_RNDN);
        mpfr_fma(r24374, r24373, r24373, r24372, MPFR_RNDN);
        mpfr_fma(r24375, r24371, r24372, r24374, MPFR_RNDN);
        ;
        mpfr_div(r24377, r24376, r24372, MPFR_RNDN);
        mpfr_neg(r24378, r24377, MPFR_RNDN);
        mpfr_pow(r24379, r24375, r24378, MPFR_RNDN);
        mpfr_set_d(r24380, B, MPFR_RNDN);
        mpfr_sin(r24381, r24380, MPFR_RNDN);
        mpfr_div(r24382, r24376, r24381, MPFR_RNDN);
        mpfr_mul(r24383, r24373, r24382, MPFR_RNDN);
        mpfr_neg(r24384, r24371, MPFR_RNDN);
        mpfr_tan(r24385, r24380, MPFR_RNDN);
        mpfr_div(r24386, r24384, r24385, MPFR_RNDN);
        mpfr_fma(r24387, r24379, r24383, r24386, MPFR_RNDN);
        return mpfr_get_d(r24387, MPFR_RNDN);
}

static mpfr_t r24388, r24389, r24390, r24391, r24392, r24393, r24394, r24395, r24396, r24397, r24398, r24399, r24400, r24401, r24402, r24403, r24404;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r24388);
        mpfr_init_set_str(r24389, "2", 10, MPFR_RNDN);
        mpfr_init(r24390);
        mpfr_init(r24391);
        mpfr_init(r24392);
        mpfr_init_set_str(r24393, "1", 10, MPFR_RNDN);
        mpfr_init(r24394);
        mpfr_init(r24395);
        mpfr_init(r24396);
        mpfr_init(r24397);
        mpfr_init(r24398);
        mpfr_init(r24399);
        mpfr_init(r24400);
        mpfr_init(r24401);
        mpfr_init(r24402);
        mpfr_init(r24403);
        mpfr_init(r24404);
}

double f_dm(double F, double B, double x) {
        mpfr_set_d(r24388, x, MPFR_RNDN);
        ;
        mpfr_set_d(r24390, F, MPFR_RNDN);
        mpfr_fma(r24391, r24390, r24390, r24389, MPFR_RNDN);
        mpfr_fma(r24392, r24388, r24389, r24391, MPFR_RNDN);
        ;
        mpfr_div(r24394, r24393, r24389, MPFR_RNDN);
        mpfr_neg(r24395, r24394, MPFR_RNDN);
        mpfr_pow(r24396, r24392, r24395, MPFR_RNDN);
        mpfr_set_d(r24397, B, MPFR_RNDN);
        mpfr_sin(r24398, r24397, MPFR_RNDN);
        mpfr_div(r24399, r24393, r24398, MPFR_RNDN);
        mpfr_mul(r24400, r24390, r24399, MPFR_RNDN);
        mpfr_neg(r24401, r24388, MPFR_RNDN);
        mpfr_tan(r24402, r24397, MPFR_RNDN);
        mpfr_div(r24403, r24401, r24402, MPFR_RNDN);
        mpfr_fma(r24404, r24396, r24400, r24403, MPFR_RNDN);
        return mpfr_get_d(r24404, MPFR_RNDN);
}