本文整理汇总了C++中MIP::IP方法的典型用法代码示例。如果您正苦于以下问题：C++ MIP::IP方法的具体用法？C++ MIP::IP怎么用？C++ MIP::IP使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类MIP的用法示例。
在下文中一共展示了MIP::IP方法的2个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: GenerateMatrix
static void GenerateMatrix (const FEL & bfel, const MIP & mip,
MAT & mat, LocalHeap & lh)
{
// must get the right elements, otherwise an exception is thrown.
const CompoundFiniteElement & cfel =
dynamic_cast<const CompoundFiniteElement&> (bfel);
// a scalar H1 element
const ScalarFiniteElement<2> & fel_u =
dynamic_cast<const ScalarFiniteElement<2>&> (cfel[0]);
const ScalarFiniteElement<2> & fel_p =
dynamic_cast<const ScalarFiniteElement<2>&> (cfel[2]);
int nd_u = fel_u.GetNDof();
int nd_p = fel_p.GetNDof();
// transformation of derivatives from reference element to general element:
FlatMatrixFixWidth<2> gradu(nd_u, lh);
fel_u.CalcMappedDShape (mip, gradu);
// the shape functions of the pressure
FlatVector<> vecp(nd_p, lh);
fel_p.CalcShape (mip.IP(), vecp);
mat = 0;
// the first nd_u shape functions belong to u_x, the next nd_u belong to u_y:
mat.Rows(0,2).Cols(cfel.GetRange(0)) = Trans (gradu);
mat.Rows(2,4).Cols(cfel.GetRange(1)) = Trans (gradu);
// ... and finally nd_p shape functions for the pressure:
mat.Row(4).Range(cfel.GetRange(2)) = vecp;
}示例2: GenerateMatrix
static void GenerateMatrix (const AFEL & bfel, const MIP & sip,
MAT & mat, LocalHeap & lh)
{
HeapReset hr(lh);
const HDivDivFiniteElement<2> & fel =
dynamic_cast<const HDivDivFiniteElement<2>&> (bfel);
int nd = fel.GetNDof();
Mat<3,2> jac = sip.GetJacobian();
double det = fabs (sip.GetJacobiDet());
FlatMatrix<> shape(nd, 3, lh);
fel.CalcShape (sip.IP(), shape);
Mat<3,9> trans;
for (int i = 0; i < 3; i++)
{
Mat<2> sigma_ref;
sigma_ref = 0.0;
switch (i)
{
case 0: sigma_ref(0,0) = 1.0; break;
case 1: sigma_ref(1,1) = 1.0; break;
case 2: sigma_ref(0,1) = sigma_ref(1,0) = 1.0; break;
}
auto hm = jac * sigma_ref;
auto sigma = hm * Trans(jac);
sigma *= (1.0 / sqr(det));
trans ( i, 0 ) = sigma(0,0);
trans ( i, 1 ) = sigma(0,1);
trans ( i, 2 ) = sigma(0,2);
trans ( i, 3 ) = sigma(1,0);
trans ( i, 4 ) = sigma(1,1);
trans ( i, 5 ) = sigma(1,2);
trans ( i, 6 ) = sigma(2,0);
trans ( i, 7 ) = sigma(2,1);
trans ( i, 8 ) = sigma(2,2);
}
mat = Trans(trans) * Trans (shape);
}