diff --git a/problems/quasiconvexity-circle-fung.parset b/problems/quasiconvexity-circle-fung.parset
new file mode 100644
index 0000000000000000000000000000000000000000..1e931a74f97fd1480dd9d70db190720b89300a43
--- /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 0000000000000000000000000000000000000000..c38393367235b63bd5f1ce8023c6dd3ccef1d3a8
--- /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 0000000000000000000000000000000000000000..396f44d46a7933004ffa26b2f8c025e274bc21a8
--- /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 0000000000000000000000000000000000000000..49ef62cb826e635273b6803db6806661a97cab05
--- /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 0000000000000000000000000000000000000000..ad598ee6d4c8a85c68a9221de4f2df40695c6ae9
--- /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