| 149 | } |
| 150 | |
| 151 | NOINLINE static void fill_first_quadrant(size_t n, double * restrict res) |
| 152 | { |
| 153 | const double hsqt2 = 0.707106781186547524400844362104849; |
| 154 | size_t quart = n>>2; |
| 155 | if ((n&7)==0) |
| 156 | res[quart] = res[quart+1] = hsqt2; |
| 157 | for (size_t i=2, j=2*quart-2; i<quart; i+=2, j-=2) |
| 158 | { |
| 159 | res[j ] = res[i+1]; |
| 160 | res[j+1] = res[i ]; |
| 161 | } |
| 162 | } |
| 163 | |
| 164 | NOINLINE static void fill_first_half(size_t n, double * restrict res) |
| 165 | { |
no outgoing calls
no test coverage detected