Add odd shaped arrays to rfft tests; these fail!

Author: egpbos

test/basic_interface_rfft.cpp

@@ -71,6 +71,50 @@ TYPED_TEST(TransformAndInvert_realFFT, realFFT_nD_n_equals_1_xarray) {
   assert_results(a, a_fourier, should_be_a);
 }
 
+// odd data sizes
+
+TYPED_TEST(TransformAndInvert_realFFT, realFFT_1D_oddsize_xarray) {
+  xt::xarray<TypeParam, xt::layout_type::row_major> a = generate_data<TypeParam, 1>(5);
+  auto a_fourier = xt::fftw::rfft(a);
+  // std::cout << "fourier transform of input before ifft (which is destructive!): " << a_fourier << std::endl;
+  auto should_be_a = xt::fftw::irfft(a_fourier);
+  assert_results(a, a_fourier, should_be_a);
+}
+
+TYPED_TEST(TransformAndInvert_realFFT, realFFT_2D_oddsize_xarray) {
+  xt::xarray<TypeParam, xt::layout_type::row_major> a = generate_data<TypeParam, 2>(5);
+  auto a_fourier = xt::fftw::rfft2(a);
+  // std::cout << "fourier transform of input before ifft (which is destructive!): " << a_fourier << std::endl;
+  auto should_be_a = xt::fftw::irfft2(a_fourier);
+  assert_results(a, a_fourier, should_be_a);
+}
+
+
+TYPED_TEST(TransformAndInvert_realFFT, realFFT_3D_oddsize_xarray) {
+  xt::xarray<TypeParam> a = generate_data<TypeParam, 3>(5);
+  auto a_fourier = xt::fftw::rfft3(a);
+  // std::cout << "fourier transform of input before ifft (which is destructive!): " << a_fourier << std::endl;
+  auto should_be_a = xt::fftw::irfft3(a_fourier);
+  assert_results(a, a_fourier, should_be_a);
+}
+
+
+TYPED_TEST(TransformAndInvert_realFFT, realFFT_nD_n_equals_4_oddsize_xarray) {
+  xt::xarray<TypeParam> a = generate_data<TypeParam, 4>(5);
+  auto a_fourier = xt::fftw::rfftn<4>(a);
+  // std::cout << "fourier transform of input before ifft (which is destructive!): " << a_fourier << std::endl;
+  auto should_be_a = xt::fftw::irfftn<4>(a_fourier);
+  assert_results(a, a_fourier, should_be_a);
+}
+
+TYPED_TEST(TransformAndInvert_realFFT, realFFT_nD_n_equals_1_oddsize_xarray) {
+  xt::xarray<TypeParam> a = generate_data<TypeParam, 1>(5);
+  auto a_fourier = xt::fftw::rfftn<1>(a);
+  // std::cout << "fourier transform of input before ifft (which is destructive!): " << a_fourier << std::endl;
+  auto should_be_a = xt::fftw::irfftn<1>(a_fourier);
+  assert_results(a, a_fourier, should_be_a);
+}
+
 
 /*
 ////