本文整理汇总了C++中dataPtr_Type::bulk方法的典型用法代码示例。如果您正苦于以下问题：C++ dataPtr_Type::bulk方法的具体用法？C++ dataPtr_Type::bulk怎么用？C++ dataPtr_Type::bulk使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类dataPtr_Type的用法示例。
在下文中一共展示了dataPtr_Type::bulk方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: dk_loc
void NeoHookeanMaterialNonLinear<Mesh>::computeStiffness ( const vector_Type& sol,
Real /*factor*/,
const dataPtr_Type& dataMaterial,
const mapMarkerVolumesPtr_Type mapsMarkerVolumes,
const displayerPtr_Type& displayer )
{
this->M_stiff.reset (new vector_Type (*this->M_localMap) );
displayer->leaderPrint (" \n******************************************************************\n ");
displayer->leaderPrint (" Non-Linear S- Computing the Neo-Hookean nonlinear stiffness vector" );
displayer->leaderPrint (" \n******************************************************************\n ");
UInt totalDof = this->M_FESpace->dof().numTotalDof();
UInt dim = this->M_FESpace->dim();
VectorElemental dk_loc ( this->M_FESpace->fe().nbFEDof(), nDimensions );
//vector_Type disp(sol);
vector_Type dRep (sol, Repeated);
mapIterator_Type it;
for ( it = (*mapsMarkerVolumes).begin(); it != (*mapsMarkerVolumes).end(); it++ )
{
//Given the marker pointed by the iterator, let's extract the material parameters
UInt marker = it->first;
Real mu = dataMaterial->mu (marker);
Real bulk = dataMaterial->bulk (marker);
for ( UInt j (0); j < it->second.size(); j++ )
{
this->M_FESpace->fe().updateFirstDerivQuadPt ( * (it->second[j]) );
UInt eleID = this->M_FESpace->fe().currentLocalId();
for ( UInt iNode = 0 ; iNode < ( UInt ) this->M_FESpace->fe().nbFEDof() ; iNode++ )
{
UInt iloc = this->M_FESpace->fe().patternFirst ( iNode );
for ( UInt iComp = 0; iComp < nDimensions; ++iComp )
{
UInt ig = this->M_FESpace->dof().localToGlobalMap ( eleID, iloc ) + iComp * dim + this->M_offset;
dk_loc[ iloc + iComp * this->M_FESpace->fe().nbFEDof() ] = dRep[ig];
}
}
this->M_elvecK->zero();
computeKinematicsVariables ( dk_loc );
//! Stiffness for non-linear terms of the Neo-Hookean model
/*!
The results of the integrals are stored at each step into elvecK, until to build K matrix of the bilinear form
*/
//! Volumetric part
/*!
Source term Pvol: int { bulk /2* (J1^2 - J1 + log(J1) ) * 1/J1 * (CofF1 : \nabla v) }
*/
AssemblyElementalStructure::source_Pvol ( 0.5 * bulk, (*M_CofFk), (*M_Jack),
*this->M_elvecK, this->M_FESpace->fe() );
//! Isochoric part
/*!
Source term P1iso_NH
*/
AssemblyElementalStructure::source_P1iso_NH ( mu, (*M_CofFk) , (*M_Fk), (*M_Jack), (*M_trCisok) ,
*this->M_elvecK, this->M_FESpace->fe() );
for ( UInt ic = 0; ic < nDimensions; ++ic )
{
/*!
M_elvecK is assemble into *vec_stiff vector that is recall
from updateSystem(matrix_ptrtype& mat_stiff, vector_ptr_type& vec_stiff)
*/
assembleVector ( *this->M_stiff,
*this->M_elvecK,
this->M_FESpace->fe(),
this->M_FESpace->dof(),
ic, this->M_offset + ic * totalDof );
}
}
}
this->M_stiff->globalAssemble();
}