void
MCLinOp::smooth (MultiFab& solnL,
const MultiFab& rhsL,
int level,
MCBC_Mode bc_mode)
{
for (int phaseflag = 0; phaseflag < numphase; phaseflag++)
{
applyBC(solnL, level, bc_mode);
Fsmooth(solnL, rhsL, level, phaseflag);
}
}
void
LinOp::smooth (MultiFab& solnL,
const MultiFab& rhsL,
int level,
LinOp::BC_Mode bc_mode)
{
for (int redBlackFlag = 0; redBlackFlag < 2; redBlackFlag++)
{
applyBC(solnL, 0, 1, level, bc_mode);
Fsmooth(solnL, rhsL, level, redBlackFlag);
}
}
