Add further quasiconvexity problem files

problems/quasiconvexity-circle-fung.parset

+#############################################
+#  Grid parameters
+#############################################
+
+structuredGrid = false
+path = /home/sander/dune/dune-elasticity/grids/
+gridFile = circle2ndorder.msh
+
+#structuredGrid = true
+#lower = 0 0
+#upper = 1 1
+#elements = 2 2
+
+# Number of grid levels
+numLevels = 6
+
+#############################################
+#  Solver parameters
+#############################################
+
+# Number of homotopy steps for the Dirichlet boundary conditions
+numHomotopySteps = 1
+
+# Tolerance of the trust region solver
+tolerance = 1e-6
+
+# Max number of steps of the trust region solver
+maxTrustRegionSteps = 500
+
+# Initial trust-region radius
+initialTrustRegionRadius = 0.01
+
+# Number of multigrid iterations per trust-region step
+numIt = 200
+
+# Number of presmoothing steps
+nu1 = 3
+
+# Number of postsmoothing steps
+nu2 = 3
+
+# Number of coarse grid corrections
+mu = 1
+
+# Number of base solver iterations
+baseIt = 100
+
+# Tolerance of the multigrid solver
+mgTolerance = 1e-7
+
+# Tolerance of the base grid solver
+baseTolerance = 1e-8
+
+############################
+#   Material parameters
+############################
+
+energy = voss
+
+[materialParameters]
+
+c = 4.25
+A1 = 56.55
+A2 = 86.27
+A3 = -16.76
+A4 = -76.08
+A5 = -9.07
+A6 = -22.29
+
+[]
+
+#############################################
+#  Boundary values
+#############################################
+
+problem = neff-quasiconvexity
+
+###  Python predicate specifying all Dirichlet grid vertices
+# x is the vertex coordinate
+dirichletVerticesPredicate = "True"
+dirichletValues = affine
+
+# Initial iterate given as a Python function
+initialIteratePython = shear-rotate
+
+# Setting initialIterateFile overrides the x0 parameter
+#initialIterateFile = quasiconvexity-test-result-1.075-x0_1.1_8levels.data

problems/quasiconvexity-circle-teichmüller.parset

+#############################################
+#  Grid parameters
+#############################################
+
+structuredGrid = false
+path = /home/sander/dune/dune-elasticity/grids/
+gridFile = circle2ndorder.msh
+
+#structuredGrid = true
+#lower = 0 0
+#upper = 1 1
+#elements = 2 2
+
+# Number of grid levels
+numLevels = 8
+
+#############################################
+#  Solver parameters
+#############################################
+
+# Number of homotopy steps for the Dirichlet boundary conditions
+numHomotopySteps = 1
+
+# Tolerance of the trust region solver
+tolerance = 1e-6
+
+# Max number of steps of the trust region solver
+maxTrustRegionSteps = 3000
+
+# Initial trust-region radius
+initialTrustRegionRadius = 0.01
+
+# Number of multigrid iterations per trust-region step
+numIt = 200
+
+# Number of presmoothing steps
+nu1 = 3
+
+# Number of postsmoothing steps
+nu2 = 3
+
+# Number of coarse grid corrections
+mu = 1
+
+# Number of base solver iterations
+baseIt = 100
+
+# Tolerance of the multigrid solver
+mgTolerance = 1e-7
+
+# Tolerance of the base grid solver
+baseTolerance = 1e-8
+
+############################
+#   Material parameters
+############################
+
+energy = kpower
+
+[materialParameters]
+
+p = 100
+
+[]
+
+#############################################
+#  Boundary values
+#############################################
+
+problem = neff-quasiconvexity
+
+###  Python predicate specifying all Dirichlet grid vertices
+# x is the vertex coordinate
+dirichletVerticesPredicate = "True"
+dirichletValues = affine
+
+# Initial iterate given as a Python function
+initialIteratePython = teichmüller

problems/quasiconvexity-circle-voss.parset

+#############################################
+#  Grid parameters
+#############################################
+
+structuredGrid = false
+path = /home/sander/dune/dune-elasticity/grids/
+gridFile = circle2ndorder.msh
+
+#structuredGrid = true
+#lower = 0 0
+#upper = 1 1
+#elements = 2 2
+
+# Number of grid levels
+numLevels = 8
+
+#############################################
+#  Solver parameters
+#############################################
+
+# Number of homotopy steps for the Dirichlet boundary conditions
+numHomotopySteps = 1
+
+# Tolerance of the trust region solver
+tolerance = 1e-6
+
+# Max number of steps of the trust region solver
+maxTrustRegionSteps = 200
+
+# Initial trust-region radius
+initialTrustRegionRadius = 0.01
+
+# Number of multigrid iterations per trust-region step
+numIt = 200
+
+# Number of presmoothing steps
+nu1 = 3
+
+# Number of postsmoothing steps
+nu2 = 3
+
+# Number of coarse grid corrections
+mu = 1
+
+# Number of base solver iterations
+baseIt = 100
+
+# Tolerance of the multigrid solver
+mgTolerance = 1e-7
+
+# Tolerance of the base grid solver
+baseTolerance = 1e-8
+
+############################
+#   Material parameters
+############################
+
+energy = voss
+
+[materialParameters]
+
+c = -0.1
+
+[]
+
+#############################################
+#  Boundary values
+#############################################
+
+problem = neff-quasiconvexity
+
+###  Python predicate specifying all Dirichlet grid vertices
+# x is the vertex coordinate
+dirichletVerticesPredicate = "True"
+dirichletValues = affine
+
+# Initial deformation is (x[0] * sqrt(x0), x[1] / sqrt(x0))
+x0 = 10
+
+# Initial iterate given as a Python function
+#initialIteratePython = one-circle
+
+# Setting initialIterateFile overrides the x0 parameter
+#initialIterateFile = quasiconvexity-test-result-1.075-x0_1.1_8levels.data

problems/quasiconvexity-neural-networks.parset

+#############################################
+#  Grid parameters
+#############################################
+
+structuredGrid = true
+lower = 0 0
+upper = 1 1
+elements = 20 20
+
+periodic = true
+
+# Number of grid levels
+numLevels = 1
+
+#############################################
+#  Solver parameters
+#############################################
+
+# Number of homotopy steps for the Dirichlet boundary conditions
+numHomotopySteps = 1
+
+# Tolerance of the trust region solver
+tolerance = 1e-10
+
+# Max number of steps of the trust region solver
+maxTrustRegionSteps = 100
+
+# Initial trust-region radius
+initialTrustRegionRadius = 0.01
+
+# Number of multigrid iterations per trust-region step
+numIt = 200
+
+# Number of presmoothing steps
+nu1 = 3
+
+# Number of postsmoothing steps
+nu2 = 3
+
+# Number of coarse grid corrections
+mu = 1
+
+# Number of base solver iterations
+baseIt = 100
+
+# Tolerance of the multigrid solver
+mgTolerance = 1e-7
+
+# Tolerance of the base grid solver
+baseTolerance = 1e-8
+
+############################
+#   Material parameters
+############################
+
+energy = magic
+
+[materialParameters]
+
+c = 1.0
+
+[]
+
+#############################################
+#  Boundary values
+#############################################
+
+problem = neff-quasiconvexity
+
+dirichletVerticesPredicate = "True"
+
+# Initial deformation is (x[0] * sqrt(x0), x[1] / sqrt(x0))
+x0 = 2.25    # 1.5^2
+
+# Setting initialIterateFile overrides the x0 parameter
+initialIterateFile = /home/sander/patrizio/neural-networks/package/output.csv

problems/quasiconvexity-periodic-expacosh.parset

+#############################################
+#  Grid parameters
+#############################################
+
+#structuredGrid = false
+#path = /home/sander/dune/dune-elasticity/grids/
+#gridFile = circle2ndorder.msh
+
+structuredGrid = true
+lower = 0 0
+upper = 1 1
+elements = 40 40
+
+# Number of grid levels
+numLevels = 1
+
+periodic = false
+
+#############################################
+#  Solver parameters
+#############################################
+
+# Number of homotopy steps for the Dirichlet boundary conditions
+numHomotopySteps = 1
+
+# Tolerance of the trust region solver
+tolerance = 1e-8
+
+# Max number of steps of the trust region solver
+maxTrustRegionSteps = 200
+
+# Initial trust-region radius
+initialTrustRegionRadius = 0.01
+
+# Number of multigrid iterations per trust-region step
+numIt = 1
+
+# Number of presmoothing steps
+nu1 = 3
+
+# Number of postsmoothing steps
+nu2 = 3
+
+# Number of coarse grid corrections
+mu = 1
+
+# Number of base solver iterations
+baseIt = 100
+
+# Tolerance of the multigrid solver
+mgTolerance = 1e-7
+
+# Tolerance of the base grid solver
+baseTolerance = 1e-8
+
+############################
+#   Material parameters
+############################
+
+energy = expacosh
+
+[materialParameters]
+
+c1 = 0.11
+
+c2 = 0
+
+[]
+
+#############################################
+#  Boundary values
+#############################################
+
+problem = neff-quasiconvexity
+
+###  Python predicate specifying all Dirichlet grid vertices
+# x is the vertex coordinate
+dirichletVerticesPredicate = "True"
+dirichletValues = affine
+
+# Initial deformation is (x[0] * sqrt(x0), x[1] / sqrt(x0))
+#x0 = 4.38
+x0 = 12.0186
+
+# Additional initial perturbation wrt the affine deformation
+#initialIteratePython = sine-wave