[cig-commits] r5214 - short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d

[willic3 at geodynamics.org]

Author: willic3
Date: 2006-11-09 06:11:18 -0800 (Thu, 09 Nov 2006)
New Revision: 5214

Modified:
   short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/mat_6.f
Log:
Fixed silly mistake in computation of tangent matrix.


Modified: short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/mat_6.f
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/mat_6.f	2006-11-09 03:10:06 UTC (rev 5213)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/mat_6.f	2006-11-09 14:11:18 UTC (rev 5214)
@@ -713,14 +713,14 @@
       dsde=f1*(one-dgde*(c2*sdevt(1)+
      & f3*(epptdt(1)-f4*sdevt(1)+ae*sdev0(1))))
       dmat(iddmat(1,1))=c3+two*third*dsde
-      dmat(iddmat(1,2))=c3-half*third*dsde
+      dmat(iddmat(1,2))=c3-third*dsde
       dmat(iddmat(1,3))=dmat(iddmat(1,2))
 c
       dgde=rkfac*(half*epptdt(2)+ae*sdev0(2)-f4*sdevt(2))
       dsde=f1*(one-dgde*(c2*sdevt(2)+
      & f3*(epptdt(2)-f4*sdevt(2)+ae*sdev0(2))))
       dmat(iddmat(2,2))=c3+two*third*dsde
-      dmat(iddmat(2,3))=c3-half*third*dsde
+      dmat(iddmat(2,3))=c3-third*dsde
 c
       dgde=rkfac*(half*epptdt(3)+ae*sdev0(3)-f4*sdevt(3))
       dsde=f1*(one-dgde*(c2*sdevt(3)+