Inheritance diagram for OpValsDiffVals:

Collaboration diagram for OpValsDiffVals:

Public Member Functions
	OpValsDiffVals (VectorDouble &vals, MatrixDouble &diff_vals, bool check_grads)

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

Public Attributes
VectorDouble &	vAls

MatrixDouble &	diffVals

const bool	checkGradients

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

Detailed Description

Examples: scalar_check_approximation.cpp.

Definition at line 98 of file scalar_check_approximation.cpp.

Constructor & Destructor Documentation

◆ OpValsDiffVals()

OpValsDiffVals::OpValsDiffVals	(	VectorDouble &	vals,
		MatrixDouble &	diff_vals,
		bool	check_grads )

inline

Examples: scalar_check_approximation.cpp.

Definition at line 102 of file scalar_check_approximation.cpp.

103 : DomainEleOp("FIELD1", OPROW), vAls(vals), diffVals(diff_vals),

104 checkGradients(check_grads) {}

DomainEleOp

DomainEle::UserDataOperator DomainEleOp

Finire element operator type.

Definition child_and_parent.cpp:36

OpValsDiffVals::checkGradients

const bool checkGradients

Definition scalar_check_approximation.cpp:101

OpValsDiffVals::vAls

VectorDouble & vAls

Definition scalar_check_approximation.cpp:99

OpValsDiffVals::diffVals

MatrixDouble & diffVals

Definition scalar_check_approximation.cpp:100

Member Function Documentation

◆ doWork()

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

inline

Examples: scalar_check_approximation.cpp.

Definition at line 108 of file scalar_check_approximation.cpp.

                                                        {
    MoFEMFunctionBegin;
    const int nb_gauss_pts = getGaussPts().size2();
    if (type == MBVERTEX) {
      vAls.resize(nb_gauss_pts, false);
      diffVals.resize(SPACE_DIM, nb_gauss_pts, false);
      vAls.clear();
      diffVals.clear();
    }
 
    const int nb_dofs = data.getIndices().size();
    if (nb_dofs) {
 
      MOFEM_LOG("AT", Sev::noisy)
          << "Type  " << moab::CN::EntityTypeName(type) << " side " << side;
      MOFEM_LOG("AT", Sev::noisy) << data.getN();
      MOFEM_LOG("AT", Sev::noisy) << data.getDiffN();
 
      auto t_vals = getFTensor0FromVec(vAls);
      auto t_base_fun = data.getFTensor0N();
      for (int gg = 0; gg != nb_gauss_pts; gg++) {
        auto t_data = data.getFTensor0FieldData();
        for (int bb = 0; bb != nb_dofs; bb++) {
          t_vals += t_base_fun * t_data;
          ++t_base_fun;
          ++t_data;
        }
        const double v = t_vals;
        if (!std::isnormal(v))
          SETERRQ(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID, "Not a number");
        ++t_vals;
      }
 
      if (checkGradients) {
        auto t_diff_vals = getFTensor1FromMat<SPACE_DIM>(diffVals);
        auto t_diff_base_fun = data.getFTensor1DiffN<SPACE_DIM>();
        for (int gg = 0; gg != nb_gauss_pts; gg++) {
          auto t_data = data.getFTensor0FieldData();
          for (int bb = 0; bb != nb_dofs; bb++) {
            t_diff_vals(i) += t_diff_base_fun(i) * t_data;
            ++t_diff_base_fun;
            ++t_data;
          }
          for (int d = 0; d != SPACE_DIM; ++d)
            if (!std::isnormal(t_diff_vals(d)))
              SETERRQ(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID, "Not a number");
          ++t_diff_vals;
        }
      }
    }
    MoFEMFunctionReturn(0);
  }