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

[Stephane Popinet]

The following commit has been merged in the upstream branch:
commit cd43cdfaff13ecb23f28a04f5073795c6504e11e
Author: Stephane Popinet <popinet at users.sf.net>
Date:   Mon Jan 30 20:46:30 2006 +1100

    VariableCurvature uses Youngs divergence
    
    Rather than gfs_center_gradient(). This improves the spurious currents slightly.
    
    darcs-hash:20060130094630-d4795-887c800389dc62fc4484b59d3802925497592eb0.gz

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

-- 
Gerris Flow Solver