[cig-commits] r9348 - in short/3D/PyLith/trunk: libsrc/bc pylith/bc

[brad at geodynamics.org]

Author: brad
Date: 2008-02-14 14:15:29 -0800 (Thu, 14 Feb 2008)
New Revision: 9348

Modified:
   short/3D/PyLith/trunk/libsrc/bc/DirichletBoundary.cc
   short/3D/PyLith/trunk/libsrc/bc/Neumann.cc
   short/3D/PyLith/trunk/pylith/bc/DirichletBoundary.py
Log:
Add coordinates to boundary meshes for Neumann and DirichletBoundary. Fixed some bugs.

Modified: short/3D/PyLith/trunk/libsrc/bc/DirichletBoundary.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/bc/DirichletBoundary.cc	2008-02-14 22:13:16 UTC (rev 9347)
+++ short/3D/PyLith/trunk/libsrc/bc/DirichletBoundary.cc	2008-02-14 22:15:29 UTC (rev 9348)
@@ -64,6 +64,8 @@
 	<< "condition '" << _label << "'.";
     throw std::runtime_error(msg.str());
   } // if
+  _boundaryMesh->setRealSection("coordinates", 
+				mesh->getRealSection("coordinates"));
 
   // Get values for degrees of freedom
   char** valueNames = (numFixedDOF > 0) ? new char*[numFixedDOF] : 0;

Modified: short/3D/PyLith/trunk/libsrc/bc/Neumann.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/bc/Neumann.cc	2008-02-14 22:13:16 UTC (rev 9347)
+++ short/3D/PyLith/trunk/libsrc/bc/Neumann.cc	2008-02-14 22:15:29 UTC (rev 9348)
@@ -65,6 +65,8 @@
 	<< "boundary condition '" << _label << "'.";
     throw std::runtime_error(msg.str());
   } // if
+  _boundaryMesh->setRealSection("coordinates", 
+				mesh->getRealSection("coordinates"));
 
   //_boundaryMesh->view("TRACTION BOUNDARY MESH");
 

Modified: short/3D/PyLith/trunk/pylith/bc/DirichletBoundary.py
===================================================================
--- short/3D/PyLith/trunk/pylith/bc/DirichletBoundary.py	2008-02-14 22:13:16 UTC (rev 9347)
+++ short/3D/PyLith/trunk/pylith/bc/DirichletBoundary.py	2008-02-14 22:15:29 UTC (rev 9348)
@@ -149,9 +149,8 @@
     self.boundaryMesh.initialize(self.mesh.coordsys)
     self.cppHandle.boundaryMesh(self.boundaryMesh.cppHandle)
 
-    if None != self.output:
-      self.output.initialize()
-      self.output.writeInfo()
+    self.output.initialize()
+    self.output.writeInfo()
 
     self._logger.eventEnd(logEvent)    
     return