[SCM] Gerris Flow Solver branch, upstream, updated. b3aa46814a06c9cb2912790b23916ffb44f1f203

[Stephane Popinet]

The following commit has been merged in the upstream branch:
commit 533eac01cc2ad65d209cfbc0fd85aeaefee7076f
Author: Stephane Popinet <popinet at users.sf.net>
Date:   Fri Feb 3 14:26:15 2006 +1100

    For curvature estimation, normals are computed using simple gradients
    
    ...rather than Youngs gradients. This stabilises the surface tension,
    improves the spurious currents and it is simpler.
    
    darcs-hash:20060203032615-d4795-a36b81dce283e403f2f05d2365236671baeccd47.gz

diff --git a/src/levelset.c b/src/levelset.c
index fc00c08..4fbab62 100644
--- a/src/levelset.c
+++ b/src/levelset.c
@@ -208,7 +208,8 @@ static void normal (FttCell * cell, gpointer * data)
   GtsVector n = { 0., 0., 0. };
   FttComponent c;
 
-  gfs_youngs_normal (cell, d, (FttVector *) n);
+  for (c = 0; c < FTT_DIMENSION; c++)
+    n[c] = gfs_center_gradient (cell, c, d->i);
   gts_vector_normalize (n);
   for (c = 0; c < FTT_DIMENSION; c++)
     GFS_VARIABLE (cell, nv[c]->i) = n[c];

-- 
Gerris Flow Solver