v0.10.0
Public Member Functions | Public Attributes | Protected Member Functions | Protected Attributes | List of all members
OpK Struct Reference

assemble the stiffness matrix More...

#include <users_modules/linear_isotropic_hardening/src/LinearIsotropicHardening.hpp>

Inheritance diagram for OpK:
[legend]
Collaboration diagram for OpK:
[legend]

Public Member Functions

 OpK (CommonData &common_data, bool symm=true)
 
MoFEMErrorCode doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)
 Do calculations for give operator. More...
 
 OpK (bool symm=true)
 
MoFEMErrorCode doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)
 Do calculations for give operator. More...
 
- Public Member Functions inherited from MoFEM::VolumeElementForcesAndSourcesCoreBase::UserDataOperator
int getNumNodes ()
 get element number of nodes More...
 
const EntityHandlegetConn ()
 get element connectivity More...
 
double getVolume () const
 element volume (linear geometry) More...
 
doublegetVolume ()
 element volume (linear geometry) More...
 
FTensor::Tensor2< double *, 3, 3 > & getJac ()
 get element Jacobian More...
 
FTensor::Tensor2< double *, 3, 3 > & getInvJac ()
 get element inverse Jacobian More...
 
double getMeasure () const
 get measure of element More...
 
doublegetMeasure ()
 get measure of element More...
 
VectorDoublegetCoords ()
 nodal coordinates More...
 
MatrixDoublegetCoordsAtGaussPts ()
 Gauss points and weight, matrix (nb. of points x 3) More...
 
MatrixDoublegetHoCoordsAtGaussPts ()
 coordinate at Gauss points (if hierarchical approximation of element geometry) More...
 
MatrixDoublegetHoGaussPtsJac ()
 
MatrixDoublegetHoGaussPtsInvJac ()
 
VectorDoublegetHoGaussPtsDetJac ()
 
auto getFTenosr0HoMeasure ()
 
auto getFTensor1CoordsAtGaussPts ()
 Get coordinates at integration points assuming linear geometry. More...
 
auto getFTensor1HoCoordsAtGaussPts ()
 Get coordinates at integration points taking geometry from field. More...
 
auto getFTensor2HoGaussPtsJac ()
 
auto getFTensor2HoGaussPtsInvJac ()
 
VolumeElementForcesAndSourcesCoreBasegetVolumeFE () const
 return pointer to Generic Volume Finite Element object More...
 
MoFEMErrorCode getDivergenceOfHDivBaseFunctions (int side, EntityType type, DataForcesAndSourcesCore::EntData &data, int gg, VectorDouble &div)
 Get divergence of base functions at integration point. More...
 
MoFEMErrorCode getCurlOfHCurlBaseFunctions (int side, EntityType type, DataForcesAndSourcesCore::EntData &data, int gg, MatrixDouble &curl)
 Get curl of base functions at integration point. More...
 
- Public Member Functions inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
 UserDataOperator (const FieldSpace space, const char type=OPLAST, const bool symm=true)
 
 UserDataOperator (const std::string &field_name, const char type, const bool symm=true)
 
 UserDataOperator (const std::string &row_field_name, const std::string &col_field_name, const char type, const bool symm=true)
 
boost::shared_ptr< const NumeredEntFiniteElementgetNumeredEntFiniteElementPtr () const
 Return raw pointer to NumeredEntFiniteElement. More...
 
EntityHandle getFEEntityHandle () const
 Return finite element entity handle. More...
 
boost::weak_ptr< SideNumber > getSideNumberPtr (const int side_number, const EntityType type)
 Get the side number pointer. More...
 
EntityHandle getSideEntity (const int side_number, const EntityType type)
 Get the side entity. More...
 
int getNumberOfNodesOnElement ()
 Get the number of nodes on finite element. More...
 
MoFEMErrorCode getProblemRowIndices (const std::string filed_name, const EntityType type, const int side, VectorInt &indices) const
 Get row indices. More...
 
MoFEMErrorCode getProblemColIndices (const std::string filed_name, const EntityType type, const int side, VectorInt &indices) const
 Get col indices. More...
 
const FEMethodgetFEMethod () const
 Return raw pointer to Finite Element Method object. More...
 
int getOpType () const
 Get operator types. More...
 
void setOpType (const OpType type)
 Set operator type. More...
 
void addOpType (const OpType type)
 Add operator type. More...
 
int getNinTheLoop () const
 get number of finite element in the loop More...
 
int getLoopSize () const
 get size of elements in the loop More...
 
const std::string & getFEName () const
 Get name of the element. More...
 
const PetscData::SwitchesgetDataCtx () const
 
const KspMethod::KSPContext getKSPCtx () const
 
const SnesMethod::SNESContext getSNESCtx () const
 
const TSMethod::TSContext getTSCtx () const
 
Vec getKSPf () const
 
Mat getKSPA () const
 
Mat getKSPB () const
 
Vec getSNESf () const
 
Vec getSNESx () const
 
Mat getSNESA () const
 
Mat getSNESB () const
 
DEPRECATED Vec getSnesF () const
 
DEPRECATED Vec getSnesX () const
 
DEPRECATED Mat getSnesA () const
 
DEPRECATED Mat getSnesB () const
 
Vec getTSu () const
 
Vec getTSu_t () const
 
Vec getTSu_tt () const
 
Vec getTSf () const
 
Mat getTSA () const
 
Mat getTSB () const
 
int getTSstep () const
 
double getTStime () const
 
double getTSa () const
 
double getTSaa () const
 
MatrixDoublegetGaussPts ()
 matrix of integration (Gauss) points for Volume Element More...
 
auto getFTensor0IntegrationWeight ()
 Get integration weights. More...
 
DEPRECATED MoFEMErrorCode getPorblemRowIndices (const std::string filed_name, const EntityType type, const int side, VectorInt &indices) const
 
- Public Member Functions inherited from MoFEM::DataOperator
 DataOperator (const bool symm=true)
 
virtual ~DataOperator ()=default
 
virtual MoFEMErrorCode opLhs (DataForcesAndSourcesCore &row_data, DataForcesAndSourcesCore &col_data)
 
virtual MoFEMErrorCode doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)
 Operator for linear form, usually to calculate values on right hand side. More...
 
virtual MoFEMErrorCode opRhs (DataForcesAndSourcesCore &data, const bool error_if_no_base=false)
 
bool getSymm () const
 Get if operator uses symmetry of DOFs or not. More...
 
void setSymm ()
 set if operator is executed taking in account symmetry More...
 
void unSetSymm ()
 unset if operator is executed for non symmetric problem More...
 

Public Attributes

MatrixDouble K
 
CommonDatacommonData
 
FTensor::Ddg< double, 3, 3 > tD
 
double yOung
 
double pOisson
 
- Public Attributes inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
char opType
 
std::string rowFieldName
 
std::string colFieldName
 
FieldSpace sPace
 
- Public Attributes inherited from MoFEM::DataOperator
bool sYmm
 If true assume that matrix is symmetric structure. More...
 
std::array< bool, MBMAXTYPE > doEntities
 If true operator is executed for entity. More...
 
booldoVertices
 \deprectaed If false skip vertices More...
 
booldoEdges
 \deprectaed If false skip edges More...
 
booldoQuads
 \deprectaed More...
 
booldoTris
 \deprectaed More...
 
booldoTets
 \deprectaed More...
 
booldoPrisms
 \deprectaed More...
 

Protected Member Functions

MoFEMErrorCode iNtegrate (DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)
 
MoFEMErrorCode aSsemble (DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)
 
MoFEMErrorCode iNtegrate (DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)
 Integrate B^T D B operator. More...
 
MoFEMErrorCode aSsemble (DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)
 Assemble local entity block matrix. More...
 
- Protected Member Functions inherited from MoFEM::VolumeElementForcesAndSourcesCoreBase::UserDataOperator
MoFEMErrorCode setPtrFE (ForcesAndSourcesCore *ptr)
 
- Protected Member Functions inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
ForcesAndSourcesCoregetPtrFE () const
 
ForcesAndSourcesCoregetSidePtrFE () const
 

Protected Attributes

int nbRows
 number of dofs on rows More...
 
int nbCols
 number if dof on column More...
 
int nbIntegrationPts
 number of integration points More...
 
bool isDiag
 true if this block is on diagonal More...
 
- Protected Attributes inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
ForcesAndSourcesCoreptrFE
 

Additional Inherited Members

- Public Types inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
enum  OpType { OPROW = 1 << 0, OPCOL = 1 << 1, OPROWCOL = 1 << 2, OPLAST = 1 << 3 }
 Controls loop over entities on element. More...
 

Detailed Description

assemble the stiffness matrix

Examples
dynamic_elastic.cpp, eigen_elastic.cpp, PoissonOperators.hpp, and simple_elasticity.cpp.

Definition at line 59 of file LinearIsotropicHardening.hpp.

Constructor & Destructor Documentation

◆ OpK() [1/2]

OpK::OpK ( CommonData common_data,
bool  symm = true 
)

Definition at line 66 of file LinearIsotropicHardening.hpp.

68  symm),
69  commonData(common_data) {}

◆ OpK() [2/2]

OpK::OpK ( bool  symm = true)

Definition at line 42 of file simple_elasticity.cpp.

44  symm) {
45 
46  // Evaluation of the elastic stiffness tensor, D
47 
48  // hardcoded choice of elastic parameters
49  pOisson = 0.1;
50  yOung = 10;
51 
52  // coefficient used in intermediate calculation
53  const double coefficient = yOung / ((1 + pOisson) * (1 - 2 * pOisson));
54 
59 
60  tD(i, j, k, l) = 0.;
61 
62  tD(0, 0, 0, 0) = 1 - pOisson;
63  tD(1, 1, 1, 1) = 1 - pOisson;
64  tD(2, 2, 2, 2) = 1 - pOisson;
65 
66  tD(0, 1, 0, 1) = 0.5 * (1 - 2 * pOisson);
67  tD(0, 2, 0, 2) = 0.5 * (1 - 2 * pOisson);
68  tD(1, 2, 1, 2) = 0.5 * (1 - 2 * pOisson);
69 
70  tD(0, 0, 1, 1) = pOisson;
71  tD(1, 1, 0, 0) = pOisson;
72  tD(0, 0, 2, 2) = pOisson;
73  tD(2, 2, 0, 0) = pOisson;
74  tD(1, 1, 2, 2) = pOisson;
75  tD(2, 2, 1, 1) = pOisson;
76 
77  tD(i, j, k, l) *= coefficient;
78 
79  }

Member Function Documentation

◆ aSsemble() [1/2]

MoFEMErrorCode OpK::aSsemble ( DataForcesAndSourcesCore::EntData row_data,
DataForcesAndSourcesCore::EntData col_data 
)
protected

Definition at line 208 of file LinearIsotropicHardening.hpp.

209  {
211  // get pointer to first global index on row
212  const int *row_indices = &*row_data.getIndices().data().begin();
213  // get pointer to first global index on column
214  const int *col_indices = &*col_data.getIndices().data().begin();
215  // Mat B = getFEMethod()->ksp_B != PETSC_NULL ? getFEMethod()->ksp_B
216  // : getFEMethod()->snes_B;
217  // assemble local matrix
218  CHKERR MatSetValues(getFEMethod()->snes_B, nbRows, row_indices, nbCols, col_indices,
219  &*K.data().begin(), ADD_VALUES);
220 
221  if (!isDiag && sYmm) {
222  // if not diagonal term and since global matrix is symmetric assemble
223  // transpose term.
224  K = trans(K);
225  CHKERR MatSetValues(getFEMethod()->snes_B, nbCols, col_indices, nbRows, row_indices,
226  &*K.data().begin(), ADD_VALUES);
227  }
229  }

◆ aSsemble() [2/2]

MoFEMErrorCode OpK::aSsemble ( DataForcesAndSourcesCore::EntData row_data,
DataForcesAndSourcesCore::EntData col_data 
)
protected

Assemble local entity block matrix.

Parameters
row_datarow data (consist base functions on row entity)
col_datacolumn data (consist base functions on column entity)
Returns
error code

Definition at line 220 of file simple_elasticity.cpp.

221  {
223  // get pointer to first global index on row
224  const int *row_indices = &*row_data.getIndices().data().begin();
225  // get pointer to first global index on column
226  const int *col_indices = &*col_data.getIndices().data().begin();
227  Mat B = getFEMethod()->ksp_B != PETSC_NULL ? getFEMethod()->ksp_B
228  : getFEMethod()->snes_B;
229  // assemble local matrix
230  CHKERR MatSetValues(B, nbRows, row_indices, nbCols, col_indices,
231  &*K.data().begin(), ADD_VALUES);
232 
233  if (!isDiag && sYmm) {
234  // if not diagonal term and since global matrix is symmetric assemble
235  // transpose term.
236  K = trans(K);
237  CHKERR MatSetValues(B, nbCols, col_indices, nbRows, row_indices,
238  &*K.data().begin(), ADD_VALUES);
239  }
241  }

◆ doWork() [1/2]

MoFEMErrorCode OpK::doWork ( int  row_side,
int  col_side,
EntityType  row_type,
EntityType  col_type,
DataForcesAndSourcesCore::EntData row_data,
DataForcesAndSourcesCore::EntData col_data 
)

Do calculations for give operator.

Parameters
row_siderow side number (local number) of entity on element
col_sidecolumn side number (local number) of entity on element
row_typetype of row entity MBVERTEX, MBEDGE, MBTRI or MBTET
col_typetype of column entity MBVERTEX, MBEDGE, MBTRI or MBTET
row_datadata for row
col_datadata for column
Returns
error code

Definition at line 81 of file LinearIsotropicHardening.hpp.

84  {
85 
87 
88  // get number of dofs on row
89  nbRows = row_data.getIndices().size();
90  // if no dofs on row, exit that work, nothing to do here
91  if (!nbRows)
93 
94  // get number of dofs on column
95  nbCols = col_data.getIndices().size();
96  // if no dofs on columns, exit nothing to do here
97  if (!nbCols)
99 
100  // K_ij matrix will have 3 times the number of degrees of freedom of the
101  // i-th entity set (nbRows)
102  // and 3 times the number of degrees of freedom of the j-th entity set
103  // (nbCols)
104  K.resize(nbRows, nbCols, false);
105  K.clear();
106 
107  // get number of integration points
108  nbIntegrationPts = getGaussPts().size2();
109  // check if entity block is on matrix diagonal
110  if (row_side == col_side && row_type == col_type) {
111  isDiag = true;
112  } else {
113  isDiag = false;
114  }
115 
116  // integrate local matrix for entity block
117  CHKERR iNtegrate(row_data, col_data);
118 
119  // assemble local matrix
120  CHKERR aSsemble(row_data, col_data);
121 
123  }

◆ doWork() [2/2]

MoFEMErrorCode OpK::doWork ( int  row_side,
int  col_side,
EntityType  row_type,
EntityType  col_type,
DataForcesAndSourcesCore::EntData row_data,
DataForcesAndSourcesCore::EntData col_data 
)
virtual

Do calculations for give operator.

Parameters
row_siderow side number (local number) of entity on element
col_sidecolumn side number (local number) of entity on element
row_typetype of row entity MBVERTEX, MBEDGE, MBTRI or MBTET
col_typetype of column entity MBVERTEX, MBEDGE, MBTRI or MBTET
row_datadata for row
col_datadata for column
Returns
error code

Reimplemented from MoFEM::DataOperator.

Definition at line 91 of file simple_elasticity.cpp.

94  {
95 
97 
98  // get number of dofs on row
99  nbRows = row_data.getIndices().size();
100  // if no dofs on row, exit that work, nothing to do here
101  if (!nbRows)
103 
104  // get number of dofs on column
105  nbCols = col_data.getIndices().size();
106  // if no dofs on Columbia, exit nothing to do here
107  if (!nbCols)
109 
110  // K_ij matrix will have 3 times the number of degrees of freedom of the
111  // i-th entity set (nbRows)
112  // and 3 times the number of degrees of freedom of the j-th entity set
113  // (nbCols)
114  K.resize(nbRows, nbCols, false);
115  K.clear();
116 
117  // get number of integration points
118  nbIntegrationPts = getGaussPts().size2();
119  // check if entity block is on matrix diagonal
120  if (row_side == col_side && row_type == col_type) {
121  isDiag = true;
122  } else {
123  isDiag = false;
124  }
125 
126  // integrate local matrix for entity block
127  CHKERR iNtegrate(row_data, col_data);
128 
129  // assemble local matrix
130  CHKERR aSsemble(row_data, col_data);
131 
133  }

◆ iNtegrate() [1/2]

MoFEMErrorCode OpK::iNtegrate ( DataForcesAndSourcesCore::EntData row_data,
DataForcesAndSourcesCore::EntData col_data 
)
protected

Definition at line 133 of file LinearIsotropicHardening.hpp.

134  {
136 
137  // get sub-block (3x3) of local stiffens matrix, here represented by second
138  // order tensor
139  auto get_tensor2 = [](MatrixDouble &m, const int r, const int c) {
141  &m(r + 0, c + 0), &m(r + 0, c + 1), &m(r + 0, c + 2),
142  &m(r + 1, c + 0), &m(r + 1, c + 1), &m(r + 1, c + 2),
143  &m(r + 2, c + 0), &m(r + 2, c + 1), &m(r + 2, c + 2));
144  };
145 
150 
151  int nb_base_functions = row_data.getN().size2();
153 
155  // get element volume
156  double vol = getVolume();
157 
158  // get intergrayion weights
159  auto t_w = getFTensor0IntegrationWeight();
160 
161  // get derivatives of base functions on rows
162  auto t_row_diff_base = row_data.getFTensor1DiffN<3>();
163  // iterate over integration points
164  for (int gg = 0; gg != nbIntegrationPts; ++gg) {
165 
166  // calculate scalar weight times element volume
167  const double a = t_w * vol;
168 
169  // iterate over row base functions
170  int rr = 0;
171  for (; rr != nbRows / 3; ++rr) {
172 
173  // get sub matrix for the row
174  auto t_m = get_tensor2(K, 3 * rr, 0);
175 
176  // get derivatives of base functions for columns
177  auto t_col_diff_base = col_data.getFTensor1DiffN<3>(gg, 0);
178 
179  // iterate column base functions
180  for (int cc = 0; cc != nbCols / 3;++cc) {
181 
182  // integrate block local stiffens matrix
183  t_m(i, k) += a * (MT_1(i, j, k, l) *
184  (t_row_diff_base(j) * t_col_diff_base(l)));
185 
186  // move to next column base function
187  ++t_col_diff_base;
188 
189  // move to next block of local stiffens matrix
190  ++t_m;
191  }
192 
193  // move to next row base function
194  ++t_row_diff_base;
195  }
196  for (; rr != nb_base_functions; rr++) {
197  ++t_row_diff_base;
198  }
199  // move to next integration weight
200  ++t_w;
201  ++MT_1;
202  }
203 
205  }

◆ iNtegrate() [2/2]

MoFEMErrorCode OpK::iNtegrate ( DataForcesAndSourcesCore::EntData row_data,
DataForcesAndSourcesCore::EntData col_data 
)
protected

Integrate B^T D B operator.

Parameters
row_datarow data (consist base functions on row entity)
col_datacolumn data (consist base functions on column entity)
Returns
error code

Definition at line 148 of file simple_elasticity.cpp.

149  {
151 
152  // get sub-block (3x3) of local stiffens matrix, here represented by second
153  // order tensor
154  auto get_tensor2 = [](MatrixDouble &m, const int r, const int c) {
156  &m(r + 0, c + 0), &m(r + 0, c + 1), &m(r + 0, c + 2),
157  &m(r + 1, c + 0), &m(r + 1, c + 1), &m(r + 1, c + 2),
158  &m(r + 2, c + 0), &m(r + 2, c + 1), &m(r + 2, c + 2));
159  };
160 
165 
166  // get element volume
167  double vol = getVolume();
168 
169  // get intergration weights
170  auto t_w = getFTensor0IntegrationWeight();
171 
172  // get derivatives of base functions on rows
173  auto t_row_diff_base = row_data.getFTensor1DiffN<3>();
174  // iterate over integration points
175  for (int gg = 0; gg != nbIntegrationPts; ++gg) {
176 
177  // calculate scalar weight times element volume
178  const double a = t_w * vol;
179 
180  // iterate over row base functions
181  for (int rr = 0; rr != nbRows / 3; ++rr) {
182 
183  // get sub matrix for the row
184  auto t_m = get_tensor2(K, 3 * rr, 0);
185 
186  // get derivatives of base functions for columns
187  auto t_col_diff_base = col_data.getFTensor1DiffN<3>(gg, 0);
188 
189  // iterate column base functions
190  for (int cc = 0; cc != nbCols / 3;++cc) {
191 
192  // integrate block local stiffens matrix
193  t_m(i, k) +=
194  a * (tD(i, j, k, l) * (t_row_diff_base(j) * t_col_diff_base(l)));
195 
196  // move to next column base function
197  ++t_col_diff_base;
198 
199  // move to next block of local stiffens matrix
200  ++t_m;
201  }
202 
203  // move to next row base function
204  ++t_row_diff_base;
205  }
206 
207  // move to next integration weight
208  ++t_w;
209  }
210 
212  }

Member Data Documentation

◆ commonData

CommonData& OpK::commonData

Definition at line 64 of file LinearIsotropicHardening.hpp.

◆ isDiag

bool OpK::isDiag
protected

true if this block is on diagonal

Definition at line 129 of file LinearIsotropicHardening.hpp.

◆ K

MatrixDouble OpK::K

Definition at line 62 of file LinearIsotropicHardening.hpp.

◆ nbCols

int OpK::nbCols
protected

number if dof on column

Definition at line 127 of file LinearIsotropicHardening.hpp.

◆ nbIntegrationPts

int OpK::nbIntegrationPts
protected

number of integration points

Definition at line 128 of file LinearIsotropicHardening.hpp.

◆ nbRows

int OpK::nbRows
protected

number of dofs on rows

Definition at line 126 of file LinearIsotropicHardening.hpp.

◆ pOisson

double OpK::pOisson

Definition at line 40 of file simple_elasticity.cpp.

◆ tD

FTensor::Ddg<double, 3, 3> OpK::tD

Definition at line 35 of file simple_elasticity.cpp.

◆ yOung

double OpK::yOung

Definition at line 38 of file simple_elasticity.cpp.


The documentation for this struct was generated from the following files:
MoFEMFunctionReturnHot
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:516
FTensor::Tensor2
Definition: Tensor2_value.hpp:17
OpK::K
MatrixDouble K
Definition: LinearIsotropicHardening.hpp:62
MoFEM::ForcesAndSourcesCore::UserDataOperator::getFTensor0IntegrationWeight
auto getFTensor0IntegrationWeight()
Get integration weights.
Definition: ForcesAndSourcesCore.hpp:1020
OpK::tD
FTensor::Ddg< double, 3, 3 > tD
Definition: simple_elasticity.cpp:35
OpK::nbRows
int nbRows
number of dofs on rows
Definition: LinearIsotropicHardening.hpp:126
OpK::yOung
double yOung
Definition: simple_elasticity.cpp:38
OpK::nbCols
int nbCols
number if dof on column
Definition: LinearIsotropicHardening.hpp:127
MoFEM::DataForcesAndSourcesCore::EntData::getIndices
const VectorInt & getIndices() const
Get global indices of dofs on entity.
Definition: DataStructures.hpp:1178
MoFEM::Types::MatrixDouble
ublas::matrix< double, ublas::row_major, DoubleAllocator > MatrixDouble
Definition: Types.hpp:76
OpK::nbIntegrationPts
int nbIntegrationPts
number of integration points
Definition: LinearIsotropicHardening.hpp:128
m
static Index< 'm', 3 > m
Definition: BasicFeTools.hpp:77
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:604
TENSOR4_MAT_PTR
#define TENSOR4_MAT_PTR(MAT)
Definition: definitions.h:654
UserDataOperator
ForcesAndSourcesCore::UserDataOperator UserDataOperator
Definition: HookeElement.hpp:88
OpK::aSsemble
MoFEMErrorCode aSsemble(DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)
Definition: LinearIsotropicHardening.hpp:208
FTensor::Tensor4
Definition: Tensor4_value.hpp:19
OpK::iNtegrate
MoFEMErrorCode iNtegrate(DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)
Definition: LinearIsotropicHardening.hpp:133
Poisson2DTraditionalDDOperators::MatSetValues
constexpr auto MatSetValues
Definition: poisson_dd_H.hpp:20
MoFEM::KspMethod::ksp_B
Mat & ksp_B
Definition: LoopMethods.hpp:114
FTensor::Index< 'i', 3 >
OpK::pOisson
double pOisson
Definition: simple_elasticity.cpp:40
OpK::isDiag
bool isDiag
true if this block is on diagonal
Definition: LinearIsotropicHardening.hpp:129
MoFEM::ForcesAndSourcesCore::UserDataOperator::OPROWCOL
@ OPROWCOL
Definition: ForcesAndSourcesCore.hpp:101
MoFEM::ForcesAndSourcesCore::UserDataOperator::getFEMethod
const FEMethod * getFEMethod() const
Return raw pointer to Finite Element Method object.
Definition: ForcesAndSourcesCore.hpp:902
OpK::commonData
CommonData & commonData
Definition: LinearIsotropicHardening.hpp:64
ElectroPhysiology::B
const double B
Definition: ElecPhysOperators2D.hpp:27
CommonData::consistentMatrix
MatrixDouble consistentMatrix
Definition: PlDataAtGaussPts.hpp:44
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
MoFEM::VolumeElementForcesAndSourcesCoreBase::UserDataOperator::getVolume
double getVolume() const
element volume (linear geometry)
Definition: VolumeElementForcesAndSourcesCore.hpp:403
ElectroPhysiology::c
const double c
Definition: ElecPhysOperators.hpp:31
i
FTensor::Index< 'i', 3 > i
Definition: matrix_function.cpp:18
MoFEM::SnesMethod::snes_B
Mat & snes_B
preconditioner of jacobian matrix
Definition: LoopMethods.hpp:152
k
FTensor::Index< 'k', 3 > k
Definition: matrix_function.cpp:20
MoFEM::DataOperator::sYmm
bool sYmm
If true assume that matrix is symmetric structure.
Definition: DataOperators.hpp:72
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:485
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:415
l
FTensor::Index< 'l', 3 > l
Definition: matrix_function.cpp:21
ReactionDiffusionEquation::r
const double r
rate factor
Definition: reaction_diffusion.cpp:33
MoFEM::DataForcesAndSourcesCore::EntData::getFTensor1DiffN
FTensor::Tensor1< double *, Tensor_Dim > getFTensor1DiffN(const FieldApproximationBase base)
Get derivatives of base functions.
Definition: DataStructures.cpp:415
MoFEM::ForcesAndSourcesCore::UserDataOperator::getGaussPts
MatrixDouble & getGaussPts()
matrix of integration (Gauss) points for Volume Element
Definition: ForcesAndSourcesCore.hpp:1016