Dalotia Integration

examples/configs/checkerboard_MN.cfg

@@ -8,51 +8,50 @@
 %
 OUTPUT_DIR = ../result
 % Output file
-OUTPUT_FILE = checkerboard_MN_Harm
+OUTPUT_FILE = checkerboard_MN_Normal_newton
 % Log directory
 LOG_DIR = ../result/logs
 % Mesh File
-MESH_FILE = meshes/checkerboard.su2
-%MESH_FILE = meshes/checkerboard_debug.su2
+MESH_FILE = ../meshes/checkerboard.su2
 %
 % ---- Problem specifications ----
 %
 PROBLEM = CHECKERBOARD
-SPATIAL_DIM = 3
+SPATIAL_DIM = 2
 SOURCE_MAGNITUDE = 1.0
 %
 % ---- Solver specifications ----
 %
 SPHERICAL_BASIS = SPHERICAL_HARMONICS
 %
 CFL_NUMBER = 0.45
-TIME_FINAL = 2
-SOLVER = MN_SOLVER
-RECONS_ORDER = 1
+TIME_FINAL = 3.2
+SOLVER = MN_SOLVER_NORMALIZED
 MAX_MOMENT_SOLVER = 3
+RECONS_ORDER = 2
+TIME_INTEGRATION_ORDER = 2
 %
 % ---- Entropy settings ----
 %
 ENTROPY_FUNCTIONAL = MAXWELL_BOLTZMANN
-ENTROPY_OPTIMIZER = NEWTON
-%NEURAL_MODEL = 7
+ENTROPY_OPTIMIZER = PART_REGULARIZED_NEWTON
 %
 % ----- Newton Solver Specifications ----
 %
-NEWTON_FAST_MODE = NO
-NEWTON_ITER = 100000
+REGULARIZER_GAMMA=0.001
+NEWTON_ITER = 1000
 NEWTON_EPSILON = 1e-6
-NEWTON_STEP_SIZE = 0.2
-NEWTON_LINE_SEARCH_ITER = 100000
+NEWTON_STEP_SIZE = 1
+NEWTON_LINE_SEARCH_ITER = 1000
 %
 % ---- Boundary Conditions ----
 %
 BC_DIRICHLET = ( void )
 %
 % ---- Quadrature ----
 %
-QUAD_TYPE = GAUSS_LEGENDRE_TENSORIZED
-QUAD_ORDER = 30
+QUAD_TYPE = GAUSS_LEGENDRE_TENSORIZED_2D
+QUAD_ORDER = 6
 %
 % ----- Output ---- 
 %
@@ -61,4 +60,4 @@ VOLUME_OUTPUT_FREQUENCY = 0
 SCREEN_OUTPUT = (ITER, MASS, RMS_FLUX, VTK_OUTPUT, CSV_OUTPUT)
 SCREEN_OUTPUT_FREQUENCY = 20
 HISTORY_OUTPUT = (ITER, MASS, RMS_FLUX, VTK_OUTPUT, CSV_OUTPUT)
-HISTORY_OUTPUT_FREQUENCY = 0
+HISTORY_OUTPUT_FREQUENCY = 1

examples/configs/checkerboard_MN_neural.cfg

@@ -6,43 +6,38 @@
 %
 % ---- File specifications ----
 %
-OUTPUT_DIR = ../../result
+OUTPUT_DIR = ../result
 % Output file
-OUTPUT_FILE = M1_neural_checker
+OUTPUT_FILE = checkerboard_MN_neural
 % Log directory
-LOG_DIR = ../../result/logs
+LOG_DIR = ../result/logs
 % Mesh File
 MESH_FILE = ../meshes/checkerboard.su2
-%MESH_FILE = ../meshes/checkerboard_debug.su2
 %
 % ---- Problem specifications ----
 %
 PROBLEM = CHECKERBOARD
 SPATIAL_DIM = 2
+SOURCE_MAGNITUDE = 1.0
 %
 % ---- Solver specifications ----
 %
-SPHERICAL_BASIS = SPHERICAL_MONOMIALS
+SPHERICAL_BASIS = SPHERICAL_HARMONICS
 %
-CFL_NUMBER = 0.7
-TIME_FINAL = 3
-SOLVER = MN_SOLVER
+CFL_NUMBER = 0.45
+TIME_FINAL = 3.2
+SOLVER = MN_SOLVER_NORMALIZED
+MAX_MOMENT_SOLVER = 3
 RECONS_ORDER = 2
-MAX_MOMENT_SOLVER = 1
+TIME_INTEGRATION_ORDER = 2
 %
 % ---- Entropy settings ----
 %
 ENTROPY_FUNCTIONAL = MAXWELL_BOLTZMANN
 ENTROPY_OPTIMIZER = ML
 NEURAL_MODEL_MK = 11
-%
-% ----- Newton Solver Specifications ----
-%
-NEWTON_FAST_MODE = NO
-NEWTON_ITER = 100000
-NEWTON_EPSILON = 0.0001
-NEWTON_STEP_SIZE = 0.7
-NEWTON_LINE_SEARCH_ITER = 10000
+NEURAL_MODEL_GAMMA = 3
+NEURAL_MODEL_ENFORCE_ROTATION_SYMMETRY = NO
 %
 % ---- Boundary Conditions ----
 %
@@ -51,13 +46,13 @@ BC_DIRICHLET = ( void )
 % ---- Quadrature ----
 %
 QUAD_TYPE = GAUSS_LEGENDRE_TENSORIZED_2D
-QUAD_ORDER = 12
+QUAD_ORDER = 6
 %
 % ----- Output ---- 
 %
 VOLUME_OUTPUT = (MINIMAL, MOMENTS, DUAL_MOMENTS)
-VOLUME_OUTPUT_FREQUENCY = 1
+VOLUME_OUTPUT_FREQUENCY = 0
 SCREEN_OUTPUT = (ITER, MASS, RMS_FLUX, VTK_OUTPUT, CSV_OUTPUT)
-SCREEN_OUTPUT_FREQUENCY = 1
+SCREEN_OUTPUT_FREQUENCY = 20
 HISTORY_OUTPUT = (ITER, MASS, RMS_FLUX, VTK_OUTPUT, CSV_OUTPUT)
 HISTORY_OUTPUT_FREQUENCY = 1

examples/configs/checkerboard_SN.cfg

@@ -10,12 +10,13 @@
 OUTPUT_DIR = ../result
 OUTPUT_FILE = checkerboard_SN
 LOG_DIR = ../result/logs
-MESH_FILE = meshes/checkerboard.su2
+% Mesh File
+MESH_FILE = ../meshes/checkerboard.su2
 %
 % --- Problem definition ---
 %
 PROBLEM = CHECKERBOARD
-TIME_FINAL = 2
+TIME_FINAL = 3.2
 SPATIAL_DIM = 3
 SOURCE_MAGNITUDE = 1.0
 %
@@ -26,7 +27,8 @@ SOLVER = SN_SOLVER
 % CFL number
 CFL_NUMBER = 0.45
 % Reconstruction order
-RECONS_ORDER = 1
+RECONS_ORDER = 2
+TIME_INTEGRATION_ORDER = 2
 %
 % ---- Boundary Conditions ----
 %
@@ -35,7 +37,7 @@ BC_DIRICHLET = ( void )
 % ----- Quadrature Specification ---
 %
 QUAD_TYPE = GAUSS_LEGENDRE_TENSORIZED
-QUAD_ORDER = 6
+QUAD_ORDER = 12
 %
 % ----- Output ---- 
 %

examples/configs/hohlraum_MN.cfg

@@ -8,51 +8,51 @@
 %
 OUTPUT_DIR = ../result
 % Output file
-OUTPUT_FILE = hohlraum
+OUTPUT_FILE = hohlraum_MN_newton
 % Log directory
 LOG_DIR =  ../result/logs
 % Mesh File
-MESH_FILE = meshes/hohlraum.su2
+MESH_FILE = ../meshes/hohlraum.su2
 %
 % ---- Problem specifications ----
 %
 PROBLEM = HOHLRAUM
-SPATIAL_DIM = 3
-SOURCE_MAGNITUDE = 1.0
+SPATIAL_DIM = 2
 %
 % ---- Solver specifications ----
 %
 SPHERICAL_BASIS = SPHERICAL_HARMONICS
 %
-CFL_NUMBER = 0.1
-TIME_FINAL = 2
-SOLVER = MN_SOLVER
+CFL_NUMBER = 0.45
+TIME_FINAL = 2.6
+SOLVER = MN_SOLVER_NORMALIZED
 RECONS_ORDER = 2
+TIME_INTEGRATION_ORDER = 2
 MAX_MOMENT_SOLVER = 3
 %
 % ---- Entropy settings ----
 %
 ENTROPY_FUNCTIONAL = MAXWELL_BOLTZMANN
 ENTROPY_OPTIMIZER = PART_REGULARIZED_NEWTON
-REGULARIZER_GAMMA =  0.00005
+REGULARIZER_GAMMA =  0.001
 %NEURAL_MODEL = 7
 %
 % ----- Newton Solver Specifications ----
 %
 NEWTON_FAST_MODE = NO
 NEWTON_ITER = 100000
-NEWTON_EPSILON = 1e-5
-NEWTON_STEP_SIZE = 0.2
-NEWTON_LINE_SEARCH_ITER = 100000
+NEWTON_EPSILON = 1e-4
+NEWTON_STEP_SIZE = 1
+NEWTON_LINE_SEARCH_ITER = 100
 %
 % ---- Boundary Conditions ----
 %
 BC_DIRICHLET = ( void )
 %
 % ---- Quadrature ----
 %
-QUAD_TYPE = GAUSS_LEGENDRE_TENSORIZED
-QUAD_ORDER = 20
+QUAD_TYPE = GAUSS_LEGENDRE_TENSORIZED_2D
+QUAD_ORDER = 6
 %
 % ----- Output ---- 
 %

examples/configs/hohlraum_MN_neural.cfg

@@ -8,34 +8,34 @@
 %
 OUTPUT_DIR = ../result
 % Output file
-OUTPUT_FILE = hohlraum_neural
+OUTPUT_FILE = hohlraum_MN_neural
 % Log directory
 LOG_DIR = ../result/logs
 % Mesh File
-MESH_FILE = meshes/hohlraum.su2
+MESH_FILE = ../meshes/hohlraum.su2
 %
 % ---- Problem specifications ----
 %
 PROBLEM = HOHLRAUM
 SPATIAL_DIM = 2
-SOURCE_MAGNITUDE = 1.0
 %
 % ---- Solver specifications ----
 %
 SPHERICAL_BASIS = SPHERICAL_MONOMIALS
 %
-CFL_NUMBER = 0.1
-TIME_FINAL = 2
+CFL_NUMBER = 0.45
+TIME_FINAL = 2.6
 SOLVER = MN_SOLVER_NORMALIZED
 RECONS_ORDER = 2
+TIME_INTEGRATION_ORDER = 2
 MAX_MOMENT_SOLVER = 3
 %
 % ---- Entropy settings ----
 %
 ENTROPY_FUNCTIONAL = MAXWELL_BOLTZMANN
 ENTROPY_OPTIMIZER = ML
 NEURAL_MODEL_MK = 11
-NEURAL_MODEL_GAMMA = 0
+NEURAL_MODEL_GAMMA = 3
 NEURAL_MODEL_ENFORCE_ROTATION_SYMMETRY = NO
 %
 % ---- Boundary Conditions ----
@@ -45,7 +45,7 @@ BC_DIRICHLET = ( void )
 % ---- Quadrature ----
 %
 QUAD_TYPE = GAUSS_LEGENDRE_TENSORIZED_2D
-QUAD_ORDER = 10
+QUAD_ORDER = 6
 %
 % ----- Output ---- 
 %

examples/example_runs.sh

@@ -0,0 +1,5 @@
+# Test case: Hohlraum
+# Newton-based Solver 
+../build/KiT-RT configs/hohlraum_MN.cfg
+# Neural Network-based Solver
+../build/KiT-RT configs/hohlraum_MN_neural.cfg