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

char *name = "2-ancestry mixing, negative discriminant";

double f_if(float g, float h) {
        float r27164 = 2;
        float r27165 = atan2(1.0, 0.0);
        float r27166 = r27164 * r27165;
        float r27167 = 3;
        float r27168 = r27166 / r27167;
        float r27169 = g;
        float r27170 = -r27169;
        float r27171 = h;
        float r27172 = r27170 / r27171;
        float r27173 = acos(r27172);
        float r27174 = r27173 / r27167;
        float r27175 = r27168 + r27174;
        float r27176 = cos(r27175);
        float r27177 = r27164 * r27176;
        return r27177;
}

double f_id(double g, double h) {
        double r27178 = 2;
        double r27179 = atan2(1.0, 0.0);
        double r27180 = r27178 * r27179;
        double r27181 = 3;
        double r27182 = r27180 / r27181;
        double r27183 = g;
        double r27184 = -r27183;
        double r27185 = h;
        double r27186 = r27184 / r27185;
        double r27187 = acos(r27186);
        double r27188 = r27187 / r27181;
        double r27189 = r27182 + r27188;
        double r27190 = cos(r27189);
        double r27191 = r27178 * r27190;
        return r27191;
}


double f_of(float g, float h) {
        float r27192 = 2;
        float r27193 = g;
        float r27194 = h;
        float r27195 = r27193 / r27194;
        float r27196 = acos(r27195);
        float r27197 = 3;
        float r27198 = r27196 / r27197;
        float r27199 = cos(r27198);
        float r27200 = 1;
        float r27201 = r27192 + r27200;
        float r27202 = atan2(1.0, 0.0);
        float r27203 = r27202 / r27197;
        float r27204 = r27201 * r27203;
        float r27205 = cos(r27204);
        float r27206 = r27199 * r27205;
        float r27207 = sin(r27198);
        float r27208 = sin(r27204);
        float r27209 = r27207 * r27208;
        float r27210 = r27206 + r27209;
        float r27211 = r27192 * r27210;
        return r27211;
}

double f_od(double g, double h) {
        double r27212 = 2;
        double r27213 = g;
        double r27214 = h;
        double r27215 = r27213 / r27214;
        double r27216 = acos(r27215);
        double r27217 = 3;
        double r27218 = r27216 / r27217;
        double r27219 = cos(r27218);
        double r27220 = 1;
        double r27221 = r27212 + r27220;
        double r27222 = atan2(1.0, 0.0);
        double r27223 = r27222 / r27217;
        double r27224 = r27221 * r27223;
        double r27225 = cos(r27224);
        double r27226 = r27219 * r27225;
        double r27227 = sin(r27218);
        double r27228 = sin(r27224);
        double r27229 = r27227 * r27228;
        double r27230 = r27226 + r27229;
        double r27231 = r27212 * r27230;
        return r27231;
}

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 r27232, r27233, r27234, r27235, r27236, r27237, r27238, r27239, r27240, r27241, r27242, r27243, r27244, r27245;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r27232, "2", 10, MPFR_RNDN);
        mpfr_init(r27233);
        mpfr_init(r27234);
        mpfr_init_set_str(r27235, "3", 10, MPFR_RNDN);
        mpfr_init(r27236);
        mpfr_init(r27237);
        mpfr_init(r27238);
        mpfr_init(r27239);
        mpfr_init(r27240);
        mpfr_init(r27241);
        mpfr_init(r27242);
        mpfr_init(r27243);
        mpfr_init(r27244);
        mpfr_init(r27245);
}

double f_im(double g, double h) {
        ;
        mpfr_const_pi(r27233, MPFR_RNDN);
        mpfr_mul(r27234, r27232, r27233, MPFR_RNDN);
        ;
        mpfr_div(r27236, r27234, r27235, MPFR_RNDN);
        mpfr_set_d(r27237, g, MPFR_RNDN);
        mpfr_neg(r27238, r27237, MPFR_RNDN);
        mpfr_set_d(r27239, h, MPFR_RNDN);
        mpfr_div(r27240, r27238, r27239, MPFR_RNDN);
        mpfr_acos(r27241, r27240, MPFR_RNDN);
        mpfr_div(r27242, r27241, r27235, MPFR_RNDN);
        mpfr_add(r27243, r27236, r27242, MPFR_RNDN);
        mpfr_cos(r27244, r27243, MPFR_RNDN);
        mpfr_mul(r27245, r27232, r27244, MPFR_RNDN);
        return mpfr_get_d(r27245, MPFR_RNDN);
}

static mpfr_t r27246, r27247, r27248, r27249, r27250, r27251, r27252, r27253, r27254, r27255, r27256, r27257, r27258, r27259, r27260, r27261, r27262, r27263, r27264, r27265;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r27246, "2", 10, MPFR_RNDN);
        mpfr_init(r27247);
        mpfr_init(r27248);
        mpfr_init(r27249);
        mpfr_init(r27250);
        mpfr_init_set_str(r27251, "3", 10, MPFR_RNDN);
        mpfr_init(r27252);
        mpfr_init(r27253);
        mpfr_init_set_str(r27254, "1", 10, MPFR_RNDN);
        mpfr_init(r27255);
        mpfr_init(r27256);
        mpfr_init(r27257);
        mpfr_init(r27258);
        mpfr_init(r27259);
        mpfr_init(r27260);
        mpfr_init(r27261);
        mpfr_init(r27262);
        mpfr_init(r27263);
        mpfr_init(r27264);
        mpfr_init(r27265);
}

double f_fm(double g, double h) {
        ;
        mpfr_set_d(r27247, g, MPFR_RNDN);
        mpfr_set_d(r27248, h, MPFR_RNDN);
        mpfr_div(r27249, r27247, r27248, MPFR_RNDN);
        mpfr_acos(r27250, r27249, MPFR_RNDN);
        ;
        mpfr_div(r27252, r27250, r27251, MPFR_RNDN);
        mpfr_cos(r27253, r27252, MPFR_RNDN);
        ;
        mpfr_add(r27255, r27246, r27254, MPFR_RNDN);
        mpfr_const_pi(r27256, MPFR_RNDN);
        mpfr_div(r27257, r27256, r27251, MPFR_RNDN);
        mpfr_mul(r27258, r27255, r27257, MPFR_RNDN);
        mpfr_cos(r27259, r27258, MPFR_RNDN);
        mpfr_mul(r27260, r27253, r27259, MPFR_RNDN);
        mpfr_sin(r27261, r27252, MPFR_RNDN);
        mpfr_sin(r27262, r27258, MPFR_RNDN);
        mpfr_mul(r27263, r27261, r27262, MPFR_RNDN);
        mpfr_add(r27264, r27260, r27263, MPFR_RNDN);
        mpfr_mul(r27265, r27246, r27264, MPFR_RNDN);
        return mpfr_get_d(r27265, MPFR_RNDN);
}

static mpfr_t r27266, r27267, r27268, r27269, r27270, r27271, r27272, r27273, r27274, r27275, r27276, r27277, r27278, r27279, r27280, r27281, r27282, r27283, r27284, r27285;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r27266, "2", 10, MPFR_RNDN);
        mpfr_init(r27267);
        mpfr_init(r27268);
        mpfr_init(r27269);
        mpfr_init(r27270);
        mpfr_init_set_str(r27271, "3", 10, MPFR_RNDN);
        mpfr_init(r27272);
        mpfr_init(r27273);
        mpfr_init_set_str(r27274, "1", 10, MPFR_RNDN);
        mpfr_init(r27275);
        mpfr_init(r27276);
        mpfr_init(r27277);
        mpfr_init(r27278);
        mpfr_init(r27279);
        mpfr_init(r27280);
        mpfr_init(r27281);
        mpfr_init(r27282);
        mpfr_init(r27283);
        mpfr_init(r27284);
        mpfr_init(r27285);
}

double f_dm(double g, double h) {
        ;
        mpfr_set_d(r27267, g, MPFR_RNDN);
        mpfr_set_d(r27268, h, MPFR_RNDN);
        mpfr_div(r27269, r27267, r27268, MPFR_RNDN);
        mpfr_acos(r27270, r27269, MPFR_RNDN);
        ;
        mpfr_div(r27272, r27270, r27271, MPFR_RNDN);
        mpfr_cos(r27273, r27272, MPFR_RNDN);
        ;
        mpfr_add(r27275, r27266, r27274, MPFR_RNDN);
        mpfr_const_pi(r27276, MPFR_RNDN);
        mpfr_div(r27277, r27276, r27271, MPFR_RNDN);
        mpfr_mul(r27278, r27275, r27277, MPFR_RNDN);
        mpfr_cos(r27279, r27278, MPFR_RNDN);
        mpfr_mul(r27280, r27273, r27279, MPFR_RNDN);
        mpfr_sin(r27281, r27272, MPFR_RNDN);
        mpfr_sin(r27282, r27278, MPFR_RNDN);
        mpfr_mul(r27283, r27281, r27282, MPFR_RNDN);
        mpfr_add(r27284, r27280, r27283, MPFR_RNDN);
        mpfr_mul(r27285, r27266, r27284, MPFR_RNDN);
        return mpfr_get_d(r27285, MPFR_RNDN);
}