[cig-commits] r16228 - short/3D/PyLith/trunk/libsrc/materials

[brad at geodynamics.org]

Author: brad
Date: 2010-02-03 16:31:49 -0800 (Wed, 03 Feb 2010)
New Revision: 16228

Added:
   short/3D/PyLith/trunk/libsrc/materials/ElasticIsotropic3D.icc
Log:
More optimization for lumped solve. Also added defines to control level of logging detail.

Added: short/3D/PyLith/trunk/libsrc/materials/ElasticIsotropic3D.icc
===================================================================
--- short/3D/PyLith/trunk/libsrc/materials/ElasticIsotropic3D.icc	                        (rev 0)
+++ short/3D/PyLith/trunk/libsrc/materials/ElasticIsotropic3D.icc	2010-02-04 00:31:49 UTC (rev 16228)
@@ -0,0 +1,36 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+//                           Brad T. Aagaard
+//                        U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+#if !defined(pylith_materials_elasticisotropic3d_hh)
+#error "ElasticIsotropic3D.icc can only be included from ElasticIsotropic3D.hh"
+#endif
+
+#include <cassert> // USES assert()
+
+// Compute density at location from properties.
+inline
+void
+pylith::materials::ElasticIsotropic3D::_calcDensity(double* const density,
+						    const double* properties,
+						    const int numProperties,
+						    const double* stateVars,
+						    const int numStateVars)
+{ // _calcDensity
+  assert(0 != density);
+  assert(0 != properties);
+  assert(_numPropsQuadPt == numProperties);
+  assert(0 == numStateVars);
+
+  density[0] = properties[p_density];
+} // _calcDensity
+
+// End of file