StdRDOperators::OpError Struct Reference

#include <users_modules/softmech/chemo_mech/src/rd_stdOperators.hpp>

Inheritance diagram for StdRDOperators::OpError:

Collaboration diagram for StdRDOperators::OpError:

Public Types
typedef boost::function< double(const double, const double, const double)>	FVal
typedef boost::function< FTensor::Tensor1< double, 3 > const double, const double, const double)>	FGrad

Public Member Functions
	OpError (FVal exact_value, FVal exact_lap, FGrad exact_grad, boost::shared_ptr< PreviousData > &prev_data, std::map< int, BlockData > &block_map, double &err)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Public Attributes
double &	eRror

Private Attributes
FVal	exactVal
FVal	exactLap
FGrad	exactGrad
boost::shared_ptr< PreviousData >	prevData
std::map< int, BlockData >	setOfBlock
FTensor::Number< 0 >	NX
FTensor::Number< 1 >	NY

Detailed Description

Definition at line 551 of file rd_stdOperators.hpp.

Member Typedef Documentation

FGrad

typedef boost::function<FTensor::Tensor1<double, 3> const double, const double, const double)> StdRDOperators::OpError::FGrad

Definition at line 556 of file rd_stdOperators.hpp.

FVal

typedef boost::function<double(const double, const double, const double)> StdRDOperators::OpError::FVal

Definition at line 553 of file rd_stdOperators.hpp.

Constructor & Destructor Documentation

OpError()

StdRDOperators::OpError::OpError	(	FVal	exact_value,
		FVal	exact_lap,
		FGrad	exact_grad,
		boost::shared_ptr< PreviousData > &	prev_data,
		std::map< int, BlockData > &	block_map,
		double &	err
	)

Definition at line 558 of file rd_stdOperators.hpp.

      : OpEle("ERROR", OpEle::OPROW), exactVal(exact_value),
        exactLap(exact_lap), exactGrad(exact_grad), prevData(prev_data),
        setOfBlock(block_map), eRror(err) {}

Member Function Documentation

doWork()

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

Definition at line 564 of file rd_stdOperators.hpp.

                                                                  {
    MoFEMFunctionBegin;
    const int nb_dofs = data.getFieldData().size();
    // cout << "nb_error_dofs : " << nb_dofs << endl;
    if (nb_dofs) {
 
      auto find_block_data = [&]() {
        EntityHandle fe_ent = getFEEntityHandle();
        BlockData *block_raw_ptr = nullptr;
        for (auto &m : setOfBlock) {
          if (m.second.block_ents.find(fe_ent) != m.second.block_ents.end()) {
            block_raw_ptr = &m.second;
            break;
          }
        }
        return block_raw_ptr;
      };
 
      auto block_data_ptr = find_block_data();
      if (!block_data_ptr)
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "Block not found");
 
      auto &block_data = *block_data_ptr;
 
      auto t_value = getFTensor0FromVec(prevData->values);
      auto t_grad = getFTensor1FromMat<2>(prevData->grads);
      // cout << "t_grad : " << t_grad << endl;
      // auto t_grad = getFTensor1FromMat<3>(prevData->grads);
      data.getFieldData().clear();
      const double vol = getMeasure();
      const int nb_integration_pts = getGaussPts().size2();
      auto t_w = getFTensor0IntegrationWeight();
      double dt;
      CHKERR TSGetTimeStep(getFEMethod()->ts, &dt);
      double ct = getFEMethod()->ts_t - dt;
      auto t_coords = getFTensor1CoordsAtGaussPts();
 
      FTensor::Index<'j', 2> j;
      // cout << "B0 : " << block_data.B0 << endl;
 
      for (int gg = 0; gg != nb_integration_pts; ++gg) {
        const double a = vol * t_w;
 
        double mass_exact = exactVal(t_coords(NX), t_coords(NY), ct);
        double flux_lap =
            -block_data.B0 * exactLap(t_coords(NX), t_coords(NY), ct);
        auto flux_exact = exactGrad(t_coords(NX), t_coords(NY), ct);
 
        // cout << "grad_exact : " << flux_exact << endl;
        // cout << "grad_value : " << t_grad << endl;
        // cout << "--------------------" << endl;
 
        double flux_error =
            pow(block_data.B0, 2) * (pow(flux_exact(0) - t_grad(0), 2) +
                                     pow(flux_exact(1) - t_grad(1), 2));
 
        double local_error = pow(mass_exact - t_value, 2); // + flux_error;
 
        data.getFieldData()[0] += a * local_error;
        eRror += a * local_error;
 
        ++t_w;
        ++t_value;
        ++t_grad;
        ++t_coords;
      }
 
      data.getFieldDofs()[0]->getFieldData() = data.getFieldData()[0];
    }
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

eRror

double& StdRDOperators::OpError::eRror

Definition at line 557 of file rd_stdOperators.hpp.

exactGrad

FGrad StdRDOperators::OpError::exactGrad

private

Definition at line 639 of file rd_stdOperators.hpp.

exactLap

FVal StdRDOperators::OpError::exactLap

private

Definition at line 638 of file rd_stdOperators.hpp.

exactVal

FVal StdRDOperators::OpError::exactVal

private

Definition at line 637 of file rd_stdOperators.hpp.

NX

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

private

Definition at line 643 of file rd_stdOperators.hpp.

NY

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

private

Definition at line 644 of file rd_stdOperators.hpp.

prevData

boost::shared_ptr<PreviousData> StdRDOperators::OpError::prevData

private

Definition at line 640 of file rd_stdOperators.hpp.

setOfBlock

std::map<int, BlockData> StdRDOperators::OpError::setOfBlock

private

Definition at line 641 of file rd_stdOperators.hpp.

---

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

• users_modules/softmech/chemo_mech/src/rd_stdOperators.hpp