Update cellsolver.parset

inputs/cellsolver.parset

@@ -5,65 +5,43 @@
 # --------------------------------------------------------
 
 
-
-
-#path for logfile
-####outputPath = "../../outputs/output.txt"
-
-
-### Remove/Comment this when running via Python-Script:
-#outputPath = "../../outputs"
-
-
-
+### choose different path for logfile
+#outputPath = "../../outputs/output.txt"
 #outputPath = "/home/klaus/Desktop/DUNE/dune-microstructure/outputs/output.txt"
 #outputPath = "/home/klaus/Desktop/DUNE/dune-microstructure/outputs"
-
+### Remove/Comment this when running via Python-Script:
+outputPath = "../../outputs"
 
 #############################################
 #  Cell Domain
 #############################################
 # Domain 1: (-1/2, 1/2)^3  , Domain 2 : [0,1)^2 x (-1/2, 1/2)
-
 cellDomain=1
 
 #############################################
 #  Grid parameters
 #############################################
 
-
-#######################################################################
+#----------------------------------------------------
 ## numLevels : Number of Levels on which solution is computed. starting with a 2x2x2 cube mesh.
 ## {start,finish} computes on all grid from 2^(start) to 2^finish refinement
-########################################################################
+#----------------------------------------------------
 
 #numLevels =  1 1   # computes all levels from first to second entry
-#numLevels =  2 2   # computes all levels from first to second entry
+numLevels =  2 2   # computes all levels from first to second entry
 #numLevels =  1 3   # computes all levels from first to second entry
 #numLevels = 3 3     # computes all levels from first to second entry
 #numLevels = 4 4    # computes all levels from first to second entry
-numLevels = 5 5    # computes all levels from first to second entry
+#numLevels = 5 5    # computes all levels from first to second entry
 #numLevels = 6 6    # computes all levels from first to second entry
 #numLevels = 1 6
 
+########################################################################################
 
-#Elements_Cell = 20 20 20	               # number elements in each direction (y1 y2 x3)
-#nElements_Cell = 30 30 30
-#nElements_Cell = 30 30 30
-#nElements_Cell = 50 50 50
-#nElements_Cell = 100 100 2
-#nElements_Cell = 100 100 100  // does not work
-#nElements_Cell = 10 10 10
-#nElements_Cell = 2 2 2
-#nElements_Cell = 4 4 4
-#nElements_Cell = 8 8 8
-#nElements_Cell = 16 16 16
-#nElements_Cell = 32 32 32
-#nElements_Cell = 64 64 64
-
+# --- Choose Gamma value:
 
+gamma=1.0  #(default)
 #gamma=50.0
-gamma=1.0
 #gamma=0.01
 #gamma=3.0
 #gamma=0.5
@@ -71,93 +49,82 @@ gamma=1.0
 #gamma=2.5
 
 #############################################
-#  Material parameters
+#  Material / Prestrain parameters and ratios
 #############################################
-beta = 2.0    # ratio between material parameters mu1 & mu2 .... beta = 1.0 corresponds to homogeneous case
 
-mu1=1.0
-lambda1=0.0
-#lambda1=1.0
+#-- stiffness-ratio (ratio between material parameters mu1 & mu2 .... beta = 1.0 corresponds to homogeneous case)
+beta = 2.0
 
 
+$--- strength of prestrain
 rho1 = 1.0
-#alpha = 5.0    # ratio between prestrain parameters rho1 & rho2
-alpha = 2.0     # ratio between prestrain parameters rho1 & rho2
 
+#--- prestrain-contrast (ratio between prestrain parameters rho1 & rho2)
+#alpha = 5.0
+alpha = 2.0
 
-theta = 0.125
-#theta = 0.5
 
-#theta = 0.25
-#theta = 0.3    # volume fraction   #default = 1.0/4.0
-#theta = 0.25   # volume fraction
-#theta = 0.75   # volume fraction
+#--- Lame-Parameters
+mu1=1.0
 
+lambda1=0.0
+#lambda1=1.0
 
 
+# ---volume fraction  (default value = 1.0/4.0)
+#theta = 1.0/4.0
+#theta = 0.3
+#theta = 0.75
+theta = 0.125
+#theta = 0.5
 
 
 
+#--- choose composite-Implementation:
 material_prestrain_imp= "parametrized_Laminate"
 #material_prestrain_imp= "analytical_Example"
 #material_prestrain_imp="isotropic_bilayer"
 #material_prestrain_imp= "circle_fiber"    #TEST
 
 
-# -- write Material and/or prestrain functions:
-
-write_materialFunctions = true
-write_prestrainFunctions = true  # VTK norm of B ,
-
-write_VTK = true
-
-
+# --- (Optional output) write Material / prestrain / Corrector functions to .vtk-Files:
+#write_materialFunctions = true
+#write_prestrainFunctions = true  # VTK norm of B ,
+#write_VTK = true
 
-# Prestrain Types:
 
-#1 Isotropic Pressure
-#    Func2Tensor B1_ = [this] (const Domain& x) {              // ISOTROPIC PRESSURE
-#                if (abs(x[0]) > (theta/2)  && x[2] > 0)
-#                    return MatrixRT{{p1, 0.0 , 0.0}, {0.0, p1, 0.0}, {0.0, 0.0, p1}};
-#                if (abs(x[0]) < (theta/2)  && x[2] < 0)
-#                    return MatrixRT{{p2, 0.0 , 0.0}, {0.0, p2, 0.0}, {0.0, 0.0, p2}};
-#                else
-#                    return MatrixRT{{0.0, 0.0 , 0.0}, {0.0, 0.0, 0.0}, {0.0, 0.0, 0.0}};
-#            };
+# --- (Optional output) L2Error, compute integral mean:
+#write_L2Error = true
+#write_IntegralMean = true
 
 
 #############################################
 #  Assembly options
 #############################################
-
 set_IntegralZero = true
 #set_IntegralZero = false
 
 #arbitraryLocalIndex = 7
 #arbitraryElementNumber = 3
-
 #arbitraryLocalIndex = 0
 #arbitraryElementNumber = 0
+#############################################
 
 
 #############################################
-#  Solver Type
-
+#  Solver Type: #1: CG - SOLVER (default), #2: GMRES - SOLVER, #3: QR - SOLVER
+#############################################
 Solvertype = 1
+Solver_verbosity = 0  #(default = 2)  degree of information for solver output
+
 
-Solver_verbosity = 0
 
 
+
+# --- Write corrector-coefficients to log-File (default=false):
 #write_corrector_phi1 = false
 #write_corrector_phi2 = false
 #write_corrector_phi3 = false
 #write_corrector_phi1 = true
 #write_corrector_phi2 = true
 #write_corrector_phi3 = true
-
-#write_L2Error = true
-#write_IntegralMean = true
-
-
-#############################################
-#  Define Analytic Solutions