Extend the parameters to allow measuring the discretization error

[[Imported from SVN: r10026]]

src/harmonicmaps-stereographic-measurement.parset

@@ -5,10 +5,10 @@
 structuredGrid = true
 lower = -5 -5
 upper =  5  5
-elements = 30 30
+elements = 10 10
 
 # Number of grid levels
-numLevels = 1
+numLevels = 2
 
 #############################################
 #  Solver parameters
@@ -18,7 +18,7 @@ numLevels = 1
 tolerance = 1e-12
 
 # Max number of steps of the trust region solver
-maxTrustRegionSteps = 10
+maxTrustRegionSteps = 0
 
 # Initial trust-region radius
 initialTrustRegionRadius = 1
@@ -53,3 +53,14 @@ 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)]"
+
+##########################################
+#   Discretization error measurements
+##########################################
+
+# none / analytical / gridfunction
+discretizationErrorMode = analytical
+
+# Again, the inverse stereographic projection
+referenceSolution = "[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)]"
+