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

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

double f_if(float re, float im, float base) {
        float r20672 = im;
        float r20673 = re;
        float r20674 = atan2(r20672, r20673);
        float r20675 = base;
        float r20676 = log(r20675);
        float r20677 = r20674 * r20676;
        float r20678 = r20673 * r20673;
        float r20679 = r20672 * r20672;
        float r20680 = r20678 + r20679;
        float r20681 = sqrt(r20680);
        float r20682 = log(r20681);
        float r20683 = 0;
        float r20684 = r20682 * r20683;
        float r20685 = r20677 - r20684;
        float r20686 = r20676 * r20676;
        float r20687 = r20683 * r20683;
        float r20688 = r20686 + r20687;
        float r20689 = r20685 / r20688;
        return r20689;
}

double f_id(double re, double im, double base) {
        double r20690 = im;
        double r20691 = re;
        double r20692 = atan2(r20690, r20691);
        double r20693 = base;
        double r20694 = log(r20693);
        double r20695 = r20692 * r20694;
        double r20696 = r20691 * r20691;
        double r20697 = r20690 * r20690;
        double r20698 = r20696 + r20697;
        double r20699 = sqrt(r20698);
        double r20700 = log(r20699);
        double r20701 = 0;
        double r20702 = r20700 * r20701;
        double r20703 = r20695 - r20702;
        double r20704 = r20694 * r20694;
        double r20705 = r20701 * r20701;
        double r20706 = r20704 + r20705;
        double r20707 = r20703 / r20706;
        return r20707;
}


double f_of(float re, float im, float base) {
        float r20708 = im;
        float r20709 = re;
        float r20710 = atan2(r20708, r20709);
        float r20711 = 0;
        float r20712 = r20710 - r20711;
        float r20713 = 1;
        float r20714 = base;
        float r20715 = log(r20714);
        float r20716 = r20713 / r20715;
        float r20717 = r20712 * r20716;
        return r20717;
}

double f_od(double re, double im, double base) {
        double r20718 = im;
        double r20719 = re;
        double r20720 = atan2(r20718, r20719);
        double r20721 = 0;
        double r20722 = r20720 - r20721;
        double r20723 = 1;
        double r20724 = base;
        double r20725 = log(r20724);
        double r20726 = r20723 / r20725;
        double r20727 = r20722 * r20726;
        return r20727;
}

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 r20728, r20729, r20730, r20731, r20732, r20733, r20734, r20735, r20736, r20737, r20738, r20739, r20740, r20741, r20742, r20743, r20744, r20745;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r20728);
        mpfr_init(r20729);
        mpfr_init(r20730);
        mpfr_init(r20731);
        mpfr_init(r20732);
        mpfr_init(r20733);
        mpfr_init(r20734);
        mpfr_init(r20735);
        mpfr_init(r20736);
        mpfr_init(r20737);
        mpfr_init(r20738);
        mpfr_init_set_str(r20739, "0", 10, MPFR_RNDN);
        mpfr_init(r20740);
        mpfr_init(r20741);
        mpfr_init(r20742);
        mpfr_init(r20743);
        mpfr_init(r20744);
        mpfr_init(r20745);
}

double f_im(double re, double im, double base) {
        mpfr_set_d(r20728, im, MPFR_RNDN);
        mpfr_set_d(r20729, re, MPFR_RNDN);
        mpfr_atan2(r20730, r20728, r20729, MPFR_RNDN);
        mpfr_set_d(r20731, base, MPFR_RNDN);
        mpfr_log(r20732, r20731, MPFR_RNDN);
        mpfr_mul(r20733, r20730, r20732, MPFR_RNDN);
        mpfr_mul(r20734, r20729, r20729, MPFR_RNDN);
        mpfr_mul(r20735, r20728, r20728, MPFR_RNDN);
        mpfr_add(r20736, r20734, r20735, MPFR_RNDN);
        mpfr_sqrt(r20737, r20736, MPFR_RNDN);
        mpfr_log(r20738, r20737, MPFR_RNDN);
        ;
        mpfr_mul(r20740, r20738, r20739, MPFR_RNDN);
        mpfr_sub(r20741, r20733, r20740, MPFR_RNDN);
        mpfr_mul(r20742, r20732, r20732, MPFR_RNDN);
        mpfr_mul(r20743, r20739, r20739, MPFR_RNDN);
        mpfr_add(r20744, r20742, r20743, MPFR_RNDN);
        mpfr_div(r20745, r20741, r20744, MPFR_RNDN);
        return mpfr_get_d(r20745, MPFR_RNDN);
}

static mpfr_t r20746, r20747, r20748, r20749, r20750, r20751, r20752, r20753, r20754, r20755;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20746);
        mpfr_init(r20747);
        mpfr_init(r20748);
        mpfr_init_set_str(r20749, "0", 10, MPFR_RNDN);
        mpfr_init(r20750);
        mpfr_init_set_str(r20751, "1", 10, MPFR_RNDN);
        mpfr_init(r20752);
        mpfr_init(r20753);
        mpfr_init(r20754);
        mpfr_init(r20755);
}

double f_fm(double re, double im, double base) {
        mpfr_set_d(r20746, im, MPFR_RNDN);
        mpfr_set_d(r20747, re, MPFR_RNDN);
        mpfr_atan2(r20748, r20746, r20747, MPFR_RNDN);
        ;
        mpfr_sub(r20750, r20748, r20749, MPFR_RNDN);
        ;
        mpfr_set_d(r20752, base, MPFR_RNDN);
        mpfr_log(r20753, r20752, MPFR_RNDN);
        mpfr_div(r20754, r20751, r20753, MPFR_RNDN);
        mpfr_mul(r20755, r20750, r20754, MPFR_RNDN);
        return mpfr_get_d(r20755, MPFR_RNDN);
}

static mpfr_t r20756, r20757, r20758, r20759, r20760, r20761, r20762, r20763, r20764, r20765;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20756);
        mpfr_init(r20757);
        mpfr_init(r20758);
        mpfr_init_set_str(r20759, "0", 10, MPFR_RNDN);
        mpfr_init(r20760);
        mpfr_init_set_str(r20761, "1", 10, MPFR_RNDN);
        mpfr_init(r20762);
        mpfr_init(r20763);
        mpfr_init(r20764);
        mpfr_init(r20765);
}

double f_dm(double re, double im, double base) {
        mpfr_set_d(r20756, im, MPFR_RNDN);
        mpfr_set_d(r20757, re, MPFR_RNDN);
        mpfr_atan2(r20758, r20756, r20757, MPFR_RNDN);
        ;
        mpfr_sub(r20760, r20758, r20759, MPFR_RNDN);
        ;
        mpfr_set_d(r20762, base, MPFR_RNDN);
        mpfr_log(r20763, r20762, MPFR_RNDN);
        mpfr_div(r20764, r20761, r20763, MPFR_RNDN);
        mpfr_mul(r20765, r20760, r20764, MPFR_RNDN);
        return mpfr_get_d(r20765, MPFR_RNDN);
}