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

[Stephane Popinet]

The following commit has been merged in the upstream branch:
commit e67284e479b47b1b1eed204664b246a38dae47e4
Author: Stephane Popinet <popinet at users.sf.net>
Date:   Thu Feb 2 21:39:56 2006 +1100

    Youngs-averaging of normal divergence was unstable in some cases
    
    darcs-hash:20060202103956-d4795-6e31ee0d81989bbcfff5be60caf3bead9d46dfb5.gz

diff --git a/src/variable.c b/src/variable.c
index c87e010..eb2f568 100644
--- a/src/variable.c
+++ b/src/variable.c
@@ -524,12 +524,9 @@ static void curvature (FttCell * cell, gpointer * data)
   GfsVariableCurvature * k = GFS_VARIABLE_CURVATURE (v);
   gdouble kappa = 0.;
   FttComponent c;
-  GtsVector n = { 0., 0., 0. };
-  
-  for (c = 0; c < FTT_DIMENSION; c++) {
-    gfs_youngs_normal (cell, nv[c], (FttVector *) n);
-    kappa += n[c];
-  }
+
+  for (c = 0; c < FTT_DIMENSION; c++)
+    kappa += gfs_center_gradient (cell, c, nv[c]->i);
   GFS_VARIABLE (cell, v->i) = (k->a*k->sigma*kappa/ftt_cell_size (cell) +
 			       (1. - k->a)*GFS_VARIABLE (cell, v->i));
 }

-- 
Gerris Flow Solver