#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 r8154 = x;
        float r8155 = 1;
        float r8156 = B;
        float r8157 = tan(r8156);
        float r8158 = r8155 / r8157;
        float r8159 = r8154 * r8158;
        float r8160 = -r8159;
        float r8161 = F;
        float r8162 = sin(r8156);
        float r8163 = r8161 / r8162;
        float r8164 = r8161 * r8161;
        float r8165 = 2;
        float r8166 = r8164 + r8165;
        float r8167 = r8165 * r8154;
        float r8168 = r8166 + r8167;
        float r8169 = r8155 / r8165;
        float r8170 = -r8169;
        float r8171 = pow(r8168, r8170);
        float r8172 = r8163 * r8171;
        float r8173 = r8160 + r8172;
        return r8173;
}

double f_id(double F, double B, double x) {
        double r8174 = x;
        double r8175 = 1;
        double r8176 = B;
        double r8177 = tan(r8176);
        double r8178 = r8175 / r8177;
        double r8179 = r8174 * r8178;
        double r8180 = -r8179;
        double r8181 = F;
        double r8182 = sin(r8176);
        double r8183 = r8181 / r8182;
        double r8184 = r8181 * r8181;
        double r8185 = 2;
        double r8186 = r8184 + r8185;
        double r8187 = r8185 * r8174;
        double r8188 = r8186 + r8187;
        double r8189 = r8175 / r8185;
        double r8190 = -r8189;
        double r8191 = pow(r8188, r8190);
        double r8192 = r8183 * r8191;
        double r8193 = r8180 + r8192;
        return r8193;
}


double f_of(float F, float B, float x) {
        float r8194 = F;
        float r8195 = -27714782.050513014;
        bool r8196 = r8194 <= r8195;
        float r8197 = 1;
        float r8198 = B;
        float r8199 = sin(r8198);
        float r8200 = 2;
        float r8201 = pow(r8194, r8200);
        float r8202 = r8199 * r8201;
        float r8203 = r8197 / r8202;
        float r8204 = r8197 / r8199;
        float r8205 = r8203 - r8204;
        float r8206 = x;
        float r8207 = tan(r8198);
        float r8208 = r8206 / r8207;
        float r8209 = r8205 - r8208;
        float r8210 = 77357060.20485163;
        bool r8211 = r8194 <= r8210;
        float r8212 = r8194 * r8194;
        float r8213 = r8200 * r8206;
        float r8214 = r8213 + r8200;
        float r8215 = r8212 + r8214;
        float r8216 = r8197 / r8200;
        float r8217 = pow(r8215, r8216);
        float r8218 = r8217 * r8199;
        float r8219 = r8194 / r8218;
        float r8220 = r8219 - r8208;
        float r8221 = r8204 - r8208;
        float r8222 = r8204 / r8212;
        float r8223 = r8221 - r8222;
        float r8224 = r8211 ? r8220 : r8223;
        float r8225 = r8196 ? r8209 : r8224;
        return r8225;
}

double f_od(double F, double B, double x) {
        double r8226 = F;
        double r8227 = -27714782.050513014;
        bool r8228 = r8226 <= r8227;
        double r8229 = 1;
        double r8230 = B;
        double r8231 = sin(r8230);
        double r8232 = 2;
        double r8233 = pow(r8226, r8232);
        double r8234 = r8231 * r8233;
        double r8235 = r8229 / r8234;
        double r8236 = r8229 / r8231;
        double r8237 = r8235 - r8236;
        double r8238 = x;
        double r8239 = tan(r8230);
        double r8240 = r8238 / r8239;
        double r8241 = r8237 - r8240;
        double r8242 = 77357060.20485163;
        bool r8243 = r8226 <= r8242;
        double r8244 = r8226 * r8226;
        double r8245 = r8232 * r8238;
        double r8246 = r8245 + r8232;
        double r8247 = r8244 + r8246;
        double r8248 = r8229 / r8232;
        double r8249 = pow(r8247, r8248);
        double r8250 = r8249 * r8231;
        double r8251 = r8226 / r8250;
        double r8252 = r8251 - r8240;
        double r8253 = r8236 - r8240;
        double r8254 = r8236 / r8244;
        double r8255 = r8253 - r8254;
        double r8256 = r8243 ? r8252 : r8255;
        double r8257 = r8228 ? r8241 : r8256;
        return r8257;
}

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 r8258, r8259, r8260, r8261, r8262, r8263, r8264, r8265, r8266, r8267, r8268, r8269, r8270, r8271, r8272, r8273, r8274, r8275, r8276, r8277;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r8258);
        mpfr_init_set_str(r8259, "1", 10, MPFR_RNDN);
        mpfr_init(r8260);
        mpfr_init(r8261);
        mpfr_init(r8262);
        mpfr_init(r8263);
        mpfr_init(r8264);
        mpfr_init(r8265);
        mpfr_init(r8266);
        mpfr_init(r8267);
        mpfr_init(r8268);
        mpfr_init_set_str(r8269, "2", 10, MPFR_RNDN);
        mpfr_init(r8270);
        mpfr_init(r8271);
        mpfr_init(r8272);
        mpfr_init(r8273);
        mpfr_init(r8274);
        mpfr_init(r8275);
        mpfr_init(r8276);
        mpfr_init(r8277);
}

double f_im(double F, double B, double x) {
        mpfr_set_d(r8258, x, MPFR_RNDN);
        ;
        mpfr_set_d(r8260, B, MPFR_RNDN);
        mpfr_tan(r8261, r8260, MPFR_RNDN);
        mpfr_div(r8262, r8259, r8261, MPFR_RNDN);
        mpfr_mul(r8263, r8258, r8262, MPFR_RNDN);
        mpfr_neg(r8264, r8263, MPFR_RNDN);
        mpfr_set_d(r8265, F, MPFR_RNDN);
        mpfr_sin(r8266, r8260, MPFR_RNDN);
        mpfr_div(r8267, r8265, r8266, MPFR_RNDN);
        mpfr_mul(r8268, r8265, r8265, MPFR_RNDN);
        ;
        mpfr_add(r8270, r8268, r8269, MPFR_RNDN);
        mpfr_mul(r8271, r8269, r8258, MPFR_RNDN);
        mpfr_add(r8272, r8270, r8271, MPFR_RNDN);
        mpfr_div(r8273, r8259, r8269, MPFR_RNDN);
        mpfr_neg(r8274, r8273, MPFR_RNDN);
        mpfr_pow(r8275, r8272, r8274, MPFR_RNDN);
        mpfr_mul(r8276, r8267, r8275, MPFR_RNDN);
        mpfr_add(r8277, r8264, r8276, MPFR_RNDN);
        return mpfr_get_d(r8277, MPFR_RNDN);
}

static mpfr_t r8278, r8279, r8280, r8281, r8282, r8283, r8284, r8285, r8286, r8287, r8288, r8289, r8290, r8291, r8292, r8293, r8294, r8295, r8296, r8297, r8298, r8299, r8300, r8301, r8302, r8303, r8304, r8305, r8306, r8307, r8308, r8309;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r8278);
        mpfr_init_set_str(r8279, "-27714782.050513014", 10, MPFR_RNDN);
        mpfr_init(r8280);
        mpfr_init_set_str(r8281, "1", 10, MPFR_RNDN);
        mpfr_init(r8282);
        mpfr_init(r8283);
        mpfr_init_set_str(r8284, "2", 10, MPFR_RNDN);
        mpfr_init(r8285);
        mpfr_init(r8286);
        mpfr_init(r8287);
        mpfr_init(r8288);
        mpfr_init(r8289);
        mpfr_init(r8290);
        mpfr_init(r8291);
        mpfr_init(r8292);
        mpfr_init(r8293);
        mpfr_init_set_str(r8294, "77357060.20485163", 10, MPFR_RNDN);
        mpfr_init(r8295);
        mpfr_init(r8296);
        mpfr_init(r8297);
        mpfr_init(r8298);
        mpfr_init(r8299);
        mpfr_init(r8300);
        mpfr_init(r8301);
        mpfr_init(r8302);
        mpfr_init(r8303);
        mpfr_init(r8304);
        mpfr_init(r8305);
        mpfr_init(r8306);
        mpfr_init(r8307);
        mpfr_init(r8308);
        mpfr_init(r8309);
}

double f_fm(double F, double B, double x) {
        mpfr_set_d(r8278, F, MPFR_RNDN);
        ;
        mpfr_set_si(r8280, mpfr_cmp(r8278, r8279) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r8282, B, MPFR_RNDN);
        mpfr_sin(r8283, r8282, MPFR_RNDN);
        ;
        mpfr_pow(r8285, r8278, r8284, MPFR_RNDN);
        mpfr_mul(r8286, r8283, r8285, MPFR_RNDN);
        mpfr_div(r8287, r8281, r8286, MPFR_RNDN);
        mpfr_div(r8288, r8281, r8283, MPFR_RNDN);
        mpfr_sub(r8289, r8287, r8288, MPFR_RNDN);
        mpfr_set_d(r8290, x, MPFR_RNDN);
        mpfr_tan(r8291, r8282, MPFR_RNDN);
        mpfr_div(r8292, r8290, r8291, MPFR_RNDN);
        mpfr_sub(r8293, r8289, r8292, MPFR_RNDN);
        ;
        mpfr_set_si(r8295, mpfr_cmp(r8278, r8294) <= 0, MPFR_RNDN);
        mpfr_mul(r8296, r8278, r8278, MPFR_RNDN);
        mpfr_mul(r8297, r8284, r8290, MPFR_RNDN);
        mpfr_add(r8298, r8297, r8284, MPFR_RNDN);
        mpfr_add(r8299, r8296, r8298, MPFR_RNDN);
        mpfr_div(r8300, r8281, r8284, MPFR_RNDN);
        mpfr_pow(r8301, r8299, r8300, MPFR_RNDN);
        mpfr_mul(r8302, r8301, r8283, MPFR_RNDN);
        mpfr_div(r8303, r8278, r8302, MPFR_RNDN);
        mpfr_sub(r8304, r8303, r8292, MPFR_RNDN);
        mpfr_sub(r8305, r8288, r8292, MPFR_RNDN);
        mpfr_div(r8306, r8288, r8296, MPFR_RNDN);
        mpfr_sub(r8307, r8305, r8306, MPFR_RNDN);
        if (mpfr_get_si(r8295, MPFR_RNDN)) { mpfr_set(r8308, r8304, MPFR_RNDN); } else { mpfr_set(r8308, r8307, MPFR_RNDN); };
        if (mpfr_get_si(r8280, MPFR_RNDN)) { mpfr_set(r8309, r8293, MPFR_RNDN); } else { mpfr_set(r8309, r8308, MPFR_RNDN); };
        return mpfr_get_d(r8309, MPFR_RNDN);
}

static mpfr_t r8310, r8311, r8312, r8313, r8314, r8315, r8316, r8317, r8318, r8319, r8320, r8321, r8322, r8323, r8324, r8325, r8326, r8327, r8328, r8329, r8330, r8331, r8332, r8333, r8334, r8335, r8336, r8337, r8338, r8339, r8340, r8341;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r8310);
        mpfr_init_set_str(r8311, "-27714782.050513014", 10, MPFR_RNDN);
        mpfr_init(r8312);
        mpfr_init_set_str(r8313, "1", 10, MPFR_RNDN);
        mpfr_init(r8314);
        mpfr_init(r8315);
        mpfr_init_set_str(r8316, "2", 10, MPFR_RNDN);
        mpfr_init(r8317);
        mpfr_init(r8318);
        mpfr_init(r8319);
        mpfr_init(r8320);
        mpfr_init(r8321);
        mpfr_init(r8322);
        mpfr_init(r8323);
        mpfr_init(r8324);
        mpfr_init(r8325);
        mpfr_init_set_str(r8326, "77357060.20485163", 10, MPFR_RNDN);
        mpfr_init(r8327);
        mpfr_init(r8328);
        mpfr_init(r8329);
        mpfr_init(r8330);
        mpfr_init(r8331);
        mpfr_init(r8332);
        mpfr_init(r8333);
        mpfr_init(r8334);
        mpfr_init(r8335);
        mpfr_init(r8336);
        mpfr_init(r8337);
        mpfr_init(r8338);
        mpfr_init(r8339);
        mpfr_init(r8340);
        mpfr_init(r8341);
}

double f_dm(double F, double B, double x) {
        mpfr_set_d(r8310, F, MPFR_RNDN);
        ;
        mpfr_set_si(r8312, mpfr_cmp(r8310, r8311) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r8314, B, MPFR_RNDN);
        mpfr_sin(r8315, r8314, MPFR_RNDN);
        ;
        mpfr_pow(r8317, r8310, r8316, MPFR_RNDN);
        mpfr_mul(r8318, r8315, r8317, MPFR_RNDN);
        mpfr_div(r8319, r8313, r8318, MPFR_RNDN);
        mpfr_div(r8320, r8313, r8315, MPFR_RNDN);
        mpfr_sub(r8321, r8319, r8320, MPFR_RNDN);
        mpfr_set_d(r8322, x, MPFR_RNDN);
        mpfr_tan(r8323, r8314, MPFR_RNDN);
        mpfr_div(r8324, r8322, r8323, MPFR_RNDN);
        mpfr_sub(r8325, r8321, r8324, MPFR_RNDN);
        ;
        mpfr_set_si(r8327, mpfr_cmp(r8310, r8326) <= 0, MPFR_RNDN);
        mpfr_mul(r8328, r8310, r8310, MPFR_RNDN);
        mpfr_mul(r8329, r8316, r8322, MPFR_RNDN);
        mpfr_add(r8330, r8329, r8316, MPFR_RNDN);
        mpfr_add(r8331, r8328, r8330, MPFR_RNDN);
        mpfr_div(r8332, r8313, r8316, MPFR_RNDN);
        mpfr_pow(r8333, r8331, r8332, MPFR_RNDN);
        mpfr_mul(r8334, r8333, r8315, MPFR_RNDN);
        mpfr_div(r8335, r8310, r8334, MPFR_RNDN);
        mpfr_sub(r8336, r8335, r8324, MPFR_RNDN);
        mpfr_sub(r8337, r8320, r8324, MPFR_RNDN);
        mpfr_div(r8338, r8320, r8328, MPFR_RNDN);
        mpfr_sub(r8339, r8337, r8338, MPFR_RNDN);
        if (mpfr_get_si(r8327, MPFR_RNDN)) { mpfr_set(r8340, r8336, MPFR_RNDN); } else { mpfr_set(r8340, r8339, MPFR_RNDN); };
        if (mpfr_get_si(r8312, MPFR_RNDN)) { mpfr_set(r8341, r8325, MPFR_RNDN); } else { mpfr_set(r8341, r8340, MPFR_RNDN); };
        return mpfr_get_d(r8341, MPFR_RNDN);
}