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

[brad at geodynamics.org]

Author: brad
Date: 2010-05-13 23:09:06 -0700 (Thu, 13 May 2010)
New Revision: 16703

Modified:
   short/3D/PyLith/trunk/libsrc/materials/DruckerPrager3D.cc
Log:
Hooked up tests for checking whether materials will generate symmetric Jacobian. Added check to see if constraints permit using matrix with block size equal to spatial dimension.

Modified: short/3D/PyLith/trunk/libsrc/materials/DruckerPrager3D.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/materials/DruckerPrager3D.cc	2010-05-14 06:08:15 UTC (rev 16702)
+++ short/3D/PyLith/trunk/libsrc/materials/DruckerPrager3D.cc	2010-05-14 06:09:06 UTC (rev 16703)
@@ -209,7 +209,7 @@
     throw std::runtime_error(msg.str());
   } // if
 
-  if (fabs(frictionAngle != dilatationAngle) > 1.0e-6)
+  if (fabs(frictionAngle - dilatationAngle) > 1.0e-6)
     _isJacobianSymmetric = false;
 
   const double mu = density * vs*vs;