add graduation_width in mesh_config

Author: kivvix

demos/FiniteVolume/heat.cpp

@@ -127,7 +127,6 @@ int main(int argc, char* argv[])
     {
         auto config = samurai::mesh_config<dim>().min_level(min_level).max_level(max_level);
         mesh        = {config, box};
-        // mesh = {box, min_level, max_level};
         u.resize();
         // Initial solution
         if (init_sol == "dirac")

demos/FiniteVolume/linear_convection.cpp

@@ -76,8 +76,6 @@ int main(int argc, char* argv[])
     app.add_flag("--implicit", implicit_scheme, "Implicit scheme instead of explicit")->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");
@@ -100,7 +98,8 @@ int main(int argc, char* argv[])
 
     if (restart_file.empty())
     {
-        mesh = {box, min_level, max_level, periodic};
+        auto config = samurai::mesh_config<dim>().min_level(min_level).max_level(max_level).periodic(periodic);
+        mesh        = {config, box};
         // Initial solution
         u = samurai::make_scalar_field<double>("u",
                                                mesh,
@@ -145,7 +144,7 @@ int main(int argc, char* argv[])
 
     if (dt == 0)
     {
-        double dx             = mesh.cell_length(max_level);
+        double dx             = mesh.cell_length(mesh.max_level());
         auto a                = xt::abs(velocity);
         double sum_velocities = xt::sum(xt::abs(velocity))();
         dt                    = cfl * dx / sum_velocities;

demos/FiniteVolume/linear_convection_obstacle.cpp

@@ -62,8 +62,6 @@ 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");
@@ -77,7 +75,9 @@ int main(int argc, char* argv[])
     samurai::DomainBuilder<dim> domain({-1., -1.}, {1., 1.});
     domain.remove({0.0, 0.0}, {0.4, 0.4});
 
-    Mesh mesh(domain, min_level, max_level);
+    auto config = samurai::mesh_config<dim>().min_level(min_level).max_level(max_level);
+    Mesh mesh   = {config, domain};
+    // Mesh mesh(domain, min_level, max_level);
 
     // Initial solution
     auto u = samurai::make_scalar_field<double>("u",
@@ -117,7 +117,7 @@ int main(int argc, char* argv[])
 
     if (dt == 0)
     {
-        double dx             = mesh.cell_length(max_level);
+        double dx             = mesh.cell_length(mesh.max_level());
         auto a                = xt::abs(constant_velocity);
         double sum_velocities = xt::sum(xt::abs(constant_velocity))();
         dt                    = cfl * dx / sum_velocities;

include/samurai/arguments.hpp

@@ -9,8 +9,9 @@ namespace samurai
     namespace args
     {
         // Mesh arguments
-        static std::size_t min_level = std::numeric_limits<std::size_t>::max();
-        static std::size_t max_level = std::numeric_limits<std::size_t>::max();
+        static std::size_t min_level        = std::numeric_limits<std::size_t>::max();
+        static std::size_t max_level        = std::numeric_limits<std::size_t>::max();
+        static std::size_t graduation_width = std::numeric_limits<std::size_t>::max();
 
         static bool timers = false;
 #ifdef SAMURAI_WITH_MPI
@@ -30,6 +31,7 @@ namespace samurai
     {
         app.add_option("--min-level", args::min_level, "The minimum level of the mesh")->group("SAMURAI");
         app.add_option("--max-level", args::max_level, "The maximum level of the mesh")->group("SAMURAI");
+        app.add_option("--graduation-width", args::graduation_width, "The graduation width of the mesh")->group("SAMURAI");
 
 #ifdef SAMURAI_WITH_MPI
         app.add_flag("--dont-redirect-output", args::dont_redirect_output, "Redirect the output for all ranks different of 0")