本文整理汇总了C++中SparseMatrix::MultTranspose方法的典型用法代码示例。如果您正苦于以下问题:C++ SparseMatrix::MultTranspose方法的具体用法?C++ SparseMatrix::MultTranspose怎么用?C++ SparseMatrix::MultTranspose使用的例子?那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类SparseMatrix
的用法示例。
在下文中一共展示了SparseMatrix::MultTranspose方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: UpdateAndInterpolate
void FiniteElementSpace::UpdateAndInterpolate(int num_grid_fns, ...)
{
if (mesh->GetState() == Mesh::NORMAL)
{
MFEM_ABORT("Mesh must be in two-level state, please call "
"Mesh::UseTwoLevelState before refining.");
}
FiniteElementSpace *cfes = SaveUpdate();
// obtain the (transpose of) interpolation matrix between mesh levels
SparseMatrix *R = GlobalRestrictionMatrix(cfes, 0);
delete cfes;
// interpolate the grid functions
std::va_list vl;
va_start(vl, num_grid_fns);
for (int i = 0; i < num_grid_fns; i++)
{
GridFunction* gf = va_arg(vl, GridFunction*);
if (gf->FESpace() != this)
{
MFEM_ABORT("Cannot interpolate: grid function is not based "
"on this space.");
}
Vector coarse_gf = *gf;
gf->Update();
R->MultTranspose(coarse_gf, *gf);
}
va_end(vl);
delete R;
mesh->SetState(Mesh::TWO_LEVEL_FINE);
}