[cig-commits] r16958 - in short/3D/PyLith/trunk: . doc/userguide libsrc/faults libsrc/topology

[brad at geodynamics.org]

Author: brad
Date: 2010-06-09 11:35:34 -0700 (Wed, 09 Jun 2010)
New Revision: 16958

Removed:
   short/3D/PyLith/trunk/doc/userguide/pylith_cover.pdf
Modified:
   short/3D/PyLith/trunk/TODO
   short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc
   short/3D/PyLith/trunk/libsrc/topology/Field.cc
Log:
Only create fibration for Lagrange multipliers if necessary.

Modified: short/3D/PyLith/trunk/TODO
===================================================================
--- short/3D/PyLith/trunk/TODO	2010-06-09 17:32:09 UTC (rev 16957)
+++ short/3D/PyLith/trunk/TODO	2010-06-09 18:35:34 UTC (rev 16958)
@@ -6,10 +6,10 @@
 
   need cygwin binbot
 
+  Add ML to petsc on Darwin and cygwin binbots (requires fortran compiler)
+
   need Darwin 10.6 binary
 
-  Add ML to petsc on all binbots (requires fortran compiler)
-
 top-level README [add fault preconditioner when done]
 release announcement [add fault preconditioner when done]
 
@@ -68,9 +68,6 @@
 
     Add time step information
 
-  Savage-Prescott
-    change .cfg file to output velocities, update post-processing
-
 CODE
 
   * Better preconditioning [Matt/Brad]

Deleted: short/3D/PyLith/trunk/doc/userguide/pylith_cover.pdf
===================================================================
(Binary files differ)

Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc	2010-06-09 17:32:09 UTC (rev 16957)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc	2010-06-09 18:35:34 UTC (rev 16958)
@@ -148,6 +148,10 @@
   const int spaceDim = field->mesh().dimension();
   const int fibrationLagrange = spaceDim;
 
+  // Add space for Lagrange multipliers if it does not yet exist.
+  if (spaceDim == section->getNumSpaces())
+    section->addSpace();
+
   const int numVertices = _cohesiveVertices.size();
   for (int iVertex=0; iVertex < numVertices; ++iVertex) {
     const int v_lagrange = _cohesiveVertices[iVertex].lagrange;

Modified: short/3D/PyLith/trunk/libsrc/topology/Field.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/topology/Field.cc	2010-06-09 17:32:09 UTC (rev 16957)
+++ short/3D/PyLith/trunk/libsrc/topology/Field.cc	2010-06-09 18:35:34 UTC (rev 16958)
@@ -742,7 +742,7 @@
 } // scatterVectorToSection
 
 // ----------------------------------------------------------------------
-// Setup split field with all entries set to a default space of 0.
+// Setup split field with all one space per spatial dimension.
 template<typename mesh_type>
 void 
 pylith::topology::Field<mesh_type>::splitDefault(void)
@@ -751,7 +751,7 @@
   const int spaceDim = _mesh.dimension();
   for (int iDim=0; iDim < spaceDim; ++iDim)
     _section->addSpace(); // displacements
-  _section->addSpace(); // Lagrange multipliers
+  //_section->addSpace(); // Lagrange multipliers
 
   const chart_type& chart = _section->getChart();