Evaluate error. More...

#include <tutorials/cor-2to5/src/PoissonOperators.hpp>

Inheritance diagram for PoissonExample::OpError:

Collaboration diagram for PoissonExample::OpError:

Public Types
typedef boost::function< double(const double, const double, const double)>	UVal

typedef boost::function< FTensor::Tensor1< double, 3 > const double, const double, const double)>	GVal

Public Member Functions
	OpError (UVal u_value, GVal g_value, boost::shared_ptr< VectorDouble > &field_vals, boost::shared_ptr< MatrixDouble > &grad_vals, Vec global_error)

MoFEMErrorCode	doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &row_data)

Public Member Functions inherited from PoissonExample::OpBaseRhs< VolumeElementForcesAndSourcesCore::UserDataOperator >
	OpBaseRhs (const std::string field_name)

MoFEMErrorCode	doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &row_data)
	This function is called by finite element. More...

virtual MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)=0
	Class dedicated to integrate operator. More...

virtual MoFEMErrorCode	aSsemble (EntitiesFieldData::EntData &data)=0
	Class dedicated to assemble operator to global system vector. More...

Private Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)
	Integrate error. More...

MoFEMErrorCode	aSsemble (EntitiesFieldData::EntData &data)
	Assemble error. More...

Private Attributes
Vec	globalError
	ghost vector with global (integrated over volume) error More...

FTensor::Number< 0 >	NX

FTensor::Number< 1 >	NY

FTensor::Number< 2 >	NZ

FTensor::Index< 'i', 3 >	i

UVal	uValue
	function with exact solution More...

GVal	gValue
	function with exact solution for gradient More...

boost::shared_ptr< VectorDouble >	fieldVals

boost::shared_ptr< MatrixDouble >	gradVals

Additional Inherited Members
Protected Attributes inherited from PoissonExample::OpBaseRhs< VolumeElementForcesAndSourcesCore::UserDataOperator >
int	nbRows
	< error code More...

int	nbIntegrationPts
	number of integration points More...

Detailed Description

Evaluate error.

Definition at line 494 of file PoissonOperators.hpp.

Member Typedef Documentation

◆ GVal

typedef boost::function<FTensor::Tensor1<double, 3> const double, const double, const double)> PoissonExample::OpError::GVal

Definition at line 501 of file PoissonOperators.hpp.

◆ UVal

typedef boost::function<double(const double, const double, const double)> PoissonExample::OpError::UVal

Definition at line 498 of file PoissonOperators.hpp.

Constructor & Destructor Documentation

◆ OpError()

PoissonExample::OpError::OpError	(	UVal	u_value,
		GVal	g_value,
		boost::shared_ptr< VectorDouble > &	field_vals,
		boost::shared_ptr< MatrixDouble > &	grad_vals,
		Vec	global_error
	)

inline

Definition at line 503 of file PoissonOperators.hpp.

       : OpBaseRhs<VolumeElementForcesAndSourcesCore::UserDataOperator>("ERROR"),
         globalError(global_error), uValue(u_value), gValue(g_value),
         fieldVals(field_vals), gradVals(grad_vals) {}

Member Function Documentation

◆ aSsemble()

MoFEMErrorCode PoissonExample::OpError::aSsemble ( EntitiesFieldData::EntData & data )

inlineprivate

Assemble error.

Definition at line 578 of file PoissonOperators.hpp.

                                                           {
     MoFEMFunctionBegin;
     // set error on mesh
     data.getFieldDofs()[0]->getFieldData() = sqrt(data.getFieldData()[0]);
     // assemble vector to global error
     CHKERR VecSetValue(globalError, 0, data.getFieldData()[0], ADD_VALUES);
     MoFEMFunctionReturn(0);
   }

◆ doWork()

MoFEMErrorCode PoissonExample::OpError::doWork	(	int	row_side,
		EntityType	row_type,
		EntitiesFieldData::EntData &	row_data
	)

inline

Definition at line 510 of file PoissonOperators.hpp.

                                                             {
     MoFEMFunctionBegin;
     nbRows = row_data.getFieldData().size();
     if (!nbRows)
       MoFEMFunctionReturnHot(0);
     nbIntegrationPts = getGaussPts().size2();
     CHKERR iNtegrate(row_data);
     CHKERR aSsemble(row_data);
     MoFEMFunctionReturn(0);
   }

◆ iNtegrate()

MoFEMErrorCode PoissonExample::OpError::iNtegrate ( EntitiesFieldData::EntData & data )

inlineprivate

Integrate error.

Definition at line 538 of file PoissonOperators.hpp.

                                                            {
     MoFEMFunctionBegin;
     // clear field dofs
     data.getFieldData().clear();
     // get volume of element
     const double vol = getVolume();
     // get integration weight
     auto t_w = getFTensor0IntegrationWeight();
     // get solution at integration point
     auto t_u = getFTensor0FromVec(*fieldVals);
     // get solution at integration point
     auto t_grad = getFTensor1FromMat<3>(*gradVals);
     // get coordinates at integration point
     auto t_coords = getFTensor1CoordsAtGaussPts();
     // keep exact gradient and error or gradient
     FTensor::Tensor1<double, 3> t_exact_grad, t_error_grad;
     // integrate over
     for (int gg = 0; gg != nbIntegrationPts; gg++) {
       double alpha = vol * t_w;
       // evaluate exact value
       double exact_u = uValue(t_coords(NX), t_coords(NY), t_coords(NZ));
       // evaluate exact gradient
       t_exact_grad = gValue(t_coords(NX), t_coords(NY), t_coords(NZ));
       // calculate gradient errro
       t_error_grad(i) = t_grad(i) - t_exact_grad(i);
       // error
       double error = pow(t_u - exact_u, 2) + t_error_grad(i) * t_error_grad(i);
       // iterate over base functions
       data.getFieldData()[0] += alpha * error;
       ++t_w;      // move to next integration point
       ++t_u;      // next value of function at integration point
       ++t_grad;   // next gradient at integration point
       ++t_coords; // next coordinate at integration point
     }
     MoFEMFunctionReturn(0);
   }

Member Data Documentation

◆ fieldVals

boost::shared_ptr<VectorDouble> PoissonExample::OpError::fieldVals

private

Definition at line 532 of file PoissonOperators.hpp.

◆ globalError

Vec PoissonExample::OpError::globalError

private

ghost vector with global (integrated over volume) error

Definition at line 523 of file PoissonOperators.hpp.

◆ gradVals

boost::shared_ptr<MatrixDouble> PoissonExample::OpError::gradVals

private

Definition at line 533 of file PoissonOperators.hpp.

◆ gValue

GVal PoissonExample::OpError::gValue

private

function with exact solution for gradient

Definition at line 530 of file PoissonOperators.hpp.

◆ i

FTensor::Index<'i', 3> PoissonExample::OpError::i

private

Definition at line 528 of file PoissonOperators.hpp.

◆ NX

FTensor::Number<0> PoissonExample::OpError::NX

private

Definition at line 525 of file PoissonOperators.hpp.

◆ NY

FTensor::Number<1> PoissonExample::OpError::NY

private

Definition at line 526 of file PoissonOperators.hpp.

◆ NZ

FTensor::Number<2> PoissonExample::OpError::NZ

private

Definition at line 527 of file PoissonOperators.hpp.

◆ uValue

UVal PoissonExample::OpError::uValue

private

function with exact solution

Definition at line 529 of file PoissonOperators.hpp.

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

tutorials/cor-2to5/src/PoissonOperators.hpp

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ GVal

◆ UVal

Constructor & Destructor Documentation

◆ OpError()

Member Function Documentation

◆ aSsemble()

◆ doWork()

◆ iNtegrate()

Member Data Documentation

◆ fieldVals

◆ globalError

◆ gradVals

◆ gValue

◆ i

◆ NX

◆ NY

◆ NZ

◆ uValue