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

[brad at geodynamics.org]

Author: brad
Date: 2010-07-15 09:47:40 -0700 (Thu, 15 Jul 2010)
New Revision: 17044

Modified:
   short/3D/PyLith/branches/v1.5-stable/libsrc/faults/FaultCohesiveLagrange.cc
Log:
Fixed comment.

Modified: short/3D/PyLith/branches/v1.5-stable/libsrc/faults/FaultCohesiveLagrange.cc
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/libsrc/faults/FaultCohesiveLagrange.cc	2010-07-14 15:21:24 UTC (rev 17043)
+++ short/3D/PyLith/branches/v1.5-stable/libsrc/faults/FaultCohesiveLagrange.cc	2010-07-15 16:47:40 UTC (rev 17044)
@@ -587,7 +587,7 @@
   /** We have J = [A C^T]
    *              [C   0]
    *
-   * We want to approximate -( C A^(-1) C^T)^(-1)
+   * We want to approximate -( C A^(-1) C^T )
    *
    * Consider Lagrange vertex L that constrains the relative
    * displacement between vertex N on the negative side of the fault
@@ -730,7 +730,7 @@
     // Compute -[C] [Adiag]^(-1) [C]^T
     //   C_{ij}          = orientationVertex[i*spaceDim+j]
     //   C^T_{ij}        = orientationVertex[j*spaceDim+i]
-    //   Adiag^{-1}_{ii} = jacobianInvVertexN[i] + jacobianInvVertexP[i] (BRAD: Are you sure its not a minus sign here?)
+    //   Adiag^{-1}_{ii} = jacobianInvVertexN[i] + jacobianInvVertexP[i]
     //  \sum_{j} C_{ij} Adiag^{-1}_{jj} C^T_{ji}
     precondVertexL = 0.0;
     for (int kDim=0; kDim < spaceDim; ++kDim) {