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

char *name = "Falkner and Boettcher, Equation (20:1,3)";

double f_if(float v, float t) {
        float r27194 = 1;
        float r27195 = 5;
        float r27196 = v;
        float r27197 = r27196 * r27196;
        float r27198 = r27195 * r27197;
        float r27199 = r27194 - r27198;
        float r27200 = atan2(1.0, 0.0);
        float r27201 = t;
        float r27202 = r27200 * r27201;
        float r27203 = 2;
        float r27204 = 3;
        float r27205 = r27204 * r27197;
        float r27206 = r27194 - r27205;
        float r27207 = r27203 * r27206;
        float r27208 = sqrt(r27207);
        float r27209 = r27202 * r27208;
        float r27210 = r27194 - r27197;
        float r27211 = r27209 * r27210;
        float r27212 = r27199 / r27211;
        return r27212;
}

double f_id(double v, double t) {
        double r27213 = 1;
        double r27214 = 5;
        double r27215 = v;
        double r27216 = r27215 * r27215;
        double r27217 = r27214 * r27216;
        double r27218 = r27213 - r27217;
        double r27219 = atan2(1.0, 0.0);
        double r27220 = t;
        double r27221 = r27219 * r27220;
        double r27222 = 2;
        double r27223 = 3;
        double r27224 = r27223 * r27216;
        double r27225 = r27213 - r27224;
        double r27226 = r27222 * r27225;
        double r27227 = sqrt(r27226);
        double r27228 = r27221 * r27227;
        double r27229 = r27213 - r27216;
        double r27230 = r27228 * r27229;
        double r27231 = r27218 / r27230;
        return r27231;
}


double f_of(float v, float t) {
        float r27232 = 5;
        float r27233 = -r27232;
        float r27234 = v;
        float r27235 = r27234 * r27234;
        float r27236 = 1;
        float r27237 = fma(r27233, r27235, r27236);
        float r27238 = 3;
        float r27239 = -r27238;
        float r27240 = r27235 * r27239;
        float r27241 = r27235 * r27238;
        float r27242 = fma(r27240, r27241, r27236);
        float r27243 = 2;
        float r27244 = r27242 * r27243;
        float r27245 = sqrt(r27244);
        float r27246 = r27237 / r27245;
        float r27247 = pow(r27234, r27238);
        float r27248 = -r27234;
        float r27249 = fma(r27247, r27248, r27236);
        float r27250 = t;
        float r27251 = atan2(1.0, 0.0);
        float r27252 = r27250 * r27251;
        float r27253 = r27249 * r27252;
        float r27254 = r27246 / r27253;
        float r27255 = r27238 * r27235;
        float r27256 = r27236 + r27255;
        float r27257 = sqrt(r27256);
        float r27258 = r27236 + r27235;
        float r27259 = r27257 * r27258;
        float r27260 = r27254 * r27259;
        return r27260;
}

double f_od(double v, double t) {
        double r27261 = 5;
        double r27262 = -r27261;
        double r27263 = v;
        double r27264 = r27263 * r27263;
        double r27265 = 1;
        double r27266 = fma(r27262, r27264, r27265);
        double r27267 = 3;
        double r27268 = -r27267;
        double r27269 = r27264 * r27268;
        double r27270 = r27264 * r27267;
        double r27271 = fma(r27269, r27270, r27265);
        double r27272 = 2;
        double r27273 = r27271 * r27272;
        double r27274 = sqrt(r27273);
        double r27275 = r27266 / r27274;
        double r27276 = pow(r27263, r27267);
        double r27277 = -r27263;
        double r27278 = fma(r27276, r27277, r27265);
        double r27279 = t;
        double r27280 = atan2(1.0, 0.0);
        double r27281 = r27279 * r27280;
        double r27282 = r27278 * r27281;
        double r27283 = r27275 / r27282;
        double r27284 = r27267 * r27264;
        double r27285 = r27265 + r27284;
        double r27286 = sqrt(r27285);
        double r27287 = r27265 + r27264;
        double r27288 = r27286 * r27287;
        double r27289 = r27283 * r27288;
        return r27289;
}

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 r27290, r27291, r27292, r27293, r27294, r27295, r27296, r27297, r27298, r27299, r27300, r27301, r27302, r27303, r27304, r27305, r27306, r27307, r27308;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27290, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27291, "5", 10, MPFR_RNDN);
        mpfr_init(r27292);
        mpfr_init(r27293);
        mpfr_init(r27294);
        mpfr_init(r27295);
        mpfr_init(r27296);
        mpfr_init(r27297);
        mpfr_init(r27298);
        mpfr_init_set_str(r27299, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r27300, "3", 10, MPFR_RNDN);
        mpfr_init(r27301);
        mpfr_init(r27302);
        mpfr_init(r27303);
        mpfr_init(r27304);
        mpfr_init(r27305);
        mpfr_init(r27306);
        mpfr_init(r27307);
        mpfr_init(r27308);
}

double f_im(double v, double t) {
        ;
        ;
        mpfr_set_d(r27292, v, MPFR_RNDN);
        mpfr_mul(r27293, r27292, r27292, MPFR_RNDN);
        mpfr_mul(r27294, r27291, r27293, MPFR_RNDN);
        mpfr_sub(r27295, r27290, r27294, MPFR_RNDN);
        mpfr_const_pi(r27296, MPFR_RNDN);
        mpfr_set_d(r27297, t, MPFR_RNDN);
        mpfr_mul(r27298, r27296, r27297, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r27301, r27300, r27293, MPFR_RNDN);
        mpfr_sub(r27302, r27290, r27301, MPFR_RNDN);
        mpfr_mul(r27303, r27299, r27302, MPFR_RNDN);
        mpfr_sqrt(r27304, r27303, MPFR_RNDN);
        mpfr_mul(r27305, r27298, r27304, MPFR_RNDN);
        mpfr_sub(r27306, r27290, r27293, MPFR_RNDN);
        mpfr_mul(r27307, r27305, r27306, MPFR_RNDN);
        mpfr_div(r27308, r27295, r27307, MPFR_RNDN);
        return mpfr_get_d(r27308, MPFR_RNDN);
}

static mpfr_t r27309, r27310, r27311, r27312, r27313, r27314, r27315, r27316, r27317, r27318, r27319, r27320, r27321, r27322, r27323, r27324, r27325, r27326, r27327, r27328, r27329, r27330, r27331, r27332, r27333, r27334, r27335, r27336, r27337;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27309, "5", 10, MPFR_RNDN);
        mpfr_init(r27310);
        mpfr_init(r27311);
        mpfr_init(r27312);
        mpfr_init_set_str(r27313, "1", 10, MPFR_RNDN);
        mpfr_init(r27314);
        mpfr_init_set_str(r27315, "3", 10, MPFR_RNDN);
        mpfr_init(r27316);
        mpfr_init(r27317);
        mpfr_init(r27318);
        mpfr_init(r27319);
        mpfr_init_set_str(r27320, "2", 10, MPFR_RNDN);
        mpfr_init(r27321);
        mpfr_init(r27322);
        mpfr_init(r27323);
        mpfr_init(r27324);
        mpfr_init(r27325);
        mpfr_init(r27326);
        mpfr_init(r27327);
        mpfr_init(r27328);
        mpfr_init(r27329);
        mpfr_init(r27330);
        mpfr_init(r27331);
        mpfr_init(r27332);
        mpfr_init(r27333);
        mpfr_init(r27334);
        mpfr_init(r27335);
        mpfr_init(r27336);
        mpfr_init(r27337);
}

double f_fm(double v, double t) {
        ;
        mpfr_neg(r27310, r27309, MPFR_RNDN);
        mpfr_set_d(r27311, v, MPFR_RNDN);
        mpfr_mul(r27312, r27311, r27311, MPFR_RNDN);
        ;
        mpfr_fma(r27314, r27310, r27312, r27313, MPFR_RNDN);
        ;
        mpfr_neg(r27316, r27315, MPFR_RNDN);
        mpfr_mul(r27317, r27312, r27316, MPFR_RNDN);
        mpfr_mul(r27318, r27312, r27315, MPFR_RNDN);
        mpfr_fma(r27319, r27317, r27318, r27313, MPFR_RNDN);
        ;
        mpfr_mul(r27321, r27319, r27320, MPFR_RNDN);
        mpfr_sqrt(r27322, r27321, MPFR_RNDN);
        mpfr_div(r27323, r27314, r27322, MPFR_RNDN);
        mpfr_pow(r27324, r27311, r27315, MPFR_RNDN);
        mpfr_neg(r27325, r27311, MPFR_RNDN);
        mpfr_fma(r27326, r27324, r27325, r27313, MPFR_RNDN);
        mpfr_set_d(r27327, t, MPFR_RNDN);
        mpfr_const_pi(r27328, MPFR_RNDN);
        mpfr_mul(r27329, r27327, r27328, MPFR_RNDN);
        mpfr_mul(r27330, r27326, r27329, MPFR_RNDN);
        mpfr_div(r27331, r27323, r27330, MPFR_RNDN);
        mpfr_mul(r27332, r27315, r27312, MPFR_RNDN);
        mpfr_add(r27333, r27313, r27332, MPFR_RNDN);
        mpfr_sqrt(r27334, r27333, MPFR_RNDN);
        mpfr_add(r27335, r27313, r27312, MPFR_RNDN);
        mpfr_mul(r27336, r27334, r27335, MPFR_RNDN);
        mpfr_mul(r27337, r27331, r27336, MPFR_RNDN);
        return mpfr_get_d(r27337, MPFR_RNDN);
}

static mpfr_t r27338, r27339, r27340, r27341, r27342, r27343, r27344, r27345, r27346, r27347, r27348, r27349, r27350, r27351, r27352, r27353, r27354, r27355, r27356, r27357, r27358, r27359, r27360, r27361, r27362, r27363, r27364, r27365, r27366;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27338, "5", 10, MPFR_RNDN);
        mpfr_init(r27339);
        mpfr_init(r27340);
        mpfr_init(r27341);
        mpfr_init_set_str(r27342, "1", 10, MPFR_RNDN);
        mpfr_init(r27343);
        mpfr_init_set_str(r27344, "3", 10, MPFR_RNDN);
        mpfr_init(r27345);
        mpfr_init(r27346);
        mpfr_init(r27347);
        mpfr_init(r27348);
        mpfr_init_set_str(r27349, "2", 10, MPFR_RNDN);
        mpfr_init(r27350);
        mpfr_init(r27351);
        mpfr_init(r27352);
        mpfr_init(r27353);
        mpfr_init(r27354);
        mpfr_init(r27355);
        mpfr_init(r27356);
        mpfr_init(r27357);
        mpfr_init(r27358);
        mpfr_init(r27359);
        mpfr_init(r27360);
        mpfr_init(r27361);
        mpfr_init(r27362);
        mpfr_init(r27363);
        mpfr_init(r27364);
        mpfr_init(r27365);
        mpfr_init(r27366);
}

double f_dm(double v, double t) {
        ;
        mpfr_neg(r27339, r27338, MPFR_RNDN);
        mpfr_set_d(r27340, v, MPFR_RNDN);
        mpfr_mul(r27341, r27340, r27340, MPFR_RNDN);
        ;
        mpfr_fma(r27343, r27339, r27341, r27342, MPFR_RNDN);
        ;
        mpfr_neg(r27345, r27344, MPFR_RNDN);
        mpfr_mul(r27346, r27341, r27345, MPFR_RNDN);
        mpfr_mul(r27347, r27341, r27344, MPFR_RNDN);
        mpfr_fma(r27348, r27346, r27347, r27342, MPFR_RNDN);
        ;
        mpfr_mul(r27350, r27348, r27349, MPFR_RNDN);
        mpfr_sqrt(r27351, r27350, MPFR_RNDN);
        mpfr_div(r27352, r27343, r27351, MPFR_RNDN);
        mpfr_pow(r27353, r27340, r27344, MPFR_RNDN);
        mpfr_neg(r27354, r27340, MPFR_RNDN);
        mpfr_fma(r27355, r27353, r27354, r27342, MPFR_RNDN);
        mpfr_set_d(r27356, t, MPFR_RNDN);
        mpfr_const_pi(r27357, MPFR_RNDN);
        mpfr_mul(r27358, r27356, r27357, MPFR_RNDN);
        mpfr_mul(r27359, r27355, r27358, MPFR_RNDN);
        mpfr_div(r27360, r27352, r27359, MPFR_RNDN);
        mpfr_mul(r27361, r27344, r27341, MPFR_RNDN);
        mpfr_add(r27362, r27342, r27361, MPFR_RNDN);
        mpfr_sqrt(r27363, r27362, MPFR_RNDN);
        mpfr_add(r27364, r27342, r27341, MPFR_RNDN);
        mpfr_mul(r27365, r27363, r27364, MPFR_RNDN);
        mpfr_mul(r27366, r27360, r27365, MPFR_RNDN);
        return mpfr_get_d(r27366, MPFR_RNDN);
}