[cig-commits] [commit] knepley/upgrade-petsc-interface: Visitor: Now create an index for closure of coordinates and fields (0b9b0aa)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/pylith

On branch  : knepley/upgrade-petsc-interface
Link       : https://github.com/geodynamics/pylith/compare/f4faf9408c04ec442346ce0c70a530c12202ae29...0b9b0aaf3cbf745e5b76f3177a1ddddb5bf75fe7

>---------------------------------------------------------------

commit 0b9b0aaf3cbf745e5b76f3177a1ddddb5bf75fe7
Author: Matthew G. Knepley <knepley at gmail.com>
Date:   Sat Nov 23 11:57:53 2013 -0600

    Visitor: Now create an index for closure of coordinates and fields


>---------------------------------------------------------------

0b9b0aaf3cbf745e5b76f3177a1ddddb5bf75fe7
 libsrc/pylith/topology/CoordsVisitor.icc | 1 +
 libsrc/pylith/topology/VisitorMesh.icc   | 1 +
 2 files changed, 2 insertions(+)

diff --git a/libsrc/pylith/topology/CoordsVisitor.icc b/libsrc/pylith/topology/CoordsVisitor.icc
index 7fe6926..1bd7898 100644
--- a/libsrc/pylith/topology/CoordsVisitor.icc
+++ b/libsrc/pylith/topology/CoordsVisitor.icc
@@ -53,6 +53,7 @@ pylith::topology::CoordsVisitor::initialize(void)
   PetscInt       size;
   PetscErrorCode err;
   err = DMPlexGetCoordinateSection(_dm, &_section);PYLITH_CHECK_ERROR(err);assert(_section);
+  err = DMPlexCreateClosureIndex(_dm, _section);PYLITH_CHECK_ERROR(err);
   err = DMGetCoordinatesLocal(_dm, &_localVec);PYLITH_CHECK_ERROR(err);assert(_localVec);
   err = VecGetLocalSize(_localVec, &size);PYLITH_CHECK_ERROR(err);
   err = VecGetArray(_localVec, &_localArray);PYLITH_CHECK_ERROR(err);assert(!size || _localArray);
diff --git a/libsrc/pylith/topology/VisitorMesh.icc b/libsrc/pylith/topology/VisitorMesh.icc
index 7261e60..2ee889f 100644
--- a/libsrc/pylith/topology/VisitorMesh.icc
+++ b/libsrc/pylith/topology/VisitorMesh.icc
@@ -59,6 +59,7 @@ pylith::topology::VecVisitorMesh::initialize(const Field& field)
   _localVec = field.localVector();assert(_localVec);
 
   PetscErrorCode err = VecGetArray(_localVec, &_localArray);PYLITH_CHECK_ERROR(err);
+  err = DMPlexCreateClosureIndex(_dm, _section);PYLITH_CHECK_ERROR(err);
 } // initialize
 
 // ----------------------------------------------------------------------