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

[brad at geodynamics.org]

Author: brad
Date: 2007-04-24 11:40:36 -0700 (Tue, 24 Apr 2007)
New Revision: 6652

Modified:
   short/3D/PyLith/trunk/libsrc/faults/BruneSlipFn.cc
Log:
Switched to specifying parameters and slip field on section over fault mesh (should be slightly more efficient).

Modified: short/3D/PyLith/trunk/libsrc/faults/BruneSlipFn.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/BruneSlipFn.cc	2007-04-24 18:19:33 UTC (rev 6651)
+++ short/3D/PyLith/trunk/libsrc/faults/BruneSlipFn.cc	2007-04-24 18:40:36 UTC (rev 6652)
@@ -70,7 +70,8 @@
   const ALE::Obj<Mesh::label_sequence>& vertices = faultMesh->depthStratum(0);
 
   // Create sections for fields
-  delete _parameters; _parameters = new feassemble::ParameterManager(mesh);
+  delete _parameters; 
+  _parameters = new feassemble::ParameterManager(faultMesh);
   if (0 == _parameters)
     throw std::runtime_error("Could not create manager for parameters of "
 			     "Brune slip time function.");
@@ -104,9 +105,9 @@
     slipTime->setFiberDimension(*v_iter, 1);
     peakRate->setFiberDimension(*v_iter, 1);
   } // for
-  mesh->allocate(finalSlip);
-  mesh->allocate(slipTime);
-  mesh->allocate(peakRate);
+  faultMesh->allocate(finalSlip);
+  faultMesh->allocate(slipTime);
+  faultMesh->allocate(peakRate);
 
   // Open databases and set query values
   _dbFinalSlip->open();
@@ -182,7 +183,7 @@
        v_iter != vEnd;
        ++v_iter)
     _slipField->setFiberDimension(*v_iter, 3);
-  mesh->allocate(_slipField);
+  faultMesh->allocate(_slipField);
 } // initialize
 
 // ----------------------------------------------------------------------