CalculateOnIntegrationPoints

From KratosWiki
(Difference between revisions)
Jump to: navigation, search
Line 1: Line 1:
double LinearElement::CalculateIntegrationWeight(double GaussPointWeight, double DetJ0)
+
  virtual void CalculateOnIntegrationPoints(const Variable<double>& rVariable, Vector& Output, const ProcessInfo&  rCurrentProcessInfo){}
+
 
{
+
  virtual void CalculateOnIntegrationPoints(const Variable<array_1d<double,3> >& rVariable, std::vector< array_1d<double,3> >& Output, const ProcessInfo& rCurrentProcessInfo){}
+
 
//to permorm the integration over the reference domain we need to include
+
    virtual void CalculateOnIntegrationPoints(const Variable<Vector >& rVariable, std::vector< Vector >& Output, const ProcessInfo& rCurrentProcessInfo){}
+
 
  // the thickness in 2D
+
  virtual void CalculateOnIntegrationPoints(const Variable<Matrix >& rVariable, std::vector< Matrix >& Output, const ProcessInfo& rCurrentProcessInfo){}
+
unsigned int dimension = GetGeometry().WorkingSpaceDimension();
+
   
+
double weight = GaussPointWeight;
+
+
     
+
   
+
weight *= DetJ0;
+
   
+
if (dimension == 2) weight *= GetProperties()[THICKNESS];
+
+
return weight;
+

Revision as of 16:10, 3 December 2007

 virtual void CalculateOnIntegrationPoints(const Variable<double>& rVariable, Vector& Output, const ProcessInfo&  rCurrentProcessInfo){}
	  virtual void CalculateOnIntegrationPoints(const Variable<array_1d<double,3> >& rVariable, std::vector<  array_1d<double,3> >& Output, const ProcessInfo& rCurrentProcessInfo){}
	  virtual void CalculateOnIntegrationPoints(const Variable<Vector >& rVariable, std::vector< Vector >& Output,  const ProcessInfo& rCurrentProcessInfo){}
	  virtual void CalculateOnIntegrationPoints(const Variable<Matrix >& rVariable, std::vector< Matrix >& Output,  const ProcessInfo& rCurrentProcessInfo){}
Personal tools
Categories