amr block levels

[Ramses][AMR] AMR levels
Get AMR levels for all blocks in each patch.
This is done by using the ComputeAMRLvel node
implemented by @tdavidcl.

updates

Author: Akos299

src/shammodels/ramses/include/shammodels/ramses/modules/SolverStorage.hpp

@@ -77,6 +77,7 @@ namespace shammodels::basegodunov {
         std::shared_ptr<shamrock::solvergraph::FieldRefs<TgridVec>> refs_block_max;
 
         std::shared_ptr<shamrock::solvergraph::ScalarsEdge<TgridVec>> level0_size;
+        std::shared_ptr<shamrock::solvergraph::Field<TgridUint>> amr_block_levels;
 
         std::shared_ptr<shamrock::solvergraph::Indexes<u32>> block_counts;
         std::shared_ptr<shamrock::solvergraph::Indexes<u32>> block_counts_with_ghost;

src/shammodels/ramses/src/Solver.cpp

@@ -24,6 +24,7 @@
 #include "shammodels/ramses/SolverConfig.hpp"
 #include "shammodels/ramses/modules/AMRGridRefinementHandler.hpp"
 #include "shammodels/ramses/modules/BlockNeighToCellNeigh.hpp"
+#include "shammodels/ramses/modules/ComputeAMRLevel.hpp"
 #include "shammodels/ramses/modules/ComputeCFL.hpp"
 #include "shammodels/ramses/modules/ComputeCellAABB.hpp"
 #include "shammodels/ramses/modules/ComputeLevel0CellSize.hpp"
@@ -40,6 +41,7 @@
 #include "shammodels/ramses/modules/InterpolateToFace.hpp"
 #include "shammodels/ramses/modules/NodeComputeFlux.hpp"
 #include "shammodels/ramses/modules/SlopeLimitedGradient.hpp"
+#include "shammodels/ramses/modules/SolverStorage.hpp"
 #include "shammodels/ramses/modules/TimeIntegrator.hpp"
 #include "shammodels/ramses/modules/TransformGhostLayer.hpp"
 #include "shammodels/ramses/solvegraph/OrientedAMRGraphEdge.hpp"
@@ -412,6 +414,10 @@ void shammodels::basegodunov::Solver<Tvec, TgridVec>::init_solver_graph() {
     // get blocks at level0 sizes for all patches
     storage.level0_size = std::make_shared<shamrock::solvergraph::ScalarsEdge<TgridVec>>(
         "level0_amr", "level0_amr");
+    // get blocks amr leves
+    using TgridUint_         = typename SolverStorage<Tvec, TgridVec, u_morton>::TgridUint;
+    storage.amr_block_levels = std::make_shared<shamrock::solvergraph::Field<TgridUint_>>(
+        1, "amr_block_levels", "L_{\\rm amr}");
 
     storage.grad_rho = std::make_shared<shamrock::solvergraph::Field<Tvec>>(
         AMRBlock::block_size, "grad_rho", "\\nabla \\rho");
@@ -921,6 +927,18 @@ void shammodels::basegodunov::Solver<Tvec, TgridVec>::init_solver_graph() {
             std::make_shared<decltype(node_level0_sizes)>(std::move(node_level0_sizes)));
     }
 
+    { // compute block amr level in patch
+        modules::ComputeAMRLevel<TgridVec> node_amr_level{};
+        node_amr_level.set_edges(
+            storage.block_counts,
+            storage.level0_size,
+            storage.refs_block_min,
+            storage.refs_block_max,
+            storage.amr_block_levels);
+        solver_sequence.push_back(
+            std::make_shared<decltype(node_amr_level)>(std::move(node_amr_level)));
+    }
+
     if (solver_config.should_compute_rho_mean()) {
         modules::NodeComputeMass<Tvec, TgridVec> node{AMRBlock::block_size};
         node.set_edges(