[cig-commits] r16483 - short/3D/PyLith/trunk/unittests/libtests/faults/data

[brad at geodynamics.org]

Author: brad
Date: 2010-04-01 12:52:46 -0700 (Thu, 01 Apr 2010)
New Revision: 16483

Modified:
   short/3D/PyLith/trunk/unittests/libtests/faults/data/cohesivedyn.py
Log:
A little work on adding sensitivity solve stuff.

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/data/cohesivedyn.py
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/data/cohesivedyn.py	2010-04-01 18:53:41 UTC (rev 16482)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/data/cohesivedyn.py	2010-04-01 19:52:46 UTC (rev 16483)
@@ -64,7 +64,18 @@
                                  [8.1, 9.1, -10.1]])
         area = numpy.array([1.0/3.0, 1.0/3.0, 1.0/3.0])
     
+        jacobianN = numpy.array(
+            [[1.0, 10.0, 10.1, 10.2, 10.3, 10.4, 10.5, 10.6, 10.7],
+             [13.1, 1.0, 13.2, 13.3, 13.4, 13.5, 13.6, 13.7, 13.8], 
+             [16.2, 16.3, 1.0, 16.4, 16.5, 16.6, 16.7, 16.8, 16.9],
+             [19.3, 19.4, 19.5, 1.0, 19.6, 19.7, 19.8, 19.9, 20.0],
+             [22.4, 22.5, 22.6, 22.7, 1.0, 22.8, 22.9, 23.0, 23.1],
+             [25.5, 25.6, 25.7, 25.8, 25.9, 1.0, 26.0, 26.1, 26.2],
+             [28.6, 28.7, 28.8, 28.9, 29.0, 29.1, 1.0, 29.2, 29.3],
+             [31.7, 31.8, 31.9, 32.0, 32.1, 32.2, 32.3, 1.0, 32.4],
+             [34.8, 34.9, 35.0, 35.1, 35.2, 35.3, 35.4, 35.5, 1.0]])
 
+
     elif cell == "hex8":
         fieldT = numpy.array([[5.4, 7.4, 9.4],
                               [5.6, 7.6, 9.6],