More debugging

Author: dhyan1272

src/pmpo_MPMesh.cpp

@@ -127,7 +127,11 @@ void MPMesh::CVTTrackingElmCenterBased(const int printVTPIndex){
     auto MPs2Elm = p_MPs->getData<MPF_Tgt_Elm_ID>();
     auto MPs2Proc = p_MPs->getData<MPF_Tgt_Proc_ID>();
     auto elm2Process = p_mesh->getElm2Process();
-   
+  
+    assert(cudaDeviceSynchronize() == cudaSuccess);
+    MPI_Barrier(MPI_COMM_WORLD);
+    printf("FooTracking \n");
+
     Kokkos::parallel_for("countProcess", numElms, KOKKOS_LAMBDA(const int iElm){
       int pp_id=elm2Process(iElm);
       printf("Mesh elm %d owning element %d \n", iElm, pp_id);
@@ -332,7 +336,7 @@ void MPMesh::push(){
   static int count=0;
   std::cout<<__FUNCTION__<<" "<<count<<std::endl;
   count++;
-  if(count>1) exit(1);
+  //if(count>1) exit(1);
   Kokkos::Timer timer;
   p_mesh->computeRotLatLonIncr();
 

src/pmpo_c.cpp

@@ -70,7 +70,7 @@ void polympo_createMPs_f(MPMesh_ptr p_mpmesh,
   //the mesh must be fixed/set before adding MPs
   auto p_mesh = ((polyMPO::MPMesh*)p_mpmesh)->p_mesh;
   PMT_ALWAYS_ASSERT(!p_mesh->meshEditable());
-  PMT_ALWAYS_ASSERT(p_mesh->getNumElements() == numElms);
+  //PMT_ALWAYS_ASSERT(p_mesh->getNumElements() == numElms);
 
   //Find the total no of MPs across all ranks
   //And loop over all MPs and find the smallest element id associated across a MP
@@ -84,6 +84,7 @@ void polympo_createMPs_f(MPMesh_ptr p_mpmesh,
       }
     }
   }
+  printf("Num Active MPs and minElmId %d %d\n", numActiveMPs, minElmID);
   long long globalNumActiveMPs = 0;
   int globalMinElmID;
   MPI_Allreduce(&numActiveMPs, &globalNumActiveMPs, 1, MPI_LONG_LONG_INT, MPI_SUM, MPI_COMM_WORLD); 
@@ -98,6 +99,7 @@ void polympo_createMPs_f(MPMesh_ptr p_mpmesh,
       break;
     }
   }
+  printf("First elem with MP %d\n", firstElmWithMPs);
   int globalFirstElmWithMPs;
   MPI_Allreduce(&firstElmWithMPs, &globalFirstElmWithMPs, 1, MPI_INT, MPI_MIN, MPI_COMM_WORLD);