Make field profile keys simpler.

Author: skang67

miniapps/tools/reconstruction.cpp

@@ -32,7 +32,7 @@ int main(int argc, char* argv[])
    int order_im = 3; // itermediate order, only used for L2 projection method
    int lref = order_im+1;
 
-   std::string field_profile = "1 + kx x + ky y";
+   std::string field_profile = "plane";
    real_t field_kx = 2.0;
    real_t field_ky = 4.0;
    bool use_ea = false;
@@ -245,21 +245,21 @@ std::unordered_map<std::string, profile_t> GetFieldProfiles()
 {
    std::unordered_map<std::string, profile_t> field_profiles;
    // plane profile
-   field_profiles["1 + kx x + ky y"] =
+   field_profiles["plane"] =
       [](const Vector &x, const Vector &k)
       {
          return 1.0 + x*k;
       };
    // sinusoidal profile
-   field_profiles["sin(2pi kx x) sin(2pi ky y)"] =
+   field_profiles["sinusoidal"] =
       [](const Vector &x, const Vector &k)
       {
          real_t result = 1.0;
          for(int i=0; i < x.Size(); i++) result *= std::sin(2.0*M_PI*k(i)*x(i));
          return result;
       };
    // exponential-sinusoidal profile
-   field_profiles["exp(r) cos(kx x) sin(ky y)"] =
+   field_profiles["exponential"] =
       [](const Vector &x, const Vector &k)
       {
          real_t result = 1.0;

miniapps/tools/run_convergence_study.py

@@ -18,7 +18,8 @@ def run_reconstruction(refinement_level, ho_value, method):
     try:
         # Run the executable
         result = subprocess.run(
-            ['./reconstruction', '-r', str(refinement_level), '-f', '2',
+            ['./reconstruction', '-r', str(refinement_level),
+             '-f', 'exponential',
              '-ho', str(ho_value), '-m', method],
             capture_output=True,
             text=True,
@@ -231,10 +232,10 @@ def main():
     cli_args = parser.parse_args()
 
     # Refinement levels to test
-    refinement_levels = [0, 1, 2, 3, 4]
+    refinement_levels = [0, 1, 2, 3]
 
     # ho values to test
-    ho_values = [1, 7]
+    ho_values = [1, 3]
 
     # Error data filename
     error_file = 'convergence_errors.json'