Adds history_interaction_test.cpp and its path to CMakelists.txt

.clang-format

@@ -19,7 +19,7 @@ AlwaysBreakAfterReturnType: None
 AlwaysBreakBeforeMultilineStrings: true
 AlwaysBreakTemplateDeclarations: true
 BinPackArguments: true
-BinPackParameters: true
+BinPackParameters: false
 BraceWrapping:   
   AfterClass:      false
   AfterControlStatement: false

CMakeLists.txt

@@ -18,7 +18,9 @@ enable_cxx_compiler_flag_if_supported("-Wextra")
 enable_cxx_compiler_flag_if_supported("-fdiagnostics-color=auto")
 enable_cxx_compiler_flag_if_supported("-pedantic")
 
-set(CPP_FEATURES cxx_auto_type cxx_constexpr cxx_delegating_constructors cxx_lambdas cxx_range_for cxx_strong_enums)
+# Require a C++14 standard
+set(CMAKE_CXX_STANDARD 14)
+set(CMAKE_CXX_STANDARD_REQUIRED ON)
 
 add_library(quest "")
 target_compile_features(quest PRIVATE ${CPP_FEATURES})

examples/input.cfg

@@ -1,10 +1,10 @@
 [parameters]
-  num_particles = 2
-  total_time = 20
-  timestep = 0.5e-2
-  interpolation_order = 3
+  num_particles = 2 
+  total_time = 1e-10
+  timestep = 1e-15
+  interpolation_order = 5
 [constants]
-  c0   = 299.792458 # um/ps
-  hbar = 0.65821193 # meV*ps
-  mu0  = 2.0133545e-04 # meV ps^2/(e^2 um)
-  laser_frequency = 2278.9013 # rad/ps
+  c0   = 2.99792458e8 # m/s
+  hbar = 6.62607004e-34 # m^2 kg/s
+  e0  = 8.85418782e-12 #s^4/m^3 A^2/kg
+  laser_frequency = 6.283185307e9 #rad/s = 1GHz

examples/pulse.cfg

@@ -1 +1 @@
-15589.2260227 5 227.89013 2278.9013 0 0 7.6015964 1 0 0
+0000 .5e-10 .2 6.283185307e9 0 1 0 0 0 1 

mp.txt

@@ -0,0 +1,3 @@
+0 1 1 1 2 3 4 5 6
+3 2 1 1 1 1 6 5 4
+4 4 0 2 3 5 1 2 3

src/CMakeLists.txt

@@ -1,22 +1,19 @@
+include("${CMAKE_CURRENT_LIST_DIR}/integrator/CMakeLists.txt")
 include("${CMAKE_CURRENT_LIST_DIR}/interactions/CMakeLists.txt")
 
 target_sources(quest
   PRIVATE
     "${CMAKE_CURRENT_LIST_DIR}/common.h"
     "${CMAKE_CURRENT_LIST_DIR}/configuration.cpp"
     "${CMAKE_CURRENT_LIST_DIR}/configuration.h"
-    "${CMAKE_CURRENT_LIST_DIR}/history.cpp"
-    "${CMAKE_CURRENT_LIST_DIR}/history.h"
-    "${CMAKE_CURRENT_LIST_DIR}/integrator.cpp"
-    "${CMAKE_CURRENT_LIST_DIR}/integrator.h"
     "${CMAKE_CURRENT_LIST_DIR}/lagrange_set.cpp"
     "${CMAKE_CURRENT_LIST_DIR}/lagrange_set.h"
     "${CMAKE_CURRENT_LIST_DIR}/math_utils.cpp"
     "${CMAKE_CURRENT_LIST_DIR}/math_utils.h"
     "${CMAKE_CURRENT_LIST_DIR}/pulse.cpp"
     "${CMAKE_CURRENT_LIST_DIR}/pulse.h"
-    "${CMAKE_CURRENT_LIST_DIR}/quantum_dot.cpp"
-    "${CMAKE_CURRENT_LIST_DIR}/quantum_dot.h"
+    "${CMAKE_CURRENT_LIST_DIR}/magnetic_particle.cpp"
+    "${CMAKE_CURRENT_LIST_DIR}/magnetic_particle.h"
 )
 
 add_executable(quest-bin "")

src/configuration.cpp

@@ -24,7 +24,7 @@ po::variables_map parse_configs(int argc, char *argv[]) {
   constants_description.add_options()
     ("constants.c0",   po::value<double>(&config.c0)->required(), "speed of light in vacuum")
     ("constants.hbar", po::value<double>(&config.hbar)->required(), "reduced Planck constant")
-    ("constants.mu0",  po::value<double>(&config.mu0)->required(), "vacuum permeability")
+    ("constants.e0",  po::value<double>(&config.e0)->required(), "vacuum permittivity")
     ("constants.laser_frequency",  po::value<double>(&config.omega)->required(), "(angular) frequency of incident pulse")
   ;
 

src/configuration.h

@@ -12,7 +12,7 @@
 struct Configuration {
   int num_particles;
   int num_timesteps;
-  double c0, hbar, mu0, omega;
+  double c0, hbar, e0, omega;
   double dt, total_time;
   int interpolation_order;
   std::string qd_path, pulse_path;