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

char *name = "math.log/1 on complex, real part";

double f_if(float re, float im) {
        float r20922 = re;
        float r20923 = r20922 * r20922;
        float r20924 = im;
        float r20925 = r20924 * r20924;
        float r20926 = r20923 + r20925;
        float r20927 = sqrt(r20926);
        float r20928 = log(r20927);
        return r20928;
}

double f_id(double re, double im) {
        double r20929 = re;
        double r20930 = r20929 * r20929;
        double r20931 = im;
        double r20932 = r20931 * r20931;
        double r20933 = r20930 + r20932;
        double r20934 = sqrt(r20933);
        double r20935 = log(r20934);
        return r20935;
}


double f_of(float re, float im) {
        float r20936 = re;
        float r20937 = -2.130257577176994e+140;
        bool r20938 = r20936 <= r20937;
        float r20939 = -r20936;
        float r20940 = log(r20939);
        float r20941 = 1.1347647040549314e-245;
        bool r20942 = r20936 <= r20941;
        float r20943 = r20936 * r20936;
        float r20944 = im;
        float r20945 = r20944 * r20944;
        float r20946 = r20943 + r20945;
        float r20947 = sqrt(r20946);
        float r20948 = log(r20947);
        float r20949 = 6.107346154734439e-166;
        bool r20950 = r20936 <= r20949;
        float r20951 = log(r20936);
        float r20952 = 1.7956000393992656e+141;
        bool r20953 = r20936 <= r20952;
        float r20954 = r20953 ? r20948 : r20951;
        float r20955 = r20950 ? r20951 : r20954;
        float r20956 = r20942 ? r20948 : r20955;
        float r20957 = r20938 ? r20940 : r20956;
        return r20957;
}

double f_od(double re, double im) {
        double r20958 = re;
        double r20959 = -2.130257577176994e+140;
        bool r20960 = r20958 <= r20959;
        double r20961 = -r20958;
        double r20962 = log(r20961);
        double r20963 = 1.1347647040549314e-245;
        bool r20964 = r20958 <= r20963;
        double r20965 = r20958 * r20958;
        double r20966 = im;
        double r20967 = r20966 * r20966;
        double r20968 = r20965 + r20967;
        double r20969 = sqrt(r20968);
        double r20970 = log(r20969);
        double r20971 = 6.107346154734439e-166;
        bool r20972 = r20958 <= r20971;
        double r20973 = log(r20958);
        double r20974 = 1.7956000393992656e+141;
        bool r20975 = r20958 <= r20974;
        double r20976 = r20975 ? r20970 : r20973;
        double r20977 = r20972 ? r20973 : r20976;
        double r20978 = r20964 ? r20970 : r20977;
        double r20979 = r20960 ? r20962 : r20978;
        return r20979;
}

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 r20980, r20981, r20982, r20983, r20984, r20985, r20986;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r20980);
        mpfr_init(r20981);
        mpfr_init(r20982);
        mpfr_init(r20983);
        mpfr_init(r20984);
        mpfr_init(r20985);
        mpfr_init(r20986);
}

double f_im(double re, double im) {
        mpfr_set_d(r20980, re, MPFR_RNDN);
        mpfr_mul(r20981, r20980, r20980, MPFR_RNDN);
        mpfr_set_d(r20982, im, MPFR_RNDN);
        mpfr_mul(r20983, r20982, r20982, MPFR_RNDN);
        mpfr_add(r20984, r20981, r20983, MPFR_RNDN);
        mpfr_sqrt(r20985, r20984, MPFR_RNDN);
        mpfr_log(r20986, r20985, MPFR_RNDN);
        return mpfr_get_d(r20986, MPFR_RNDN);
}

static mpfr_t r20987, r20988, r20989, r20990, r20991, r20992, r20993, r20994, r20995, r20996, r20997, r20998, r20999, r21000, r21001, r21002, r21003, r21004, r21005, r21006, r21007, r21008;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20987);
        mpfr_init_set_str(r20988, "-2.130257577176994e+140", 10, MPFR_RNDN);
        mpfr_init(r20989);
        mpfr_init(r20990);
        mpfr_init(r20991);
        mpfr_init_set_str(r20992, "1.1347647040549314e-245", 10, MPFR_RNDN);
        mpfr_init(r20993);
        mpfr_init(r20994);
        mpfr_init(r20995);
        mpfr_init(r20996);
        mpfr_init(r20997);
        mpfr_init(r20998);
        mpfr_init(r20999);
        mpfr_init_set_str(r21000, "6.107346154734439e-166", 10, MPFR_RNDN);
        mpfr_init(r21001);
        mpfr_init(r21002);
        mpfr_init_set_str(r21003, "1.7956000393992656e+141", 10, MPFR_RNDN);
        mpfr_init(r21004);
        mpfr_init(r21005);
        mpfr_init(r21006);
        mpfr_init(r21007);
        mpfr_init(r21008);
}

double f_fm(double re, double im) {
        mpfr_set_d(r20987, re, MPFR_RNDN);
        ;
        mpfr_set_si(r20989, mpfr_cmp(r20987, r20988) <= 0, MPFR_RNDN);
        mpfr_neg(r20990, r20987, MPFR_RNDN);
        mpfr_log(r20991, r20990, MPFR_RNDN);
        ;
        mpfr_set_si(r20993, mpfr_cmp(r20987, r20992) <= 0, MPFR_RNDN);
        mpfr_mul(r20994, r20987, r20987, MPFR_RNDN);
        mpfr_set_d(r20995, im, MPFR_RNDN);
        mpfr_mul(r20996, r20995, r20995, MPFR_RNDN);
        mpfr_add(r20997, r20994, r20996, MPFR_RNDN);
        mpfr_sqrt(r20998, r20997, MPFR_RNDN);
        mpfr_log(r20999, r20998, MPFR_RNDN);
        ;
        mpfr_set_si(r21001, mpfr_cmp(r20987, r21000) <= 0, MPFR_RNDN);
        mpfr_log(r21002, r20987, MPFR_RNDN);
        ;
        mpfr_set_si(r21004, mpfr_cmp(r20987, r21003) <= 0, MPFR_RNDN);
        if (mpfr_get_si(r21004, MPFR_RNDN)) { mpfr_set(r21005, r20999, MPFR_RNDN); } else { mpfr_set(r21005, r21002, MPFR_RNDN); };
        if (mpfr_get_si(r21001, MPFR_RNDN)) { mpfr_set(r21006, r21002, MPFR_RNDN); } else { mpfr_set(r21006, r21005, MPFR_RNDN); };
        if (mpfr_get_si(r20993, MPFR_RNDN)) { mpfr_set(r21007, r20999, MPFR_RNDN); } else { mpfr_set(r21007, r21006, MPFR_RNDN); };
        if (mpfr_get_si(r20989, MPFR_RNDN)) { mpfr_set(r21008, r20991, MPFR_RNDN); } else { mpfr_set(r21008, r21007, MPFR_RNDN); };
        return mpfr_get_d(r21008, MPFR_RNDN);
}

static mpfr_t r21009, r21010, r21011, r21012, r21013, r21014, r21015, r21016, r21017, r21018, r21019, r21020, r21021, r21022, r21023, r21024, r21025, r21026, r21027, r21028, r21029, r21030;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r21009);
        mpfr_init_set_str(r21010, "-2.130257577176994e+140", 10, MPFR_RNDN);
        mpfr_init(r21011);
        mpfr_init(r21012);
        mpfr_init(r21013);
        mpfr_init_set_str(r21014, "1.1347647040549314e-245", 10, MPFR_RNDN);
        mpfr_init(r21015);
        mpfr_init(r21016);
        mpfr_init(r21017);
        mpfr_init(r21018);
        mpfr_init(r21019);
        mpfr_init(r21020);
        mpfr_init(r21021);
        mpfr_init_set_str(r21022, "6.107346154734439e-166", 10, MPFR_RNDN);
        mpfr_init(r21023);
        mpfr_init(r21024);
        mpfr_init_set_str(r21025, "1.7956000393992656e+141", 10, MPFR_RNDN);
        mpfr_init(r21026);
        mpfr_init(r21027);
        mpfr_init(r21028);
        mpfr_init(r21029);
        mpfr_init(r21030);
}

double f_dm(double re, double im) {
        mpfr_set_d(r21009, re, MPFR_RNDN);
        ;
        mpfr_set_si(r21011, mpfr_cmp(r21009, r21010) <= 0, MPFR_RNDN);
        mpfr_neg(r21012, r21009, MPFR_RNDN);
        mpfr_log(r21013, r21012, MPFR_RNDN);
        ;
        mpfr_set_si(r21015, mpfr_cmp(r21009, r21014) <= 0, MPFR_RNDN);
        mpfr_mul(r21016, r21009, r21009, MPFR_RNDN);
        mpfr_set_d(r21017, im, MPFR_RNDN);
        mpfr_mul(r21018, r21017, r21017, MPFR_RNDN);
        mpfr_add(r21019, r21016, r21018, MPFR_RNDN);
        mpfr_sqrt(r21020, r21019, MPFR_RNDN);
        mpfr_log(r21021, r21020, MPFR_RNDN);
        ;
        mpfr_set_si(r21023, mpfr_cmp(r21009, r21022) <= 0, MPFR_RNDN);
        mpfr_log(r21024, r21009, MPFR_RNDN);
        ;
        mpfr_set_si(r21026, mpfr_cmp(r21009, r21025) <= 0, MPFR_RNDN);
        if (mpfr_get_si(r21026, MPFR_RNDN)) { mpfr_set(r21027, r21021, MPFR_RNDN); } else { mpfr_set(r21027, r21024, MPFR_RNDN); };
        if (mpfr_get_si(r21023, MPFR_RNDN)) { mpfr_set(r21028, r21024, MPFR_RNDN); } else { mpfr_set(r21028, r21027, MPFR_RNDN); };
        if (mpfr_get_si(r21015, MPFR_RNDN)) { mpfr_set(r21029, r21021, MPFR_RNDN); } else { mpfr_set(r21029, r21028, MPFR_RNDN); };
        if (mpfr_get_si(r21011, MPFR_RNDN)) { mpfr_set(r21030, r21013, MPFR_RNDN); } else { mpfr_set(r21030, r21029, MPFR_RNDN); };
        return mpfr_get_d(r21030, MPFR_RNDN);
}

