diff --git a/inputs/cellsolver.parset b/inputs/cellsolver.parset
index eff49c5decd3cbcd2a02260b30ca4440fdd3228d..6285f6a43e6aff9477b18e50c58e458a03832ac4 100644
--- a/inputs/cellsolver.parset
+++ b/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