StdRDOperators::OpAssembleStiffLhs< DIM > Struct Template Reference

#include <users_modules/softmech/chemo_mech/src/rd_stdOperators.hpp>

Inheritance diagram for StdRDOperators::OpAssembleStiffLhs< DIM >:

Collaboration diagram for StdRDOperators::OpAssembleStiffLhs< DIM >:

Public Member Functions
	OpAssembleStiffLhs (std::string fieldu, boost::shared_ptr< PreviousData > &data, std::map< int, BlockData > &block_map)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data)

Private Attributes
boost::shared_ptr< PreviousData >	commonData
std::map< int, BlockData >	setOfBlock
MatrixDouble	mat
MatrixDouble	transMat

Detailed Description

template<int DIM>
struct StdRDOperators::OpAssembleStiffLhs< DIM >

Definition at line 395 of file rd_stdOperators.hpp.

Constructor & Destructor Documentation

OpAssembleStiffLhs()

template<int DIM>

StdRDOperators::OpAssembleStiffLhs< DIM >::OpAssembleStiffLhs	(	std::string	fieldu,
		boost::shared_ptr< PreviousData > &	data,
		std::map< int, BlockData > &	block_map
	)

Definition at line 397 of file rd_stdOperators.hpp.

      : OpEle(fieldu, fieldu, OpEle::OPROWCOL), commonData(data),
        setOfBlock(block_map) {
    sYmm = true;
  }

Member Function Documentation

doWork()

template<int DIM>

MoFEMErrorCode StdRDOperators::OpAssembleStiffLhs< DIM >::doWork	(	int	row_side,
		int	col_side,
		EntityType	row_type,
		EntityType	col_type,
		EntData &	row_data,
		EntData &	col_data
	)

Definition at line 403 of file rd_stdOperators.hpp.

                                           {
    MoFEMFunctionBegin;
 
    const int nb_row_dofs = row_data.getIndices().size();
    const int nb_col_dofs = col_data.getIndices().size();
    // cerr << "In doWork() : (row, col) = (" << nb_row_dofs << ", " <<
    // nb_col_dofs << ")" << endl;
    if (nb_row_dofs && nb_col_dofs) {
      auto find_block_data = [&]() {
        EntityHandle fe_ent = getFEEntityHandle();
        BlockData *block_raw_ptr = nullptr;
        for (auto &m : setOfBlock) {
          if (m.second.block_ents.find(fe_ent) != m.second.block_ents.end()) {
            block_raw_ptr = &m.second;
            break;
          }
        }
        return block_raw_ptr;
      };
 
      auto block_data_ptr = find_block_data();
      if (!block_data_ptr)
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "Block not found");
 
      auto &block_data = *block_data_ptr;
 
      mat.resize(nb_row_dofs, nb_col_dofs, false);
      mat.clear();
      const int nb_integration_pts = getGaussPts().size2();
      auto t_row_base = row_data.getFTensor0N();
 
      auto t_val = getFTensor0FromVec(commonData->values);
      auto t_grad = getFTensor1FromMat<DIM>(commonData->grads);
 
      auto t_row_diff_base = row_data.getFTensor1DiffN<DIM>();
      auto t_w = getFTensor0IntegrationWeight();
      const double ts_a = getFEMethod()->ts_a;
      const double vol = getMeasure();
 
      FTensor::Index<'i', DIM> i;
 
      // cout << "B0 : " << block_data.B0 << endl;
 
      for (int gg = 0; gg != nb_integration_pts; ++gg) {
        const double a = vol * t_w;
 
        for (int rr = 0; rr != nb_row_dofs; ++rr) {
          auto t_col_base = col_data.getFTensor0N(gg, 0);
          auto t_col_diff_base = col_data.getFTensor1DiffN<DIM>(gg, 0);
          for (int cc = 0; cc != nb_col_dofs; ++cc) {
 
            mat(rr, cc) +=
                a * (t_row_base * t_col_base * ts_a +
                     (block_data.B0 + B * t_val) * t_row_diff_base(i) *
                         t_col_diff_base(i) +
                     B * t_col_base * t_grad(i) * t_row_diff_base(i));
 
            ++t_col_base;
            ++t_col_diff_base;
          }
          ++t_row_base;
          ++t_row_diff_base;
        }
        ++t_w;
        ++t_val;
        ++t_grad;
      }
      CHKERR MatSetValues(getFEMethod()->ts_B, row_data, col_data, &mat(0, 0),
                          ADD_VALUES);
      if (row_side != col_side || row_type != col_type) {
        transMat.resize(nb_col_dofs, nb_row_dofs, false);
        noalias(transMat) = trans(mat);
        CHKERR MatSetValues(getFEMethod()->ts_B, col_data, row_data,
                            &transMat(0, 0), ADD_VALUES);
      }
    }
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

commonData

template<int DIM>

boost::shared_ptr<PreviousData> StdRDOperators::OpAssembleStiffLhs< DIM >::commonData

private

Definition at line 485 of file rd_stdOperators.hpp.

mat

template<int DIM>

MatrixDouble StdRDOperators::OpAssembleStiffLhs< DIM >::mat

private

Definition at line 487 of file rd_stdOperators.hpp.

setOfBlock

template<int DIM>

std::map<int, BlockData> StdRDOperators::OpAssembleStiffLhs< DIM >::setOfBlock

private

Definition at line 486 of file rd_stdOperators.hpp.

transMat

template<int DIM>

MatrixDouble StdRDOperators::OpAssembleStiffLhs< DIM >::transMat

private

Definition at line 487 of file rd_stdOperators.hpp.

---

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

• users_modules/softmech/chemo_mech/src/rd_stdOperators.hpp