FFT: Fix for Intel limitation (AMReX-Codes#4222)

It turns out Intel does not support batched 2D FFT.
JustinRayAngus · Nov 12, 2024 · 2aa22d1 · 2aa22d1
1 parent d6c9caa
commit 2aa22d1
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 3 deletions.
diff --git a/Src/FFT/AMReX_FFT_OpenBCSolver.H b/Src/FFT/AMReX_FFT_OpenBCSolver.H
@@ -3,8 +3,6 @@
 
 #include <AMReX_FFT_R2C.H>
 
-#include <AMReX_VisMF.H>
-
 namespace amrex::FFT
 {
 
@@ -54,6 +52,9 @@ OpenBCSolver<T>::OpenBCSolver (Box const& domain, Info const& info)
 {
 #if (AMREX_SPACEDIM == 3)
     if (m_info.batch_mode) {
+#if defined(AMREX_USE_SYCL)
+        amrex::Abort("FFT::OpenBCSolver: batch mode not support because Intel does not support batched 2d fft.");
+#endif
         auto gdom = make_grown_domain(domain,m_info);
         gdom.enclosedCells(2);
         gdom.setSmall(2, 0);

diff --git a/Src/FFT/AMReX_FFT_R2C.H b/Src/FFT/AMReX_FFT_R2C.H
@@ -246,7 +246,7 @@ R2C<T,D,S>::R2C (Box const& domain, Info const& info)
     int myproc = ParallelContext::MyProcSub();
     int nprocs = std::min(ParallelContext::NProcsSub(), m_info.nprocs);
 
-#if (AMREX_SPACEDIM == 3)
+#if (AMREX_SPACEDIM == 3) && !defined(AMREX_USE_SYCL) /* Intel does not support batched 2d fft */
     if (S == DomainStrategy::slab && (m_real_domain.length(1) > 1)) {
         m_slab_decomp = true;
     }