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

char *name = "powComplex, imaginary part";

double f_if(float x_re, float x_im, float y_re, float y_im) {
        float r21688 = x_re;
        float r21689 = r21688 * r21688;
        float r21690 = x_im;
        float r21691 = r21690 * r21690;
        float r21692 = r21689 + r21691;
        float r21693 = sqrt(r21692);
        float r21694 = log(r21693);
        float r21695 = y_re;
        float r21696 = r21694 * r21695;
        float r21697 = atan2(r21690, r21688);
        float r21698 = y_im;
        float r21699 = r21697 * r21698;
        float r21700 = r21696 - r21699;
        float r21701 = exp(r21700);
        float r21702 = r21694 * r21698;
        float r21703 = r21697 * r21695;
        float r21704 = r21702 + r21703;
        float r21705 = sin(r21704);
        float r21706 = r21701 * r21705;
        return r21706;
}

double f_id(double x_re, double x_im, double y_re, double y_im) {
        double r21707 = x_re;
        double r21708 = r21707 * r21707;
        double r21709 = x_im;
        double r21710 = r21709 * r21709;
        double r21711 = r21708 + r21710;
        double r21712 = sqrt(r21711);
        double r21713 = log(r21712);
        double r21714 = y_re;
        double r21715 = r21713 * r21714;
        double r21716 = atan2(r21709, r21707);
        double r21717 = y_im;
        double r21718 = r21716 * r21717;
        double r21719 = r21715 - r21718;
        double r21720 = exp(r21719);
        double r21721 = r21713 * r21717;
        double r21722 = r21716 * r21714;
        double r21723 = r21721 + r21722;
        double r21724 = sin(r21723);
        double r21725 = r21720 * r21724;
        return r21725;
}


double f_of(float x_re, float x_im, float y_re, float y_im) {
        float r21726 = x_re;
        float r21727 = r21726 * r21726;
        float r21728 = x_im;
        float r21729 = r21728 * r21728;
        float r21730 = r21727 + r21729;
        float r21731 = sqrt(r21730);
        float r21732 = log(r21731);
        float r21733 = y_re;
        float r21734 = r21732 * r21733;
        float r21735 = atan2(r21728, r21726);
        float r21736 = y_im;
        float r21737 = r21735 * r21736;
        float r21738 = r21734 - r21737;
        float r21739 = exp(r21738);
        float r21740 = log(r21726);
        float r21741 = r21740 * r21736;
        float r21742 = r21735 * r21733;
        float r21743 = r21741 + r21742;
        float r21744 = sin(r21743);
        float r21745 = r21739 * r21744;
        float r21746 = -0.0;
        bool r21747 = r21745 <= r21746;
        float r21748 = +inf.0;
        bool r21749 = r21745 <= r21748;
        float r21750 = cbrt(r21735);
        float r21751 = r21750 * r21750;
        float r21752 = cbrt(r21750);
        float r21753 = r21752 * r21752;
        float r21754 = r21753 * r21752;
        float r21755 = r21754 * r21736;
        float r21756 = r21751 * r21755;
        float r21757 = r21734 - r21756;
        float r21758 = exp(r21757);
        float r21759 = r21732 * r21736;
        float r21760 = r21759 + r21742;
        float r21761 = sin(r21760);
        float r21762 = r21758 * r21761;
        float r21763 = -1;
        float r21764 = r21763 / r21726;
        float r21765 = log(r21764);
        float r21766 = r21765 * r21736;
        float r21767 = r21742 - r21766;
        float r21768 = sin(r21767);
        float r21769 = r21739 * r21768;
        float r21770 = r21749 ? r21762 : r21769;
        float r21771 = r21747 ? r21745 : r21770;
        return r21771;
}

double f_od(double x_re, double x_im, double y_re, double y_im) {
        double r21772 = x_re;
        double r21773 = r21772 * r21772;
        double r21774 = x_im;
        double r21775 = r21774 * r21774;
        double r21776 = r21773 + r21775;
        double r21777 = sqrt(r21776);
        double r21778 = log(r21777);
        double r21779 = y_re;
        double r21780 = r21778 * r21779;
        double r21781 = atan2(r21774, r21772);
        double r21782 = y_im;
        double r21783 = r21781 * r21782;
        double r21784 = r21780 - r21783;
        double r21785 = exp(r21784);
        double r21786 = log(r21772);
        double r21787 = r21786 * r21782;
        double r21788 = r21781 * r21779;
        double r21789 = r21787 + r21788;
        double r21790 = sin(r21789);
        double r21791 = r21785 * r21790;
        double r21792 = -0.0;
        bool r21793 = r21791 <= r21792;
        double r21794 = +inf.0;
        bool r21795 = r21791 <= r21794;
        double r21796 = cbrt(r21781);
        double r21797 = r21796 * r21796;
        double r21798 = cbrt(r21796);
        double r21799 = r21798 * r21798;
        double r21800 = r21799 * r21798;
        double r21801 = r21800 * r21782;
        double r21802 = r21797 * r21801;
        double r21803 = r21780 - r21802;
        double r21804 = exp(r21803);
        double r21805 = r21778 * r21782;
        double r21806 = r21805 + r21788;
        double r21807 = sin(r21806);
        double r21808 = r21804 * r21807;
        double r21809 = -1;
        double r21810 = r21809 / r21772;
        double r21811 = log(r21810);
        double r21812 = r21811 * r21782;
        double r21813 = r21788 - r21812;
        double r21814 = sin(r21813);
        double r21815 = r21785 * r21814;
        double r21816 = r21795 ? r21808 : r21815;
        double r21817 = r21793 ? r21791 : r21816;
        return r21817;
}

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 r21818, r21819, r21820, r21821, r21822, r21823, r21824, r21825, r21826, r21827, r21828, r21829, r21830, r21831, r21832, r21833, r21834, r21835, r21836;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2128);
        mpfr_init(r21818);
        mpfr_init(r21819);
        mpfr_init(r21820);
        mpfr_init(r21821);
        mpfr_init(r21822);
        mpfr_init(r21823);
        mpfr_init(r21824);
        mpfr_init(r21825);
        mpfr_init(r21826);
        mpfr_init(r21827);
        mpfr_init(r21828);
        mpfr_init(r21829);
        mpfr_init(r21830);
        mpfr_init(r21831);
        mpfr_init(r21832);
        mpfr_init(r21833);
        mpfr_init(r21834);
        mpfr_init(r21835);
        mpfr_init(r21836);
}

double f_im(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r21818, x_re, MPFR_RNDN);
        mpfr_mul(r21819, r21818, r21818, MPFR_RNDN);
        mpfr_set_d(r21820, x_im, MPFR_RNDN);
        mpfr_mul(r21821, r21820, r21820, MPFR_RNDN);
        mpfr_add(r21822, r21819, r21821, MPFR_RNDN);
        mpfr_sqrt(r21823, r21822, MPFR_RNDN);
        mpfr_log(r21824, r21823, MPFR_RNDN);
        mpfr_set_d(r21825, y_re, MPFR_RNDN);
        mpfr_mul(r21826, r21824, r21825, MPFR_RNDN);
        mpfr_atan2(r21827, r21820, r21818, MPFR_RNDN);
        mpfr_set_d(r21828, y_im, MPFR_RNDN);
        mpfr_mul(r21829, r21827, r21828, MPFR_RNDN);
        mpfr_sub(r21830, r21826, r21829, MPFR_RNDN);
        mpfr_exp(r21831, r21830, MPFR_RNDN);
        mpfr_mul(r21832, r21824, r21828, MPFR_RNDN);
        mpfr_mul(r21833, r21827, r21825, MPFR_RNDN);
        mpfr_add(r21834, r21832, r21833, MPFR_RNDN);
        mpfr_sin(r21835, r21834, MPFR_RNDN);
        mpfr_mul(r21836, r21831, r21835, MPFR_RNDN);
        return mpfr_get_d(r21836, MPFR_RNDN);
}

static mpfr_t r21837, r21838, r21839, r21840, r21841, r21842, r21843, r21844, r21845, r21846, r21847, r21848, r21849, r21850, r21851, r21852, r21853, r21854, r21855, r21856, r21857, r21858, r21859, r21860, r21861, r21862, r21863, r21864, r21865, r21866, r21867, r21868, r21869, r21870, r21871, r21872, r21873, r21874, r21875, r21876, r21877, r21878, r21879, r21880, r21881, r21882;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2128);
        mpfr_init(r21837);
        mpfr_init(r21838);
        mpfr_init(r21839);
        mpfr_init(r21840);
        mpfr_init(r21841);
        mpfr_init(r21842);
        mpfr_init(r21843);
        mpfr_init(r21844);
        mpfr_init(r21845);
        mpfr_init(r21846);
        mpfr_init(r21847);
        mpfr_init(r21848);
        mpfr_init(r21849);
        mpfr_init(r21850);
        mpfr_init(r21851);
        mpfr_init(r21852);
        mpfr_init(r21853);
        mpfr_init(r21854);
        mpfr_init(r21855);
        mpfr_init(r21856);
        mpfr_init_set_str(r21857, "-0.0", 10, MPFR_RNDN);
        mpfr_init(r21858);
        mpfr_init_set_str(r21859, "+inf.0", 10, MPFR_RNDN);
        mpfr_init(r21860);
        mpfr_init(r21861);
        mpfr_init(r21862);
        mpfr_init(r21863);
        mpfr_init(r21864);
        mpfr_init(r21865);
        mpfr_init(r21866);
        mpfr_init(r21867);
        mpfr_init(r21868);
        mpfr_init(r21869);
        mpfr_init(r21870);
        mpfr_init(r21871);
        mpfr_init(r21872);
        mpfr_init(r21873);
        mpfr_init_set_str(r21874, "-1", 10, MPFR_RNDN);
        mpfr_init(r21875);
        mpfr_init(r21876);
        mpfr_init(r21877);
        mpfr_init(r21878);
        mpfr_init(r21879);
        mpfr_init(r21880);
        mpfr_init(r21881);
        mpfr_init(r21882);
}

double f_fm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r21837, x_re, MPFR_RNDN);
        mpfr_mul(r21838, r21837, r21837, MPFR_RNDN);
        mpfr_set_d(r21839, x_im, MPFR_RNDN);
        mpfr_mul(r21840, r21839, r21839, MPFR_RNDN);
        mpfr_add(r21841, r21838, r21840, MPFR_RNDN);
        mpfr_sqrt(r21842, r21841, MPFR_RNDN);
        mpfr_log(r21843, r21842, MPFR_RNDN);
        mpfr_set_d(r21844, y_re, MPFR_RNDN);
        mpfr_mul(r21845, r21843, r21844, MPFR_RNDN);
        mpfr_atan2(r21846, r21839, r21837, MPFR_RNDN);
        mpfr_set_d(r21847, y_im, MPFR_RNDN);
        mpfr_mul(r21848, r21846, r21847, MPFR_RNDN);
        mpfr_sub(r21849, r21845, r21848, MPFR_RNDN);
        mpfr_exp(r21850, r21849, MPFR_RNDN);
        mpfr_log(r21851, r21837, MPFR_RNDN);
        mpfr_mul(r21852, r21851, r21847, MPFR_RNDN);
        mpfr_mul(r21853, r21846, r21844, MPFR_RNDN);
        mpfr_add(r21854, r21852, r21853, MPFR_RNDN);
        mpfr_sin(r21855, r21854, MPFR_RNDN);
        mpfr_mul(r21856, r21850, r21855, MPFR_RNDN);
        ;
        mpfr_set_si(r21858, mpfr_cmp(r21856, r21857) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r21860, mpfr_cmp(r21856, r21859) <= 0, MPFR_RNDN);
        mpfr_cbrt(r21861, r21846, MPFR_RNDN);
        mpfr_mul(r21862, r21861, r21861, MPFR_RNDN);
        mpfr_cbrt(r21863, r21861, MPFR_RNDN);
        mpfr_mul(r21864, r21863, r21863, MPFR_RNDN);
        mpfr_mul(r21865, r21864, r21863, MPFR_RNDN);
        mpfr_mul(r21866, r21865, r21847, MPFR_RNDN);
        mpfr_mul(r21867, r21862, r21866, MPFR_RNDN);
        mpfr_sub(r21868, r21845, r21867, MPFR_RNDN);
        mpfr_exp(r21869, r21868, MPFR_RNDN);
        mpfr_mul(r21870, r21843, r21847, MPFR_RNDN);
        mpfr_add(r21871, r21870, r21853, MPFR_RNDN);
        mpfr_sin(r21872, r21871, MPFR_RNDN);
        mpfr_mul(r21873, r21869, r21872, MPFR_RNDN);
        ;
        mpfr_div(r21875, r21874, r21837, MPFR_RNDN);
        mpfr_log(r21876, r21875, MPFR_RNDN);
        mpfr_mul(r21877, r21876, r21847, MPFR_RNDN);
        mpfr_sub(r21878, r21853, r21877, MPFR_RNDN);
        mpfr_sin(r21879, r21878, MPFR_RNDN);
        mpfr_mul(r21880, r21850, r21879, MPFR_RNDN);
        if (mpfr_get_si(r21860, MPFR_RNDN)) { mpfr_set(r21881, r21873, MPFR_RNDN); } else { mpfr_set(r21881, r21880, MPFR_RNDN); };
        if (mpfr_get_si(r21858, MPFR_RNDN)) { mpfr_set(r21882, r21856, MPFR_RNDN); } else { mpfr_set(r21882, r21881, MPFR_RNDN); };
        return mpfr_get_d(r21882, MPFR_RNDN);
}

static mpfr_t r21883, r21884, r21885, r21886, r21887, r21888, r21889, r21890, r21891, r21892, r21893, r21894, r21895, r21896, r21897, r21898, r21899, r21900, r21901, r21902, r21903, r21904, r21905, r21906, r21907, r21908, r21909, r21910, r21911, r21912, r21913, r21914, r21915, r21916, r21917, r21918, r21919, r21920, r21921, r21922, r21923, r21924, r21925, r21926, r21927, r21928;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2128);
        mpfr_init(r21883);
        mpfr_init(r21884);
        mpfr_init(r21885);
        mpfr_init(r21886);
        mpfr_init(r21887);
        mpfr_init(r21888);
        mpfr_init(r21889);
        mpfr_init(r21890);
        mpfr_init(r21891);
        mpfr_init(r21892);
        mpfr_init(r21893);
        mpfr_init(r21894);
        mpfr_init(r21895);
        mpfr_init(r21896);
        mpfr_init(r21897);
        mpfr_init(r21898);
        mpfr_init(r21899);
        mpfr_init(r21900);
        mpfr_init(r21901);
        mpfr_init(r21902);
        mpfr_init_set_str(r21903, "-0.0", 10, MPFR_RNDN);
        mpfr_init(r21904);
        mpfr_init_set_str(r21905, "+inf.0", 10, MPFR_RNDN);
        mpfr_init(r21906);
        mpfr_init(r21907);
        mpfr_init(r21908);
        mpfr_init(r21909);
        mpfr_init(r21910);
        mpfr_init(r21911);
        mpfr_init(r21912);
        mpfr_init(r21913);
        mpfr_init(r21914);
        mpfr_init(r21915);
        mpfr_init(r21916);
        mpfr_init(r21917);
        mpfr_init(r21918);
        mpfr_init(r21919);
        mpfr_init_set_str(r21920, "-1", 10, MPFR_RNDN);
        mpfr_init(r21921);
        mpfr_init(r21922);
        mpfr_init(r21923);
        mpfr_init(r21924);
        mpfr_init(r21925);
        mpfr_init(r21926);
        mpfr_init(r21927);
        mpfr_init(r21928);
}

double f_dm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r21883, x_re, MPFR_RNDN);
        mpfr_mul(r21884, r21883, r21883, MPFR_RNDN);
        mpfr_set_d(r21885, x_im, MPFR_RNDN);
        mpfr_mul(r21886, r21885, r21885, MPFR_RNDN);
        mpfr_add(r21887, r21884, r21886, MPFR_RNDN);
        mpfr_sqrt(r21888, r21887, MPFR_RNDN);
        mpfr_log(r21889, r21888, MPFR_RNDN);
        mpfr_set_d(r21890, y_re, MPFR_RNDN);
        mpfr_mul(r21891, r21889, r21890, MPFR_RNDN);
        mpfr_atan2(r21892, r21885, r21883, MPFR_RNDN);
        mpfr_set_d(r21893, y_im, MPFR_RNDN);
        mpfr_mul(r21894, r21892, r21893, MPFR_RNDN);
        mpfr_sub(r21895, r21891, r21894, MPFR_RNDN);
        mpfr_exp(r21896, r21895, MPFR_RNDN);
        mpfr_log(r21897, r21883, MPFR_RNDN);
        mpfr_mul(r21898, r21897, r21893, MPFR_RNDN);
        mpfr_mul(r21899, r21892, r21890, MPFR_RNDN);
        mpfr_add(r21900, r21898, r21899, MPFR_RNDN);
        mpfr_sin(r21901, r21900, MPFR_RNDN);
        mpfr_mul(r21902, r21896, r21901, MPFR_RNDN);
        ;
        mpfr_set_si(r21904, mpfr_cmp(r21902, r21903) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r21906, mpfr_cmp(r21902, r21905) <= 0, MPFR_RNDN);
        mpfr_cbrt(r21907, r21892, MPFR_RNDN);
        mpfr_mul(r21908, r21907, r21907, MPFR_RNDN);
        mpfr_cbrt(r21909, r21907, MPFR_RNDN);
        mpfr_mul(r21910, r21909, r21909, MPFR_RNDN);
        mpfr_mul(r21911, r21910, r21909, MPFR_RNDN);
        mpfr_mul(r21912, r21911, r21893, MPFR_RNDN);
        mpfr_mul(r21913, r21908, r21912, MPFR_RNDN);
        mpfr_sub(r21914, r21891, r21913, MPFR_RNDN);
        mpfr_exp(r21915, r21914, MPFR_RNDN);
        mpfr_mul(r21916, r21889, r21893, MPFR_RNDN);
        mpfr_add(r21917, r21916, r21899, MPFR_RNDN);
        mpfr_sin(r21918, r21917, MPFR_RNDN);
        mpfr_mul(r21919, r21915, r21918, MPFR_RNDN);
        ;
        mpfr_div(r21921, r21920, r21883, MPFR_RNDN);
        mpfr_log(r21922, r21921, MPFR_RNDN);
        mpfr_mul(r21923, r21922, r21893, MPFR_RNDN);
        mpfr_sub(r21924, r21899, r21923, MPFR_RNDN);
        mpfr_sin(r21925, r21924, MPFR_RNDN);
        mpfr_mul(r21926, r21896, r21925, MPFR_RNDN);
        if (mpfr_get_si(r21906, MPFR_RNDN)) { mpfr_set(r21927, r21919, MPFR_RNDN); } else { mpfr_set(r21927, r21926, MPFR_RNDN); };
        if (mpfr_get_si(r21904, MPFR_RNDN)) { mpfr_set(r21928, r21902, MPFR_RNDN); } else { mpfr_set(r21928, r21927, MPFR_RNDN); };
        return mpfr_get_d(r21928, MPFR_RNDN);
}