Some more clean up

Author: dhyan1272

src/pmpo_MPMesh.cpp

@@ -2,7 +2,6 @@
 #include "pmpo_utils.hpp"
 #include "pmpo_MPMesh.hpp"
 #include "pmpo_wachspressBasis.hpp"
-#include <unistd.h>
 namespace polyMPO{
 
 void printVTP_mesh(MPMesh& mpMesh, int printVTPIndex=-1);
@@ -129,7 +128,6 @@ void MPMesh::CVTTrackingElmCenterBased(const int printVTPIndex){
     auto elm2Process = p_mesh->getElm2Process();
     auto elm2global = p_mesh->getElmGlobal();
 
-    //Since Mesh is static print pnly for 1 time step 
     if(printVTPIndex>=0) {
       printVTP_mesh(printVTPIndex);
     }
@@ -167,10 +165,6 @@ void MPMesh::CVTTrackingElmCenterBased(const int printVTPIndex){
                         minDistSq = neighborDistSq;
                     }
                 }
-                if(abs(mpTgtPos(mp,0)+5.724649188625998e+06) < 1e-10)
-                  printf("Track: %.15e %.15e %.15e => %.15e %.15e %.15e \n", mpTgtPos(mp,0), mpTgtPos(mp,1), mpTgtPos(mp,2), 
-                                                                             mpPositions(mp,0), mpPositions(mp,1), mpPositions(mp,2) );
-
                 if(closestElm<0){
                     MPs2Elm(mp) = iElm;
                     if (elm2Process.size() > 0)

src/pmpo_MPMesh_assembly.hpp

@@ -103,7 +103,7 @@ void MPMesh::resetPreComputeFlag(){
   isPreComputed = false;
 }
 
-//Method 1
+//Method 1 for coefficients
 void MPMesh::computeMatricesAndSolve(){
   Kokkos::Timer timer;
   //Mesh Information
@@ -209,7 +209,7 @@ void MPMesh::computeMatricesAndSolve(){
   pumipic::RecordTime("PolyMPO_Calculate_MLS_Coeff", timer.seconds());
 }
 
-//Method 2
+//Method 2 for coefficients
 void MPMesh::subAssemblyCoeffs(int dim1, int dim2, double* m11, double* m12, double* m13, double* m14, 
                                                    double* m22, double* m23, double* m24, 
                                                    double* m33, double* m34, 
@@ -299,21 +299,12 @@ void MPMesh::subAssemblyCoeffs(int dim1, int dim2, double* m11, double* m12, dou
 
 }
 
-//Method 2
-void MPMesh::solveMatrixAndRegularize( int dim1, double* m11, double* m12, double* m13, double* m14, 
+//Method 2 for coefficients Solve matrix
+void MPMesh::solveMatrixAndRegularize(int dim1, double* m11, double* m12, double* m13, double* m14, 
                                        double* m22, double* m23, double* m24, 
                                        double* m33, double* m34,
                                        double* m44){
 
-  MPI_Comm comm = p_MPs->getMPIComm(); 
-  int comm_rank;
-  MPI_Comm_rank(comm, &comm_rank);
-
-  static int count=0;
-  if(!comm_rank)
-    std::cout<<__FUNCTION__<<count<<std::endl;
-  count ++; 
-
   auto dual_triangle_area=p_mesh->getMeshField<MeshF_DualTriangleArea>();
 
   kkViewHostU<const double*> m11_h(m11, dim1);

src/pmpo_c.cpp

@@ -265,10 +265,6 @@ void polympo_getMPTgtElmID_f(MPMesh_ptr p_mpmesh,
   kkIntViewHostU arrayHost(elmIDs,numMPs);
   polyMPO::IntView mpTgtElmIDCopy("mpTgtElmIDNewValue",numMPs);
 
-  int rank;
-  auto mpi_comm=p_MPs->getMPIComm();
-  MPI_Comm_rank(mpi_comm, &rank);
-   
   auto setTgtElmId = PS_LAMBDA(const int&, const int& mp, const int& mask){
     if(mask){
         mpTgtElmIDCopy(mpAppID(mp)) = mpTgtElmID(mp)+elmIDoffset;
@@ -294,10 +290,6 @@ void polympo_getMPCurElmID_f(MPMesh_ptr p_mpmesh,
   kkIntViewHostU arrayHost(elmIDs,numMPs);
   polyMPO::IntView mpCurElmIDCopy("mpCurElmIDNewValue",numMPs);
 
-  int rank;
-  auto mpi_comm=p_MPs->getMPIComm();
-  MPI_Comm_rank(mpi_comm, &rank);
-   
   auto getElmId = PS_LAMBDA(const int&, const int& mp, const int& mask){
     if(mask){
         mpCurElmIDCopy(mpAppID(mp)) = mpCurElmID(mp)+elmIDoffset;