[cig-commits] commit: Grow a box in fix_viscosity so that viscosities that are not in a

[Mercurial]

changeset:   247:70f53147e2fb
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Thu May 03 11:58:46 2012 -0700
files:       src/FACStokes/fix_viscosity.C
description:
Grow a box in fix_viscosity so that viscosities that are not in a
"node" have ghost values set.


diff -r f089bf615532 -r 70f53147e2fb src/FACStokes/fix_viscosity.C
--- a/src/FACStokes/fix_viscosity.C	Thu May 05 13:22:05 2011 -0700
+++ b/src/FACStokes/fix_viscosity.C	Thu May 03 11:58:46 2012 -0700
@@ -105,7 +105,10 @@ void SAMRAI::FACStokes::fix_viscosity()
               for(int axis=0;axis<3;++axis)
                 {
                   const int axis2((axis+1)%3), axis3((axis+2)%3);
-                  for(pdat::EdgeIterator ni(edge_viscosity.getBox(),axis); ni; ni++)
+                  hier::Box pbox=patch->getBox();
+                  pbox.grow(axis,edge_viscosity.getGhostCellWidth()[axis]);
+
+                  for(pdat::EdgeIterator ni(pbox,axis); ni; ni++)
                     {
                       pdat::EdgeIndex e=ni();
                       pdat::CellIndex c(e);