#include <src/boundary_conditions/EssentialDisplacementCubitBcData.hpp>

Inheritance diagram for MoFEM::OpEssentialRhsImpl< DisplacementCubitBcData, 1, FIELD_DIM, A, I, OpBase >:

Collaboration diagram for MoFEM::OpEssentialRhsImpl< DisplacementCubitBcData, 1, FIELD_DIM, A, I, OpBase >:

Public Types
using	OpSource = typename FormsIntegrators< OpBase >::template Assembly< A >::template LinearForm< I >::template OpSource< 1, FIELD_DIM >

Public Member Functions
	OpEssentialRhsImpl (const std::string field_name, boost::shared_ptr< DisplacementCubitBcData > bc_data, boost::shared_ptr< Range > ents_ptr, std::vector< boost::shared_ptr< ScalingMethod > > smv)

Private Member Functions
FTensor::Tensor1< double, FIELD_DIM >	rotFunction (double x, double y, double z)

Private Attributes
FTensor::Tensor1< double, FIELD_DIM >	tVal

FTensor::Tensor1< double, 3 >	tAngles

FTensor::Tensor1< double, 3 >	tOffset

VecOfTimeScalingMethods	vecOfTimeScalingMethods

VectorFun< FIELD_DIM >	dispFunction

Detailed Description

template<int FIELD_DIM, AssemblyType A, IntegrationType I, typename OpBase>
struct MoFEM::OpEssentialRhsImpl< DisplacementCubitBcData, 1, FIELD_DIM, A, I, OpBase >

Definition at line 167 of file EssentialDisplacementCubitBcData.hpp.

Member Typedef Documentation

◆ OpSource

template<int FIELD_DIM, AssemblyType A, IntegrationType I, typename OpBase >

using MoFEM::OpEssentialRhsImpl< DisplacementCubitBcData, 1, FIELD_DIM, A, I, OpBase >::OpSource = typename FormsIntegrators<OpBase>::template Assembly< A>::template LinearForm<I>::template OpSource<1, FIELD_DIM>

Definition at line 171 of file EssentialDisplacementCubitBcData.hpp.

Constructor & Destructor Documentation

◆ OpEssentialRhsImpl()

template<int FIELD_DIM, AssemblyType A, IntegrationType I, typename OpBase >

MoFEM::OpEssentialRhsImpl< DisplacementCubitBcData, 1, FIELD_DIM, A, I, OpBase >::OpEssentialRhsImpl	(	const std::string	field_name,
		boost::shared_ptr< DisplacementCubitBcData >	bc_data,
		boost::shared_ptr< Range >	ents_ptr,
		std::vector< boost::shared_ptr< ScalingMethod > >	smv
	)

Definition at line 198 of file EssentialDisplacementCubitBcData.hpp.

    : OpSource(
          field_name, [this](double, double, double) { return tVal; },
          ents_ptr),
      vecOfTimeScalingMethods(smv) {
  static_assert(FIELD_DIM > 1, "Is not implemented for scalar field");
 
  FTensor::Index<'i', FIELD_DIM> i;
  tVal(i) = 0;
  tAngles(i) = 0;
  tOffset(i) = 0;
 
  if (bc_data->data.flag1 == 1)
    tVal(0) = -bc_data->data.value1;
  if (bc_data->data.flag2 == 1 && FIELD_DIM > 1)
    tVal(1) = -bc_data->data.value2;
  if (bc_data->data.flag3 == 1 && FIELD_DIM > 2)
    tVal(2) = -bc_data->data.value3;
  if (bc_data->data.flag4 == 1 && FIELD_DIM > 2)
    tAngles(0) = -bc_data->data.value4;
  if (bc_data->data.flag5 == 1 && FIELD_DIM > 2)
    tAngles(1) = -bc_data->data.value5;
  if (bc_data->data.flag6 == 1 && FIELD_DIM > 1)
    tAngles(2) = -bc_data->data.value6;
 
  if (auto ext_bc_data =
          dynamic_cast<DisplacementCubitBcDataWithRotation const *>(
              bc_data.get())) {
    for (int a = 0; a != 3; ++a)
      tOffset(a) = ext_bc_data->rotOffset[a];
  }
 
  this->timeScalingFun = [this](const double t) {
    double s = 1;
    for (auto &o : vecOfTimeScalingMethods) {
      s *= o->getScale(t);
    }
    return s;
  };
  if (bc_data->data.flag4 || bc_data->data.flag5 || bc_data->data.flag6) {
    this->dispFunction = [this](double x, double y, double z) {
      return this->rotFunction(x, y, z);
    };
  } else {
    // use default set at construction
  }
}

Member Function Documentation

◆ rotFunction()

template<int FIELD_DIM, AssemblyType A, IntegrationType I, typename OpBase >

FTensor::Tensor1< double, FIELD_DIM > MoFEM::OpEssentialRhsImpl< DisplacementCubitBcData, 1, FIELD_DIM, A, I, OpBase >::rotFunction	(	double	x,
		double	y,
		double	z
	)

inlineprivate

Definition at line 186 of file EssentialDisplacementCubitBcData.hpp.

                                                            {
    FTensor::Index<'i', FIELD_DIM> i;
    auto rot = DisplacementCubitBcDataWithRotation::GetRotDisp(
        tAngles, FTensor::Tensor1<double, 3>{x, y, z}, tOffset);
    FTensor::Tensor1<double, FIELD_DIM> t_ret;
    t_ret(i) = tVal(i) + rot(i);
    return t_ret;
  };