CarpetLib: Correct error in omit_prolongation_points_when_restricting

Correct the logic error when omit_prolongation_points_when_restricting
is used.  It remains off by default.

darcs-hash:20060504231609-dae7b-d2b89530463b9c67e120426f20e86e6f21211a3d.gz

1 files changed, 12 insertions, 8 deletions

diff --git a/Carpet/CarpetLib/src/dh.cc b/Carpet/CarpetLib/src/dh.cc
index d54f7dee8..937f27410 100644
--- a/Carpet/CarpetLib/src/dh.cc
+++ b/Carpet/CarpetLib/src/dh.cc
@@ -415,15 +415,19 @@ void dh::setup_refinement_restriction_boxes (dh::dboxes & box,
           // TODO: this needs to remove what is sent from other
           // processors as well, not only what is sent from processor
           // c
-          for (iblistvect::const_iterator slvi = box.send_ref_bnd_fine.begin();
-               slvi != box.send_ref_bnd_fine.end(); ++ slvi)
-          {
-            const iblist& sendlist = * slvi;
-            for (iblist::const_iterator sli = sendlist.begin();
-                 sli != sendlist.end(); ++sli)
+          for (int ccc=0; ccc<h.components(rl); ++ccc) {
+            const dh::dboxes& box2 = boxes.at(ml).at(rl).at(ccc);
+            for (iblistvect::const_iterator slvi =
+                   box2.send_ref_bnd_fine.begin();
+                 slvi != box2.send_ref_bnd_fine.end(); ++ slvi)
             {
-              const ibbox& send = * sli;
-              recvs -= send.expand(pss,pss);
+              const iblist& sendlist = * slvi;
+              for (iblist::const_iterator sli = sendlist.begin();
+                   sli != sendlist.end(); ++sli)
+              {
+                const ibbox& send = * sli;
+                recvs -= send.expand(pss,pss);
+              }
             }
           }
         }