[SCM] Gerris Flow Solver branch, upstream, updated. b3aa46814a06c9cb2912790b23916ffb44f1f203
Stephane Popinet
popinet at users.sf.net
Fri May 15 02:53:34 UTC 2009
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
More information about the debian-science-commits
mailing list