[cig-commits] r16533 - short/3D/PyLith/tags/v1.4.3/pylith/topology

[brad at geodynamics.org]

Author: brad
Date: 2010-04-12 15:17:38 -0700 (Mon, 12 Apr 2010)
New Revision: 16533

Modified:
   short/3D/PyLith/tags/v1.4.3/pylith/topology/MeshGenerator.py
Log:
Updated interface for changes to adjustTopology().

Modified: short/3D/PyLith/tags/v1.4.3/pylith/topology/MeshGenerator.py
===================================================================
--- short/3D/PyLith/tags/v1.4.3/pylith/topology/MeshGenerator.py	2010-04-12 20:12:47 UTC (rev 16532)
+++ short/3D/PyLith/tags/v1.4.3/pylith/topology/MeshGenerator.py	2010-04-12 22:17:38 UTC (rev 16533)
@@ -95,20 +95,22 @@
     #mesh.view("===== MESH BEFORE ADJUSTING TOPOLOGY =====")
 
     if not interfaces is None:
-      firstFaultVertex = 0
-      firstFaultCell   = 0
+      firstFaultVertex    = 0
+      firstLagrangeVertex = 0
+      firstFaultCell      = 0
       for interface in interfaces:
         self._info.log("Counting vertices for fault '%s'." % interface.label())
         nvertices = interface.numVertices(mesh)
-        firstFaultCell += nvertices
+        firstLagrangeVertex += nvertices
+        firstFaultCell      += nvertices
         if interface.useLagrangeConstraints():
           firstFaultCell += nvertices
       for interface in interfaces:
         nvertices = interface.numVertices(mesh)
         self._info.log("Adjusting topology for fault '%s' with %d vertices." % \
                          (interface.label(), nvertices))
-        firstFaultVertex, firstFaultCell = \
-            interface.adjustTopology(mesh, firstFaultVertex, firstFaultCell)
+        firstFaultVertex, firstLagrangeVertex, firstFaultCell = \
+            interface.adjustTopology(mesh, firstFaultVertex, firstLagrangeVertex, firstFaultCell)
         
     #mesh.view("===== MESH AFTER ADJUSTING TOPOLOGY =====")
     #self._info.deactivate()