Evaluate function values and gradients at Gauss Pts. More...

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

Inheritance diagram for MinimalSurfaceEquation::MinimalSurfaceElement::OpCalculateCoefficientsAtGaussPts:

Collaboration diagram for MinimalSurfaceEquation::MinimalSurfaceElement::OpCalculateCoefficientsAtGaussPts:

Public Member Functions
	OpCalculateCoefficientsAtGaussPts (const string field_name, CommonData &common_data, bool rhs_and_not_lhs)

PetscErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Public Attributes
CommonData &	commonData

bool	rhsAndNotLhs

Detailed Description

Evaluate function values and gradients at Gauss Pts.

Examples: minimal_surface_area.cpp.

Definition at line 228 of file MinimalSurfaceElement.hpp.

Constructor & Destructor Documentation

◆ OpCalculateCoefficientsAtGaussPts()

MinimalSurfaceEquation::MinimalSurfaceElement::OpCalculateCoefficientsAtGaussPts::OpCalculateCoefficientsAtGaussPts	(	const string	field_name,
		CommonData &	common_data,
		bool	rhs_and_not_lhs
	)

inline

Definition at line 232 of file MinimalSurfaceElement.hpp.

        : FaceElementForcesAndSourcesCore::UserDataOperator(
              field_name, UserDataOperator::OPCOL),
          commonData(common_data), rhsAndNotLhs(rhs_and_not_lhs) {}

Member Function Documentation

◆ doWork()

PetscErrorCode MinimalSurfaceEquation::MinimalSurfaceElement::OpCalculateCoefficientsAtGaussPts::doWork	(	int	side,
		EntityType	type,
		DataForcesAndSourcesCore::EntData &	data
	)

inline

Definition at line 239 of file MinimalSurfaceElement.hpp.

                                                                 {
      MoFEMFunctionBegin;
 
      // Element loops over entities, it start from vertices
      if (type == MBVERTEX) {
        int nb_gauss_pts = data.getN().size1();
        commonData.normGradU2.resize(nb_gauss_pts, false);
        for (int gg = 0; gg != nb_gauss_pts; gg++) {
          ublas::matrix_row<MatrixDouble> grad_at_gauss_pt(commonData.gradU,
                                                           gg);
          commonData.normGradU2[gg] =
              inner_prod(grad_at_gauss_pt, grad_at_gauss_pt);
        }
        commonData.aN.resize(nb_gauss_pts, false);
        for (int gg = 0; gg != nb_gauss_pts; gg++) {
          commonData.aN[gg] = 1. / sqrt(1 + commonData.normGradU2[gg]);
        }
        if (rhsAndNotLhs) {
          // needed at right hand side when residual is calculated
          commonData.aNbyGradU.resize(nb_gauss_pts, 2, false);
          for (int gg = 0; gg != nb_gauss_pts; gg++) {
            ublas::matrix_row<MatrixDouble> grad_at_gauss_pt(commonData.gradU,
                                                             gg);
            for (int rr = 0; rr != 2; rr++) {
              commonData.aNbyGradU(gg, rr) =
                  commonData.aN[gg] * grad_at_gauss_pt[rr];
            }
          }
        } else {
          // needed at left hand side when matrix is calculated
          commonData.aNpow3byGradU.resize(nb_gauss_pts, 2, false);
          for (int gg = 0; gg != nb_gauss_pts; gg++) {
            ublas::matrix_row<MatrixDouble> grad_at_gauss_pt(commonData.gradU,
                                                             gg);
            for (int rr = 0; rr != 2; rr++) {
              commonData.aNpow3byGradU(gg, rr) =
                  pow(commonData.aN[gg], 3) * grad_at_gauss_pt[rr];
            }
          }
        }
      }
 
      // doVerticesRow = true;
      // doEdgesRow = false;
      // doQuadsRow = false;
      // doTrisRow = false;
      // doTetsRow = false;
      // doPrismsRow = false;
 
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

◆ commonData

CommonData& MinimalSurfaceEquation::MinimalSurfaceElement::OpCalculateCoefficientsAtGaussPts::commonData

Definition at line 230 of file MinimalSurfaceElement.hpp.

◆ rhsAndNotLhs

bool MinimalSurfaceEquation::MinimalSurfaceElement::OpCalculateCoefficientsAtGaussPts::rhsAndNotLhs

Definition at line 231 of file MinimalSurfaceElement.hpp.

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

users_modules/minimal_surface_equation/src/MinimalSurfaceElement.hpp

Public Member Functions

Public Attributes