KratosMultiphysics
KRATOS Multiphysics (Kratos) is a framework for building parallel, multi-disciplinary simulation software, aiming at modularity, extensibility, and high performance. Kratos is written in C++, and counts with an extensive Python interface.
helmholtz_bulk_shape_element.h
Go to the documentation of this file.
1 // | / |
2 // ' / __| _` | __| _ \ __|
3 // . \ | ( | | ( |\__ `
4 // _|\_\_| \__,_|\__|\___/ ____/
5 // Multi-Physics
6 //
7 // License: BSD License
8 // license: OptimizationApplication/license.txt
9 //
10 // Main authors: Reza Najian Asl
11 //
12 
13 #if !defined(KRATOS_HELMHOLTZ_BULK_SHAPE_ELEMENT_H_INCLUDED )
14 #define KRATOS_HELMHOLTZ_BULK_SHAPE_ELEMENT_H_INCLUDED
15 
16 // System includes
17 
18 // External includes
19 
20 // Project includes
21 #include "includes/define.h"
22 #include "includes/element.h"
24 #include "utilities/geometry_utilities.h"
26 #include "includes/variables.h"
28 
29 
30 namespace Kratos
31 {
32 
35 
39 
43 
47 
51 
53 
55 class KRATOS_API(OPTIMIZATION_APPLICATION) HelmholtzBulkShapeElement
56  : public Element
57 {
58 public:
61 
64 
68 
70  HelmholtzBulkShapeElement(IndexType NewId, GeometryType::Pointer pGeometry);
71  HelmholtzBulkShapeElement(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
72 
74  virtual ~HelmholtzBulkShapeElement();
75 
76 
80 
81 
85 
86  Element::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, PropertiesType::Pointer pProperties) const override;
87 
88  Element::Pointer Create(IndexType NewId, GeometryType::Pointer pGeom, PropertiesType::Pointer pProperties) const override;
89 
90  void CalculateLocalSystem(MatrixType& rLeftHandSideMatrix, VectorType& rRightHandSideVector, const ProcessInfo& rCurrentProcessInfo) override;
91 
92  void CalculateLeftHandSide(MatrixType& rLeftHandSideMatrix, const ProcessInfo& rCurrentProcessInfo) override;
93 
94  void CalculateRightHandSide(VectorType& rRightHandSideVector, const ProcessInfo& rCurrentProcessInfo) override;
95 
96  void EquationIdVector(EquationIdVectorType& rResult, const ProcessInfo& rCurrentProcessInfo) const override;
97 
98  void GetDofList(DofsVectorType& ElementalDofList, const ProcessInfo& CurrentProcessInfo) const override;
99 
100  void GetValuesVector(VectorType &rValues, int Step = 0) const override;
101 
102  void Calculate(const Variable<Matrix>& rVariable, Matrix& rOutput, const ProcessInfo& rCurrentProcessInfo) override;
103 
104  void Calculate(const Variable<double>& rVariable, double& rOutput, const ProcessInfo& rCurrentProcessInfo) override;
105 
106  int Check(const ProcessInfo& rCurrentProcessInfo) const override;
107 
111 
112 
116 
120 
124 
125 
127 
128 protected:
131 
132 
136 
137 
141 
142 
146 
147 
151 
152 
156 
157 
161 
162  // Protected default constructor necessary for serialization
164  {
165  }
166 
168 
169 private:
172 
173 
177 
178 
182  friend class Serializer;
183 
184  void save(Serializer& rSerializer) const override
185  {
187  }
188 
189  void load(Serializer& rSerializer) override
190  {
192  }
193 
197 
198 
202  MatrixType SetAndModifyConstitutiveLaw(const int Dimension, const int PointNumber) const;
203  MatrixType CalculateBMatrix(const int Dimension, const int PointNumber) const;
204  void CalculateBulkMassMatrix(MatrixType& rMassMatrix,const ProcessInfo& rCurrentProcessInfo) const;
205  void CalculateBulkStiffnessMatrix(MatrixType& rStiffnessMatrix,const ProcessInfo& rCurrentProcessInfo) const;
206 
210 
211 
215 
216 
220 
222  //HelmholtzBulkShapeElement& operator=(const HelmholtzBulkShapeElement& rOther);
223 
225  //HelmholtzBulkShapeElement(const HelmholtzBulkShapeElement& rOther);
226 
227 
229 
230 }; // Class HelmholtzBulkShapeElement
231 
233 
236 
237 
241 
242 
244 /* inline std::istream& operator >> (std::istream& rIStream,
245  HelmholtzBulkShapeElement& rThis);
246 */
248 /* inline std::ostream& operator << (std::ostream& rOStream,
249  const HelmholtzBulkShapeElement& rThis)
250  {
251  rThis.PrintInfo(rOStream);
252  rOStream << std::endl;
253  rThis.PrintData(rOStream);
254 
255  return rOStream;
256  }*/
258 
259 } // namespace Kratos.
260 
261 #endif // KRATOS_HELMHOLTZ_BULK_SHAPE_ELEMENT_H_INCLUDED defined
262 
263 
Base class for all Elements.
Definition: element.h:60
std::vector< DofType::Pointer > DofsVectorType
Definition: element.h:100
std::vector< std::size_t > EquationIdVectorType
Definition: element.h:98
std::size_t IndexType
Definition: flags.h:74
Short class definition.
Definition: helmholtz_bulk_shape_element.h:57
KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(HelmholtzBulkShapeElement)
Counted pointer of HelmholtzBulkShapeElement.
HelmholtzBulkShapeElement()
Definition: helmholtz_bulk_shape_element.h:163
PointerVector is a container like stl vector but using a vector to store pointers to its data.
Definition: pointer_vector.h:72
ProcessInfo holds the current value of different solution parameters.
Definition: process_info.h:59
The serialization consists in storing the state of an object into a storage format like data file or ...
Definition: serializer.h:123
Variable class contains all information needed to store and retrive data from a data container.
Definition: variable.h:63
#define KRATOS_SERIALIZE_SAVE_BASE_CLASS(Serializer, BaseType)
Definition: define.h:812
#define KRATOS_SERIALIZE_LOAD_BASE_CLASS(Serializer, BaseType)
Definition: define.h:815
Matrix MatrixType
Definition: geometrical_transformation_utilities.h:55
Modeler::Pointer Create(const std::string &ModelerName, Model &rModel, const Parameters ModelParameters)
Checks if the modeler is registered.
Definition: modeler_factory.cpp:30
TDataType Calculate(GeometryType &dummy, const Variable< TDataType > &rVariable)
Definition: add_geometries_to_python.cpp:103
REF: G. R. Cowper, GAUSSIAN QUADRATURE FORMULAS FOR TRIANGLES.
Definition: mesh_condition.cpp:21
ProcessInfo
Definition: edgebased_PureConvection.py:116
def load(f)
Definition: ode_solve.py:307