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

fixed: make sure to clear out old blocks in DomainDecomposition

Browse Source 
if not, we end up with bad data in adaptive simulations since
the DomainDecomposition object is reused
This commit is contained in:
Arne Morten Kvarving 2016-11-17 10:30:40 +01:00
parent 18d3d976e2
commit ffb7db391e
1 changed files with 1 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
src/ASM/DomainDecomposition.C
View File

@ -897,6 +897,7 @@ bool DomainDecomposition::setup(const ProcessAdm& adm, const SIMbase& sim)
if (sim.getSolParams() && sim.getSolParams()->getNoBlocks() > 1) { if (sim.getSolParams() && sim.getSolParams()->getNoBlocks() > 1) {
IFEM::cout << "\tEstablishing local block equation numbers" << std::endl; IFEM::cout << "\tEstablishing local block equation numbers" << std::endl;
const LinSolParams& solParams = *sim.getSolParams(); const LinSolParams& solParams = *sim.getSolParams();
blocks.clear();
blocks.resize(solParams.getNoBlocks()+1); blocks.resize(solParams.getNoBlocks()+1);
// Find local equations for each block // Find local equations for each block