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

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

double f_if(float F, float l) {
        float r23366 = atan2(1.0, 0.0);
        float r23367 = l;
        float r23368 = r23366 * r23367;
        float r23369 = 1;
        float r23370 = F;
        float r23371 = r23370 * r23370;
        float r23372 = r23369 / r23371;
        float r23373 = tan(r23368);
        float r23374 = r23372 * r23373;
        float r23375 = r23368 - r23374;
        return r23375;
}

double f_id(double F, double l) {
        double r23376 = atan2(1.0, 0.0);
        double r23377 = l;
        double r23378 = r23376 * r23377;
        double r23379 = 1;
        double r23380 = F;
        double r23381 = r23380 * r23380;
        double r23382 = r23379 / r23381;
        double r23383 = tan(r23378);
        double r23384 = r23382 * r23383;
        double r23385 = r23378 - r23384;
        return r23385;
}


double f_of(float F, float l) {
        float r23386 = l;
        float r23387 = -3.724350483518799e+27;
        bool r23388 = r23386 <= r23387;
        float r23389 = atan2(1.0, 0.0);
        float r23390 = r23389 * r23386;
        float r23391 = F;
        float r23392 = r23389 / r23386;
        float r23393 = sin(r23392);
        float r23394 = cos(r23392);
        float r23395 = r23393 / r23394;
        float r23396 = r23391 * r23395;
        float r23397 = cbrt(r23396);
        float r23398 = tan(r23390);
        float r23399 = r23398 / r23391;
        float r23400 = cbrt(r23399);
        float r23401 = r23400 / r23391;
        float r23402 = r23397 * r23401;
        float r23403 = r23402 * r23400;
        float r23404 = r23390 - r23403;
        float r23405 = 1.1221113420603007e+89;
        bool r23406 = r23386 <= r23405;
        float r23407 = r23386 * r23389;
        float r23408 = tan(r23407);
        float r23409 = r23408 / r23391;
        float r23410 = r23409 / r23391;
        float r23411 = r23407 - r23410;
        float r23412 = r23406 ? r23411 : r23404;
        float r23413 = r23388 ? r23404 : r23412;
        return r23413;
}

double f_od(double F, double l) {
        double r23414 = l;
        double r23415 = -3.724350483518799e+27;
        bool r23416 = r23414 <= r23415;
        double r23417 = atan2(1.0, 0.0);
        double r23418 = r23417 * r23414;
        double r23419 = F;
        double r23420 = r23417 / r23414;
        double r23421 = sin(r23420);
        double r23422 = cos(r23420);
        double r23423 = r23421 / r23422;
        double r23424 = r23419 * r23423;
        double r23425 = cbrt(r23424);
        double r23426 = tan(r23418);
        double r23427 = r23426 / r23419;
        double r23428 = cbrt(r23427);
        double r23429 = r23428 / r23419;
        double r23430 = r23425 * r23429;
        double r23431 = r23430 * r23428;
        double r23432 = r23418 - r23431;
        double r23433 = 1.1221113420603007e+89;
        bool r23434 = r23414 <= r23433;
        double r23435 = r23414 * r23417;
        double r23436 = tan(r23435);
        double r23437 = r23436 / r23419;
        double r23438 = r23437 / r23419;
        double r23439 = r23435 - r23438;
        double r23440 = r23434 ? r23439 : r23432;
        double r23441 = r23416 ? r23432 : r23440;
        return r23441;
}

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 r23442, r23443, r23444, r23445, r23446, r23447, r23448, r23449, r23450, r23451;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2448);
        mpfr_init(r23442);
        mpfr_init(r23443);
        mpfr_init(r23444);
        mpfr_init_set_str(r23445, "1", 10, MPFR_RNDN);
        mpfr_init(r23446);
        mpfr_init(r23447);
        mpfr_init(r23448);
        mpfr_init(r23449);
        mpfr_init(r23450);
        mpfr_init(r23451);
}

double f_im(double F, double l) {
        mpfr_const_pi(r23442, MPFR_RNDN);
        mpfr_set_d(r23443, l, MPFR_RNDN);
        mpfr_mul(r23444, r23442, r23443, MPFR_RNDN);
        ;
        mpfr_set_d(r23446, F, MPFR_RNDN);
        mpfr_mul(r23447, r23446, r23446, MPFR_RNDN);
        mpfr_div(r23448, r23445, r23447, MPFR_RNDN);
        mpfr_tan(r23449, r23444, MPFR_RNDN);
        mpfr_mul(r23450, r23448, r23449, MPFR_RNDN);
        mpfr_sub(r23451, r23444, r23450, MPFR_RNDN);
        return mpfr_get_d(r23451, MPFR_RNDN);
}

static mpfr_t r23452, r23453, r23454, r23455, r23456, r23457, r23458, r23459, r23460, r23461, r23462, r23463, r23464, r23465, r23466, r23467, r23468, r23469, r23470, r23471, r23472, r23473, r23474, r23475, r23476, r23477, r23478, r23479;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r23452);
        mpfr_init_set_str(r23453, "-3.724350483518799e+27", 10, MPFR_RNDN);
        mpfr_init(r23454);
        mpfr_init(r23455);
        mpfr_init(r23456);
        mpfr_init(r23457);
        mpfr_init(r23458);
        mpfr_init(r23459);
        mpfr_init(r23460);
        mpfr_init(r23461);
        mpfr_init(r23462);
        mpfr_init(r23463);
        mpfr_init(r23464);
        mpfr_init(r23465);
        mpfr_init(r23466);
        mpfr_init(r23467);
        mpfr_init(r23468);
        mpfr_init(r23469);
        mpfr_init(r23470);
        mpfr_init_set_str(r23471, "1.1221113420603007e+89", 10, MPFR_RNDN);
        mpfr_init(r23472);
        mpfr_init(r23473);
        mpfr_init(r23474);
        mpfr_init(r23475);
        mpfr_init(r23476);
        mpfr_init(r23477);
        mpfr_init(r23478);
        mpfr_init(r23479);
}

double f_fm(double F, double l) {
        mpfr_set_d(r23452, l, MPFR_RNDN);
        ;
        mpfr_set_si(r23454, mpfr_cmp(r23452, r23453) <= 0, MPFR_RNDN);
        mpfr_const_pi(r23455, MPFR_RNDN);
        mpfr_mul(r23456, r23455, r23452, MPFR_RNDN);
        mpfr_set_d(r23457, F, MPFR_RNDN);
        mpfr_div(r23458, r23455, r23452, MPFR_RNDN);
        mpfr_sin(r23459, r23458, MPFR_RNDN);
        mpfr_cos(r23460, r23458, MPFR_RNDN);
        mpfr_div(r23461, r23459, r23460, MPFR_RNDN);
        mpfr_mul(r23462, r23457, r23461, MPFR_RNDN);
        mpfr_cbrt(r23463, r23462, MPFR_RNDN);
        mpfr_tan(r23464, r23456, MPFR_RNDN);
        mpfr_div(r23465, r23464, r23457, MPFR_RNDN);
        mpfr_cbrt(r23466, r23465, MPFR_RNDN);
        mpfr_div(r23467, r23466, r23457, MPFR_RNDN);
        mpfr_mul(r23468, r23463, r23467, MPFR_RNDN);
        mpfr_mul(r23469, r23468, r23466, MPFR_RNDN);
        mpfr_sub(r23470, r23456, r23469, MPFR_RNDN);
        ;
        mpfr_set_si(r23472, mpfr_cmp(r23452, r23471) <= 0, MPFR_RNDN);
        mpfr_mul(r23473, r23452, r23455, MPFR_RNDN);
        mpfr_tan(r23474, r23473, MPFR_RNDN);
        mpfr_div(r23475, r23474, r23457, MPFR_RNDN);
        mpfr_div(r23476, r23475, r23457, MPFR_RNDN);
        mpfr_sub(r23477, r23473, r23476, MPFR_RNDN);
        if (mpfr_get_si(r23472, MPFR_RNDN)) { mpfr_set(r23478, r23477, MPFR_RNDN); } else { mpfr_set(r23478, r23470, MPFR_RNDN); };
        if (mpfr_get_si(r23454, MPFR_RNDN)) { mpfr_set(r23479, r23470, MPFR_RNDN); } else { mpfr_set(r23479, r23478, MPFR_RNDN); };
        return mpfr_get_d(r23479, MPFR_RNDN);
}

static mpfr_t r23480, r23481, r23482, r23483, r23484, r23485, r23486, r23487, r23488, r23489, r23490, r23491, r23492, r23493, r23494, r23495, r23496, r23497, r23498, r23499, r23500, r23501, r23502, r23503, r23504, r23505, r23506, r23507;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r23480);
        mpfr_init_set_str(r23481, "-3.724350483518799e+27", 10, MPFR_RNDN);
        mpfr_init(r23482);
        mpfr_init(r23483);
        mpfr_init(r23484);
        mpfr_init(r23485);
        mpfr_init(r23486);
        mpfr_init(r23487);
        mpfr_init(r23488);
        mpfr_init(r23489);
        mpfr_init(r23490);
        mpfr_init(r23491);
        mpfr_init(r23492);
        mpfr_init(r23493);
        mpfr_init(r23494);
        mpfr_init(r23495);
        mpfr_init(r23496);
        mpfr_init(r23497);
        mpfr_init(r23498);
        mpfr_init_set_str(r23499, "1.1221113420603007e+89", 10, MPFR_RNDN);
        mpfr_init(r23500);
        mpfr_init(r23501);
        mpfr_init(r23502);
        mpfr_init(r23503);
        mpfr_init(r23504);
        mpfr_init(r23505);
        mpfr_init(r23506);
        mpfr_init(r23507);
}

double f_dm(double F, double l) {
        mpfr_set_d(r23480, l, MPFR_RNDN);
        ;
        mpfr_set_si(r23482, mpfr_cmp(r23480, r23481) <= 0, MPFR_RNDN);
        mpfr_const_pi(r23483, MPFR_RNDN);
        mpfr_mul(r23484, r23483, r23480, MPFR_RNDN);
        mpfr_set_d(r23485, F, MPFR_RNDN);
        mpfr_div(r23486, r23483, r23480, MPFR_RNDN);
        mpfr_sin(r23487, r23486, MPFR_RNDN);
        mpfr_cos(r23488, r23486, MPFR_RNDN);
        mpfr_div(r23489, r23487, r23488, MPFR_RNDN);
        mpfr_mul(r23490, r23485, r23489, MPFR_RNDN);
        mpfr_cbrt(r23491, r23490, MPFR_RNDN);
        mpfr_tan(r23492, r23484, MPFR_RNDN);
        mpfr_div(r23493, r23492, r23485, MPFR_RNDN);
        mpfr_cbrt(r23494, r23493, MPFR_RNDN);
        mpfr_div(r23495, r23494, r23485, MPFR_RNDN);
        mpfr_mul(r23496, r23491, r23495, MPFR_RNDN);
        mpfr_mul(r23497, r23496, r23494, MPFR_RNDN);
        mpfr_sub(r23498, r23484, r23497, MPFR_RNDN);
        ;
        mpfr_set_si(r23500, mpfr_cmp(r23480, r23499) <= 0, MPFR_RNDN);
        mpfr_mul(r23501, r23480, r23483, MPFR_RNDN);
        mpfr_tan(r23502, r23501, MPFR_RNDN);
        mpfr_div(r23503, r23502, r23485, MPFR_RNDN);
        mpfr_div(r23504, r23503, r23485, MPFR_RNDN);
        mpfr_sub(r23505, r23501, r23504, MPFR_RNDN);
        if (mpfr_get_si(r23500, MPFR_RNDN)) { mpfr_set(r23506, r23505, MPFR_RNDN); } else { mpfr_set(r23506, r23498, MPFR_RNDN); };
        if (mpfr_get_si(r23482, MPFR_RNDN)) { mpfr_set(r23507, r23498, MPFR_RNDN); } else { mpfr_set(r23507, r23506, MPFR_RNDN); };
        return mpfr_get_d(r23507, MPFR_RNDN);
}