From f6248ffc6ea971af1fd4b016e5835ca449293bfd Mon Sep 17 00:00:00 2001
From: ScSteffen <steffen.schotthoefer@outlook.de>
Date: Thu, 18 Jan 2024 14:44:47 -0500
Subject: [PATCH] fixed bug in max absorption computation

---
 examples_hpc/lattice/mesh/lattice_rectangular.geo | 4 ++--
 examples_hpc/lattice/mesh/make_lattice_mesh.sh    | 2 +-
 src/problems/lattice.cpp                          | 1 -
 3 files changed, 3 insertions(+), 4 deletions(-)

diff --git a/examples_hpc/lattice/mesh/lattice_rectangular.geo b/examples_hpc/lattice/mesh/lattice_rectangular.geo
index 9721aecd..376415fa 100644
--- 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
index ea100c9c..8e559b4c 100644
--- 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
index d70eb6ec..644b7615 100644
--- 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();