[cig-commits] r8047 - short/3D/PyLith/trunk/libsrc/materials

[willic3 at geodynamics.org]

Author: willic3
Date: 2007-09-27 17:30:01 -0700 (Thu, 27 Sep 2007)
New Revision: 8047

Modified:
   short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc
Log:
Finished putting in logging.


Modified: short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc	2007-09-27 23:48:37 UTC (rev 8046)
+++ short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc	2007-09-28 00:30:01 UTC (rev 8047)
@@ -532,6 +532,7 @@
      (*parameters)[_MaxwellIsotropic3D::pidStrainT][1] +
      (*parameters)[_MaxwellIsotropic3D::pidStrainT][2])/3.0;
   
+  PetscLogFlopsNoCheck(6);
   // The code below should probably be in a separate function since it
   // is used more than once.  I should also probably cover the possibility
   // that Maxwell time is zero (although this should never happen).
@@ -549,13 +550,17 @@
       fraction *= _dt/maxwelltime;
       dq += fSign*fraction/factorial;
     } // for
-  } else
+    PetscLogFlopsNoCheck(1 + 7 * numTerms);
+  } else {
     dq = maxwelltime*(1.0-exp(-_dt/maxwelltime))/_dt;
+    PetscLogFlopsNoCheck(7);
+  } // else
 
   const double expFac = exp(-_dt/maxwelltime);
   double devStrainTpdt = 0.0;
   double devStrainT = 0.0;
   double visStrain = 0.0;
+  PetscLogFlopsNoCheck(3);
   for (int iComp=0; iComp < _MaxwellIsotropic3D::tensorSize; ++iComp) {
     devStrainTpdt = totalStrain[iComp] - diag[iComp]*meanStrainTpdt;
     devStrainT = (*parameters)[_MaxwellIsotropic3D::pidStrainT][iComp] -
@@ -566,6 +571,7 @@
     (*parameters)[_MaxwellIsotropic3D::pidVisStrain][iComp] = visStrain;
     (*parameters)[_MaxwellIsotropic3D::pidStrainT][iComp] = totalStrain[iComp];
   } // for
+  PetscLogFlopsNoCheck(8 * MaxwellIsotropic3D::tensorSize);
 
   _needNewJacobian = false;