[cig-commits] r15208 - short/3D/PyLith/trunk/libsrc/topology

[brad at geodynamics.org]

Author: brad
Date: 2009-06-12 08:49:24 -0700 (Fri, 12 Jun 2009)
New Revision: 15208

Modified:
   short/3D/PyLith/trunk/libsrc/topology/Field.cc
Log:
Fixed bug in traversing chart.

Modified: short/3D/PyLith/trunk/libsrc/topology/Field.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/topology/Field.cc	2009-06-12 06:35:10 UTC (rev 15207)
+++ short/3D/PyLith/trunk/libsrc/topology/Field.cc	2009-06-12 15:49:24 UTC (rev 15208)
@@ -474,10 +474,12 @@
     for (typename chart_type::const_iterator c_iter = chartBegin;
 	 c_iter != chartEnd;
 	 ++c_iter) {
-      assert(fiberDim == field._section->getFiberDimension(*c_iter));
-      assert(fiberDim == _section->getFiberDimension(*c_iter));
-      field._section->restrictPoint(*c_iter, &values[0], values.size());
-      _section->updatePointAllAdd(*c_iter, &values[0]);
+      if (0 != field._section->getFiberDimension(*c_iter)) {
+	assert(fiberDim == field._section->getFiberDimension(*c_iter));
+	assert(fiberDim == _section->getFiberDimension(*c_iter));
+	field._section->restrictPoint(*c_iter, &values[0], values.size());
+	_section->updatePointAllAdd(*c_iter, &values[0]);
+      } // if
     } // for
   } // if