patchkernel: initialize partitioning in the constructor

examples/POD_example_00002.cpp

@@ -93,9 +93,6 @@ void run(int rank, int nProcs)
         int dumpBlock = (nProcs > 1) ? rank : -1;
         std::string filename = "./data/test.0.mesh";
         IBinaryArchive binaryReader(filename, dumpBlock);
-#if BITPIT_ENABLE_MPI	
-        meshr->setCommunicator(MPI_COMM_WORLD);
-#endif
         meshr->restore(binaryReader.getStream());
         binaryReader.close();
     }

src/PABLO/LocalTree.cpp

@@ -102,9 +102,10 @@ namespace bitpit {
     };
 
     /*! Get max depth reached in local tree
+     *  If the tree is empty a negative number is returned.
      * \return Max depth in local partition of the octree.
      */
-    uint8_t
+    int8_t
     LocalTree::getLocalMaxDepth() const{
         return m_localMaxDepth;
     };
@@ -364,14 +365,14 @@ namespace bitpit {
         m_lastGhostBros.clear();
         m_firstGhostBros.clear();
 
-        m_localMaxDepth = 0;
-
         clearConnectivity();
         intervector().swap(m_intersections);
 
         std::fill(m_periodic.begin(), m_periodic.end(), false);
 
-        if (createRoot) {
+        if (m_dim > 0 && createRoot) {
+            m_localMaxDepth = 0;
+
             m_octants.push_back(Octant(m_dim));
 
             Octant firstDesc(TreeConstants::MAX_LEVEL,0,0,0);
@@ -380,6 +381,8 @@ namespace bitpit {
             Octant lastDesc(m_dim,TreeConstants::MAX_LEVEL,TreeConstants::MAX_LENGTH-1,TreeConstants::MAX_LENGTH-1,(m_dim-2)*(TreeConstants::MAX_LENGTH-1));
             m_lastDescMorton = lastDesc.computeMorton();
         } else {
+            m_localMaxDepth = -1;
+
             Octant octDesc(m_dim,TreeConstants::MAX_LEVEL,pow(2,TreeConstants::MAX_LEVEL),pow(2,TreeConstants::MAX_LEVEL),(m_dim > 2 ? pow(2,TreeConstants::MAX_LEVEL) : 0));
             m_lastDescMorton  = octDesc.computeMorton();
             m_firstDescMorton = PABLO::INVALID_MORTON;
@@ -824,6 +827,9 @@ namespace bitpit {
 
         m_sizeOctants = m_octants.size();
 
+        // Update maximum depth
+        updateLocalMaxDepth();
+
         // Set final last desc
         setFirstDescMorton();
         setLastDescMorton();

src/PABLO/LocalTree.hpp

@@ -122,7 +122,7 @@ class LocalTree{
 	uint64_t			 	m_lastDescMorton;		/**< Morton number of last (Morton order) most refined octant possible in local partition */
 	uint32_t 				m_sizeGhosts;			/**< Size of vector of ghost octants */
 	uint32_t 				m_sizeOctants;			/**< Size of vector of local octants */
-	uint8_t					m_localMaxDepth;		/**< Reached max depth in local tree */
+	int8_t					m_localMaxDepth;		/**< Reached max depth in local tree */
 	uint8_t 				m_balanceCodim;			/**<Maximum codimension of the entity for 2:1 balancing (1 = 2:1 balance through faces (default);
 	 	 	 	 	 	 	 	 	 	 	 	 	 	 2 = 2:1 balance through edges and faces;
 	 	 	 	 	 	 	 	 	 	 	 	 	 	 3 = 2:1 balance through nodes, edges and faces)*/
@@ -157,7 +157,7 @@ class LocalTree{
 	uint64_t		getLastDescMorton() const;
 	uint32_t 		getNumGhosts() const;
 	uint32_t 		getNumOctants() const;
-	uint8_t 		getLocalMaxDepth() const;
+	int8_t 			getLocalMaxDepth() const;
 	int8_t 			getMarker(int32_t idx) const;
 	uint8_t 		getLevel(int32_t idx) const;
 	uint64_t 		computeMorton(int32_t idx) const;