| 41 | |
| 42 | template<int rank> |
| 43 | array fft(const array &in, double norm) { |
| 44 | switch (rank) { |
| 45 | case 1: return fftNorm(in, norm); |
| 46 | case 2: return fft2Norm(in, norm); |
| 47 | case 3: return fft3Norm(in, norm); |
| 48 | default: return in; |
| 49 | } |
| 50 | } |
| 51 | |
| 52 | #define MY_ASSERT_NEAR(aa, bb, cc) ASSERT_NEAR(abs(aa), abs(bb), (cc)) |
| 53 |