double f(double x_re, double x_im, double y_re, double y_im) {
double r1056531 = x_re;
double r1056532 = r1056531 * r1056531;
double r1056533 = x_im;
double r1056534 = r1056533 * r1056533;
double r1056535 = r1056532 + r1056534;
double r1056536 = sqrt(r1056535);
double r1056537 = log(r1056536);
double r1056538 = y_re;
double r1056539 = r1056537 * r1056538;
double r1056540 = atan2(r1056533, r1056531);
double r1056541 = y_im;
double r1056542 = r1056540 * r1056541;
double r1056543 = r1056539 - r1056542;
double r1056544 = exp(r1056543);
double r1056545 = r1056537 * r1056541;
double r1056546 = r1056540 * r1056538;
double r1056547 = r1056545 + r1056546;
double r1056548 = sin(r1056547);
double r1056549 = r1056544 * r1056548;
return r1056549;
}