double f(double x_re, double x_im, double y_re, double y_im) {
double r1027417 = x_re;
double r1027418 = r1027417 * r1027417;
double r1027419 = x_im;
double r1027420 = r1027419 * r1027419;
double r1027421 = r1027418 + r1027420;
double r1027422 = sqrt(r1027421);
double r1027423 = log(r1027422);
double r1027424 = y_re;
double r1027425 = r1027423 * r1027424;
double r1027426 = atan2(r1027419, r1027417);
double r1027427 = y_im;
double r1027428 = r1027426 * r1027427;
double r1027429 = r1027425 - r1027428;
double r1027430 = exp(r1027429);
double r1027431 = r1027423 * r1027427;
double r1027432 = r1027426 * r1027424;
double r1027433 = r1027431 + r1027432;
double r1027434 = cos(r1027433);
double r1027435 = r1027430 * r1027434;
return r1027435;
}