add automatic convergence testing

Amazingkivas · Feb 21, 2024 · 5f11954 · 5f11954
1 parent 55f9ffb
commit 5f11954
Show file tree

Hide file tree

Showing 5 changed files with 777 additions and 3 deletions.
diff --git a/sln/CMakeLists.txt b/sln/CMakeLists.txt
@@ -2,11 +2,12 @@ cmake_minimum_required(VERSION 2.8)
 
 project(FDTD)
 
-include_directories(../include ../gtest)
+include_directories(../include ../gtest ../test)
 
 
 # BUILD
 add_subdirectory(FDTD)
 add_subdirectory(FDTD_test)
+add_subdirectory(test_FDTD_method)
 add_subdirectory(gtest)
 add_subdirectory(sample)
diff --git a/sln/FDTD_test/CMakeLists.txt b/sln/FDTD_test/CMakeLists.txt
@@ -1,8 +1,8 @@
 set(target FDTD_test)
 
-file(GLOB hdrs "*.h*" "../../include/*.h" "../../gtest/*.h")
+file(GLOB hdrs "*.h*" "../../include/*.h")
 #file(GLOB hdrs "*.h*")
-file(GLOB srcs "*.cpp" "../../src/*.cpp" "../../test/*.cpp")
+file(GLOB srcs "*.cpp" "../../src/*.cpp")
 
 add_executable(${target} ${srcs} ${hdrs})
 target_link_libraries(${target} gtest)

diff --git a/sln/test_FDTD_method/CmakeLists.txt b/sln/test_FDTD_method/CmakeLists.txt
@@ -0,0 +1,8 @@
+set(target test_FDTD_method)
+
+file(GLOB hdrs "*.h*" "../../include/*.h" "../../gtest/*.h")
+#file(GLOB hdrs "*.h*")
+file(GLOB srcs "*.cpp" "../../src/*.cpp" "../../test/*.cpp")
+
+add_executable(${target} ${srcs} ${hdrs})
+target_link_libraries(${target} gtest)
diff --git a/src/FDTD.cpp b/src/FDTD.cpp
@@ -38,6 +38,19 @@ double& Field::operator() (int i, int j, int k)
 
 FDTD::FDTD(Parameters _parameters, double _dt) : parameters(_parameters), dt(_dt)
 {
+    if (parameters.ax >= parameters.bx ||
+        parameters.ay >= parameters.by ||
+        parameters.az >= parameters.bz)
+    {
+        throw std::exception("ERROR: invalid parameters");
+    }
+    if (parameters.Ni <= 0 ||
+        parameters.Nj <= 0 ||
+        parameters.Nk <= 0 ||
+        dt <= 0)
+    {
+        throw std::exception("ERROR: invalid parameters");
+    }
     Ex = Ey = Ez = Bx = By = Bz = Field(parameters.Ni, parameters.Nj, parameters.Nk);
 }
 
@@ -61,6 +74,10 @@ Field& FDTD::get_field(Component this_field)
 
 void FDTD::shifted_update_field(const int time)
 {
+    if (time < 0)
+    {
+        throw std::exception("ERROR: Invalid update field argument");
+    }
     double dx = parameters.dx;
     double dy = parameters.dy;
     double dz = parameters.dz;