[cig-commits] commit: Handle the case where the interface sits exactly on the boundary

[Mercurial]

changeset:   130:44124aed6828
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Wed Mar 28 07:10:11 2012 -0700
files:       compute_v_on_interface/compute_dv_dtt.cxx
description:
Handle the case where the interface sits exactly on the boundary
between cells.


diff -r e6db0bbfac61 -r 44124aed6828 compute_v_on_interface/compute_dv_dtt.cxx
--- a/compute_v_on_interface/compute_dv_dtt.cxx	Wed Mar 28 06:33:20 2012 -0700
+++ b/compute_v_on_interface/compute_dv_dtt.cxx	Wed Mar 28 07:10:11 2012 -0700
@@ -247,6 +247,8 @@ void compute_dv_dtt(const double zx[Nx+1
       std::vector<Valid> values[2];
       for(auto &V: valid)
         {
+          if(V.sign==0)
+            continue;
           const int pm=(V.sign==-1 ? 0 : 1);
           switch(values[pm].size())
             {
@@ -297,6 +299,8 @@ void compute_dv_dtt(const double zx[Nx+1
           std::vector<Valid> derivs[2];
           for(auto &V: d_valid(d0))
             {
+              if(V.sign==0)
+                continue;
               if(V.valid)
                 {
                   const int pm=(V.sign==-1 ? 0 : 1);