Commit c943ef74 authored by Alberto's avatar Alberto
Browse files

fix cocoon for Chebyshev

parent c4303818
......@@ -210,20 +210,29 @@ public:
assert(block1.size() > 0);
SizeType site = block1[0];
if (site > 1 && site < numberOfSites - 2)
bool doBorderIfBorder = false;
if (site > 1 && site < numberOfSites - 2) {
this->common().cocoon(block1, direction, doBorderIfBorder); // in-situ
return;
}
if (direction == ProgramGlobals::DirectionEnum::EXPAND_SYSTEM) {
if (site == 1) return;
if (site == 1) {
this->common().cocoon(block1, direction, doBorderIfBorder); // in-situ
return;
}
} else {
if (site == numberOfSites - 2) return;
if (site == numberOfSites - 2) {
this->common().cocoon(block1, direction, doBorderIfBorder); // in-situ
return;
}
}
SizeType x = (site == 1) ? 0 : numberOfSites - 1;
BlockType block(1, x);
evolveInternal(Eg, direction, block, loopNumber);
bool doBorderIfBorder = true;
doBorderIfBorder = true;
this->common().cocoon(block1, direction, doBorderIfBorder); // in-situ
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment