From 93dc6f2406e2ca5f77d9640602a46df02a47dda1 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 27 Jan 2015 12:55:31 +0000
Subject: [PATCH] New parameter file for measuring discretization errors

The plan is to minimize the harmonic energy starting from the inverse
stereographic projection.  That projection already is the minimizer
of the energy, so it should be possible to measure convergence orders
against a known solution.

[[Imported from SVN: r10019]]
---
 ...monicmaps-stereographic-measurement.parset | 58 +++++++++++++++++++
 1 file changed, 58 insertions(+)
 create mode 100644 src/harmonicmaps-stereographic-measurement.parset

diff --git a/src/harmonicmaps-stereographic-measurement.parset b/src/harmonicmaps-stereographic-measurement.parset
new file mode 100644
index 00000000..c5f9e277
--- /dev/null
+++ b/src/harmonicmaps-stereographic-measurement.parset
@@ -0,0 +1,58 @@
+#############################################
+#  Grid parameters
+#############################################
+
+structuredGrid = true
+lower = -5 -5
+upper =  5  5
+elements = 30 30
+
+# Number of grid levels
+numLevels = 1
+
+#############################################
+#  Solver parameters
+#############################################
+
+# Tolerance of the trust region solver
+tolerance = 1e-12
+
+# Max number of steps of the trust region solver
+maxTrustRegionSteps = 10
+
+# Initial trust-region radius
+initialTrustRegionRadius = 1
+
+# 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-10
+
+# Tolerance of the base grid solver
+baseTolerance = 1e-8
+
+# Measure convergence
+instrumented = false
+
+############################
+#   Problem specifications
+############################
+
+# Type of energy we are minimizing
+energy = harmonic
+
+# Inverse stereographic projection
+initialIterate = "[2*x[0] / (x[0]*x[0]+x[1]*x[1]+1), 2*x[1] / (x[0]*x[0]+x[1]*x[1]+1), (x[0]*x[0]+x[1]*x[1]-1)/ (x[0]*x[0]+x[1]*x[1]+1)]"
-- 
GitLab