fixed bug in max absorption computation

KiT-RT · Jan 18, 2024 · f6248ff · f6248ff
1 parent fee759f
commit f6248ff
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 4 deletions.
diff --git a/examples_hpc/lattice/mesh/lattice_rectangular.geo b/examples_hpc/lattice/mesh/lattice_rectangular.geo
@@ -1,6 +1,6 @@
 cl_fine = 0.1;
-n_recombine = 20;
-n_prog = 1.1;
+n_recombine = 40;
+n_prog = 1.05;
 Point(1) = {-3.5, -3.5, 0, cl_fine};
 Point(2) = {3.5, -3.5, 0, cl_fine};
 Point(3) = {-3.5, 3.5, 0, cl_fine};

diff --git a/examples_hpc/lattice/mesh/make_lattice_mesh.sh b/examples_hpc/lattice/mesh/make_lattice_mesh.sh
@@ -1 +1 @@
-gmsh lattice_rectangular.geo -2 -format su2 -save_all -o lattice_n20.su2
+gmsh lattice_rectangular.geo -2 -format su2 -save_all -o lattice_n40.su2
diff --git a/src/problems/lattice.cpp b/src/problems/lattice.cpp
@@ -146,7 +146,6 @@ void Lattice_SN::ComputeCurrentAbsorptionLattice( const Vector& scalarFlux ) {
 }
 // TODO all absorption qois can be refactored in one function
 void Lattice_SN::ComputeMaxAbsorptionLattice( const Vector& scalarFlux ) {
-    _curMaxAbsorptionLattice  = 0.0;
     unsigned nCells           = _mesh->GetNumCells();
     auto cellMids             = _mesh->GetCellMidPoints();
     std::vector<double> areas = _mesh->GetCellAreas();
Original file line number	Diff line number	Diff line change
		@@ -1 +1 @@
		gmsh lattice_rectangular.geo -2 -format su2 -save_all -o lattice_n20.su2
		gmsh lattice_rectangular.geo -2 -format su2 -save_all -o lattice_n40.su2