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