template<int DIM_0, int DIM_1>
struct OpCalculateFStab< DIM_0, DIM_1 >
- Examples
- mofem/tutorials/adv-4/dynamic_first_order_con_law.cpp.
Definition at line 103 of file dynamic_first_order_con_law.cpp.
◆ OpCalculateFStab()
template<int DIM_0, int DIM_1> 
 
Definition at line 104 of file dynamic_first_order_con_law.cpp.
  110      : ForcesAndSourcesCore::UserDataOperator(
NOSPACE, OPLAST),
 
boost::shared_ptr< MatrixDouble > gradxPtr
boost::shared_ptr< MatrixDouble > defGradStabPtr
boost::shared_ptr< MatrixDouble > gradVelPtr
boost::shared_ptr< MatrixDouble > defGradPtr
boost::shared_ptr< MatrixDouble > defGradDotPtr
 
 
◆ doWork()
template<int DIM_0, int DIM_1> 
 
- Examples
- mofem/tutorials/adv-4/dynamic_first_order_con_law.cpp.
Definition at line 115 of file dynamic_first_order_con_law.cpp.
  116                                                               {
  118    
  121 
  122    
  123    const size_t nb_gauss_pts = getGaussPts().size2();
  124 
  127 
  128    
  129    auto t_F = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*
defGradPtr);
 
  130    auto t_Fstab = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*
defGradStabPtr);
 
  131    auto t_F_dot = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*
defGradDotPtr);
 
  132 
  133    
  136    auto t_gradx = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*
gradxPtr);
 
  137    auto t_gradVel = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*
gradVelPtr);
 
  138 
  139    for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
  140      
  141      t_Fstab(
i, 
j) = t_F(
i, 
j) + tau_F * (t_gradVel(
i, 
j) - t_F_dot(
i, 
j)) +
 
  142                      xi_F * (t_gradx(
i, 
j) - t_F(
i, 
j));
 
  143 
  144      ++t_F;
  145      ++t_Fstab;
  146      ++t_gradVel;
  147      ++t_F_dot;
  148 
  149      ++t_gradx;
  150    }
  151 
  153  }
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'j', 3 > j
 
 
◆ defGradDotPtr
template<int DIM_0, int DIM_1> 
 
 
◆ defGradPtr
template<int DIM_0, int DIM_1> 
 
 
◆ defGradStabPtr
template<int DIM_0, int DIM_1> 
 
 
◆ gradVelPtr
template<int DIM_0, int DIM_1> 
 
 
◆ gradxPtr
template<int DIM_0, int DIM_1> 
 
 
◆ tauFPtr
template<int DIM_0, int DIM_1> 
 
 
◆ xiF
template<int DIM_0, int DIM_1> 
 
 
The documentation for this struct was generated from the following file: