[cig-commits] commit: Cleanup

Mercurial hg at geodynamics.org
Wed Nov 9 00:58:17 PST 2011


changeset:   910:b73798246e20
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Wed Nov 09 00:58:10 2011 -0800
files:       Rheology/src/ConstitutiveMatrix.cxx Rheology/src/ConstitutiveMatrixCartesian.cxx
description:
Cleanup


diff -r 1a2bed9b406a -r b73798246e20 Rheology/src/ConstitutiveMatrix.cxx
--- a/Rheology/src/ConstitutiveMatrix.cxx	Wed Nov 09 00:57:57 2011 -0800
+++ b/Rheology/src/ConstitutiveMatrix.cxx	Wed Nov 09 00:58:10 2011 -0800
@@ -289,11 +289,13 @@ void ConstitutiveMatrix_Assemble(
    MaterialPoint*          materialPoint;
 
    /* Big fat assumption!
-    * Because of Rheology framework vs PIC IP mapping change, ConstitutiveMatrix assumes that
-    * we are using a one-to-one mapping.
-    * This is because the Rheology stuff was made based on operating on MaterialPoints
-    * rather than IntegrationPoints even though its assembling. However it does require xi
-    * which we are passing from the integration point for speed (rather than re-computing it)
+    * Because of Rheology framework vs PIC IP mapping change,
+    * ConstitutiveMatrix assumes that we are using a one-to-one
+    * mapping.  This is because the Rheology stuff was made based on
+    * operating on MaterialPoints rather than IntegrationPoints even
+    * though its assembling. However it does require xi which we are
+    * passing from the integration point for speed (rather than
+    * re-computing it)
     *
     * _Init() firewalls this assumption.
     *
diff -r 1a2bed9b406a -r b73798246e20 Rheology/src/ConstitutiveMatrixCartesian.cxx
--- a/Rheology/src/ConstitutiveMatrixCartesian.cxx	Wed Nov 09 00:57:57 2011 -0800
+++ b/Rheology/src/ConstitutiveMatrixCartesian.cxx	Wed Nov 09 00:58:10 2011 -0800
@@ -360,6 +360,12 @@ void _ConstitutiveMatrixCartesian_Assemb
 
      eta = self->matrixData[2][2];
 
+       // std::cout << "eta "
+       //           << lElement_I << " "
+       //           << cParticle_I << " "
+       //           << eta << " "
+       //           << "\n";
+
      /* Turn D Matrix into D~ Matrix by multiplying in the weight and the detJac (this is a shortcut for speed) */
      ConstitutiveMatrix_MultiplyByValue( constitutiveMatrix, detJac * particle->weight );
 



More information about the CIG-COMMITS mailing list