feat: add dynamic configuration for mesh

demos/FiniteVolume/AMR_Burgers_Hat.cpp

@@ -197,7 +197,6 @@ int main(int argc, char* argv[])
     auto& app = samurai::initialize("Finite volume example for the Burgers equation in 2d using AMR", argc, argv);
 
     constexpr std::size_t dim = 1; // cppcheck-suppress unreadVariable
-    using Config              = samurai::amr::Config<dim>;
 
     // Simulation parameters
     double left_box  = -3;
@@ -208,10 +207,7 @@ int main(int argc, char* argv[])
     std::string restart_file;
 
     // AMR parameters
-    std::size_t start_level = 7;
-    std::size_t min_level   = 2;
-    std::size_t max_level   = 7;
-    bool correction         = false;
+    bool correction = false;
 
     // Output parameters
     fs::path path        = fs::current_path();
@@ -224,9 +220,6 @@ int main(int argc, char* argv[])
     app.add_option("--Ti", t, "Initial time")->capture_default_str()->group("Simulation parameters");
     app.add_option("--Tf", Tf, "Final time")->capture_default_str()->group("Simulation parameters");
     app.add_option("--restart-file", restart_file, "Restart file")->capture_default_str()->group("Simulation parameters");
-    app.add_option("--start-level", start_level, "Start level of AMR")->capture_default_str()->group("AMR parameters");
-    app.add_option("--min-level", min_level, "Minimum level of AMR")->capture_default_str()->group("AMR parameters");
-    app.add_option("--max-level", max_level, "Maximum level of AMR")->capture_default_str()->group("AMR parameters");
     app.add_option("--with-correction", correction, "Apply flux correction at the interface of two refinement levels")
         ->capture_default_str()
         ->group("AMR parameters");
@@ -236,12 +229,13 @@ int main(int argc, char* argv[])
     SAMURAI_PARSE(argc, argv);
 
     const samurai::Box<double, dim> box({left_box}, {right_box});
-    samurai::amr::Mesh<Config> mesh;
-    auto phi = samurai::make_scalar_field<double>("phi", mesh);
+    auto config = samurai::mesh_config<dim>().min_level(2).max_level(7).max_stencil_size(2).start_level(7).disable_minimal_ghost_width();
+    auto mesh   = samurai::amr::make_empty_Mesh(config);
+    auto phi    = samurai::make_scalar_field<double>("phi", mesh);
 
     if (restart_file.empty())
     {
-        mesh = {box, start_level, min_level, max_level};
+        mesh = samurai::amr::make_Mesh(config, box);
         init_solution(phi);
     }
     else
@@ -256,7 +250,7 @@ int main(int argc, char* argv[])
     auto tag = samurai::make_scalar_field<int>("tag", mesh);
     const xt::xtensor_fixed<int, xt::xshape<2, 1>> stencil_grad{{1}, {-1}};
 
-    const double dx      = mesh.cell_length(max_level);
+    const double dx      = mesh.min_cell_length();
     double dt            = 0.99 * dx;
     const double dt_save = Tf / static_cast<double>(nfiles);
 

demos/FiniteVolume/BZ/bz_2d.cpp

@@ -289,8 +289,8 @@ int main()
     const double regularity = 1.;    // Regularity guess for multiresolution
     const double epsilon_MR = 2.e-4; // Threshold used by multiresolution
 
-    using Config = samurai::MRConfig<dim, 2>;
-    samurai::MRMesh<Config> mesh(box, min_level, max_level); // Constructing mesh from the box
+    auto config = samurai::mesh_config<dim>().min_level(min_level).max_level(max_level).max_stencil_radius(2);
+    auto mesh   = samurai::make_MRMesh(config, box); // Constructing mesh from the box
 
     using mesh_id_t = typename decltype(mesh)::mesh_id_t;
 

demos/FiniteVolume/advection_1d.cpp

@@ -62,7 +62,6 @@ int main(int argc, char* argv[])
     auto& app = samurai::initialize("Finite volume example for the advection equation in 1d using multiresolution", argc, argv);
 
     constexpr std::size_t dim = 1; // cppcheck-suppress unreadVariable
-    using Config              = samurai::MRConfig<dim>;
 
     // Simulation parameters
     double left_box  = -2;
@@ -74,10 +73,6 @@ int main(int argc, char* argv[])
     double t         = 0.;
     std::string restart_file;
 
-    // Multiresolution parameters
-    std::size_t min_level = 6;
-    std::size_t max_level = 12;
-
     // Output parameters
     fs::path path        = fs::current_path();
     std::string filename = "FV_advection_1d";
@@ -91,29 +86,29 @@ int main(int argc, char* argv[])
     app.add_option("--Ti", t, "Initial time")->capture_default_str()->group("Simulation parameters");
     app.add_option("--Tf", Tf, "Final time")->capture_default_str()->group("Simulation parameters");
     app.add_option("--restart-file", restart_file, "Restart file")->capture_default_str()->group("Simulation parameters");
-    app.add_option("--min-level", min_level, "Minimum level of the multiresolution")->capture_default_str()->group("Multiresolution");
-    app.add_option("--max-level", max_level, "Maximum level of the multiresolution")->capture_default_str()->group("Multiresolution");
     app.add_option("--path", path, "Output path")->capture_default_str()->group("Output");
     app.add_option("--filename", filename, "File name prefix")->capture_default_str()->group("Output");
     app.add_option("--nfiles", nfiles, "Number of output files")->capture_default_str()->group("Output");
 
     SAMURAI_PARSE(argc, argv);
 
     const samurai::Box<double, dim> box({left_box}, {right_box});
-    samurai::MRMesh<Config> mesh;
-    auto u = samurai::make_scalar_field<double>("u", mesh);
+
+    auto config = samurai::mesh_config<dim>().min_level(6).max_level(12).periodic(is_periodic).max_stencil_size(2).disable_minimal_ghost_width();
+    auto mesh = samurai::make_empty_MRMesh(config);
+    auto u    = samurai::make_scalar_field<double>("u", mesh);
 
     if (restart_file.empty())
     {
-        mesh = {box, min_level, max_level, std::array<bool, dim>{is_periodic}};
+        mesh = samurai::make_MRMesh(config, box);
         init(u);
     }
     else
     {
         samurai::load(restart_file, mesh, u);
     }
 
-    double dt            = cfl * mesh.cell_length(max_level);
+    double dt            = cfl * mesh.min_cell_length();
     const double dt_save = Tf / static_cast<double>(nfiles);
 
     if (!is_periodic)

demos/FiniteVolume/advection_2d.cpp

@@ -63,7 +63,6 @@ int main(int argc, char* argv[])
     auto& app = samurai::initialize("Finite volume example for the advection equation in 2d using multiresolution", argc, argv);
 
     constexpr std::size_t dim = 2;
-    using Config              = samurai::MRConfig<dim>;
 
     // Simulation parameters
     xt::xtensor_fixed<double, xt::xshape<dim>> min_corner = {0., 0.};
@@ -76,10 +75,6 @@ int main(int argc, char* argv[])
     double t   = 0.;
     std::string restart_file;
 
-    // Multiresolution parameters
-    std::size_t min_level = 4;
-    std::size_t max_level = 10;
-
     // Output parameters
     fs::path path        = fs::current_path();
     std::string filename = "FV_advection_2d";
@@ -92,21 +87,20 @@ int main(int argc, char* argv[])
     app.add_option("--Ti", t, "Initial time")->capture_default_str()->group("Simulation parameters");
     app.add_option("--Tf", Tf, "Final time")->capture_default_str()->group("Simulation parameters");
     app.add_option("--restart-file", restart_file, "Restart file")->capture_default_str()->group("Simulation parameters");
-    app.add_option("--min-level", min_level, "Minimum level of the multiresolution")->capture_default_str()->group("Multiresolution");
-    app.add_option("--max-level", max_level, "Maximum level of the multiresolution")->capture_default_str()->group("Multiresolution");
     app.add_option("--path", path, "Output path")->capture_default_str()->group("Output");
     app.add_option("--filename", filename, "File name prefix")->capture_default_str()->group("Output");
     app.add_option("--nfiles", nfiles, "Number of output files")->capture_default_str()->group("Output");
 
     SAMURAI_PARSE(argc, argv);
 
     const samurai::Box<double, dim> box(min_corner, max_corner);
-    samurai::MRMesh<Config> mesh;
-    auto u = samurai::make_scalar_field<double>("u", mesh);
+    auto config = samurai::mesh_config<dim>().min_level(4).max_level(10).max_stencil_size(2).disable_minimal_ghost_width();
+    auto mesh   = samurai::make_empty_MRMesh(config);
+    auto u      = samurai::make_scalar_field<double>("u", mesh);
 
     if (restart_file.empty())
     {
-        mesh = {box, min_level, max_level};
+        mesh = samurai::make_MRMesh(config, box);
         init(u);
     }
     else
@@ -115,7 +109,7 @@ int main(int argc, char* argv[])
     }
     samurai::make_bc<samurai::Dirichlet<1>>(u, 0.);
 
-    double dt            = cfl * mesh.cell_length(max_level);
+    double dt            = cfl * mesh.min_cell_length();
     const double dt_save = Tf / static_cast<double>(nfiles);
 
     auto unp1 = samurai::make_scalar_field<double>("unp1", mesh);

demos/FiniteVolume/advection_2d_user_bc.cpp

@@ -95,7 +95,6 @@ int main(int argc, char* argv[])
     auto& app = samurai::initialize("Finite volume example for the advection equation in 2d using multiresolution", argc, argv);
 
     constexpr std::size_t dim = 2;
-    using Config              = samurai::MRConfig<dim>;
 
     // Simulation parameters
     xt::xtensor_fixed<double, xt::xshape<dim>> min_corner = {0., 0.};
@@ -108,10 +107,6 @@ int main(int argc, char* argv[])
     double t   = 0.;
     std::string restart_file;
 
-    // Multiresolution parameters
-    std::size_t min_level = 4;
-    std::size_t max_level = 10;
-
     // Output parameters
     fs::path path        = fs::current_path();
     std::string filename = "FV_advection_2d";
@@ -124,21 +119,20 @@ int main(int argc, char* argv[])
     app.add_option("--Ti", t, "Initial time")->capture_default_str()->group("Simulation parameters");
     app.add_option("--Tf", Tf, "Final time")->capture_default_str()->group("Simulation parameters");
     app.add_option("--restart-file", restart_file, "Restart file")->capture_default_str()->group("Simulation parameters");
-    app.add_option("--min-level", min_level, "Minimum level of the multiresolution")->capture_default_str()->group("Multiresolution");
-    app.add_option("--max-level", max_level, "Maximum level of the multiresolution")->capture_default_str()->group("Multiresolution");
     app.add_option("--path", path, "Output path")->capture_default_str()->group("Output");
     app.add_option("--filename", filename, "File name prefix")->capture_default_str()->group("Output");
     app.add_option("--nfiles", nfiles, "Number of output files")->capture_default_str()->group("Output");
 
     SAMURAI_PARSE(argc, argv);
 
     const samurai::Box<double, dim> box(min_corner, max_corner);
-    samurai::MRMesh<Config> mesh;
-    auto u = samurai::make_scalar_field<double>("u", mesh);
+    auto config = samurai::mesh_config<dim>().min_level(4).max_level(10).max_stencil_size(2).disable_minimal_ghost_width();
+    auto mesh   = samurai::make_empty_MRMesh(config);
+    auto u      = samurai::make_scalar_field<double>("u", mesh);
 
     if (restart_file.empty())
     {
-        mesh = {box, min_level, max_level};
+        mesh = samurai::make_MRMesh(config, box);
         init(u);
     }
     else
@@ -147,7 +141,7 @@ int main(int argc, char* argv[])
     }
     samurai::make_bc<Mybc>(u, 0.);
 
-    double dt            = cfl * mesh.cell_length(max_level);
+    double dt            = cfl * mesh.min_cell_length();
     const double dt_save = Tf / static_cast<double>(nfiles);
 
     auto unp1 = samurai::make_scalar_field<double>("unp1", mesh);

demos/FiniteVolume/burgers.cpp

@@ -52,7 +52,6 @@ int main_dim(int argc, char* argv[])
 {
     auto& app = samurai::initialize("Finite volume example for the Burgers equation", argc, argv);
 
-    using Config  = samurai::MRConfig<dim, 3>;
     using Box     = samurai::Box<double, dim>;
     using point_t = typename Box::point_t;
 
@@ -91,8 +90,6 @@ int main_dim(int argc, char* argv[])
     app.add_option("--restart-file", restart_file, "Restart file")->capture_default_str()->group("Simulation parameters");
     app.add_option("--dt", dt, "Time step")->capture_default_str()->group("Simulation parameters");
     app.add_option("--cfl", cfl, "The CFL")->capture_default_str()->group("Simulation parameters");
-    app.add_option("--min-level", min_level, "Minimum level of the multiresolution")->capture_default_str()->group("Multiresolution");
-    app.add_option("--max-level", max_level, "Maximum level of the multiresolution")->capture_default_str()->group("Multiresolution");
     app.add_option("--path", path, "Output path")->capture_default_str()->group("Output");
     app.add_option("--filename", filename, "File name prefix")->capture_default_str()->group("Output");
     app.add_option("--nfiles", nfiles, "Number of output files")->capture_default_str()->group("Output");
@@ -108,7 +105,9 @@ int main_dim(int argc, char* argv[])
     box_corner1.fill(left_box);
     box_corner2.fill(right_box);
     Box box(box_corner1, box_corner2);
-    samurai::MRMesh<Config> mesh;
+
+    auto config = samurai::mesh_config<dim>().min_level(min_level).max_level(max_level).max_stencil_size(6);
+    auto mesh   = samurai::make_empty_MRMesh(config);
 
     auto u    = samurai::make_vector_field<n_comp>("u", mesh);
     auto u1   = samurai::make_vector_field<n_comp>("u1", mesh);
@@ -117,7 +116,7 @@ int main_dim(int argc, char* argv[])
 
     if (restart_file.empty())
     {
-        mesh = {box, min_level, max_level};
+        mesh = samurai::make_MRMesh(config, box);
         u.resize();
 
         // Initial solution
@@ -221,7 +220,7 @@ int main_dim(int argc, char* argv[])
     //   Time iteration   //
     //--------------------//
 
-    double dx = mesh.cell_length(max_level);
+    double dx = mesh.min_cell_length();
     dt        = cfl * dx / pow(2, dim);
 
     auto MRadaptation = samurai::make_MRAdapt(u);

demos/FiniteVolume/burgers_mra.cpp

@@ -226,7 +226,6 @@ void run_simulation(Field& u,
 int main(int argc, char* argv[])
 {
     constexpr std::size_t dim = 1;
-    using Config              = samurai::MRConfig<dim, 2, 2>;
     using Box                 = samurai::Box<double, dim>;
 
     auto& app = samurai::initialize("Finite volume example for the Burgers equation in 1d", argc, argv);
@@ -245,19 +244,13 @@ int main(int argc, char* argv[])
     // Time integration
     double cfl = 0.95;
 
-    // Multiresolution parameters
-    std::size_t min_level = 1;
-    std::size_t max_level = 7;
-
     // Output parameters
     fs::path path        = fs::current_path();
     std::string filename = "burgers_mra";
     std::size_t nfiles   = 50;
 
     app.add_option("--init-sol", init_sol, "Initial solution: hat/gaussian")->capture_default_str()->group("Simulation parameters");
     app.add_option("--cfl", cfl, "The CFL")->capture_default_str()->group("Simulation parameters");
-    app.add_option("--min-level", min_level, "Minimum level of the multiresolution")->capture_default_str()->group("Multiresolution");
-    app.add_option("--max-level", max_level, "Maximum level of the multiresolution")->capture_default_str()->group("Multiresolution");
     app.add_option("--path", path, "Output path")->capture_default_str()->group("Ouput");
     app.add_option("--filename", filename, "File name prefix")->capture_default_str()->group("Ouput");
     app.add_option("--nfiles", nfiles, "Number of output files")->capture_default_str()->group("Ouput");
@@ -270,8 +263,10 @@ int main(int argc, char* argv[])
     //--------------------//
 
     Box box({left_box}, {right_box});
-    samurai::MRMesh<Config> mesh{box, min_level, max_level};
-    samurai::MRMesh<Config> max_level_mesh{box, max_level, max_level};
+    auto config           = samurai::mesh_config<dim>().min_level(1).max_level(7).max_stencil_radius(2).graduation_width(2);
+    auto mesh             = samurai::make_MRMesh(config, box);
+    auto max_level_config = samurai::mesh_config<dim>().min_level(mesh.max_level()).max_level(mesh.max_level()).max_stencil_radius(2);
+    auto max_level_mesh   = samurai::make_MRMesh(max_level_config, box);
 
     auto u    = samurai::make_scalar_field<double>("u", mesh);
     auto unp1 = samurai::make_scalar_field<double>("unp1", mesh);

demos/FiniteVolume/burgers_os.cpp

@@ -90,9 +90,9 @@ int main(int argc, char* argv[])
     auto& app = samurai::initialize("Finite volume example for the linear convection equation", argc, argv);
 
     static constexpr std::size_t dim = 1;
-    using Config                     = samurai::MRConfig<dim, 2, 2>;
-    using Box                        = samurai::Box<double, dim>;
-    using point_t                    = typename Box::point_t;
+
+    using Box     = samurai::Box<double, dim>;
+    using point_t = typename Box::point_t;
 
     std::cout << "------------------------- Burgers 1D -------------------------" << std::endl;
 
@@ -125,17 +125,13 @@ int main(int argc, char* argv[])
     app.add_option("--Tf", Tf, "Final time")->capture_default_str()->group("Simulation parameters");
     app.add_option("--dt", dt, "Time step")->capture_default_str()->group("Simulation parameters");
     app.add_option("--cfl", cfl, "The CFL")->capture_default_str()->group("Simulation parameters");
-    app.add_option("--min-level", min_level, "Minimum level of the multiresolution")->capture_default_str()->group("Multiresolution");
-    app.add_option("--max-level", max_level, "Maximum level of the multiresolution")->capture_default_str()->group("Multiresolution");
     app.add_option("--path", path, "Output path")->capture_default_str()->group("Output");
     app.add_option("--filename", filename, "File name prefix")->capture_default_str()->group("Output");
     app.add_option("--nfiles", nfiles, "Number of output files")->capture_default_str()->group("Output");
     app.allow_extras();
 
     SAMURAI_PARSE(argc, argv);
 
-    std::cout << "  max_level = " << max_level << "   min_level = " << min_level << std::endl;
-
     //--------------------//
     // Problem definition //
     //--------------------//
@@ -148,13 +144,18 @@ int main(int argc, char* argv[])
     Box box_left(box_corner1, 0.5 * (box_corner1 + box_corner2));
     Box box_right(0.5 * (box_corner1 + box_corner2), box_corner2);
 
-    samurai::LevelCellArray<1> lca_left(max_level, box_left, {-1}, 0.05, 2);
-    samurai::LevelCellArray<1> lca_right(max_level, box_right, {-1}, 0.05, 2);
+    samurai::LevelCellArray<1> lca_left(max_level, box_left, {-1}, -1., 2);
+    samurai::LevelCellArray<1> lca_right(max_level, box_right, {-1}, -1., 2);
 
-    std::array<bool, dim> periodic;
-    periodic.fill(true);
-    samurai::MRMesh<Config> mesh{box, min_level, max_level, periodic, 0.05, 2};
+    auto config = samurai::mesh_config<dim>()
+                      .min_level(min_level)
+                      .max_level(max_level)
+                      .periodic(true)
+                      .scaling_factor(2)
+                      .max_stencil_radius(1)
+                      .graduation_width(2);
 
+    auto mesh = samurai::make_MRMesh(config, box);
     auto u    = samurai::make_scalar_field<double>("u", mesh);
     auto unp1 = samurai::make_scalar_field<double>("unp1", mesh);
 
@@ -187,7 +188,7 @@ int main(int argc, char* argv[])
     //   Time iteration   //
     //--------------------//
 
-    double dx = mesh.cell_length(max_level);
+    double dx = mesh.min_cell_length();
     dt        = cfl * dx / velocity(0);
 
     while (t != Tf)