double f(double x_re, double x_im, double y_re, double y_im) {
double r1020127 = x_re;
double r1020128 = r1020127 * r1020127;
double r1020129 = x_im;
double r1020130 = r1020129 * r1020129;
double r1020131 = r1020128 + r1020130;
double r1020132 = sqrt(r1020131);
double r1020133 = log(r1020132);
double r1020134 = y_re;
double r1020135 = r1020133 * r1020134;
double r1020136 = atan2(r1020129, r1020127);
double r1020137 = y_im;
double r1020138 = r1020136 * r1020137;
double r1020139 = r1020135 - r1020138;
double r1020140 = exp(r1020139);
double r1020141 = r1020133 * r1020137;
double r1020142 = r1020136 * r1020134;
double r1020143 = r1020141 + r1020142;
double r1020144 = sin(r1020143);
double r1020145 = r1020140 * r1020144;
return r1020145;
}