[cig-commits] r16043 - in short/3D/PyLith/branches/pylith-friction: libsrc/faults playpen/friction

[surendra at geodynamics.org]

Author: surendra
Date: 2009-11-25 23:51:12 -0800 (Wed, 25 Nov 2009)
New Revision: 16043

Modified:
   short/3D/PyLith/branches/pylith-friction/libsrc/faults/FaultCohesiveDynL.cc
   short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-axial.mesh
   short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-opening.mesh
   short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4.mesh
Log:
Updated cell-dimension in case of 2-D friction examples

Modified: short/3D/PyLith/branches/pylith-friction/libsrc/faults/FaultCohesiveDynL.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/libsrc/faults/FaultCohesiveDynL.cc	2009-11-26 04:10:35 UTC (rev 16042)
+++ short/3D/PyLith/branches/pylith-friction/libsrc/faults/FaultCohesiveDynL.cc	2009-11-26 07:51:12 UTC (rev 16043)
@@ -855,12 +855,18 @@
 
 	    // Sensitivity of slip to changes in the Lagrange multipliers
 	    // Aixjx = 1.0/Aix + 1.0/Ajx
+	    assert(jacobianVertex[0] > 0.0);
+	    assert(jacobianVertex[spaceDim+0] > 0.0);
 	    const double Aixjx = 
 	      1.0/jacobianVertex[0] + 1.0/jacobianVertex[spaceDim+0];
 	    // Aiyjy = 1.0/Aiy + 1.0/Ajy
+	    assert(jacobianVertex[1] > 0.0);
+	    assert(jacobianVertex[spaceDim+1] > 0.0);
 	    const double Aiyjy = 
 	      1.0/jacobianVertex[1] + 1.0/jacobianVertex[spaceDim+1];
 	    // Aizjz = 1.0/Aiz + 1.0/Ajz
+	    assert(jacobianVertex[2] > 0.0);
+	    assert(jacobianVertex[spaceDim+2] > 0.0);
 	    const double Aizjz = 
 	      1.0/jacobianVertex[2] + 1.0/jacobianVertex[spaceDim+2];
 	    const double Cpx = orientationVertex[0];

Modified: short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-axial.mesh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-axial.mesh	2009-11-26 04:10:35 UTC (rev 16042)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-axial.mesh	2009-11-26 07:51:12 UTC (rev 16043)
@@ -31,12 +31,12 @@
     // Vertex number (starting from zero), x-coord, y-coord
     // Use coordinate units that are consistent with the other units used.
     coordinates = {
-             0     -2.0 -1.0
-             1     -2.0  1.0
-             2      0.0 -1.0
-             3      0.0  1.0
-             4      2.0 -1.0
-             5      2.0  1.0
+             0     -1.0 -0.5
+             1     -1.0  0.5
+             2      0.0 -0.5
+             3      0.0  0.5
+             4      1.0 -0.5
+             5      1.0  0.5
     }
   }
 

Modified: short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-opening.mesh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-opening.mesh	2009-11-26 04:10:35 UTC (rev 16042)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-opening.mesh	2009-11-26 07:51:12 UTC (rev 16043)
@@ -31,12 +31,12 @@
     // Vertex number (starting from zero), x-coord, y-coord
     // Use coordinate units that are consistent with the other units used.
     coordinates = {
-             0     -2.0 -1.0
-             1     -2.0  1.0
-             2      0.0 -1.0
-             3      0.0  1.0
-             4      2.0 -1.0
-             5      2.0  1.0
+             0     -1.0 -0.5
+             1     -1.0  0.5
+             2      0.0 -0.5
+             3      0.0  0.5
+             4      1.0 -0.5
+             5      1.0  0.5
     }
   }
 

Modified: short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4.mesh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4.mesh	2009-11-26 04:10:35 UTC (rev 16042)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4.mesh	2009-11-26 07:51:12 UTC (rev 16043)
@@ -31,12 +31,12 @@
     // Vertex number (starting from zero), x-coord, y-coord
     // Use coordinate units that are consistent with the other units used.
     coordinates = {
-             0     -2.0 -1.0
-             1     -2.0  1.0
-             2      0.0 -1.0
-             3      0.0  1.0
-             4      2.0 -1.0
-             5      2.0  1.0
+             0     -1.0 -0.5
+             1     -1.0  0.5
+             2      0.0 -0.5
+             3      0.0  0.5
+             4      1.0 -0.5
+             5      1.0  0.5
     }
   }