#include <users_modules/softmech/poisson_problem/src/PoissonOperators.hpp>

Inheritance diagram for PoissonOps::OpLhsUU:

Collaboration diagram for PoissonOps::OpLhsUU:

Public Member Functions
	OpLhsUU (std::string domain_field1, std::string domain_field2)

MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data)

Private Attributes
MatrixDouble	locLhs

MatrixDouble	transLocLhs

Detailed Description

Definition at line 23 of file PoissonOperators.hpp.

Constructor & Destructor Documentation

◆ OpLhsUU()

PoissonOps::OpLhsUU::OpLhsUU	(	std::string	domain_field1,
		std::string	domain_field2
	)

Definition at line 25 of file PoissonOperators.hpp.

      : OpVolEle(domain_field1, domain_field2, OpVolEle::OPROWCOL) {
    sYmm = true;
  }

Member Function Documentation

◆ doWork()

MoFEMErrorCode PoissonOps::OpLhsUU::doWork	(	int	row_side,
		int	col_side,
		EntityType	row_type,
		EntityType	col_type,
		EntData &	row_data,
		EntData &	col_data
	)

Definition at line 30 of file PoissonOperators.hpp.

                                           {
    MoFEMFunctionBegin;
 
    const int nb_row_dofs = row_data.getIndices().size();
    const int nb_col_dofs = col_data.getIndices().size();
 
    if (nb_row_dofs && nb_col_dofs) {
      locLhs.resize(nb_row_dofs, nb_col_dofs, false);
      locLhs.clear();
      const int nb_integration_pts = getGaussPts().size2();
      auto t_w = getFTensor0IntegrationWeight();
 
      auto t_row_diff_base = row_data.getFTensor1DiffN<2>();
 
      const double vol = getMeasure();
      for (int gg = 0; gg != nb_integration_pts; ++gg) {
        const double a = t_w * vol;
        for (int rr = 0; rr != nb_row_dofs; ++rr) {
          auto t_col_diff_base = col_data.getFTensor1DiffN<2>(gg, 0);
          for (int cc = 0; cc != nb_col_dofs; ++cc) {
            locLhs(rr, cc) += a * t_row_diff_base(i) * t_col_diff_base(i);
            ++t_col_diff_base;
          } // endFor cc
          ++t_row_diff_base;
        } // endFor rr
        ++t_w;
      } // endFor gg
      CHKERR MatSetValues(getFEMethod()->ksp_B, row_data, col_data,
                          &locLhs(0, 0), ADD_VALUES);
      if (row_side != col_side || row_type != col_type) {
        transLocLhs.resize(nb_col_dofs, nb_row_dofs, false);
        noalias(transLocLhs) = trans(locLhs);
        CHKERR MatSetValues(getFEMethod()->ksp_B, col_data, row_data,
                            &transLocLhs(0, 0), ADD_VALUES);
      } // (row_side != col_side || row_type != col_type)
    }   // endIf (nb_row_dofs && nb_col_dofs)
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

◆ locLhs

MatrixDouble PoissonOps::OpLhsUU::locLhs

private

Definition at line 73 of file PoissonOperators.hpp.

◆ transLocLhs

MatrixDouble PoissonOps::OpLhsUU::transLocLhs