Added: Permutation of integration points in PlasticMaterial to ensure…

… they

always are evaluated in the same order (or else the history data is corrupted).

git-svn-id: http://svn.sintef.no/trondheim/IFEM/trunk@1944 e10b68d5-8a6e-419e-a041-bce267b0401d

Elasticity.C

@@ -644,6 +644,19 @@ const char* Elasticity::getField2Name (size_t i, const char* prefix) const
 }
 
 
+bool Elasticity::hasItgBuffers () const
+{
+  return material ? material->hasItgBuffers() : false;
+}
+
+
+void Elasticity::setItgPtMap (int ip, int jp) const
+{
+  if (material)
+    material->setItgPtMap(ip,jp);
+}
+
+
 void Elasticity::printMaxVals (std::ostream& os,
                                std::streamsize precision, size_t comp) const
 {

Elasticity.h

@@ -139,6 +139,11 @@ class Elasticity : public IntegrandBase
   //! \param[in] prefix Name prefix for all components
   virtual const char* getField2Name(size_t i, const char* prefix = 0) const;
 
+  //! \brief Returns whether this integrand has internal point buffers or not.
+  virtual bool hasItgBuffers() const;
+  //! \brief Defines the internal processing order of the integration points.
+  virtual void setItgPtMap(int ip, int jp) const;
+
   //! \brief Prints out the maximum secondary solution values.
   //! \param os Output stream to write the values to
   //! \param[in] precision Number of digits after the decimal point

MaterialBase.h

@@ -82,6 +82,11 @@ class Material
 
   //! \brief Returns whether the material model has diverged.
   virtual bool diverged() const { return false; }
+
+  //! \brief Returns whether this material has internal point buffers or not.
+  virtual bool hasItgBuffers() const { return false; }
+  //! \brief Defines the internal processing order of the integration points.
+  virtual void setItgPtMap(size_t, size_t) const {}
 };
 
 #endif