added fast functionality to set channel width

src/paper_thinstrip_full/chns.ini

@@ -12,26 +12,26 @@ CurvatureAlpha=0.001
 
 [domain]
 level = 0
-filename = grids/tper8.msh
+filename = grids/tper1.msh
 len = 1
 periodic = 01
 #YScalingFactor = 0.5
-YScaling = 0.125 #0.125
+YScaling = 1 #0.125
 
 [output]
-filename = output4
+filename = output1
 
 [Phasefield]
 eps=0.03   #0.03
 delta=0.03 #0.03
 DoubleWellDelta=0.05 #0.05
-Mobility=0.02
+Mobility=0.04
 Sigma1=1
 Sigma2=1
 Sigma3=1
 
 [Time]
-tMax = 300
+tMax = 10
 dt = 0.0015
 dtMax = 0.01
 dtMin = 0.00015

src/paper_thinstrip_full/thinstrip_equation.cc

@@ -106,7 +106,21 @@ int main(int argc, char** argv)
     ptreeparser.readINITree("chns.ini",ptree);
     ptreeparser.readOptions(argc,argv,ptree);
 
-    // read ini file
+    if(argc > 1)
+    {
+      std::cout << argv[1] << std::endl;
+      int factorsize = std::stoi( argv[1] );
+      ptree["output.filename"] = "output" + std::to_string(factorsize);
+      ptree["domain.filename"] = "grids/tper"  + std::to_string(factorsize) + ".msh";
+      ptree["domain.YScaling"] = std::to_string(1./float(factorsize));
+
+      ptree["Time.dt"] = std::to_string(1./float(factorsize) * ptree.get("Time.dt",0.1));
+      ptree["Time.dtMax"] = std::to_string(1./float(factorsize) * ptree.get("Time.dtMax",0.1));
+      ptree["Time.dtMin"] = std::to_string(1./float(factorsize) * ptree.get("Time.dtMin",0.1));
+    }
+
+    ptree["output.filename"] = "teststest";
+
 //    const int dim = ptree.get<int>("grid.dim");
     const int refinement = ptree.get<int>("domain.level");