Removed debug statements

Author: dhyan1272

src/pmpo_MPMesh.cpp

@@ -23,13 +23,7 @@ void MPMesh::calculateStrain(){
   auto solveStress = p_mesh->getMeshField<polyMPO::MeshF_SolveStress>();
 
   bool isRotated = p_mesh->getRotatedFlag();
-
-  static int strain_count=0;
-  strain_count ++; 
-  bool debug = false;
-  if (strain_count < 3) debug = true;
-  if (debug) printf("Strain count %d ===========\n", strain_count);
-
+  
   auto setMPStrainRate = PS_LAMBDA(const int& elm, const int& mp, const int& mask){
     if(mask){
 
@@ -58,7 +52,7 @@ void MPMesh::calculateStrain(){
       double v22 = 0.0;
       double uTanOverR = 0.0;
       double vTanOverR = 0.0;
-     
+ 
       for (int i = 0; i < numVtx; i++){
         int iVertex = elm2VtxConn(elm, i+1)-1;
         v11 = v11 + MPsBasisGrads(mp, i*2 + 0) * velField(iVertex, 0);
@@ -67,17 +61,11 @@ void MPMesh::calculateStrain(){
         v22 = v22 + MPsBasisGrads(mp, i*2 + 1) * velField(iVertex, 1);
         uTanOverR = uTanOverR + MPsBasis(mp, i) * tanLatVertexRotatedOverRadius(iVertex, 0) * velField(iVertex, 0);
         vTanOverR = vTanOverR + MPsBasis(mp, i) * tanLatVertexRotatedOverRadius(iVertex, 0) * velField(iVertex, 1);
-        if(debug && MPsAppID(mp)==4372){
-          printf("Velocity %.15e %.15e \n", velField(iVertex, 0), velField(iVertex, 1));
-        }
       }
+
       MPsStrainRate(mp, 0) =  v11 - vTanOverR;
       MPsStrainRate(mp, 1) =  v22;
       MPsStrainRate(mp, 2) =  0.5*(v12 + v21 + uTanOverR);
-
-      if(debug && MPsAppID(mp)==4372){
-        printf("Strain %.15e %.15e %.15e \n", MPsStrainRate(mp, 0), MPsStrainRate(mp, 1), MPsStrainRate(mp, 2));
-      }
     }
   };
   p_MPs->parallel_for(setMPStrainRate, "setMPStrainRate");
@@ -97,26 +85,17 @@ void MPMesh::calculateStress(const int constitutive_relation){
   auto MPsIcePressure = p_MPs->getData<polyMPO::MPF_IcePressure>();
   auto MPsRepPressure = p_MPs->getData<polyMPO::MPF_ReplacementPressure>();
 
-  static int stress_count=0;
-  stress_count ++;
-  bool debug = false;
-  if (stress_count < 3) debug = true;
-
   auto setMPStress = PS_LAMBDA(const int& elm, const int& mp, const int& mask){
     if(mask){
       Vec3d strain_rate (MPsStrainRate(mp, 0), MPsStrainRate(mp, 1), MPsStrainRate(mp, 2));
       Vec3d stress(MPsStress(mp, 0), MPsStress(mp, 1), MPsStress(mp, 2));
 
       if (constitutive_relation == 1)
-        constitutive_evp(strain_rate, stress, MPsIcePressure(mp, 0), MPsRepPressure(mp, 0), MPsArea(mp, 0), elasticTimeStep, dampingTimescale);
+        constitutive_evp(strain_rate, stress, MPsIcePressure(mp,0), MPsRepPressure(mp,0), MPsArea(mp,0), elasticTimeStep, dampingTimescale);
       else if(constitutive_relation == 3)
 	constitutive_linear(strain_rate, stress);
       for (int m=0 ; m<3; m++)
         MPsStress(mp, m) = stress[m]*solveStress(elm);
- 
-      if(debug && MPsAppID(mp)==4372){
-        printf("Stress %.15e %.15e %.15e \n", MPsStress(mp, 0), MPsStress(mp, 1), MPsStress(mp, 2));
-      } 
     }
   };
   p_MPs->parallel_for(setMPStress, "setMPStress");
@@ -190,7 +169,6 @@ void MPMesh::calculateStressDivergence(){
                                               2 * tanLatVertexRotatedOverRadius(vID, 0) * w_vtx * MPsStress(mp, 2));
         Kokkos::atomic_add(&divUV_edge(vID, 1), - weight_grads(mp, i*2 + 1)  * MPsStress(mp, 1)  - weight_grads(mp, i*2 + 0) *  MPsStress(mp, 2) +
                                               w_vtx * tanLatVertexRotatedOverRadius(vID, 0) * (MPsStress(mp, 0)- MPsStress(mp, 1)));
-
       }
     }
   };
@@ -216,16 +194,6 @@ void MPMesh::calculateStressDivergence(){
   });
 
   pumipic::RecordTime("Stress_Divergence_Reconstruction" + std::to_string(self), timer.seconds()); 
- 
-  static int strain_count=0;
-  strain_count ++; 
-  if (strain_count >= 3) return; 
-  
-  Kokkos::parallel_for("debug_rec", numVertices, KOKKOS_LAMBDA(const int vtx){
-    if(vtx ==45028)
-      printf("Reconstructed SD %.15e %.15e \n", stressDivergence(vtx, 0), stressDivergence(vtx, 1));
-  });
-
 }
 
 void MPMesh::calcBasis() {

src/pmpo_MPMesh_assembly.hpp

@@ -384,13 +384,6 @@ void MPMesh::assemblyVtx1(){
     communicate_and_take_halo_contributions(meshField, numVertices, numEntries, 1, 1);
   }
   pumipic::RecordTime("Communicate Field Values" + std::to_string(self), timer.seconds());
-
-  //Debugging
-  Kokkos::parallel_for("debug_rec", numVertices, KOKKOS_LAMBDA(const int iVertex){
-    if(iVertex ==45028)
-      printf("Reconstructed velocity for strain rate %.15e %.15e \n", meshField(iVertex, 0), meshField(iVertex, 1));
-  });
-
 }
 
 //Start Communication routine

src/pmpo_materialPoints.hpp

@@ -264,13 +264,6 @@ class MaterialPoints {
           MPsStress(mp,0) = stress_prev[0] * (1 - MPsStressMetric(mp, 1)) + 2 * MPsStressMetric(mp, 1) * stress_prev[2];
           MPsStress(mp,1) = stress_prev[1] + 2*MPsStressMetric(mp, 0) * stress_prev[2];
           MPsStress(mp,2) = stress_prev[2] +   MPsStressMetric(mp, 0) * (stress_prev[0] - stress_prev[1]);
-
-          if(mpAppID(mp)==4372){
-            printf("Current Position %.15e %.15e %.15e \n", curPosXYZ(mp, 0), curPosXYZ(mp, 1), curPosXYZ(mp, 2));
-            printf("Final Position   %.15e %.15e %.15e \n", tgtPosXYZ(mp, 0), tgtPosXYZ(mp, 1), tgtPosXYZ(mp, 2));
-            printf("Final Velocity   %.15e %.15e \n", velMPs(mp, 0), velMPs(mp, 1));
-              
-          }
         } 
       };
       ps::parallel_for(MPs, updateRotLatLon,"updateRotationalLatitudeLongitude"); 

src/pmpo_wachspressBasis.hpp

@@ -34,7 +34,6 @@ void sphericalInterpolation(MPMesh& mpMesh){
         double mpValue = 0.0;
         for(int i=1; i<= numVtx; i++){
           mpValue += meshField(elm2VtxConn(elm,i)-1,entry)*MPsBasis(mp,i-1);
-          //if(MPsAppID(mp) == 4372) printf("Vertex %d F %.15e \n", elm2VtxConn(elm,i)-1, meshField(elm2VtxConn(elm,i)-1,entry) );
         }
         mpField(mp,entry) = mpValue;
       }