From 1bb3394698717658e57dfb5a25a300f6f9a2aa0a Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 14 Jan 2015 20:25:45 +0000
Subject: [PATCH] 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 | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/src/harmonicmaps.parset b/src/harmonicmaps.parset
index 52ffe12d..0735c6b4 100644
--- a/src/harmonicmaps.parset
+++ b/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)]"
-- 
GitLab