v0.14.0
Public Member Functions | Public Attributes | List of all members
MoFEM::OpSetCovariantPiolaTransformOnFace Struct Reference

transform Hcurl base fluxes from reference element to physical triangle More...

#include <src/finite_elements/DataOperators.hpp>

Inheritance diagram for MoFEM::OpSetCovariantPiolaTransformOnFace:
[legend]
Collaboration diagram for MoFEM::OpSetCovariantPiolaTransformOnFace:
[legend]

Public Member Functions

 OpSetCovariantPiolaTransformOnFace (const VectorDouble &normal, const MatrixDouble &normals_at_pts, const VectorDouble &tangent0, const MatrixDouble &tangent0_at_pts, const VectorDouble &tangent1, const MatrixDouble &tangent1_at_pts)
 
MoFEMErrorCode doWork (int side, EntityType type, EntitiesFieldData::EntData &data)
 Operator for linear form, usually to calculate values on right hand side. More...
 
- Public Member Functions inherited from MoFEM::DataOperator
 DataOperator (const bool symm=true)
 
virtual ~DataOperator ()=default
 
virtual MoFEMErrorCode doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
 Operator for bi-linear form, usually to calculate values on left hand side. More...
 
virtual MoFEMErrorCode opLhs (EntitiesFieldData &row_data, EntitiesFieldData &col_data)
 
virtual MoFEMErrorCode opRhs (EntitiesFieldData &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

const VectorDouble & nOrmal
 
const MatrixDouble & normalsAtGaussPts
 
const VectorDouble & tAngent0
 
const MatrixDouble & tangent0AtGaussPt
 
const VectorDouble & tAngent1
 
const MatrixDouble & tangent1AtGaussPt
 
- Public Attributes inherited from MoFEM::DataOperator
DoWorkLhsHookFunType doWorkLhsHook
 
DoWorkRhsHookFunType doWorkRhsHook
 
bool sYmm
 If true assume that matrix is symmetric structure. More...
 
std::array< bool, MBMAXTYPE > doEntities
 If true operator is executed for entity. More...
 
bool & doVertices
 \deprectaed If false skip vertices More...
 
bool & doEdges
 \deprectaed If false skip edges More...
 
bool & doQuads
 \deprectaed More...
 
bool & doTris
 \deprectaed More...
 
bool & doTets
 \deprectaed More...
 
bool & doPrisms
 \deprectaed More...
 

Additional Inherited Members

- Public Types inherited from MoFEM::DataOperator
using DoWorkLhsHookFunType = boost::function< MoFEMErrorCode(DataOperator *op_ptr, int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)>
 
using DoWorkRhsHookFunType = boost::function< MoFEMErrorCode(DataOperator *op_ptr, int side, EntityType type, EntitiesFieldData::EntData &data)>
 

Detailed Description

transform Hcurl base fluxes from reference element to physical triangle

Deprecated:
It is used in contact elements. Contact elements should be minified to work as face element,

Definition at line 457 of file DataOperators.hpp.

Constructor & Destructor Documentation

◆ OpSetCovariantPiolaTransformOnFace()

MoFEM::OpSetCovariantPiolaTransformOnFace::OpSetCovariantPiolaTransformOnFace ( const VectorDouble &  normal,
const MatrixDouble &  normals_at_pts,
const VectorDouble &  tangent0,
const MatrixDouble &  tangent0_at_pts,
const VectorDouble &  tangent1,
const MatrixDouble &  tangent1_at_pts 
)
inline

Definition at line 466 of file DataOperators.hpp.

472  : nOrmal(normal), normalsAtGaussPts(normals_at_pts), tAngent0(tangent0),
473  tangent0AtGaussPt(tangent0_at_pts), tAngent1(tangent1),
474  tangent1AtGaussPt(tangent1_at_pts) {}

Member Function Documentation

◆ doWork()

MoFEMErrorCode MoFEM::OpSetCovariantPiolaTransformOnFace::doWork ( int  side,
EntityType  type,
EntitiesFieldData::EntData &  data 
)
virtual

Operator for linear form, usually to calculate values on right hand side.

Reimplemented from MoFEM::DataOperator.

Definition at line 585 of file DataOperators.cpp.

586  {
587  MoFEMFunctionBegin;
588 
589  const auto type_dim = moab::CN::Dimension(type);
590  if (type_dim != 1 && type_dim != 2)
591  MoFEMFunctionReturnHot(0);
592 
593  FTensor::Index<'i', 3> i;
594  FTensor::Index<'j', 3> j;
595  FTensor::Index<'k', 2> k;
596 
597  FTensor::Tensor2<FTensor::PackPtr<const double *, 3>, 3, 3> t_m(
598  &tAngent0[0], &tAngent1[0], &nOrmal[0],
599 
600  &tAngent0[1], &tAngent1[1], &nOrmal[1],
601 
602  &tAngent0[2], &tAngent1[2], &nOrmal[2]);
603  double det;
604  FTensor::Tensor2<double, 3, 3> t_inv_m;
605  CHKERR determinantTensor3by3(t_m, det);
606  CHKERR invertTensor3by3(t_m, det, t_inv_m);
607 
608  for (int b = AINSWORTH_LEGENDRE_BASE; b != LASTBASE; ++b) {
609 
610  FieldApproximationBase base = static_cast<FieldApproximationBase>(b);
611 
612  auto &baseN = data.getN(base);
613  auto &diffBaseN = data.getDiffN(base);
614 
615  int nb_dofs = baseN.size2() / 3;
616  int nb_gauss_pts = baseN.size1();
617 
618  MatrixDouble piola_n(baseN.size1(), baseN.size2());
619  MatrixDouble diff_piola_n(diffBaseN.size1(), diffBaseN.size2());
620 
621  if (nb_dofs > 0 && nb_gauss_pts > 0) {
622 
623  FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3> t_h_curl(
624  &baseN(0, HVEC0), &baseN(0, HVEC1), &baseN(0, HVEC2));
625  FTensor::Tensor2<FTensor::PackPtr<double *, 6>, 3, 2> t_diff_h_curl(
626  &diffBaseN(0, HVEC0_0), &diffBaseN(0, HVEC0_1),
627  &diffBaseN(0, HVEC1_0), &diffBaseN(0, HVEC1_1),
628  &diffBaseN(0, HVEC2_0), &diffBaseN(0, HVEC2_1));
629  FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3> t_transformed_h_curl(
630  &piola_n(0, HVEC0), &piola_n(0, HVEC1), &piola_n(0, HVEC2));
631  FTensor::Tensor2<FTensor::PackPtr<double *, 6>, 3, 2>
632  t_transformed_diff_h_curl(
633  &diff_piola_n(0, HVEC0_0), &diff_piola_n(0, HVEC0_1),
634  &diff_piola_n(0, HVEC1_0), &diff_piola_n(0, HVEC1_1),
635  &diff_piola_n(0, HVEC2_0), &diff_piola_n(0, HVEC2_1));
636 
637  int cc = 0;
638  if (normalsAtGaussPts.size1() == (unsigned int)nb_gauss_pts) {
639  // HO geometry is set, so jacobian is different at each gauss point
640  FTensor::Tensor2<FTensor::PackPtr<const double *, 3>, 3, 3> t_m_at_pts(
641  &tangent0AtGaussPt(0, 0), &tangent1AtGaussPt(0, 0),
642  &normalsAtGaussPts(0, 0), &tangent0AtGaussPt(0, 1),
643  &tangent1AtGaussPt(0, 1), &normalsAtGaussPts(0, 1),
644  &tangent0AtGaussPt(0, 2), &tangent1AtGaussPt(0, 2),
645  &normalsAtGaussPts(0, 2));
646  for (int gg = 0; gg < nb_gauss_pts; ++gg) {
647  CHKERR determinantTensor3by3(t_m_at_pts, det);
648  CHKERR invertTensor3by3(t_m_at_pts, det, t_inv_m);
649  for (int ll = 0; ll != nb_dofs; ll++) {
650  t_transformed_h_curl(i) = t_inv_m(j, i) * t_h_curl(j);
651  t_transformed_diff_h_curl(i, k) =
652  t_inv_m(j, i) * t_diff_h_curl(j, k);
653  ++t_h_curl;
654  ++t_transformed_h_curl;
655  ++t_diff_h_curl;
656  ++t_transformed_diff_h_curl;
657  ++cc;
658  }
659  ++t_m_at_pts;
660  }
661  } else {
662  for (int gg = 0; gg < nb_gauss_pts; ++gg) {
663  for (int ll = 0; ll != nb_dofs; ll++) {
664  t_transformed_h_curl(i) = t_inv_m(j, i) * t_h_curl(j);
665  t_transformed_diff_h_curl(i, k) =
666  t_inv_m(j, i) * t_diff_h_curl(j, k);
667  ++t_h_curl;
668  ++t_transformed_h_curl;
669  ++t_diff_h_curl;
670  ++t_transformed_diff_h_curl;
671  ++cc;
672  }
673  }
674  }
675  if (cc != nb_gauss_pts * nb_dofs)
676  SETERRQ(PETSC_COMM_SELF, MOFEM_IMPOSSIBLE_CASE, "Data inconsistency");
677 
678  baseN.swap(piola_n);
679  diffBaseN.swap(diff_piola_n);
680  }
681  }
682 
683  MoFEMFunctionReturn(0);
684 }

Member Data Documentation

◆ nOrmal

const VectorDouble& MoFEM::OpSetCovariantPiolaTransformOnFace::nOrmal

Definition at line 459 of file DataOperators.hpp.

◆ normalsAtGaussPts

const MatrixDouble& MoFEM::OpSetCovariantPiolaTransformOnFace::normalsAtGaussPts

Definition at line 460 of file DataOperators.hpp.

◆ tAngent0

const VectorDouble& MoFEM::OpSetCovariantPiolaTransformOnFace::tAngent0

Definition at line 461 of file DataOperators.hpp.

◆ tangent0AtGaussPt

const MatrixDouble& MoFEM::OpSetCovariantPiolaTransformOnFace::tangent0AtGaussPt

Definition at line 462 of file DataOperators.hpp.

◆ tAngent1

const VectorDouble& MoFEM::OpSetCovariantPiolaTransformOnFace::tAngent1

Definition at line 463 of file DataOperators.hpp.

◆ tangent1AtGaussPt

const MatrixDouble& MoFEM::OpSetCovariantPiolaTransformOnFace::tangent1AtGaussPt

Definition at line 464 of file DataOperators.hpp.

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:447
LASTBASE
@ LASTBASE
Definition: definitions.h:69
FTensor::Tensor1
Definition: Tensor1_value.hpp:8
MoFEM::OpSetCovariantPiolaTransformOnFace::normalsAtGaussPts
const MatrixDouble & normalsAtGaussPts
Definition: DataOperators.hpp:460
MoFEM::Types::MatrixDouble
UBlasMatrix< double > MatrixDouble
Definition: Types.hpp:77
HVEC0_1
@ HVEC0_1
Definition: definitions.h:195
MoFEM::OpSetCovariantPiolaTransformOnFace::nOrmal
const VectorDouble & nOrmal
Definition: DataOperators.hpp:459
HVEC0_0
@ HVEC0_0
Definition: definitions.h:192
HVEC1_1
@ HVEC1_1
Definition: definitions.h:196
MOFEM_IMPOSSIBLE_CASE
@ MOFEM_IMPOSSIBLE_CASE
Definition: definitions.h:35
HVEC1
@ HVEC1
Definition: definitions.h:186
FTensor::Tensor2
Definition: Tensor2_value.hpp:16
MoFEM::invertTensor3by3
MoFEMErrorCode invertTensor3by3(ublas::matrix< T, L, A > &jac_data, ublas::vector< T, A > &det_data, ublas::matrix< T, L, A > &inv_jac_data)
Calculate inverse of tensor rank 2 at integration points.
Definition: Templates.hpp:1559
HVEC2_1
@ HVEC2_1
Definition: definitions.h:197
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:535
convert.type
type
Definition: convert.py:64
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
MoFEM::OpSetCovariantPiolaTransformOnFace::tangent1AtGaussPt
const MatrixDouble & tangent1AtGaussPt
Definition: DataOperators.hpp:464
FTensor::Index< 'i', 3 >
MoFEM::determinantTensor3by3
static auto determinantTensor3by3(T &t)
Calculate the determinant of a 3x3 matrix or a tensor of rank 2.
Definition: Templates.hpp:1511
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
AINSWORTH_LEGENDRE_BASE
@ AINSWORTH_LEGENDRE_BASE
Ainsworth Cole (Legendre) approx. base .
Definition: definitions.h:60
FieldApproximationBase
FieldApproximationBase
approximation base
Definition: definitions.h:58
HVEC2_0
@ HVEC2_0
Definition: definitions.h:194
MoFEM::OpSetCovariantPiolaTransformOnFace::tangent0AtGaussPt
const MatrixDouble & tangent0AtGaussPt
Definition: DataOperators.hpp:462
MoFEM::OpSetCovariantPiolaTransformOnFace::tAngent0
const VectorDouble & tAngent0
Definition: DataOperators.hpp:461
MoFEM::OpSetCovariantPiolaTransformOnFace::tAngent1
const VectorDouble & tAngent1
Definition: DataOperators.hpp:463
HVEC2
@ HVEC2
Definition: definitions.h:186
k
FTensor::Index< 'k', 3 > k
Definition: matrix_function.cpp:20
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:416
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:346
HVEC1_0
@ HVEC1_0
Definition: definitions.h:193
HVEC0
@ HVEC0
Definition: definitions.h:186
Generated by Doxygen 1.8.17 and hosted at University of Glasgow