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

Definition at line 102 of file scalar_check_approximation.cpp.

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

104 checkGradients(check_grads) {}

Member Function Documentation

◆ doWork()

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

inline

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);
   }