109                                                        {
  111    const int nb_gauss_pts = getGaussPts().size2();
  112    if (type == MBVERTEX) {
  113      vAls.resize(nb_gauss_pts, 
false);
 
  117    }
  118 
  120    if (nb_dofs) {
  121 
  123          << "Type  " << moab::CN::EntityTypeName(type) << " side " << side;
  126 
  129      for (int gg = 0; gg != nb_gauss_pts; gg++) {
  131        for (int bb = 0; bb != nb_dofs; bb++) {
  132          t_vals += t_base_fun * t_data;
  133          ++t_base_fun;
  134          ++t_data;
  135        }
  136        const double v = t_vals;
 
  137        if (!std::isnormal(
v))
 
  139        ++t_vals;
  140      }
  141 
  143        auto t_diff_vals = getFTensor1FromMat<SPACE_DIM>(
diffVals);
 
  145        for (int gg = 0; gg != nb_gauss_pts; gg++) {
  147          for (int bb = 0; bb != nb_dofs; bb++) {
  148            t_diff_vals(
i) += t_diff_base_fun(
i) * t_data;
 
  149            ++t_diff_base_fun;
  150            ++t_data;
  151          }
  153            if (!std::isnormal(t_diff_vals(d)))
  155          ++t_diff_vals;
  156        }
  157      }
  158    }
  160  }
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
@ MOFEM_ATOM_TEST_INVALID
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define MOFEM_LOG(channel, severity)
Log.
const double v
phase velocity of light in medium (cm/ns)
static auto getFTensor0FromVec(ublas::vector< T, A > &data)
Get tensor rank 0 (scalar) form data vector.
FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim > getFTensor1DiffN(const FieldApproximationBase base)
Get derivatives of base functions.
FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > getFTensor0N(const FieldApproximationBase base)
Get base function as Tensor0.
MatrixDouble & getDiffN(const FieldApproximationBase base)
get derivatives of base functions
MatrixDouble & getN(const FieldApproximationBase base)
get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb....
FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > getFTensor0FieldData()
Return scalar files as a FTensor of rank 0.
const VectorInt & getIndices() const
Get global indices of degrees of freedom on entity.
FTensor::Index< 'i', SPACE_DIM > i