[cig-commits] r6679 - short/3D/PyLith/trunk/libsrc/feassemble

[brad at geodynamics.org]

Author: brad
Date: 2007-04-24 16:00:42 -0700 (Tue, 24 Apr 2007)
New Revision: 6679

Modified:
   short/3D/PyLith/trunk/libsrc/feassemble/ExplicitElasticity.cc
Log:
Moved consistency check of dimension of cell and space dimension out of loop.

Modified: short/3D/PyLith/trunk/libsrc/feassemble/ExplicitElasticity.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/feassemble/ExplicitElasticity.cc	2007-04-24 22:42:49 UTC (rev 6678)
+++ short/3D/PyLith/trunk/libsrc/feassemble/ExplicitElasticity.cc	2007-04-24 23:00:42 UTC (rev 6679)
@@ -97,6 +97,16 @@
   const int spaceDim = _quadrature->spaceDim();
   const int cellDim = _quadrature->cellDim();
 
+  /** :TODO:
+   *
+   * If cellDim and spaceDim are different, we need to transform
+   * displacements into cellDim, compute action, and transform result
+   * back into spaceDim. We get this information from the Jacobian and
+   * inverse of the Jacobian.
+   */
+  if (cellDim != spaceDim)
+    throw std::logic_error("Not implemented yet.");
+
   // Allocate vector for cell values (if necessary)
   _initCellVector();
 
@@ -163,16 +173,6 @@
     if (err)
       throw std::runtime_error("Logging PETSc flops failed.");
     
-    /** :TODO:
-     *
-     * If cellDim and spaceDim are different, we need to transform
-     * displacements into cellDim, compute action, and transform
-     * result back into spaceDim. Can we get this from the inverse of
-     * the Jacobian?
-     */
-    if (cellDim != spaceDim)
-      throw std::logic_error("Not implemented yet.");
-
     // Compute action for elastic terms
     if (1 == cellDim) {
       // Compute stresses