Inheritance diagram for OpError:

Collaboration diagram for OpError:

Public Member Functions
	OpError (boost::shared_ptr< VectorDouble > u_at_pts_ptr, double &l2_error)

	OpError (const std::string field_name, boost::shared_ptr< MatrixDouble > x_ptr)

MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Static Public Attributes
static SmartPetscObj< Vec >	errorVec

Private Attributes
boost::shared_ptr< VectorDouble >	uAtPtsPtr

double &	l2Error

boost::shared_ptr< MatrixDouble >	xPtr

Detailed Description

Examples: mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/dg_projection.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/higher_derivatives.cpp, mofem/tutorials/cor-2to5/src/PoissonOperators.hpp, mofem/tutorials/mix-0/mixed_poisson.cpp, mofem/tutorials/scl-0/approximaton.cpp, mofem/tutorials/scl-10/initial_diffusion.cpp, and mofem/tutorials/vec-4/approx_sphere.cpp.

Definition at line 61 of file initial_diffusion.cpp.

Constructor & Destructor Documentation

◆ OpError() [1/2]

OpError::OpError	(	boost::shared_ptr< VectorDouble >	u_at_pts_ptr,
		double &	l2_error
	)

inline

Definition at line 62 of file initial_diffusion.cpp.

      : DomainEleOp("PHOTON_FLUENCE_RATE", OPROW), uAtPtsPtr(u_at_pts_ptr),
        l2Error(l2_error) {
    // Only will be executed once on vertices
    std::fill(&doEntities[MBEDGE], &doEntities[MBMAXTYPE], false);
  }

◆ OpError() [2/2]

OpError::OpError	(	const std::string	field_name,
		boost::shared_ptr< MatrixDouble >	x_ptr
	)

inline

Definition at line 221 of file approx_sphere.cpp.

      : DomainEleOp(field_name, field_name, AssemblyDomainEleOp::OPROW),
        xPtr(x_ptr) {
 
    std::fill(&doEntities[MBEDGE], &doEntities[MBMAXTYPE], false);
  }

Member Function Documentation

◆ doWork()

MoFEMErrorCode OpError::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

inline

Examples: mofem/tutorials/vec-4/approx_sphere.cpp.

Definition at line 228 of file approx_sphere.cpp.

                                                                  {
 
    MoFEMFunctionBegin;
 
    auto t_w = getFTensor0IntegrationWeight();
    auto t_x = getFTensor1FromMat<3>(*xPtr);
    auto t_normal = getFTensor1NormalsAtGaussPts();
    auto nb_integration_pts = getGaussPts().size2();
 
    double error = 0;
 
    for (int gg = 0; gg != nb_integration_pts; gg++) {
 
      double l = std::sqrt(t_normal(i) * t_normal(i));
      error += t_w * l * std::abs((t_x(i) * t_x(i) - A * A));
 
      ++t_w;
      ++t_x;
      ++t_normal;
    }
 
    CHKERR VecSetValue(errorVec, 0, error, ADD_VALUES);
 
    MoFEMFunctionReturn(0);
  }