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

char *name = "_divideComplex, real part";

double f_if(float x_re, float x_im, float y_re, float y_im) {
        float r22762 = x_re;
        float r22763 = y_re;
        float r22764 = r22762 * r22763;
        float r22765 = x_im;
        float r22766 = y_im;
        float r22767 = r22765 * r22766;
        float r22768 = r22764 + r22767;
        float r22769 = r22763 * r22763;
        float r22770 = r22766 * r22766;
        float r22771 = r22769 + r22770;
        float r22772 = r22768 / r22771;
        return r22772;
}

double f_id(double x_re, double x_im, double y_re, double y_im) {
        double r22773 = x_re;
        double r22774 = y_re;
        double r22775 = r22773 * r22774;
        double r22776 = x_im;
        double r22777 = y_im;
        double r22778 = r22776 * r22777;
        double r22779 = r22775 + r22778;
        double r22780 = r22774 * r22774;
        double r22781 = r22777 * r22777;
        double r22782 = r22780 + r22781;
        double r22783 = r22779 / r22782;
        return r22783;
}


double f_of(float x_re, float x_im, float y_re, float y_im) {
        float r22784 = y_im;
        float r22785 = -7.893629649994627e+166;
        bool r22786 = r22784 <= r22785;
        float r22787 = x_im;
        float r22788 = -r22787;
        float r22789 = y_re;
        float r22790 = hypot(r22789, r22784);
        float r22791 = r22788 / r22790;
        float r22792 = 1.9121082623795437e+143;
        bool r22793 = r22784 <= r22792;
        float r22794 = x_re;
        float r22795 = r22794 * r22789;
        float r22796 = fma(r22784, r22787, r22795);
        float r22797 = r22796 / r22790;
        float r22798 = r22797 / r22790;
        float r22799 = r22787 / r22790;
        float r22800 = r22793 ? r22798 : r22799;
        float r22801 = r22786 ? r22791 : r22800;
        return r22801;
}

double f_od(double x_re, double x_im, double y_re, double y_im) {
        double r22802 = y_im;
        double r22803 = -7.893629649994627e+166;
        bool r22804 = r22802 <= r22803;
        double r22805 = x_im;
        double r22806 = -r22805;
        double r22807 = y_re;
        double r22808 = hypot(r22807, r22802);
        double r22809 = r22806 / r22808;
        double r22810 = 1.9121082623795437e+143;
        bool r22811 = r22802 <= r22810;
        double r22812 = x_re;
        double r22813 = r22812 * r22807;
        double r22814 = fma(r22802, r22805, r22813);
        double r22815 = r22814 / r22808;
        double r22816 = r22815 / r22808;
        double r22817 = r22805 / r22808;
        double r22818 = r22811 ? r22816 : r22817;
        double r22819 = r22804 ? r22809 : r22818;
        return r22819;
}

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 r22820, r22821, r22822, r22823, r22824, r22825, r22826, r22827, r22828, r22829, r22830;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r22820);
        mpfr_init(r22821);
        mpfr_init(r22822);
        mpfr_init(r22823);
        mpfr_init(r22824);
        mpfr_init(r22825);
        mpfr_init(r22826);
        mpfr_init(r22827);
        mpfr_init(r22828);
        mpfr_init(r22829);
        mpfr_init(r22830);
}

double f_im(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r22820, x_re, MPFR_RNDN);
        mpfr_set_d(r22821, y_re, MPFR_RNDN);
        mpfr_mul(r22822, r22820, r22821, MPFR_RNDN);
        mpfr_set_d(r22823, x_im, MPFR_RNDN);
        mpfr_set_d(r22824, y_im, MPFR_RNDN);
        mpfr_mul(r22825, r22823, r22824, MPFR_RNDN);
        mpfr_add(r22826, r22822, r22825, MPFR_RNDN);
        mpfr_mul(r22827, r22821, r22821, MPFR_RNDN);
        mpfr_mul(r22828, r22824, r22824, MPFR_RNDN);
        mpfr_add(r22829, r22827, r22828, MPFR_RNDN);
        mpfr_div(r22830, r22826, r22829, MPFR_RNDN);
        return mpfr_get_d(r22830, MPFR_RNDN);
}

static mpfr_t r22831, r22832, r22833, r22834, r22835, r22836, r22837, r22838, r22839, r22840, r22841, r22842, r22843, r22844, r22845, r22846, r22847, r22848;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r22831);
        mpfr_init_set_str(r22832, "-7.893629649994627e+166", 10, MPFR_RNDN);
        mpfr_init(r22833);
        mpfr_init(r22834);
        mpfr_init(r22835);
        mpfr_init(r22836);
        mpfr_init(r22837);
        mpfr_init(r22838);
        mpfr_init_set_str(r22839, "1.9121082623795437e+143", 10, MPFR_RNDN);
        mpfr_init(r22840);
        mpfr_init(r22841);
        mpfr_init(r22842);
        mpfr_init(r22843);
        mpfr_init(r22844);
        mpfr_init(r22845);
        mpfr_init(r22846);
        mpfr_init(r22847);
        mpfr_init(r22848);
}

double f_fm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r22831, y_im, MPFR_RNDN);
        ;
        mpfr_set_si(r22833, mpfr_cmp(r22831, r22832) <= 0, MPFR_RNDN);
        mpfr_set_d(r22834, x_im, MPFR_RNDN);
        mpfr_neg(r22835, r22834, MPFR_RNDN);
        mpfr_set_d(r22836, y_re, MPFR_RNDN);
        mpfr_hypot(r22837, r22836, r22831, MPFR_RNDN);
        mpfr_div(r22838, r22835, r22837, MPFR_RNDN);
        ;
        mpfr_set_si(r22840, mpfr_cmp(r22831, r22839) <= 0, MPFR_RNDN);
        mpfr_set_d(r22841, x_re, MPFR_RNDN);
        mpfr_mul(r22842, r22841, r22836, MPFR_RNDN);
        mpfr_fma(r22843, r22831, r22834, r22842, MPFR_RNDN);
        mpfr_div(r22844, r22843, r22837, MPFR_RNDN);
        mpfr_div(r22845, r22844, r22837, MPFR_RNDN);
        mpfr_div(r22846, r22834, r22837, MPFR_RNDN);
        if (mpfr_get_si(r22840, MPFR_RNDN)) { mpfr_set(r22847, r22845, MPFR_RNDN); } else { mpfr_set(r22847, r22846, MPFR_RNDN); };
        if (mpfr_get_si(r22833, MPFR_RNDN)) { mpfr_set(r22848, r22838, MPFR_RNDN); } else { mpfr_set(r22848, r22847, MPFR_RNDN); };
        return mpfr_get_d(r22848, MPFR_RNDN);
}

static mpfr_t r22849, r22850, r22851, r22852, r22853, r22854, r22855, r22856, r22857, r22858, r22859, r22860, r22861, r22862, r22863, r22864, r22865, r22866;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r22849);
        mpfr_init_set_str(r22850, "-7.893629649994627e+166", 10, MPFR_RNDN);
        mpfr_init(r22851);
        mpfr_init(r22852);
        mpfr_init(r22853);
        mpfr_init(r22854);
        mpfr_init(r22855);
        mpfr_init(r22856);
        mpfr_init_set_str(r22857, "1.9121082623795437e+143", 10, MPFR_RNDN);
        mpfr_init(r22858);
        mpfr_init(r22859);
        mpfr_init(r22860);
        mpfr_init(r22861);
        mpfr_init(r22862);
        mpfr_init(r22863);
        mpfr_init(r22864);
        mpfr_init(r22865);
        mpfr_init(r22866);
}

double f_dm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r22849, y_im, MPFR_RNDN);
        ;
        mpfr_set_si(r22851, mpfr_cmp(r22849, r22850) <= 0, MPFR_RNDN);
        mpfr_set_d(r22852, x_im, MPFR_RNDN);
        mpfr_neg(r22853, r22852, MPFR_RNDN);
        mpfr_set_d(r22854, y_re, MPFR_RNDN);
        mpfr_hypot(r22855, r22854, r22849, MPFR_RNDN);
        mpfr_div(r22856, r22853, r22855, MPFR_RNDN);
        ;
        mpfr_set_si(r22858, mpfr_cmp(r22849, r22857) <= 0, MPFR_RNDN);
        mpfr_set_d(r22859, x_re, MPFR_RNDN);
        mpfr_mul(r22860, r22859, r22854, MPFR_RNDN);
        mpfr_fma(r22861, r22849, r22852, r22860, MPFR_RNDN);
        mpfr_div(r22862, r22861, r22855, MPFR_RNDN);
        mpfr_div(r22863, r22862, r22855, MPFR_RNDN);
        mpfr_div(r22864, r22852, r22855, MPFR_RNDN);
        if (mpfr_get_si(r22858, MPFR_RNDN)) { mpfr_set(r22865, r22863, MPFR_RNDN); } else { mpfr_set(r22865, r22864, MPFR_RNDN); };
        if (mpfr_get_si(r22851, MPFR_RNDN)) { mpfr_set(r22866, r22856, MPFR_RNDN); } else { mpfr_set(r22866, r22865, MPFR_RNDN); };
        return mpfr_get_d(r22866, MPFR_RNDN);
}