[cig-commits] r16808 - in short/3D/PyLith/trunk/playpen: . faultpc

[brad at geodynamics.org]

Author: brad
Date: 2010-05-27 09:53:44 -0700 (Thu, 27 May 2010)
New Revision: 16808

Added:
   short/3D/PyLith/trunk/playpen/faultpc/
   short/3D/PyLith/trunk/playpen/faultpc/checkfaultpc.py
Log:
Added faultpc.

Added: short/3D/PyLith/trunk/playpen/faultpc/checkfaultpc.py
===================================================================
--- short/3D/PyLith/trunk/playpen/faultpc/checkfaultpc.py	                        (rev 0)
+++ short/3D/PyLith/trunk/playpen/faultpc/checkfaultpc.py	2010-05-27 16:53:44 UTC (rev 16808)
@@ -0,0 +1,36 @@
+import numpy
+import numpy.linalg as linalg
+
+Afull = numpy.array([[1,1, 1,1, 1,1, 1,1, 0,0, 0,0, 0,0, 0,0],
+                     [1,1, 1,1, 1,1, 1,1, 0,0, 0,0, 0,0, 0,0],
+                     [1,1, 1,1, 1,1, 0,0, 0,0, 0,0, 0,0, 0,0],
+                     [1,1, 1,1, 1,1, 0,0, 0,0, 0,0, 0,0, 0,0],
+                     [1,1, 1,1, 1,1, 1,1, 0,0, 0,0, 0,0, 0,0],
+                     [1,1, 1,1, 1,1, 1,1, 0,0, 0,0, 0,0, 0,0],
+                     [1,1, 0,0, 1,1, 1,1, 0,0, 0,0, 0,0, 0,0],
+                     [1,1, 0,0, 1,1, 1,1, 0,0, 0,0, 0,0, 0,0],
+                     [0,0, 0,0, 0,0, 0,0, 1,1, 1,1, 0,0, 1,1],
+                     [0,0, 0,0, 0,0, 0,0, 1,1, 1,1, 0,0, 1,1],
+                     [0,0, 0,0, 0,0, 0,0, 1,1, 1,1, 1,1, 1,1],
+                     [0,0, 0,0, 0,0, 0,0, 1,1, 1,1, 1,1, 1,1],
+                     [0,0, 0,0, 0,0, 0,0, 0,0, 1,1, 1,1, 1,1],
+                     [0,0, 0,0, 0,0, 0,0, 0,0, 1,1, 1,1, 1,1],
+                     [0,0, 0,0, 0,0, 0,0, 1,1, 1,1, 1,1, 1,1],
+                     [0,0, 0,0, 0,0, 0,0, 1,1, 1,1, 1,1, 1,1]],
+                    dtype=numpy.float64)
+
+#Ainv = linalg.inv(Afull)
+
+C = numpy.array([[0,0, 0,-1, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0],
+                 [0,0, -1,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0],
+                 [0,0, 0,0, 0,-1, 0,0, 0,0, 0,1, 0,0, 0,0],
+                 [0,0, 0,0, -1,0, 0,0, 0,0, 1,0, 0,0, 0,0],
+                 [0,0, 0,0, 0,0, 0,-1, 0,0, 0,0, 0,1, 0,0],
+                 [0,0, 0,0, 0,0, -1,0, 0,0, 0,0, 1,0, 0,0]],
+                dtype=numpy.float64)
+CAC = numpy.dot(numpy.dot(C, Afull), C.transpose())
+
+
+print "Afull",Afull
+#print "Afull inverse",Ainv
+print "CACfull",CAC