Initial iterate is now the inverse stereographic projection

According to Christof Melcher, this produces a field of topological
charge q=1  (or -1, not sure).

[[Imported from SVN: r10006]]

src/harmonicmaps.parset

@@ -52,3 +52,6 @@ structuredGrid = true
 lower = -10 -10
 upper =  10  10
 elements = 100 100
+
+# 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)]"