v0.14.0
Public Member Functions | List of all members
ConvectiveMassElement::OpVelocityLhs_dV_dX Struct Reference

#include <users_modules/basic_finite_elements/src/ConvectiveMassElement.hpp>

Inheritance diagram for ConvectiveMassElement::OpVelocityLhs_dV_dX:
[legend]
Collaboration diagram for ConvectiveMassElement::OpVelocityLhs_dV_dX:
[legend]

Public Member Functions

 OpVelocityLhs_dV_dX (const std::string vel_field, const std::string field_name, BlockData &data, CommonData &common_data)
 
virtual MoFEMErrorCode getJac (EntitiesFieldData::EntData &col_data, int gg)
 
- Public Member Functions inherited from ConvectiveMassElement::OpVelocityLhs_dV_dv
 OpVelocityLhs_dV_dv (const std::string vel_field, const std::string field_name, BlockData &data, CommonData &common_data)
 
- Public Member Functions inherited from ConvectiveMassElement::OpMassLhs_dM_dv
 OpMassLhs_dM_dv (const std::string vel_field, const std::string field_name, BlockData &data, CommonData &common_data, Range *forcesonlyonentities_ptr=NULL)
 
MoFEMErrorCode doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
 

Additional Inherited Members

- Public Attributes inherited from ConvectiveMassElement::OpMassLhs_dM_dv
BlockData & dAta
 
CommonData & commonData
 
Range forcesOnlyOnEntities
 
MatrixDouble k
 
MatrixDouble jac
 

Detailed Description

Definition at line 337 of file ConvectiveMassElement.hpp.

Constructor & Destructor Documentation

◆ OpVelocityLhs_dV_dX()

ConvectiveMassElement::OpVelocityLhs_dV_dX::OpVelocityLhs_dV_dX ( const std::string  vel_field,
const std::string  field_name,
BlockData &  data,
CommonData &  common_data 
)

Definition at line 1168 of file ConvectiveMassElement.cpp.

1171  : OpVelocityLhs_dV_dv(vel_field, field_name, data, common_data) {}

Member Function Documentation

◆ getJac()

MoFEMErrorCode ConvectiveMassElement::OpVelocityLhs_dV_dX::getJac ( EntitiesFieldData::EntData &  col_data,
int  gg 
)
virtual

Reimplemented from ConvectiveMassElement::OpVelocityLhs_dV_dv.

Definition at line 1173 of file ConvectiveMassElement.cpp.

1174  {
1175  MoFEMFunctionBeginHot;
1176  int nb_col = col_data.getIndices().size();
1177  jac.clear();
1178  if (!nb_col)
1179  MoFEMFunctionReturnHot(0);
1180  double *base_ptr = const_cast<double *>(&col_data.getN(gg)[0]);
1181  FTensor::Tensor0<double *> base(base_ptr, 1);
1182  double *diff_ptr =
1183  const_cast<double *>(&(col_data.getDiffN(gg, nb_col / 3)(0, 0)));
1184  FTensor::Tensor1<double *, 3> diff(diff_ptr, &diff_ptr[1], &diff_ptr[2], 3);
1185  FTensor::Tensor2<double *, 3, 3> t_jac(&jac(0, 0), &jac(0, 1), &jac(0, 2),
1186  &jac(1, 0), &jac(1, 1), &jac(1, 2),
1187  &jac(2, 0), &jac(2, 1), &jac(2, 2), 3);
1188  const int u = 3 + 3 + 9;
1189  FTensor::Tensor2<double *, 3, 3> t_mass1(
1190  &commonData.jacVel[gg](0, u + 0), &commonData.jacVel[gg](0, u + 1),
1191  &commonData.jacVel[gg](0, u + 2), &commonData.jacVel[gg](1, u + 0),
1192  &commonData.jacVel[gg](1, u + 1), &commonData.jacVel[gg](1, u + 2),
1193  &commonData.jacVel[gg](2, u + 0), &commonData.jacVel[gg](2, u + 1),
1194  &commonData.jacVel[gg](2, u + 2));
1195  const int s = 3 + 3 + 9 + 3;
1196  FTensor::Tensor3<double *, 3, 3, 3> t_mass3(
1197  // T* d000, T* d001, T* d002,
1198  // T* d010, T* d011, T* d012,
1199  // T* d020, T* d021, T* d022,
1200  // T* d100, T* d101, T* d102,
1201  // T* d110, T* d111, T* d112,
1202  // T* d120, T* d121, T* d122,
1203  // T* d200, T* d201, T* d202,
1204  // T* d210, T* d211, T* d212,
1205  // T* d220, T* d221, T* d222,
1206  &commonData.jacVel[gg](0, s + 0), &commonData.jacVel[gg](0, s + 1),
1207  &commonData.jacVel[gg](0, s + 2), &commonData.jacVel[gg](0, s + 3),
1208  &commonData.jacVel[gg](0, s + 4), &commonData.jacVel[gg](0, s + 5),
1209  &commonData.jacVel[gg](0, s + 6), &commonData.jacVel[gg](0, s + 7),
1210  &commonData.jacVel[gg](0, s + 8), &commonData.jacVel[gg](1, s + 0),
1211  &commonData.jacVel[gg](1, s + 1), &commonData.jacVel[gg](1, s + 2),
1212  &commonData.jacVel[gg](1, s + 3), &commonData.jacVel[gg](1, s + 4),
1213  &commonData.jacVel[gg](1, s + 5), &commonData.jacVel[gg](1, s + 6),
1214  &commonData.jacVel[gg](1, s + 7), &commonData.jacVel[gg](1, s + 8),
1215  &commonData.jacVel[gg](2, s + 0), &commonData.jacVel[gg](2, s + 1),
1216  &commonData.jacVel[gg](2, s + 2), &commonData.jacVel[gg](2, s + 3),
1217  &commonData.jacVel[gg](2, s + 4), &commonData.jacVel[gg](2, s + 5),
1218  &commonData.jacVel[gg](2, s + 6), &commonData.jacVel[gg](2, s + 7),
1219  &commonData.jacVel[gg](2, s + 8));
1220  FTensor::Index<'i', 3> i;
1221  FTensor::Index<'j', 3> j;
1222  FTensor::Index<'k', 3> k;
1223  for (int dd = 0; dd < nb_col / 3; dd++) {
1224  t_jac(i, j) += t_mass1(i, j) * base * getFEMethod()->ts_a;
1225  t_jac(i, j) += t_mass3(i, j, k) * diff(k);
1226  ++base_ptr;
1227  ++diff_ptr;
1228  ++t_jac;
1229  }
1230  MoFEMFunctionReturnHot(0);
1231 }
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:460
FTensor::Tensor1
Definition: Tensor1_value.hpp:8
ConvectiveMassElement::OpMassLhs_dM_dv::jac
MatrixDouble jac
Definition: ConvectiveMassElement.hpp:234
MoFEM::EntitiesFieldData::EntData::getDiffN
MatrixDouble & getDiffN(const FieldApproximationBase base)
get derivatives of base functions
Definition: EntitiesFieldData.hpp:1329
ConvectiveMassElement::OpMassLhs_dM_dv::k
MatrixDouble k
Definition: ConvectiveMassElement.hpp:234
FTensor::Tensor2< double *, 3, 3 >
FTensor::Tensor3
Definition: Tensor3_value.hpp:12
MoFEM::EntitiesFieldData::EntData::getIndices
const VectorInt & getIndices() const
Get global indices of dofs on entity.
Definition: EntitiesFieldData.hpp:1214
ConvectiveMassElement::CommonData::jacVel
std::vector< MatrixDouble > jacVel
Definition: ConvectiveMassElement.hpp:140
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
field_name
constexpr auto field_name
Definition: poisson_2d_homogeneous.cpp:13
FTensor::Index< 'i', 3 >
ConvectiveMassElement::OpMassLhs_dM_dv::commonData
CommonData & commonData
Definition: ConvectiveMassElement.hpp:227
FTensor::dd
const Tensor2_symmetric_Expr< const ddTensor0< T, Dim, i, j >, typename promote< T, double >::V, Dim, i, j > dd(const Tensor0< T * > &a, const Index< i, Dim > index1, const Index< j, Dim > index2, const Tensor1< int, Dim > &d_ijk, const Tensor1< double, Dim > &d_xyz)
Definition: ddTensor0.hpp:33
FTensor::Tensor0
Definition: Tensor0.hpp:16
MoFEM::EntitiesFieldData::EntData::getN
MatrixDouble & getN(const FieldApproximationBase base)
get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb....
Definition: EntitiesFieldData.hpp:1318
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
MoFEMFunctionBeginHot
#define MoFEMFunctionBeginHot
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:453
ConvectiveMassElement::OpVelocityLhs_dV_dv::OpVelocityLhs_dV_dv
OpVelocityLhs_dV_dv(const std::string vel_field, const std::string field_name, BlockData &data, CommonData &common_data)
Definition: ConvectiveMassElement.cpp:1059
Generated by Doxygen 1.8.17 and hosted at University of Glasgow