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

char *name = "_divideComplex, imaginary part";

double f_if(float x_re, float x_im, float y_re, float y_im) {
        float r22777 = x_im;
        float r22778 = y_re;
        float r22779 = r22777 * r22778;
        float r22780 = x_re;
        float r22781 = y_im;
        float r22782 = r22780 * r22781;
        float r22783 = r22779 - r22782;
        float r22784 = r22778 * r22778;
        float r22785 = r22781 * r22781;
        float r22786 = r22784 + r22785;
        float r22787 = r22783 / r22786;
        return r22787;
}

double f_id(double x_re, double x_im, double y_re, double y_im) {
        double r22788 = x_im;
        double r22789 = y_re;
        double r22790 = r22788 * r22789;
        double r22791 = x_re;
        double r22792 = y_im;
        double r22793 = r22791 * r22792;
        double r22794 = r22790 - r22793;
        double r22795 = r22789 * r22789;
        double r22796 = r22792 * r22792;
        double r22797 = r22795 + r22796;
        double r22798 = r22794 / r22797;
        return r22798;
}


double f_of(float x_re, float x_im, float y_re, float y_im) {
        float r22799 = y_re;
        float r22800 = -9.290821530888773e+139;
        bool r22801 = r22799 <= r22800;
        float r22802 = x_im;
        float r22803 = -r22802;
        float r22804 = y_im;
        float r22805 = hypot(r22799, r22804);
        float r22806 = r22803 / r22805;
        float r22807 = 1.3347106463364611e+187;
        bool r22808 = r22799 <= r22807;
        float r22809 = 1;
        float r22810 = r22809 / r22805;
        float r22811 = x_re;
        float r22812 = r22811 * r22804;
        float r22813 = -r22812;
        float r22814 = fma(r22799, r22802, r22813);
        float r22815 = r22814 / r22805;
        float r22816 = r22810 * r22815;
        float r22817 = r22802 / r22805;
        float r22818 = r22808 ? r22816 : r22817;
        float r22819 = r22801 ? r22806 : r22818;
        return r22819;
}

double f_od(double x_re, double x_im, double y_re, double y_im) {
        double r22820 = y_re;
        double r22821 = -9.290821530888773e+139;
        bool r22822 = r22820 <= r22821;
        double r22823 = x_im;
        double r22824 = -r22823;
        double r22825 = y_im;
        double r22826 = hypot(r22820, r22825);
        double r22827 = r22824 / r22826;
        double r22828 = 1.3347106463364611e+187;
        bool r22829 = r22820 <= r22828;
        double r22830 = 1;
        double r22831 = r22830 / r22826;
        double r22832 = x_re;
        double r22833 = r22832 * r22825;
        double r22834 = -r22833;
        double r22835 = fma(r22820, r22823, r22834);
        double r22836 = r22835 / r22826;
        double r22837 = r22831 * r22836;
        double r22838 = r22823 / r22826;
        double r22839 = r22829 ? r22837 : r22838;
        double r22840 = r22822 ? r22827 : r22839;
        return r22840;
}

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 r22841, r22842, r22843, r22844, r22845, r22846, r22847, r22848, r22849, r22850, r22851;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r22841);
        mpfr_init(r22842);
        mpfr_init(r22843);
        mpfr_init(r22844);
        mpfr_init(r22845);
        mpfr_init(r22846);
        mpfr_init(r22847);
        mpfr_init(r22848);
        mpfr_init(r22849);
        mpfr_init(r22850);
        mpfr_init(r22851);
}

double f_im(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r22841, x_im, MPFR_RNDN);
        mpfr_set_d(r22842, y_re, MPFR_RNDN);
        mpfr_mul(r22843, r22841, r22842, MPFR_RNDN);
        mpfr_set_d(r22844, x_re, MPFR_RNDN);
        mpfr_set_d(r22845, y_im, MPFR_RNDN);
        mpfr_mul(r22846, r22844, r22845, MPFR_RNDN);
        mpfr_sub(r22847, r22843, r22846, MPFR_RNDN);
        mpfr_mul(r22848, r22842, r22842, MPFR_RNDN);
        mpfr_mul(r22849, r22845, r22845, MPFR_RNDN);
        mpfr_add(r22850, r22848, r22849, MPFR_RNDN);
        mpfr_div(r22851, r22847, r22850, MPFR_RNDN);
        return mpfr_get_d(r22851, MPFR_RNDN);
}

static mpfr_t r22852, r22853, r22854, r22855, r22856, r22857, r22858, r22859, r22860, r22861, r22862, r22863, r22864, r22865, r22866, r22867, r22868, r22869, r22870, r22871, r22872;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r22852);
        mpfr_init_set_str(r22853, "-9.290821530888773e+139", 10, MPFR_RNDN);
        mpfr_init(r22854);
        mpfr_init(r22855);
        mpfr_init(r22856);
        mpfr_init(r22857);
        mpfr_init(r22858);
        mpfr_init(r22859);
        mpfr_init_set_str(r22860, "1.3347106463364611e+187", 10, MPFR_RNDN);
        mpfr_init(r22861);
        mpfr_init_set_str(r22862, "1", 10, MPFR_RNDN);
        mpfr_init(r22863);
        mpfr_init(r22864);
        mpfr_init(r22865);
        mpfr_init(r22866);
        mpfr_init(r22867);
        mpfr_init(r22868);
        mpfr_init(r22869);
        mpfr_init(r22870);
        mpfr_init(r22871);
        mpfr_init(r22872);
}

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

static mpfr_t r22873, r22874, r22875, r22876, r22877, r22878, r22879, r22880, r22881, r22882, r22883, r22884, r22885, r22886, r22887, r22888, r22889, r22890, r22891, r22892, r22893;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r22873);
        mpfr_init_set_str(r22874, "-9.290821530888773e+139", 10, MPFR_RNDN);
        mpfr_init(r22875);
        mpfr_init(r22876);
        mpfr_init(r22877);
        mpfr_init(r22878);
        mpfr_init(r22879);
        mpfr_init(r22880);
        mpfr_init_set_str(r22881, "1.3347106463364611e+187", 10, MPFR_RNDN);
        mpfr_init(r22882);
        mpfr_init_set_str(r22883, "1", 10, MPFR_RNDN);
        mpfr_init(r22884);
        mpfr_init(r22885);
        mpfr_init(r22886);
        mpfr_init(r22887);
        mpfr_init(r22888);
        mpfr_init(r22889);
        mpfr_init(r22890);
        mpfr_init(r22891);
        mpfr_init(r22892);
        mpfr_init(r22893);
}

double f_dm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r22873, y_re, MPFR_RNDN);
        ;
        mpfr_set_si(r22875, mpfr_cmp(r22873, r22874) <= 0, MPFR_RNDN);
        mpfr_set_d(r22876, x_im, MPFR_RNDN);
        mpfr_neg(r22877, r22876, MPFR_RNDN);
        mpfr_set_d(r22878, y_im, MPFR_RNDN);
        mpfr_hypot(r22879, r22873, r22878, MPFR_RNDN);
        mpfr_div(r22880, r22877, r22879, MPFR_RNDN);
        ;
        mpfr_set_si(r22882, mpfr_cmp(r22873, r22881) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r22884, r22883, r22879, MPFR_RNDN);
        mpfr_set_d(r22885, x_re, MPFR_RNDN);
        mpfr_mul(r22886, r22885, r22878, MPFR_RNDN);
        mpfr_neg(r22887, r22886, MPFR_RNDN);
        mpfr_fma(r22888, r22873, r22876, r22887, MPFR_RNDN);
        mpfr_div(r22889, r22888, r22879, MPFR_RNDN);
        mpfr_mul(r22890, r22884, r22889, MPFR_RNDN);
        mpfr_div(r22891, r22876, r22879, MPFR_RNDN);
        if (mpfr_get_si(r22882, MPFR_RNDN)) { mpfr_set(r22892, r22890, MPFR_RNDN); } else { mpfr_set(r22892, r22891, MPFR_RNDN); };
        if (mpfr_get_si(r22875, MPFR_RNDN)) { mpfr_set(r22893, r22880, MPFR_RNDN); } else { mpfr_set(r22893, r22892, MPFR_RNDN); };
        return mpfr_get_d(r22893, MPFR_RNDN);
}