small change

AMDiS/CMakeLists.txt

@@ -684,6 +684,11 @@ INSTALL(FILES ${HEADERS}
 	DESTINATION include/amdis/solver)
 list(APPEND deb_add_dirs "include/amdis/solver")
 
+FILE(GLOB HEADERS "${SOURCE_DIR}/solver/details/*.h")
+INSTALL(FILES ${HEADERS} 
+        DESTINATION include/amdis/solver/details)
+list(APPEND deb_add_dirs "include/amdis/solver/details")
+
 FILE(GLOB HEADERS "${SOURCE_DIR}/solver/itl/*.hpp")
 INSTALL(FILES ${HEADERS} 
 	DESTINATION include/amdis/solver/itl)

AMDiS/src/est/ResidualEstimator.cc

@@ -30,6 +30,7 @@
 #ifdef HAVE_PARALLEL_DOMAIN_AMDIS
 #include <mpi.h>
 #include "parallel/MeshDistributor.h"
+#include "parallel/ParallelDebug.h"
 #endif
 
 using namespace std;
@@ -184,6 +185,11 @@ namespace AMDiS {
 
     Parallel::InteriorBoundary &intBoundary = 
       Parallel::MeshDistributor::globalMeshDistributor->getIntBoundary(0);
+      
+#if (DEBUG != 0)
+    // Make sure interior boundary is correct
+    Parallel::ParallelDebug::testInteriorBoundary(*(Parallel::MeshDistributor::globalMeshDistributor));
+#endif
 
     ElInfo *elInfo = mesh->createNewElInfo();
     elInfo->setFillFlag(Mesh::FILL_COORDS);
@@ -268,6 +274,7 @@ namespace AMDiS {
       if (subBound.size() == 0)
 	continue;
 
+      // Highly possible mesh is not correct, for example: hanging node
       TEST_EXIT_DBG(subBound.size() == stdMpiDet.getRecvData(it->first).size())
 	("Should not happen: %d %d from rank %d\n", subBound.size(), stdMpiDet.getRecvData(it->first).size(), it->first);