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

[knepley at geodynamics.org]

Author: knepley
Date: 2013-01-22 17:29:49 -0800 (Tue, 22 Jan 2013)
New Revision: 21287

Modified:
   short/3D/PyLith/trunk/libsrc/pylith/topology/Field.cc
Log:
Fixes for Field copy

Modified: short/3D/PyLith/trunk/libsrc/pylith/topology/Field.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/topology/Field.cc	2013-01-23 00:37:28 UTC (rev 21286)
+++ short/3D/PyLith/trunk/libsrc/pylith/topology/Field.cc	2013-01-23 01:29:49 UTC (rev 21287)
@@ -1045,7 +1045,7 @@
   err = DMGetDefaultSection(_dm, &section);CHECK_PETSC_ERROR(err);
   assert(section);assert(_localVec);
   assert(osection);assert(ovec);
-  err = PetscSectionGetNumFields(section, &numFields);CHECK_PETSC_ERROR(err);
+  err = PetscSectionGetNumFields(osection, &numFields);CHECK_PETSC_ERROR(err);
   err = PetscSectionGetChart(section, &pStart, &pEnd);CHECK_PETSC_ERROR(err);
   err = PetscSectionGetChart(osection, &qStart, &qEnd);CHECK_PETSC_ERROR(err);
   if ((pStart != qStart) || (pEnd != qEnd)) {
@@ -1068,7 +1068,7 @@
     throw std::runtime_error(msg.str());
   }
   if (field >= 0) {
-    err = PetscSectionGetFieldComponents(section, field, &numComp);CHECK_PETSC_ERROR(err);
+    err = PetscSectionGetFieldComponents(osection, field, &numComp);CHECK_PETSC_ERROR(err);
     if (component >= numComp) {
       std::ostringstream msg;
       msg << "Invalid field component "<<component<<" should be in [0, "<<numComp<<")";