OpValsDiffVals Struct Reference

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

mofem/atom_tests/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

Definition at line 102 of file scalar_check_approximation.cpp.

      : DomainEleOp("FIELD1", OPROW), vAls(vals), diffVals(diff_vals),
        checkGradients(check_grads) {}

Member Function Documentation

doWork()

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

inline

Examples

mofem/atom_tests/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);
  }

Member Data Documentation

checkGradients

const bool OpValsDiffVals::checkGradients

Examples

mofem/atom_tests/scalar_check_approximation.cpp.

Definition at line 101 of file scalar_check_approximation.cpp.

diffVals

MatrixDouble& OpValsDiffVals::diffVals

Examples

mofem/atom_tests/scalar_check_approximation.cpp.

Definition at line 100 of file scalar_check_approximation.cpp.

i

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

Examples

mofem/atom_tests/scalar_check_approximation.cpp.

Definition at line 106 of file scalar_check_approximation.cpp.

vAls

VectorDouble& OpValsDiffVals::vAls

Examples

mofem/atom_tests/scalar_check_approximation.cpp.

Definition at line 99 of file scalar_check_approximation.cpp.

---

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

• atom_tests/scalar_check_approximation.cpp