本文整理汇总了C++中Tpetra_Vector::getLocalLength方法的典型用法代码示例。如果您正苦于以下问题：C++ Tpetra_Vector::getLocalLength方法的具体用法？C++ Tpetra_Vector::getLocalLength怎么用？C++ Tpetra_Vector::getLocalLength使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Tpetra_Vector的用法示例。
在下文中一共展示了Tpetra_Vector::getLocalLength方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1:
void
Albany::SolutionMaxValueResponseFunction::
computeMaxValueT(const Tpetra_Vector& xT, double& global_max, int& global_index)
{
//The following is needed b/c Epetra_MaxDouble comes from Trilinos Epetra package.
double Tpetra_MaxDouble = 1.0E+100;
double my_max = -Tpetra_MaxDouble;
int my_index = -1, index;
Teuchos::ArrayRCP<const ST> xT_constView = xT.get1dView();
// Loop over nodes to find max value for equation eq
int num_my_nodes = xT.getLocalLength() / neq;
for (int node=0; node<num_my_nodes; node++) {
if (interleavedOrdering) index = node*neq+eq;
else index = node + eq*num_my_nodes;
if (xT_constView[index] > my_max) {
my_max = xT_constView[index];
my_index = index;
}
}
// Check remainder (AGS: NOT SURE HOW THIS CODE GETS CALLED?)
if (num_my_nodes*neq+eq < xT.getLocalLength()) {
if (interleavedOrdering) index = num_my_nodes*neq+eq;
else index = num_my_nodes + eq*num_my_nodes;
if (xT_constView[index] > my_max) {
my_max = xT_constView[index];
my_index = index;
}
}
Teuchos::RCP<const Teuchos::Comm<int> > commT = xT.getMap()->getComm();
// Get max value across all proc's
Teuchos::reduceAll(*commT, Teuchos::REDUCE_MAX, my_max, Teuchos::ptr(&global_max));
// Compute min of all global indices equal to max value
if (my_max == global_max)
my_index = xT.getMap()->getGlobalElement(my_index);
else
my_index = xT.getGlobalLength();
Teuchos::reduceAll(*commT, Teuchos::REDUCE_MIN, my_index, Teuchos::ptr(&global_index));
}