[cig-commits] r17088 - short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble

[brad at geodynamics.org]

Author: brad
Date: 2010-08-13 19:36:33 -0700 (Fri, 13 Aug 2010)
New Revision: 17088

Modified:
   short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/IntegratorElasticity.cc
Log:
Added comments on elastic constants.

Modified: short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/IntegratorElasticity.cc
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/IntegratorElasticity.cc	2010-08-13 03:49:12 UTC (rev 17087)
+++ short/3D/PyLith/branches/v1.5-stable/libsrc/feassemble/IntegratorElasticity.cc	2010-08-14 02:36:33 UTC (rev 17088)
@@ -842,8 +842,8 @@
     const double C2211 = elasticConsts[iQuad*numConsts+3];
     const double C2222 = elasticConsts[iQuad*numConsts+4];
     const double C2212 = elasticConsts[iQuad*numConsts+5]/2.0;
-    const double C1211 = elasticConsts[iQuad*numConsts+6]/2.0;
-    const double C1222 = elasticConsts[iQuad*numConsts+7]/2.0;
+    const double C1211 = elasticConsts[iQuad*numConsts+6]/2.0; //?
+    const double C1222 = elasticConsts[iQuad*numConsts+7]/2.0; //?
     const double C1212 = elasticConsts[iQuad*numConsts+8]/2.0;
     for (int iBasis=0, iQ=iQuad*numBasis*spaceDim;
 	 iBasis < numBasis;
@@ -966,21 +966,21 @@
     const double C3312 = elasticConsts[iQuad*numConsts+15] / 2.0;
     const double C3323 = elasticConsts[iQuad*numConsts+16] / 2.0;
     const double C3313 = elasticConsts[iQuad*numConsts+17] / 2.0;
-    const double C1211 = elasticConsts[iQuad*numConsts+18] / 2.0;
-    const double C1222 = elasticConsts[iQuad*numConsts+19] / 2.0;
-    const double C1233 = elasticConsts[iQuad*numConsts+20] / 2.0;
+    const double C1211 = elasticConsts[iQuad*numConsts+18] / 2.0; //?
+    const double C1222 = elasticConsts[iQuad*numConsts+19] / 2.0; //?
+    const double C1233 = elasticConsts[iQuad*numConsts+20] / 2.0; //?
     const double C1212 = elasticConsts[iQuad*numConsts+21] / 2.0;
     const double C1223 = elasticConsts[iQuad*numConsts+22] / 2.0;
     const double C1213 = elasticConsts[iQuad*numConsts+23] / 2.0;
-    const double C2311 = elasticConsts[iQuad*numConsts+24] / 2.0;
-    const double C2322 = elasticConsts[iQuad*numConsts+25] / 2.0;
-    const double C2333 = elasticConsts[iQuad*numConsts+26] / 2.0;
+    const double C2311 = elasticConsts[iQuad*numConsts+24] / 2.0; //?
+    const double C2322 = elasticConsts[iQuad*numConsts+25] / 2.0; //?
+    const double C2333 = elasticConsts[iQuad*numConsts+26] / 2.0; //?
     const double C2312 = elasticConsts[iQuad*numConsts+27] / 2.0;
     const double C2323 = elasticConsts[iQuad*numConsts+28] / 2.0;
     const double C2313 = elasticConsts[iQuad*numConsts+29] / 2.0;
-    const double C1311 = elasticConsts[iQuad*numConsts+30] / 2.0;
-    const double C1322 = elasticConsts[iQuad*numConsts+31] / 2.0;
-    const double C1333 = elasticConsts[iQuad*numConsts+32] / 2.0;
+    const double C1311 = elasticConsts[iQuad*numConsts+30] / 2.0; //?
+    const double C1322 = elasticConsts[iQuad*numConsts+31] / 2.0; //?
+    const double C1333 = elasticConsts[iQuad*numConsts+32] / 2.0; //?
     const double C1312 = elasticConsts[iQuad*numConsts+33] / 2.0;
     const double C1323 = elasticConsts[iQuad*numConsts+34] / 2.0;
     const double C1313 = elasticConsts[iQuad*numConsts+35] / 2.0;