Merge branch 'main' of github.com:milliCoulomb/tetraSN

Author: milliCoulomb

README.md

@@ -157,6 +157,7 @@ solver_parameters:
   max_power_iterations: 500
   rays_per_face: 10
   max_ray_length: 1000
+  use_half_hemisphere: false
 
 output:
   flux_output_file: "../examples/cube/output/flux.dat"

examples/cow/cow.yaml

@@ -21,11 +21,11 @@ solver_parameters:
   max_power_iterations: 5000
   rays_per_face: 10
   max_ray_length: 1000
+  use_half_hemisphere: false
 
 output:
   flux_output_file: "../examples/cow/output/flux.txt"
   k_eff_output_file: "../examples/cow/output/k_eff.txt"
 
 logging:
-  level: "ERROR"
-  log_file: "logs/solver.log"
+  level: "ERROR"

examples/cube/cube.yaml

@@ -21,11 +21,11 @@ solver_parameters:
   max_power_iterations: 300
   rays_per_face: 10
   max_ray_length: 1000
+  use_half_hemisphere: false
 
 output:
   flux_output_file: "../examples/cube/output/flux.txt"
   k_eff_output_file: "../examples/cube/output/k_eff.txt"
 
 logging:
-  level: "RUNNING"
-  log_file: "logs/solver.log"
+  level: "RUNNING"

examples/cylinder/cylinder.yaml

@@ -21,11 +21,11 @@ solver_parameters:
   max_power_iterations: 500
   rays_per_face: 8
   max_ray_length: 5000
+  use_half_hemisphere: false
 
 output:
   flux_output_file: "../examples/sphere/output/flux.txt"
   k_eff_output_file: "../examples/sphere/output/k_eff.txt"
 
 logging:
-  level: "ERROR"
-  log_file: "logs/solver.log"
+  level: "ERROR"

examples/sphere/sphere.yaml

@@ -34,6 +34,7 @@ struct SolverParameters {
     int max_power_iterations;
     int rays_per_face;
     int max_ray_length;
+    bool use_half_hemisphere; // use only the upper hemisphere for ray tracing and symmetrize afterwards
 };
 
 // Structure to hold output settings
@@ -45,7 +46,6 @@ struct OutputSettings {
 // Structure to hold logging settings
 struct LoggingSettings {
     std::string level;
-    std::string log_file;
 };
 
 // Main structure to hold the entire input deck

include/InputDeck.hpp

@@ -43,6 +43,7 @@ InputDeck InputDeckParser::parse(const std::string& filename) {
         deck.solver_parameters.max_power_iterations = sp_node["max_power_iterations"].as<int>();
         deck.solver_parameters.rays_per_face = sp_node["rays_per_face"].as<int>();
         deck.solver_parameters.max_ray_length = sp_node["max_ray_length"].as<int>();
+        deck.solver_parameters.use_half_hemisphere = sp_node["use_half_hemisphere"].as<bool>();
 
         // Parse output settings
         YAML::Node out_node = config["output"];
@@ -57,7 +58,6 @@ InputDeck InputDeckParser::parse(const std::string& filename) {
         YAML::Node log_node = config["logging"];
         if (!log_node) throw std::runtime_error("Missing 'logging' section in YAML.");
         deck.logging.level = log_node["level"].as<std::string>();
-        deck.logging.log_file = log_node["log_file"].as<std::string>();
 
     } catch (const YAML::Exception& e) {
         std::cerr << "YAML Parsing Error: " << e.what() << std::endl;

src/InputDeck.cpp

@@ -67,7 +67,7 @@ int main(int argc, char* argv[]) {
         Field field;
 
         bool constant_dir_bool = true;
-        bool use_half_quadrature_for_constant = false;
+        bool use_half_quadrature_for_constant = input_deck.solver_parameters.use_half_hemisphere;
 
         // Initialize RayTracerManager and generate tracking data
         RayTracerManager ray_tracer_manager(mesh_handler, field, angular_quadrature, constant_dir_bool,
@@ -76,6 +76,10 @@ int main(int argc, char* argv[]) {
         int rays_per_face = input_deck.solver_parameters.rays_per_face;
         int max_ray_length = input_deck.solver_parameters.max_ray_length;
         ray_tracer_manager.generateTrackingData(rays_per_face, max_ray_length);
+        if (use_half_quadrature_for_constant) {
+            Logger::info("Using only the upper hemisphere for ray tracing, symmetrizing afterwards.");
+            ray_tracer_manager.doubleTrackingDataByReversing();
+        }
 
         Logger::running("Tracking data generated.");
 

src/main.cpp

@@ -44,14 +44,14 @@ TEST(InputDeckTest, ParseValidInputDeck) {
   max_power_iterations: 100
   rays_per_face: 8
   max_ray_length: 1000
+  use_half_hemisphere: false
 
 output:
   flux_output_file: "output/flux.dat"
   k_eff_output_file: "output/k_eff.dat"
 
 logging:
   level: "INFO"
-  log_file: "logs/solver.log"
 )";
 
     ASSERT_TRUE(createTempYAML(filename, yaml_content)) << "Failed to create temporary YAML file.";
@@ -88,7 +88,6 @@ TEST(InputDeckTest, ParseValidInputDeck) {
 
     // Verify logging settings
     EXPECT_EQ(deck.logging.level, "INFO");
-    EXPECT_EQ(deck.logging.log_file, "logs/solver.log");
 
     // Clean up
     std::remove(filename.c_str());
@@ -116,14 +115,14 @@ TEST(InputDeckTest, ParseMissingMeshSection) {
   max_power_iterations: 100
   rays_per_face: 8
   max_ray_length: 1000
+  use_half_hemisphere: false
 
 output:
   flux_output_file: "output/flux.dat"
   k_eff_output_file: "output/k_eff.dat"
 
 logging:
   level: "INFO"
-  log_file: "logs/solver.log"
 )";
 
     ASSERT_TRUE(createTempYAML(filename, yaml_content)) << "Failed to create temporary YAML file.";
@@ -158,14 +157,14 @@ TEST(InputDeckTest, ParseMissingCrossSectionsSection) {
   max_power_iterations: 100
   rays_per_face: 8
   max_ray_length: 1000
+  use_half_hemisphere: false
 
 output:
   flux_output_file: "output/flux.dat"
   k_eff_output_file: "output/k_eff.dat"
 
 logging:
   level: "INFO"
-  log_file: "logs/solver.log"
 )";
 
     ASSERT_TRUE(createTempYAML(filename, yaml_content)) << "Failed to create temporary YAML file.";
@@ -204,14 +203,14 @@ TEST(InputDeckTest, ParseInvalidAngularQuadratureData) {
   max_power_iterations: 100
   rays_per_face: 8
   max_ray_length: 1000
+  use_half_hemisphere: false
 
 output:
   flux_output_file: "output/flux.dat"
   k_eff_output_file: "output/k_eff.dat"
 
 logging:
   level: "INFO"
-  log_file: "logs/solver.log"
 )";
 
     ASSERT_TRUE(createTempYAML(filename, yaml_content)) << "Failed to create temporary YAML file.";
@@ -246,7 +245,6 @@ TEST(InputDeckTest, ParseMissingSolverParametersSection) {
 
 logging:
   level: "INFO"
-  log_file: "logs/solver.log"
 )";
 
     ASSERT_TRUE(createTempYAML(filename, yaml_content)) << "Failed to create temporary YAML file.";
@@ -285,14 +283,14 @@ TEST(InputDeckTest, ParseInvalidSolverParametersValues) {
   max_power_iterations: 100
   rays_per_face: 8
   max_ray_length: 1000
+  use_half_hemisphere: false
 
 output:
   flux_output_file: "output/flux.dat"
   k_eff_output_file: "output/k_eff.dat"
 
 logging:
   level: "INFO"
-  log_file: "logs/solver.log"
 )";
 
     ASSERT_TRUE(createTempYAML(filename, yaml_content)) << "Failed to create temporary YAML file.";
@@ -331,10 +329,10 @@ TEST(InputDeckTest, ParseMissingOutputSection) {
   max_power_iterations: 100
   rays_per_face: 8
   max_ray_length: 1000
+  use_half_hemisphere: false
 
 logging:
   level: "INFO"
-  log_file: "logs/solver.log"
 )";
 
     ASSERT_TRUE(createTempYAML(filename, yaml_content)) << "Failed to create temporary YAML file.";
@@ -373,6 +371,7 @@ TEST(InputDeckTest, ParseMissingLoggingSection) {
   max_power_iterations: 100
   rays_per_face: 8
   max_ray_length: 1000
+  use_half_hemisphere: false
 
 output:
   flux_output_file: "output/flux.dat"
@@ -458,14 +457,14 @@ TEST(InputDeckTest, ParseWithExtraSections) {
   max_power_iterations: 100
   rays_per_face: 8
   max_ray_length: 1000
+  use_half_hemisphere: false
 
 output:
   flux_output_file: "output/flux.dat"
   k_eff_output_file: "output/k_eff.dat"
 
 logging:
   level: "INFO"
-  log_file: "logs/solver.log"
 
 extra_section:
   key: "value"