v0.15.0
Loading...
Searching...
No Matches
ConvectiveMassElement::OpVelocityRhs Struct Reference

#include "users_modules/basic_finite_elements/src/ConvectiveMassElement.hpp"

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

Public Member Functions

 OpVelocityRhs (const std::string field_name, BlockData &data, CommonData &common_data)
 
MoFEMErrorCode doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &row_data)
 

Public Attributes

BlockDatadAta
 
CommonDatacommonData
 
VectorDouble nf
 

Detailed Description

Definition at line 301 of file ConvectiveMassElement.hpp.

Constructor & Destructor Documentation

◆ OpVelocityRhs()

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

Definition at line 1008 of file ConvectiveMassElement.cpp.

1010 : VolumeElementForcesAndSourcesCore::UserDataOperator(
1011 field_name, ForcesAndSourcesCore::UserDataOperator::OPROW),
1012 dAta(data), commonData(common_data) {}
constexpr auto field_name

Member Function Documentation

◆ doWork()

MoFEMErrorCode ConvectiveMassElement::OpVelocityRhs::doWork ( int row_side,
EntityType row_type,
EntitiesFieldData::EntData & row_data )

Definition at line 1014 of file ConvectiveMassElement.cpp.

1015 {
1017
1018 if (dAta.tEts.find(getNumeredEntFiniteElementPtr()->getEnt()) ==
1019 dAta.tEts.end()) {
1021 }
1022 int nb_dofs = row_data.getIndices().size();
1023 if (nb_dofs == 0)
1025
1026 auto base = row_data.getFTensor0N();
1027 int nb_base_functions = row_data.getN().size2();
1028 FTensor::Index<'i', 3> i;
1029
1030 {
1031
1032 nf.resize(nb_dofs);
1033 nf.clear();
1034
1035 for (unsigned int gg = 0; gg < row_data.getN().size1(); gg++) {
1036 FTensor::Tensor1<double *, 3> t_nf(&nf[0], &nf[1], &nf[2], 3);
1038 &commonData.valVel[gg][1],
1039 &commonData.valVel[gg][2]);
1040 int dd = 0;
1041 for (; dd < nb_dofs / 3; dd++) {
1042 t_nf(i) += base * res(i);
1043 ++base;
1044 ++t_nf;
1045 }
1046 for (; dd != nb_base_functions; dd++) {
1047 ++base;
1048 }
1049 }
1050
1051 if (row_data.getIndices().size() > 3 * row_data.getN().size2()) {
1052 SETERRQ(PETSC_COMM_SELF, 1, "data inconsistency");
1053 }
1054 CHKERR VecSetValues(getFEMethod()->ts_F, row_data.getIndices().size(),
1055 &row_data.getIndices()[0], &nf[0], ADD_VALUES);
1056 }
1058}
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define CHKERR
Inline error check.
#define MoFEMFunctionBeginHot
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
FTensor::Index< 'i', SPACE_DIM > i
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
MoFEMErrorCode VecSetValues(Vec V, const EntitiesFieldData::EntData &data, const double *ptr, InsertMode iora)
Assemble PETSc vector.
FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > getFTensor0N(const FieldApproximationBase base)
Get base function as Tensor0.
MatrixDouble & getN(const FieldApproximationBase base)
get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb....
const VectorInt & getIndices() const
Get global indices of dofs on entity.

Member Data Documentation

◆ commonData

CommonData& ConvectiveMassElement::OpVelocityRhs::commonData

Definition at line 306 of file ConvectiveMassElement.hpp.

◆ dAta

BlockData& ConvectiveMassElement::OpVelocityRhs::dAta

Definition at line 305 of file ConvectiveMassElement.hpp.

◆ nf

VectorDouble ConvectiveMassElement::OpVelocityRhs::nf

Definition at line 311 of file ConvectiveMassElement.hpp.


The documentation for this struct was generated from the following files: