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

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

double f_if(float re, float im, float base) {
        float r9061 = re;
        float r9062 = r9061 * r9061;
        float r9063 = im;
        float r9064 = r9063 * r9063;
        float r9065 = r9062 + r9064;
        float r9066 = sqrt(r9065);
        float r9067 = log(r9066);
        float r9068 = base;
        float r9069 = log(r9068);
        float r9070 = r9067 * r9069;
        float r9071 = atan2(r9063, r9061);
        float r9072 = 0;
        float r9073 = r9071 * r9072;
        float r9074 = r9070 + r9073;
        float r9075 = r9069 * r9069;
        float r9076 = r9072 * r9072;
        float r9077 = r9075 + r9076;
        float r9078 = r9074 / r9077;
        return r9078;
}

double f_id(double re, double im, double base) {
        double r9079 = re;
        double r9080 = r9079 * r9079;
        double r9081 = im;
        double r9082 = r9081 * r9081;
        double r9083 = r9080 + r9082;
        double r9084 = sqrt(r9083);
        double r9085 = log(r9084);
        double r9086 = base;
        double r9087 = log(r9086);
        double r9088 = r9085 * r9087;
        double r9089 = atan2(r9081, r9079);
        double r9090 = 0;
        double r9091 = r9089 * r9090;
        double r9092 = r9088 + r9091;
        double r9093 = r9087 * r9087;
        double r9094 = r9090 * r9090;
        double r9095 = r9093 + r9094;
        double r9096 = r9092 / r9095;
        return r9096;
}


double f_of(float re, float im, float base) {
        float r9097 = re;
        float r9098 = im;
        float r9099 = hypot(r9097, r9098);
        float r9100 = log(r9099);
        float r9101 = base;
        float r9102 = log(r9101);
        float r9103 = r9100 / r9102;
        float r9104 = 3;
        float r9105 = pow(r9103, r9104);
        float r9106 = cbrt(r9105);
        return r9106;
}

double f_od(double re, double im, double base) {
        double r9107 = re;
        double r9108 = im;
        double r9109 = hypot(r9107, r9108);
        double r9110 = log(r9109);
        double r9111 = base;
        double r9112 = log(r9111);
        double r9113 = r9110 / r9112;
        double r9114 = 3;
        double r9115 = pow(r9113, r9114);
        double r9116 = cbrt(r9115);
        return r9116;
}

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 r9117, r9118, r9119, r9120, r9121, r9122, r9123, r9124, r9125, r9126, r9127, r9128, r9129, r9130, r9131, r9132, r9133, r9134;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9117);
        mpfr_init(r9118);
        mpfr_init(r9119);
        mpfr_init(r9120);
        mpfr_init(r9121);
        mpfr_init(r9122);
        mpfr_init(r9123);
        mpfr_init(r9124);
        mpfr_init(r9125);
        mpfr_init(r9126);
        mpfr_init(r9127);
        mpfr_init_set_str(r9128, "0", 10, MPFR_RNDN);
        mpfr_init(r9129);
        mpfr_init(r9130);
        mpfr_init(r9131);
        mpfr_init(r9132);
        mpfr_init(r9133);
        mpfr_init(r9134);
}

double f_im(double re, double im, double base) {
        mpfr_set_d(r9117, re, MPFR_RNDN);
        mpfr_mul(r9118, r9117, r9117, MPFR_RNDN);
        mpfr_set_d(r9119, im, MPFR_RNDN);
        mpfr_mul(r9120, r9119, r9119, MPFR_RNDN);
        mpfr_add(r9121, r9118, r9120, MPFR_RNDN);
        mpfr_sqrt(r9122, r9121, MPFR_RNDN);
        mpfr_log(r9123, r9122, MPFR_RNDN);
        mpfr_set_d(r9124, base, MPFR_RNDN);
        mpfr_log(r9125, r9124, MPFR_RNDN);
        mpfr_mul(r9126, r9123, r9125, MPFR_RNDN);
        mpfr_atan2(r9127, r9119, r9117, MPFR_RNDN);
        ;
        mpfr_mul(r9129, r9127, r9128, MPFR_RNDN);
        mpfr_add(r9130, r9126, r9129, MPFR_RNDN);
        mpfr_mul(r9131, r9125, r9125, MPFR_RNDN);
        mpfr_mul(r9132, r9128, r9128, MPFR_RNDN);
        mpfr_add(r9133, r9131, r9132, MPFR_RNDN);
        mpfr_div(r9134, r9130, r9133, MPFR_RNDN);
        return mpfr_get_d(r9134, MPFR_RNDN);
}

static mpfr_t r9135, r9136, r9137, r9138, r9139, r9140, r9141, r9142, r9143, r9144;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9135);
        mpfr_init(r9136);
        mpfr_init(r9137);
        mpfr_init(r9138);
        mpfr_init(r9139);
        mpfr_init(r9140);
        mpfr_init(r9141);
        mpfr_init_set_str(r9142, "3", 10, MPFR_RNDN);
        mpfr_init(r9143);
        mpfr_init(r9144);
}

double f_fm(double re, double im, double base) {
        mpfr_set_d(r9135, re, MPFR_RNDN);
        mpfr_set_d(r9136, im, MPFR_RNDN);
        mpfr_hypot(r9137, r9135, r9136, MPFR_RNDN);
        mpfr_log(r9138, r9137, MPFR_RNDN);
        mpfr_set_d(r9139, base, MPFR_RNDN);
        mpfr_log(r9140, r9139, MPFR_RNDN);
        mpfr_div(r9141, r9138, r9140, MPFR_RNDN);
        ;
        mpfr_pow(r9143, r9141, r9142, MPFR_RNDN);
        mpfr_cbrt(r9144, r9143, MPFR_RNDN);
        return mpfr_get_d(r9144, MPFR_RNDN);
}

static mpfr_t r9145, r9146, r9147, r9148, r9149, r9150, r9151, r9152, r9153, r9154;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9145);
        mpfr_init(r9146);
        mpfr_init(r9147);
        mpfr_init(r9148);
        mpfr_init(r9149);
        mpfr_init(r9150);
        mpfr_init(r9151);
        mpfr_init_set_str(r9152, "3", 10, MPFR_RNDN);
        mpfr_init(r9153);
        mpfr_init(r9154);
}

double f_dm(double re, double im, double base) {
        mpfr_set_d(r9145, re, MPFR_RNDN);
        mpfr_set_d(r9146, im, MPFR_RNDN);
        mpfr_hypot(r9147, r9145, r9146, MPFR_RNDN);
        mpfr_log(r9148, r9147, MPFR_RNDN);
        mpfr_set_d(r9149, base, MPFR_RNDN);
        mpfr_log(r9150, r9149, MPFR_RNDN);
        mpfr_div(r9151, r9148, r9150, MPFR_RNDN);
        ;
        mpfr_pow(r9153, r9151, r9152, MPFR_RNDN);
        mpfr_cbrt(r9154, r9153, MPFR_RNDN);
        return mpfr_get_d(r9154, MPFR_RNDN);
}