[cig-commits] r8054 - short/3D/PyLith/trunk/libsrc/faults

[willic3 at geodynamics.org]

Author: willic3
Date: 2007-09-28 12:03:02 -0700 (Fri, 28 Sep 2007)
New Revision: 8054

Modified:
   short/3D/PyLith/trunk/libsrc/faults/BruneSlipFn.cc
Log:
Put in a very small number of missing flops.
Almost everything is Sieve operations.


Modified: short/3D/PyLith/trunk/libsrc/faults/BruneSlipFn.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/BruneSlipFn.cc	2007-09-28 18:56:39 UTC (rev 8053)
+++ short/3D/PyLith/trunk/libsrc/faults/BruneSlipFn.cc	2007-09-28 19:03:02 UTC (rev 8054)
@@ -229,6 +229,7 @@
   double_array slipValues(3);
   const vert_iterator vBegin = vertices.begin();
   const vert_iterator vEnd = vertices.end();
+  const int vSize = vertices.size();
   for (vert_iterator v_iter=vBegin; v_iter != vEnd; ++v_iter) {
     // Get values of parameters at vertex
     const int numSlipValues = finalSlip->getFiberDimension(*v_iter);
@@ -251,6 +252,7 @@
     // Update field
     _slipField->updatePoint(*v_iter, &slipValues[0]);
   } // for
+  PetscLogFlopsNoCheck(vSize * (4 + numSlipValues));
 
   return _slipField;
 } // slip
@@ -283,6 +285,7 @@
   double_array slipValues(3);
   const vert_iterator vBegin = vertices.begin();
   const vert_iterator vEnd = vertices.end();
+  const int vSize = vertices.size();
   for (vert_iterator v_iter=vBegin; v_iter != vEnd; ++v_iter) {
     // Get values of parameters at vertex
     const int numSlipValues = finalSlip->getFiberDimension(*v_iter);
@@ -306,6 +309,7 @@
     // Update field
     _slipField->updatePoint(*v_iter, &slipValues[0]);
   } // for
+  PetscLogFlopsNoCheck(vSize * (2 * numSlipValues));
 
   return _slipField;
 } // slipIncr