[cig-commits] r22351 - short/3D/PyLith/trunk/libsrc/pylith/faults

[brad at geodynamics.org]

Author: brad
Date: 2013-06-17 17:08:22 -0700 (Mon, 17 Jun 2013)
New Revision: 22351

Modified:
   short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveImpulses.cc
Log:
Added calls to complete() on output buffers with only local values since scatterSectionToVector does not seem to copy local values to output Vec as expected.

Modified: short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveImpulses.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveImpulses.cc	2013-06-18 00:06:19 UTC (rev 22350)
+++ short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveImpulses.cc	2013-06-18 00:08:22 UTC (rev 22351)
@@ -217,6 +217,7 @@
     buffer.copy(dispRel);
     buffer.label("slip");
     FaultCohesiveLagrange::globalToFault(&buffer, orientation);
+    buffer.complete();
     PYLITH_METHOD_RETURN(buffer);
 
   } else if (cohesiveDim > 0 && 0 == strcasecmp("strike_dir", name)) {
@@ -253,7 +254,12 @@
 
   } else if (0 == strcasecmp("impulse_amplitude", name)) {
     topology::Field& amplitude = _fields->get("impulse amplitude");
-    PYLITH_METHOD_RETURN(amplitude);
+    _allocateBufferScalarField();
+    topology::Field& buffer = _fields->get("buffer (scalar)");
+    buffer.copy(amplitude);
+    buffer.label("impulse_amplitude");
+    buffer.complete();
+    PYLITH_METHOD_RETURN(buffer);
 
   } else if (0 == strcasecmp("area", name)) {
     topology::Field& area = _fields->get("area");