OilfieldToolsNet/IFEM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fixed: also put multiplier in block 2 with component based blocks

Browse Source
This commit is contained in:
Arne Morten Kvarving 2021-12-04 14:38:53 +01:00
parent a2cf0932f1
commit ec6b692c3a
1 changed files with 11 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
src/ASM/DomainDecomposition.C
View File

@ -1144,16 +1144,17 @@ bool DomainDecomposition::setup(const ProcessAdm& adm, const SIMbase& sim)
dofType = cb == 1 ? 'D' : 'P'+cb-2;
std::set<int> tmp = adm.dd.getSAM()->getEquations(dofType);
blocks[i+1].localEqs.insert(tmp.begin(), tmp.end());
// HACK: multipliers always in second block
// Correct thing to do for average pressure constraint in Stokes.
if (i == 1) {
for (size_t n = 1; n <= sim.getPatch(1)->getNoNodes(); ++n) {
if (sim.getPatch(1)->isLMn(n) && sim.getPatch(1)->getLMType(n) == 'G') {
int lEq = sam->getEquation(sim.getPatch(1)->getNodeID(n), 1);
if (lEq > 0)
blocks[i+1].localEqs.insert(lEq);
}
}
}
}
// HACK: multipliers always in second block
// Correct thing to do for average pressure constraint in Stokes.
if (i == 1) {
for (size_t n = 1; n <= sim.getPatch(1)->getNoNodes(); ++n) {
if (sim.getPatch(1)->isLMn(n) && sim.getPatch(1)->getLMType(n) == 'G') {
int lEq = sam->getEquation(sim.getPatch(1)->getNodeID(n), 1);
if (lEq > 0)
blocks[i+1].localEqs.insert(lEq);
}
}
}