#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 r27039 = 2;
        float r27040 = atan2(1.0, 0.0);
        float r27041 = r27039 * r27040;
        float r27042 = 3;
        float r27043 = r27041 / r27042;
        float r27044 = g;
        float r27045 = -r27044;
        float r27046 = h;
        float r27047 = r27045 / r27046;
        float r27048 = acos(r27047);
        float r27049 = r27048 / r27042;
        float r27050 = r27043 + r27049;
        float r27051 = cos(r27050);
        float r27052 = r27039 * r27051;
        return r27052;
}

double f_id(double g, double h) {
        double r27053 = 2;
        double r27054 = atan2(1.0, 0.0);
        double r27055 = r27053 * r27054;
        double r27056 = 3;
        double r27057 = r27055 / r27056;
        double r27058 = g;
        double r27059 = -r27058;
        double r27060 = h;
        double r27061 = r27059 / r27060;
        double r27062 = acos(r27061);
        double r27063 = r27062 / r27056;
        double r27064 = r27057 + r27063;
        double r27065 = cos(r27064);
        double r27066 = r27053 * r27065;
        return r27066;
}


double f_of(float g, float h) {
        float r27067 = 2;
        float r27068 = g;
        float r27069 = h;
        float r27070 = r27068 / r27069;
        float r27071 = acos(r27070);
        float r27072 = 3;
        float r27073 = r27071 / r27072;
        float r27074 = cos(r27073);
        float r27075 = 1;
        float r27076 = r27067 + r27075;
        float r27077 = atan2(1.0, 0.0);
        float r27078 = r27077 / r27072;
        float r27079 = r27076 * r27078;
        float r27080 = cos(r27079);
        float r27081 = r27074 * r27080;
        float r27082 = sin(r27073);
        float r27083 = sin(r27079);
        float r27084 = r27082 * r27083;
        float r27085 = r27081 + r27084;
        float r27086 = r27067 * r27085;
        return r27086;
}

double f_od(double g, double h) {
        double r27087 = 2;
        double r27088 = g;
        double r27089 = h;
        double r27090 = r27088 / r27089;
        double r27091 = acos(r27090);
        double r27092 = 3;
        double r27093 = r27091 / r27092;
        double r27094 = cos(r27093);
        double r27095 = 1;
        double r27096 = r27087 + r27095;
        double r27097 = atan2(1.0, 0.0);
        double r27098 = r27097 / r27092;
        double r27099 = r27096 * r27098;
        double r27100 = cos(r27099);
        double r27101 = r27094 * r27100;
        double r27102 = sin(r27093);
        double r27103 = sin(r27099);
        double r27104 = r27102 * r27103;
        double r27105 = r27101 + r27104;
        double r27106 = r27087 * r27105;
        return r27106;
}

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 r27107, r27108, r27109, r27110, r27111, r27112, r27113, r27114, r27115, r27116, r27117, r27118, r27119, r27120;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r27107, "2", 10, MPFR_RNDN);
        mpfr_init(r27108);
        mpfr_init(r27109);
        mpfr_init_set_str(r27110, "3", 10, MPFR_RNDN);
        mpfr_init(r27111);
        mpfr_init(r27112);
        mpfr_init(r27113);
        mpfr_init(r27114);
        mpfr_init(r27115);
        mpfr_init(r27116);
        mpfr_init(r27117);
        mpfr_init(r27118);
        mpfr_init(r27119);
        mpfr_init(r27120);
}

double f_im(double g, double h) {
        ;
        mpfr_const_pi(r27108, MPFR_RNDN);
        mpfr_mul(r27109, r27107, r27108, MPFR_RNDN);
        ;
        mpfr_div(r27111, r27109, r27110, MPFR_RNDN);
        mpfr_set_d(r27112, g, MPFR_RNDN);
        mpfr_neg(r27113, r27112, MPFR_RNDN);
        mpfr_set_d(r27114, h, MPFR_RNDN);
        mpfr_div(r27115, r27113, r27114, MPFR_RNDN);
        mpfr_acos(r27116, r27115, MPFR_RNDN);
        mpfr_div(r27117, r27116, r27110, MPFR_RNDN);
        mpfr_add(r27118, r27111, r27117, MPFR_RNDN);
        mpfr_cos(r27119, r27118, MPFR_RNDN);
        mpfr_mul(r27120, r27107, r27119, MPFR_RNDN);
        return mpfr_get_d(r27120, MPFR_RNDN);
}

static mpfr_t r27121, r27122, r27123, r27124, r27125, r27126, r27127, r27128, r27129, r27130, r27131, r27132, r27133, r27134, r27135, r27136, r27137, r27138, r27139, r27140;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r27121, "2", 10, MPFR_RNDN);
        mpfr_init(r27122);
        mpfr_init(r27123);
        mpfr_init(r27124);
        mpfr_init(r27125);
        mpfr_init_set_str(r27126, "3", 10, MPFR_RNDN);
        mpfr_init(r27127);
        mpfr_init(r27128);
        mpfr_init_set_str(r27129, "1", 10, MPFR_RNDN);
        mpfr_init(r27130);
        mpfr_init(r27131);
        mpfr_init(r27132);
        mpfr_init(r27133);
        mpfr_init(r27134);
        mpfr_init(r27135);
        mpfr_init(r27136);
        mpfr_init(r27137);
        mpfr_init(r27138);
        mpfr_init(r27139);
        mpfr_init(r27140);
}

double f_fm(double g, double h) {
        ;
        mpfr_set_d(r27122, g, MPFR_RNDN);
        mpfr_set_d(r27123, h, MPFR_RNDN);
        mpfr_div(r27124, r27122, r27123, MPFR_RNDN);
        mpfr_acos(r27125, r27124, MPFR_RNDN);
        ;
        mpfr_div(r27127, r27125, r27126, MPFR_RNDN);
        mpfr_cos(r27128, r27127, MPFR_RNDN);
        ;
        mpfr_add(r27130, r27121, r27129, MPFR_RNDN);
        mpfr_const_pi(r27131, MPFR_RNDN);
        mpfr_div(r27132, r27131, r27126, MPFR_RNDN);
        mpfr_mul(r27133, r27130, r27132, MPFR_RNDN);
        mpfr_cos(r27134, r27133, MPFR_RNDN);
        mpfr_mul(r27135, r27128, r27134, MPFR_RNDN);
        mpfr_sin(r27136, r27127, MPFR_RNDN);
        mpfr_sin(r27137, r27133, MPFR_RNDN);
        mpfr_mul(r27138, r27136, r27137, MPFR_RNDN);
        mpfr_add(r27139, r27135, r27138, MPFR_RNDN);
        mpfr_mul(r27140, r27121, r27139, MPFR_RNDN);
        return mpfr_get_d(r27140, MPFR_RNDN);
}

static mpfr_t r27141, r27142, r27143, r27144, r27145, r27146, r27147, r27148, r27149, r27150, r27151, r27152, r27153, r27154, r27155, r27156, r27157, r27158, r27159, r27160;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r27141, "2", 10, MPFR_RNDN);
        mpfr_init(r27142);
        mpfr_init(r27143);
        mpfr_init(r27144);
        mpfr_init(r27145);
        mpfr_init_set_str(r27146, "3", 10, MPFR_RNDN);
        mpfr_init(r27147);
        mpfr_init(r27148);
        mpfr_init_set_str(r27149, "1", 10, MPFR_RNDN);
        mpfr_init(r27150);
        mpfr_init(r27151);
        mpfr_init(r27152);
        mpfr_init(r27153);
        mpfr_init(r27154);
        mpfr_init(r27155);
        mpfr_init(r27156);
        mpfr_init(r27157);
        mpfr_init(r27158);
        mpfr_init(r27159);
        mpfr_init(r27160);
}

double f_dm(double g, double h) {
        ;
        mpfr_set_d(r27142, g, MPFR_RNDN);
        mpfr_set_d(r27143, h, MPFR_RNDN);
        mpfr_div(r27144, r27142, r27143, MPFR_RNDN);
        mpfr_acos(r27145, r27144, MPFR_RNDN);
        ;
        mpfr_div(r27147, r27145, r27146, MPFR_RNDN);
        mpfr_cos(r27148, r27147, MPFR_RNDN);
        ;
        mpfr_add(r27150, r27141, r27149, MPFR_RNDN);
        mpfr_const_pi(r27151, MPFR_RNDN);
        mpfr_div(r27152, r27151, r27146, MPFR_RNDN);
        mpfr_mul(r27153, r27150, r27152, MPFR_RNDN);
        mpfr_cos(r27154, r27153, MPFR_RNDN);
        mpfr_mul(r27155, r27148, r27154, MPFR_RNDN);
        mpfr_sin(r27156, r27147, MPFR_RNDN);
        mpfr_sin(r27157, r27153, MPFR_RNDN);
        mpfr_mul(r27158, r27156, r27157, MPFR_RNDN);
        mpfr_add(r27159, r27155, r27158, MPFR_RNDN);
        mpfr_mul(r27160, r27141, r27159, MPFR_RNDN);
        return mpfr_get_d(r27160, MPFR_RNDN);
}