| /* PR middle-end/40669 */ |
| /* { dg-do compile } */ |
| /* { dg-options "-O2 -ffast-math" } */ |
| |
| double _Complex |
| test (int d, int t, double *x, double *y, double *z, int n, |
| double _Complex (*fnp) (double)) |
| { |
| int m = n / 2; |
| double min = y[t], max = z[t], med = x[m * d + t]; |
| double _Complex result = 0.0; |
| |
| if (n == 0) |
| return 0.0; |
| |
| if (min > med) |
| result += test (d, (t + 1) % d, x + (m + 1) * d, y, z, n - m - 1, fnp); |
| else if (max < med) |
| result += test (d, (t + 1) % d, x, y, z, m, fnp); |
| else |
| { |
| result += fnp (y[0] + x[m]); |
| result += test (d, (t + 1) % d, x + (m + 1) * d, y, z, n - m - 1, fnp); |
| } |
| return result; |
| } |