aboutsummaryrefslogtreecommitdiff
path: root/Carpet/CarpetLib/src
diff options
context:
space:
mode:
authorErik Schnetter <schnetter@aei.mpg.de>2005-01-02 17:34:00 +0000
committerErik Schnetter <schnetter@aei.mpg.de>2005-01-02 17:34:00 +0000
commit43d0a63ab79490aa39fc560b199f4956554869ca (patch)
tree2b983dcbe5d0a1ec44ce0d89848b305aacbada7c /Carpet/CarpetLib/src
parent02f6c4848a06350d26dec929a48b96acad6ec878 (diff)
CarpetLib: Add more checks to dh class
darcs-hash:20050102173453-891bb-833515dd47ce1469ebe319718ab169e5eb82c6c4.gz
Diffstat (limited to 'Carpet/CarpetLib/src')
-rw-r--r--Carpet/CarpetLib/src/dh.cc10
1 files changed, 10 insertions, 0 deletions
diff --git a/Carpet/CarpetLib/src/dh.cc b/Carpet/CarpetLib/src/dh.cc
index 446425946..fff52c35c 100644
--- a/Carpet/CarpetLib/src/dh.cc
+++ b/Carpet/CarpetLib/src/dh.cc
@@ -425,6 +425,7 @@ void dh::assert_assert_assert( dh::dboxes & box, int rl, int c, int ml )
} else { // rl!=0
const iblistvect& recv_ref_coarse = box.recv_ref_coarse;
ibset intr = box.interior;
+ ibset received;
for (iblistvect::const_iterator
lvi=recv_ref_coarse.begin();
lvi!=recv_ref_coarse.end(); ++lvi) {
@@ -436,6 +437,8 @@ void dh::assert_assert_assert( dh::dboxes & box, int rl, int c, int ml )
const int new_sz = intr.size();
// TODO
assert (new_sz + this_sz == old_sz);
+ assert ((received & *li).empty());
+ received |= *li;
}
}
// TODO
@@ -455,6 +458,7 @@ void dh::assert_assert_assert( dh::dboxes & box, int rl, int c, int ml )
const iblistvect& recv_sync = box.recv_sync;
const iblistvect& recv_ref_bnd_coarse = box.recv_ref_bnd_coarse;
ibset bnds = box.boundaries;
+ ibset received;
for (iblistvect::const_iterator lvi=recv_sync.begin();
lvi!=recv_sync.end(); ++lvi) {
for (iblist::const_iterator li=lvi->begin();
@@ -464,6 +468,8 @@ void dh::assert_assert_assert( dh::dboxes & box, int rl, int c, int ml )
bnds -= *li;
const int new_sz = bnds.size();
assert (new_sz + this_sz == old_sz);
+ assert ((received & *li).empty());
+ received |= *li;
}
}
for (iblistvect::const_iterator
@@ -480,6 +486,10 @@ void dh::assert_assert_assert( dh::dboxes & box, int rl, int c, int ml )
// prolongation went into the buffer zone.
// assert (new_sz + this_sz == old_sz);
assert (new_sz + this_sz >= old_sz);
+#if 0
+ assert ((received & *li).empty());
+ received |= *li;
+#endif
}
}
// TODO