[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