double f(double x_re, double x_im, double y_re, double y_im) {
double r1008923 = x_re;
double r1008924 = r1008923 * r1008923;
double r1008925 = x_im;
double r1008926 = r1008925 * r1008925;
double r1008927 = r1008924 + r1008926;
double r1008928 = sqrt(r1008927);
double r1008929 = log(r1008928);
double r1008930 = y_re;
double r1008931 = r1008929 * r1008930;
double r1008932 = atan2(r1008925, r1008923);
double r1008933 = y_im;
double r1008934 = r1008932 * r1008933;
double r1008935 = r1008931 - r1008934;
double r1008936 = exp(r1008935);
double r1008937 = r1008929 * r1008933;
double r1008938 = r1008932 * r1008930;
double r1008939 = r1008937 + r1008938;
double r1008940 = cos(r1008939);
double r1008941 = r1008936 * r1008940;
return r1008941;
}