[cig-commits] r3825 - short/3D/PyLith/branches/pylith-0.8/pylith3d/module

[knepley at geodynamics.org]

Author: knepley
Date: 2006-06-20 12:04:01 -0700 (Tue, 20 Jun 2006)
New Revision: 3825

Modified:
   short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc
Log:
Fixed problem with viscos VTK output


Modified: short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc	2006-06-20 17:57:51 UTC (rev 3824)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc	2006-06-20 19:04:01 UTC (rev 3825)
@@ -479,12 +479,15 @@
   ALE::Obj<ALE::Mesh::field_type::order_type::baseSequence> patches = displacement->getPatches();
   ALE::Obj<ALE::Mesh::foliation_type> boundaries = m->getBoundaries();
 
-  for(ALE::Mesh::field_type::order_type::baseSequence::iterator p_iter = patches->begin(); p_iter != patches->end(); ++p_iter) {
-    full_displacement->setPatch(displacement->getPatch(*p_iter), *p_iter);
-    full_displacement->setFiberDimensionByDepth(*p_iter, 0, 3);
+  // This is wrong if the domain changes
+  if (!full_displacement->getGlobalOrder()) {
+    for(ALE::Mesh::field_type::order_type::baseSequence::iterator p_iter = patches->begin(); p_iter != patches->end(); ++p_iter) {
+      full_displacement->setPatch(displacement->getPatch(*p_iter), *p_iter);
+      full_displacement->setFiberDimensionByDepth(*p_iter, 0, 3);
+    }
+    full_displacement->orderPatches();
+    full_displacement->createGlobalOrder();
   }
-  full_displacement->orderPatches();
-  full_displacement->createGlobalOrder();
   for(ALE::Mesh::field_type::order_type::baseSequence::iterator p_iter = patches->begin(); p_iter != patches->end(); ++p_iter) {
     ALE::Obj<ALE::Mesh::field_type::order_type::coneSequence> elements = full_displacement->getPatch(*p_iter);