MueLu: Fixing off-by-one error

trilinos · Jul 17, 2024 · cd9fc2a · cd9fc2a
1 parent 3f96163
commit cd9fc2a
Showing 1 changed file with 2 additions and 3 deletions.
diff --git a/packages/muelu/src/Transfers/Smoothed-Aggregation/MueLu_ReitzingerPFactory_def.hpp b/packages/muelu/src/Transfers/Smoothed-Aggregation/MueLu_ReitzingerPFactory_def.hpp
@@ -291,7 +291,6 @@ void ReitzingerPFactory<Scalar, LocalOrdinal, GlobalOrdinal, Node>::BuildP(Level
   // Handle empty ranks gracefully
   if (num_coarse_edges == 0) {
     D0_rowptr[0] = 0;
-    D0_rowptr[1] = 0;
   }
 
   // Count the total number of edges
@@ -344,7 +343,7 @@ void ReitzingerPFactory<Scalar, LocalOrdinal, GlobalOrdinal, Node>::BuildP(Level
         for(int i=0; i<(int)ja.size(); i++)
           fprintf(f,"%d ",(int)ja[i]);
         fclose(f);
-        
+
       }
 #endif
     D0_coarse->expertStaticFillComplete(ownedCoarseNodeMap, ownedCoarseEdgeMap);
@@ -361,7 +360,7 @@ void ReitzingerPFactory<Scalar, LocalOrdinal, GlobalOrdinal, Node>::BuildP(Level
   RCP<Matrix> Pe;
   {
     SubFactoryMonitor m2(*this, "Generate Pe (pre-fix)", coarseLevel);
-#if 0   
+#if 0
     {
       // If you're concerned about processor / rank mismatches, this debugging code might help
       int rank =  D0->getRowMap()->getComm()->getRank();