[SCM] Gerris Flow Solver branch, upstream, updated. b3aa46814a06c9cb2912790b23916ffb44f1f203
Stephane Popinet
popinet at users.sf.net
Fri May 15 02:53:36 UTC 2009
The following commit has been merged in the upstream branch:
commit 99119f787dc5da8639e6944bd4bb0a55efe4a20d
Author: Stephane Popinet <popinet at users.sf.net>
Date: Sun Feb 5 00:38:44 2006 +1100
Bug fix for VariableCurvature computation near boundaries
darcs-hash:20060204133844-d4795-f9b67b149bd7d23d5c24e61093e744fde6466415.gz
diff --git a/src/levelset.c b/src/levelset.c
index 8c5be47..e839429 100644
--- a/src/levelset.c
+++ b/src/levelset.c
@@ -299,7 +299,7 @@ static void interface_curvature (FttCell * cell, gpointer * data)
static void variable_curvature_event_half (GfsEvent * event, GfsSimulation * sim)
{
- GfsVariable * n[FTT_DIMENSION + 1];
+ GfsVariable * n[FTT_DIMENSION + 1], * kappa;
GfsDomain * domain = GFS_DOMAIN (sim);
gpointer data[2];
FttComponent c;
@@ -308,6 +308,7 @@ static void variable_curvature_event_half (GfsEvent * event, GfsSimulation * sim
n[c] = gfs_temporary_variable (domain);
gfs_variable_set_vector (n[c], c);
}
+ kappa = n[FTT_DIMENSION];
data[0] = n;
data[1] = event;
gfs_domain_cell_traverse (domain, FTT_PRE_ORDER, FTT_TRAVERSE_LEAFS, -1,
@@ -316,6 +317,7 @@ static void variable_curvature_event_half (GfsEvent * event, GfsSimulation * sim
gfs_domain_bc (domain, FTT_TRAVERSE_LEAFS, -1, n[c]);
gfs_domain_cell_traverse (domain, FTT_PRE_ORDER, FTT_TRAVERSE_LEAFS, -1,
(FttCellTraverseFunc) curvature, data);
+ gfs_domain_bc (domain, FTT_TRAVERSE_LEAFS, -1, kappa);
gfs_domain_cell_traverse (domain, FTT_PRE_ORDER, FTT_TRAVERSE_LEAFS, -1,
(FttCellTraverseFunc) interface_curvature, data);
gfs_domain_bc (domain, FTT_TRAVERSE_LEAFS, -1, GFS_VARIABLE1 (event));
--
Gerris Flow Solver
More information about the debian-science-commits
mailing list