Integrate vector on lhs,. More...

#include <users_modules/minimal_surface_equation/src/MinimalSurfaceElement.hpp>

Inheritance diagram for MinimalSurfaceEquation::MinimalSurfaceElement::OpAssmebleBcLhs:

Collaboration diagram for MinimalSurfaceEquation::MinimalSurfaceElement::OpAssmebleBcLhs:

Public Member Functions
	OpAssmebleBcLhs (const string field_name, Mat m_a)

PetscErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)

Public Attributes
Mat	mA

MatrixDouble	K

MatrixDouble	transK

Detailed Description

Integrate vector on lhs,.

\[ \mathbf{A} = \int_L \mathbf{N}^\textrm{T}\mathbf{N} \textrm{d}L \]

Examples: minimal_surface_area.cpp.

Definition at line 139 of file MinimalSurfaceElement.hpp.

Constructor & Destructor Documentation

◆ OpAssmebleBcLhs()

MinimalSurfaceEquation::MinimalSurfaceElement::OpAssmebleBcLhs::OpAssmebleBcLhs	(	const string	field_name,
		Mat	m_a
	)

inline

Definition at line 142 of file MinimalSurfaceElement.hpp.

        : EdgeElementForcesAndSourcesCore::UserDataOperator(
              field_name, UserDataOperator::OPROWCOL),
          mA(m_a) {}

Member Function Documentation

◆ doWork()

PetscErrorCode MinimalSurfaceEquation::MinimalSurfaceElement::OpAssmebleBcLhs::doWork	(	int	row_side,
		int	col_side,
		EntityType	row_type,
		EntityType	col_type,
		DataForcesAndSourcesCore::EntData &	row_data,
		DataForcesAndSourcesCore::EntData &	col_data
	)

inline

Definition at line 147 of file MinimalSurfaceElement.hpp.

                                                                     {
      MoFEMFunctionBegin;
 
      int row_nb_dofs = row_data.getIndices().size();
      if (row_nb_dofs == 0) {
        MoFEMFunctionReturnHot(0);
      }
      int col_nb_dofs = col_data.getIndices().size();
      if (col_nb_dofs == 0) {
        MoFEMFunctionReturnHot(0);
      }
      K.resize(row_nb_dofs, col_nb_dofs, false);
      K.clear();
      int nb_gauss_pts = row_data.getN().size1();
      for (int gg = 0; gg != nb_gauss_pts; gg++) {
        double val = getLength() * getGaussPts()(1, gg);
        VectorAdaptor row_n = row_data.getN(gg);
        VectorAdaptor col_n = col_data.getN(gg);
        noalias(K) += val * outer_prod(row_n, col_n);
      }
      CHKERR MatSetValues(mA, row_nb_dofs,
                          &*row_data.getIndices().data().begin(), col_nb_dofs,
                          &*col_data.getIndices().data().begin(),
                          &*K.data().begin(), ADD_VALUES);
      // Matrix is symmetric, assemble off diagonal part
      if (row_type != col_type || row_side != col_side) {
        transK.resize(col_nb_dofs, row_nb_dofs, false);
        noalias(transK) = trans(K);
        CHKERR MatSetValues(mA, col_nb_dofs,
                            &*col_data.getIndices().data().begin(), row_nb_dofs,
                            &*row_data.getIndices().data().begin(),
                            &*transK.data().begin(), ADD_VALUES);
      }
 
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

◆ K

MatrixDouble MinimalSurfaceEquation::MinimalSurfaceElement::OpAssmebleBcLhs::K

Definition at line 146 of file MinimalSurfaceElement.hpp.

◆ mA

Mat MinimalSurfaceEquation::MinimalSurfaceElement::OpAssmebleBcLhs::mA

Definition at line 141 of file MinimalSurfaceElement.hpp.

◆ transK

MatrixDouble MinimalSurfaceEquation::MinimalSurfaceElement::OpAssmebleBcLhs::transK