#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 r21551 = re;
        float r21552 = r21551 * r21551;
        float r21553 = im;
        float r21554 = r21553 * r21553;
        float r21555 = r21552 + r21554;
        float r21556 = sqrt(r21555);
        float r21557 = log(r21556);
        return r21557;
}

double f_id(double re, double im) {
        double r21558 = re;
        double r21559 = r21558 * r21558;
        double r21560 = im;
        double r21561 = r21560 * r21560;
        double r21562 = r21559 + r21561;
        double r21563 = sqrt(r21562);
        double r21564 = log(r21563);
        return r21564;
}


double f_of(float re, float im) {
        float r21565 = re;
        float r21566 = im;
        float r21567 = hypot(r21565, r21566);
        float r21568 = log(r21567);
        return r21568;
}

double f_od(double re, double im) {
        double r21569 = re;
        double r21570 = im;
        double r21571 = hypot(r21569, r21570);
        double r21572 = log(r21571);
        return r21572;
}

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 r21573, r21574, r21575, r21576, r21577, r21578, r21579;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r21573);
        mpfr_init(r21574);
        mpfr_init(r21575);
        mpfr_init(r21576);
        mpfr_init(r21577);
        mpfr_init(r21578);
        mpfr_init(r21579);
}

double f_im(double re, double im) {
        mpfr_set_d(r21573, re, MPFR_RNDN);
        mpfr_mul(r21574, r21573, r21573, MPFR_RNDN);
        mpfr_set_d(r21575, im, MPFR_RNDN);
        mpfr_mul(r21576, r21575, r21575, MPFR_RNDN);
        mpfr_add(r21577, r21574, r21576, MPFR_RNDN);
        mpfr_sqrt(r21578, r21577, MPFR_RNDN);
        mpfr_log(r21579, r21578, MPFR_RNDN);
        return mpfr_get_d(r21579, MPFR_RNDN);
}

static mpfr_t r21580, r21581, r21582, r21583;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r21580);
        mpfr_init(r21581);
        mpfr_init(r21582);
        mpfr_init(r21583);
}

double f_fm(double re, double im) {
        mpfr_set_d(r21580, re, MPFR_RNDN);
        mpfr_set_d(r21581, im, MPFR_RNDN);
        mpfr_hypot(r21582, r21580, r21581, MPFR_RNDN);
        mpfr_log(r21583, r21582, MPFR_RNDN);
        return mpfr_get_d(r21583, MPFR_RNDN);
}

static mpfr_t r21584, r21585, r21586, r21587;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r21584);
        mpfr_init(r21585);
        mpfr_init(r21586);
        mpfr_init(r21587);
}

double f_dm(double re, double im) {
        mpfr_set_d(r21584, re, MPFR_RNDN);
        mpfr_set_d(r21585, im, MPFR_RNDN);
        mpfr_hypot(r21586, r21584, r21585, MPFR_RNDN);
        mpfr_log(r21587, r21586, MPFR_RNDN);
        return mpfr_get_d(r21587, MPFR_RNDN);
}