Merge pull request #592 from MaxThevenet/crash_if_needed

Check that incompatible options are not used

1 files changed, 19 insertions, 1 deletions

diff --git a/Source/WarpX.cpp b/Source/WarpX.cpp
index 48b4bbd55..f142073d5 100644
--- a/Source/WarpX.cpp
+++ b/Source/WarpX.cpp
@@ -259,6 +259,25 @@ WarpX::WarpX ()
     // at different levels (the stencil depends on c*dt/dz)
     nci_godfrey_filter_exeybz.resize(nlevs_max);
     nci_godfrey_filter_bxbyez.resize(nlevs_max);
+
+    // Sanity checks. Must be done after calling the MultiParticleContainer
+    // constructor, as it reads additional parameters
+    // (e.g., use_fdtd_nci_corr)
+
+#ifndef WARPX_USE_PSATD
+    AMREX_ALWAYS_ASSERT_WITH_MESSAGE(
+        not ( do_pml && do_nodal ),
+        "PML + do_nodal for finite-difference not implemented"
+        );
+#endif
+    AMREX_ALWAYS_ASSERT_WITH_MESSAGE(
+        not ( do_dive_cleaning && do_nodal ),
+        "divE cleaning + do_nodal not implemented"
+        );
+#ifdef WARPX_USE_PSATD
+    AMREX_ALWAYS_ASSERT(use_fdtd_nci_corr == 0);
+    AMREX_ALWAYS_ASSERT(do_subcycling == 0);
+#endif
 }
 
 WarpX::~WarpX ()
@@ -642,7 +661,6 @@ WarpX::ReadParameters ()
        }
 
     }
-
 }
 
 // This is a virtual function.