From 1e4ad767dcb2cf413ae22b4df407a11841fd1852 Mon Sep 17 00:00:00 2001
From: Oliver Sander <oliver.sander@tu-dresden.de>
Date: Fri, 11 Sep 2020 10:15:42 +0200
Subject: [PATCH] Add further quasiconvexity problem files

---
 problems/quasiconvexity-circle-fung.parset    | 87 +++++++++++++++++++
 ...onvexity-circle-teichm\303\274ller.parset" | 78 +++++++++++++++++
 problems/quasiconvexity-circle-voss.parset    | 84 ++++++++++++++++++
 .../quasiconvexity-neural-networks.parset     | 76 ++++++++++++++++
 .../quasiconvexity-periodic-expacosh.parset   | 86 ++++++++++++++++++
 5 files changed, 411 insertions(+)
 create mode 100644 problems/quasiconvexity-circle-fung.parset
 create mode 100644 "problems/quasiconvexity-circle-teichm\303\274ller.parset"
 create mode 100644 problems/quasiconvexity-circle-voss.parset
 create mode 100644 problems/quasiconvexity-neural-networks.parset
 create mode 100644 problems/quasiconvexity-periodic-expacosh.parset

diff --git a/problems/quasiconvexity-circle-fung.parset b/problems/quasiconvexity-circle-fung.parset
new file mode 100644
index 0000000..1e931a7
--- /dev/null
+++ b/problems/quasiconvexity-circle-fung.parset
@@ -0,0 +1,87 @@
+#############################################
+#  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
diff --git "a/problems/quasiconvexity-circle-teichm\303\274ller.parset" "b/problems/quasiconvexity-circle-teichm\303\274ller.parset"
new file mode 100644
index 0000000..c383933
--- /dev/null
+++ "b/problems/quasiconvexity-circle-teichm\303\274ller.parset"
@@ -0,0 +1,78 @@
+#############################################
+#  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
diff --git a/problems/quasiconvexity-circle-voss.parset b/problems/quasiconvexity-circle-voss.parset
new file mode 100644
index 0000000..396f44d
--- /dev/null
+++ b/problems/quasiconvexity-circle-voss.parset
@@ -0,0 +1,84 @@
+#############################################
+#  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
diff --git a/problems/quasiconvexity-neural-networks.parset b/problems/quasiconvexity-neural-networks.parset
new file mode 100644
index 0000000..49ef62c
--- /dev/null
+++ b/problems/quasiconvexity-neural-networks.parset
@@ -0,0 +1,76 @@
+#############################################
+#  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
diff --git a/problems/quasiconvexity-periodic-expacosh.parset b/problems/quasiconvexity-periodic-expacosh.parset
new file mode 100644
index 0000000..ad598ee
--- /dev/null
+++ b/problems/quasiconvexity-periodic-expacosh.parset
@@ -0,0 +1,86 @@
+#############################################
+#  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
-- 
GitLab