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