本文整理汇总了C++中AssemblyContext::fixed_interior_hessian方法的典型用法代码示例。如果您正苦于以下问题：C++ AssemblyContext::fixed_interior_hessian方法的具体用法？C++ AssemblyContext::fixed_interior_hessian怎么用？C++ AssemblyContext::fixed_interior_hessian使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类AssemblyContext的用法示例。
在下文中一共展示了AssemblyContext::fixed_interior_hessian方法的2个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: rhocpU
void HeatTransferStabilizationHelper::compute_res_energy_steady_and_derivs
( AssemblyContext& context,
unsigned int qp,
const libMesh::Real rho,
const libMesh::Real Cp,
const libMesh::Real k,
libMesh::Real &res,
libMesh::Real &d_res_dT,
libMesh::Gradient &d_res_dgradT,
libMesh::Tensor &d_res_dhessT,
libMesh::Gradient &d_res_dU
) const
{
libMesh::Gradient grad_T = context.fixed_interior_gradient(this->_temp_vars.T(), qp);
libMesh::Tensor hess_T = context.fixed_interior_hessian(this->_temp_vars.T(), qp);
libMesh::RealGradient rhocpU( rho*Cp*context.fixed_interior_value(this->_flow_vars.u(), qp),
rho*Cp*context.fixed_interior_value(this->_flow_vars.v(), qp) );
if(this->_flow_vars.dim() == 3)
rhocpU(2) = rho*Cp*context.fixed_interior_value(this->_flow_vars.w(), qp);
res = rhocpU*grad_T - k*(hess_T(0,0) + hess_T(1,1) + hess_T(2,2));
d_res_dT = 0;
d_res_dgradT = rhocpU;
d_res_dhessT = 0;
d_res_dhessT(0,0) = -k;
d_res_dhessT(1,1) = -k;
d_res_dhessT(2,2) = -k;
d_res_dU = rho * Cp * grad_T;
}示例2: compute_res_energy_steady
libMesh::Real HeatTransferStabilizationHelper::compute_res_energy_steady( AssemblyContext& context,
unsigned int qp,
const libMesh::Real rho,
const libMesh::Real Cp,
const libMesh::Real k ) const
{
libMesh::Gradient grad_T = context.fixed_interior_gradient(this->_temp_vars.T(), qp);
libMesh::Tensor hess_T = context.fixed_interior_hessian(this->_temp_vars.T(), qp);
libMesh::RealGradient rhocpU( rho*Cp*context.fixed_interior_value(this->_flow_vars.u(), qp),
rho*Cp*context.fixed_interior_value(this->_flow_vars.v(), qp) );
if(this->_flow_vars.dim() == 3)
rhocpU(2) = rho*Cp*context.fixed_interior_value(this->_flow_vars.w(), qp);
return rhocpU*grad_T - k*(hess_T(0,0) + hess_T(1,1) + hess_T(2,2));
}