v0.15.0
Loading...
Searching...
No Matches
MoFEM::OpBaseTimesScalarImpl< 1, S, GAUSS, OpBase > Struct Template Reference

#include "src/finite_elements/LinearFormsIntegratorsImpl.hpp"

Inheritance diagram for MoFEM::OpBaseTimesScalarImpl< 1, S, GAUSS, OpBase >:
[legend]
Collaboration diagram for MoFEM::OpBaseTimesScalarImpl< 1, S, GAUSS, OpBase >:
[legend]

Public Member Functions

 OpBaseTimesScalarImpl (const std::string field_name, boost::shared_ptr< VectorDouble > vec, ScalarFun beta_coeff=[](double, double, double) constexpr { return 1;}, boost::shared_ptr< Range > ents_ptr=nullptr)
 
- Public Member Functions inherited from MoFEM::OpBaseImpl< A, EleOp >
 OpBaseImpl (const std::string row_field_name, const std::string col_field_name, const OpType type, boost::shared_ptr< Range > ents_ptr=nullptr)
 
MoFEMErrorCode doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data)
 Do calculations for the left hand side.
 
MoFEMErrorCode doWork (int row_side, EntityType row_type, EntData &row_data)
 Do calculations for the right hand side.
 

Protected Member Functions

MoFEMErrorCode iNtegrate (EntitiesFieldData::EntData &data)
 Class dedicated to integrate operator.
 
- Protected Member Functions inherited from MoFEM::OpBaseImpl< A, EleOp >
template<int DIM>
FTensor::Tensor1< FTensor::PackPtr< double *, DIM >, DIM > getNf ()
 
template<int DIM>
FTensor::Tensor2< FTensor::PackPtr< double *, DIM >, DIM, DIM > getLocMat (const int rr)
 
virtual MoFEMErrorCode iNtegrate (EntData &row_data, EntData &col_data)
 Integrate grad-grad operator.
 
virtual MoFEMErrorCode aSsemble (EntData &row_data, EntData &col_data, const bool trans)
 
virtual MoFEMErrorCode aSsemble (EntData &data)
 
virtual size_t getNbOfBaseFunctions (EntitiesFieldData::EntData &data)
 Get number of base functions.
 

Protected Attributes

ScalarFun betaCoeff
 
boost::shared_ptr< VectorDoublesourceVec
 
- Protected Attributes inherited from MoFEM::OpBaseImpl< A, EleOp >
int nbRows
 number of dofs on rows
 
int nbCols
 number if dof on column
 
int nbIntegrationPts
 number of integration points
 
int nbRowBaseFunctions
 number or row base functions
 
int rowSide
 row side number
 
int colSide
 column side number
 
EntityType rowType
 row type
 
EntityType colType
 column type
 
bool assembleTranspose
 
bool onlyTranspose
 
MatrixDouble locMat
 local entity block matrix
 
MatrixDouble locMatTranspose
 local entity block matrix
 
VectorDouble locF
 local entity vector
 

Additional Inherited Members

- Public Types inherited from MoFEM::OpBaseImpl< A, EleOp >
using OpType = typename EleOp::OpType
 
using EntData = EntitiesFieldData::EntData
 
using MatSetValuesHook
 
- Public Attributes inherited from MoFEM::OpBaseImpl< A, EleOp >
TimeFun timeScalingFun
 assumes that time variable is set
 
FEFun feScalingFun
 assumes that time variable is set
 
boost::shared_ptr< RangeentsPtr
 Entities on which element is run.
 
- Static Public Attributes inherited from MoFEM::OpBaseImpl< A, EleOp >
static MatSetValuesHook matSetValuesHook
 

Detailed Description

template<int S, typename OpBase>
struct MoFEM::OpBaseTimesScalarImpl< 1, S, GAUSS, OpBase >

Definition at line 123 of file LinearFormsIntegratorsImpl.hpp.

Constructor & Destructor Documentation

◆ OpBaseTimesScalarImpl()

template<int S, typename OpBase >
MoFEM::OpBaseTimesScalarImpl< 1, S, GAUSS, OpBase >::OpBaseTimesScalarImpl ( const std::string field_name,
boost::shared_ptr< VectorDouble > vec,
ScalarFun beta_coeff = [](doubledoubledouble) constexpr { return 1; },
boost::shared_ptr< Range > ents_ptr = nullptr )
inline

Definition at line 125 of file LinearFormsIntegratorsImpl.hpp.

127 { return 1; },
128 boost::shared_ptr<Range> ents_ptr = nullptr)
129 : OpBase(field_name, field_name, OpBase::OPROW, ents_ptr), sourceVec(vec),
130 betaCoeff(beta_coeff) {}
constexpr auto field_name

Member Function Documentation

◆ iNtegrate()

template<int S, typename OpBase >
MoFEMErrorCode MoFEM::OpBaseTimesScalarImpl< 1, S, GAUSS, OpBase >::iNtegrate ( EntitiesFieldData::EntData & data)
protectedvirtual

Class dedicated to integrate operator.

Parameters
dataentity data on element row
Returns
error code

Reimplemented from MoFEM::OpBaseImpl< A, EleOp >.

Definition at line 597 of file LinearFormsIntegratorsImpl.hpp.

598 {
600
601 // get element volume
602 const double vol = OpBase::getMeasure();
603 // get integration weights
604 auto t_w = OpBase::getFTensor0IntegrationWeight();
605 // get base function gradient on rows
606 auto t_row_base = row_data.getFTensor0N();
607 // get vector values
608 auto t_vec = getFTensor0FromVec<S>(*sourceVec);
609
610#ifndef NDEBUG
611 if (sourceVec->size() != OpBase::nbIntegrationPts) {
612 SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
613 "Wrong number of integration points %d != %ld",
615 }
616#endif
617
618 // get coords
619 auto t_coords = OpBase::getFTensor1CoordsAtGaussPts();
620 // loop over integration points
621 for (int gg = 0; gg != OpBase::nbIntegrationPts; gg++) {
622 // take into account Jacobian
623 const double alpha =
624 t_w * vol * betaCoeff(t_coords(0), t_coords(1), t_coords(2));
625 // loop over rows base functions
626 int rr = 0;
627 for (; rr != OpBase::nbRows; ++rr) {
628 OpBase::locF[rr] += alpha * t_row_base * t_vec;
629 ++t_row_base;
630 }
631 for (; rr < OpBase::nbRowBaseFunctions; ++rr)
632 ++t_row_base;
633 ++t_w; // move to another integration weight
634 ++t_vec;
635 ++t_coords;
636 }
638}
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
@ MOFEM_DATA_INCONSISTENCY
Definition definitions.h:31
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
static auto getFTensor0FromVec(ublas::vector< T, A > &data)
Get tensor rank 0 (scalar) form data vector.
VectorDouble locF
local entity vector
int nbRows
number of dofs on rows
int nbIntegrationPts
number of integration points
int nbRowBaseFunctions
number or row base functions

Member Data Documentation

◆ betaCoeff

template<int S, typename OpBase >
ScalarFun MoFEM::OpBaseTimesScalarImpl< 1, S, GAUSS, OpBase >::betaCoeff
protected

Definition at line 133 of file LinearFormsIntegratorsImpl.hpp.

◆ sourceVec

template<int S, typename OpBase >
boost::shared_ptr<VectorDouble> MoFEM::OpBaseTimesScalarImpl< 1, S, GAUSS, OpBase >::sourceVec
protected

Definition at line 134 of file LinearFormsIntegratorsImpl.hpp.


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