From ffb7db391edec1ee402dddddb5aa8d1c508c8be0 Mon Sep 17 00:00:00 2001
From: Arne Morten Kvarving <arne.morten.kvarving@sintef.no>
Date: Thu, 17 Nov 2016 10:30:40 +0100
Subject: [PATCH] fixed: make sure to clear out old blocks in
 DomainDecomposition

if not, we end up with bad data in adaptive simulations since
the DomainDecomposition object is reused
---
 src/ASM/DomainDecomposition.C | 1 +
 1 file changed, 1 insertion(+)

diff --git a/src/ASM/DomainDecomposition.C b/src/ASM/DomainDecomposition.C
index 915520ac..3eae5c9b 100644
--- a/src/ASM/DomainDecomposition.C
+++ b/src/ASM/DomainDecomposition.C
@@ -897,6 +897,7 @@ bool DomainDecomposition::setup(const ProcessAdm& adm, const SIMbase& sim)
   if (sim.getSolParams() && sim.getSolParams()->getNoBlocks() > 1) {
     IFEM::cout << "\tEstablishing local block equation numbers" << std::endl;
     const LinSolParams& solParams = *sim.getSolParams();
+    blocks.clear();
     blocks.resize(solParams.getNoBlocks()+1);
 
     // Find local equations for each block