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

[Stephane Popinet]

The following commit has been merged in the upstream branch:
commit 0bff000ead144532046b14c6faa95db381ef46ef
Author: Stephane Popinet <popinet at users.sf.net>
Date:   Thu Oct 19 14:16:31 2006 +1000

    Added warnings when HF-curvature fails
    
    darcs-hash:20061019041631-d4795-1f871f17e7bf85cfddf6ebc8e6ab31a7f757248b.gz

diff --git a/src/vof.c b/src/vof.c
index 3078a6e..a5e840f 100644
--- a/src/vof.c
+++ b/src/vof.c
@@ -966,15 +966,19 @@ gdouble gfs_height_curvature (FttCell * cell, GfsVariable * v)
   (&q.x)[cp] += size;
   (&q.x)[c] -= slope*size;
   //  fprintf (stderr, "\n  (%g %g) ", q.x, q.y);
-  if (!local_height (&q, &p, level, v, c, &h[0]))
+  if (!local_height (&q, &p, level, v, c, &h[0])) {
+    g_warning ("Failed to compute local height at (%g,%g)", q.x, q.y);
     return G_MAXDOUBLE;
+  }
 
   q = p;
   (&q.x)[cp] -= size;
   (&q.x)[c] += slope*size;
   //  fprintf (stderr, "\n  (%g %g) ", q.x, q.y);
-  if (!local_height (&q, &p, level, v, c, &h[1]))
+  if (!local_height (&q, &p, level, v, c, &h[1])) {
+    g_warning ("Failed to compute local height at (%g,%g)", q.x, q.y);
     return G_MAXDOUBLE;
+  }
 
   hxx = h[0] - 2*H + h[1];
   hx = (h[0] - h[1])/2.;

-- 
Gerris Flow Solver