[MOAB-dev] r1836 - MOAB/trunk

[tautges at mcs.anl.gov]

Author: tautges
Date: 2008-05-22 11:57:29 -0500 (Thu, 22 May 2008)
New Revision: 1836

Modified:
   MOAB/trunk/mbparallelcomm_test.cpp
Log:
Still better reporting of shared entities.


Modified: MOAB/trunk/mbparallelcomm_test.cpp
===================================================================
--- MOAB/trunk/mbparallelcomm_test.cpp	2008-05-22 16:51:23 UTC (rev 1835)
+++ MOAB/trunk/mbparallelcomm_test.cpp	2008-05-22 16:57:29 UTC (rev 1836)
@@ -175,7 +175,7 @@
         // now figure out which vertices are shared
       MBParallelComm *pcomm = new MBParallelComm(mbImpl);
 
-      MBRange shared_ents[4];
+      MBRange shared_ents[6];
       for (int i = 0; i < 4; i++) {
         tmp_result = pcomm->get_shared_entities(i, shared_ents[i]);
       
@@ -185,6 +185,7 @@
           PRINT_LAST_ERROR
               result = tmp_result;
         }
+        if (0 != i) shared_ents[4].merge(shared_ents[i]);
       }
       
       if (0 == rank) setime = MPI_Wtime();
@@ -201,10 +202,15 @@
                 << " succeeded." << std::endl;
 
       if (-1 == nshared) {
+        result = mbImpl->get_adjacencies(shared_ents[4], 0, false, shared_ents[5], 
+                                         MBInterface::UNION);
+        
         std::cerr << "Proc " << rank << " shared entities: " << std::endl;
         for (int i = 0; i < 4; i++)
           std::cerr << "    " << shared_ents[i].size() << " "
                     << i << "d shared entities." << std::endl;
+        std::cerr << "    (" << shared_ents[5].size() 
+                  << " verts adj to other shared ents)" << std::endl;
       }
       
       if (debug && 2 == nprocs) {