36 (
const std::string & helper_name ) :
49 unsigned int qp )
const
51 libMesh::RealGradient g( fe->get_dxidx()[qp] + fe->get_detadx()[qp],
52 fe->get_dxidy()[qp] + fe->get_detady()[qp] );
54 if( c.get_dim() == 3 )
56 g(0) += fe->get_dzetadx()[qp];
57 g(1) += fe->get_dzetady()[qp];
58 g(2) = fe->get_dxidz()[qp] + fe->get_detadz()[qp] + fe->get_dzetadz()[qp];
66 unsigned int qp )
const
68 libMesh::Real dxidx = fe->get_dxidx()[qp];
69 libMesh::Real dxidy = fe->get_dxidy()[qp];
71 libMesh::Real detadx = fe->get_detadx()[qp];
72 libMesh::Real detady = fe->get_detady()[qp];
74 libMesh::RealTensor G( dxidx*dxidx + detadx*detadx,
75 dxidx*dxidy + detadx*detady,
77 dxidy*dxidx + detady*detadx,
78 dxidy*dxidy + detady*detady,
81 if( c.get_dim() == 3 )
83 libMesh::Real dxidz = fe->get_dxidz()[qp];
85 libMesh::Real detadz = fe->get_detadz()[qp];
87 libMesh::Real dzetadx = fe->get_dzetadx()[qp];
88 libMesh::Real dzetady = fe->get_dzetady()[qp];
89 libMesh::Real dzetadz = fe->get_dzetadz()[qp];
91 G(0,0) += dzetadx*dzetadx;
92 G(0,1) += dzetadx*dzetady;
93 G(0,2) = dxidx*dxidz + detadx*detadz + dzetadx*dzetadz;
94 G(1,0) += dzetady*dzetadx;
95 G(1,1) += dzetady*dzetady;
96 G(1,2) = dxidy*dxidz + detady*detadz + dzetady*dzetadz;
97 G(2,0) = dxidz*dxidx + detadz*detadx + dzetadz*dzetadx;
98 G(2,1) = dxidz*dxidy + detadz*detady + dzetadz*dzetady;
99 G(2,2) = dxidz*dxidz + detadz*detadz + dzetadz*dzetadz;
libMesh::RealTensor compute_G(libMesh::FEBase *fe, AssemblyContext &c, unsigned int qp) const
libMesh::RealGradient compute_g(libMesh::FEBase *fe, AssemblyContext &c, unsigned int qp) const
ParameterUser base class. Utility methods for subclasses.
StabilizationHelper(const std::string &helper_name)