Shuenhoy/vtk-m
1
0
Fork 0
mirror of https://gitlab.kitware.com/vtk/vtk-m synced 2024-09-19 18:45:43 +00:00
Code Issues 2 Actions Packages Projects Releases Wiki Activity

Fix bad CopyIf call

Browse Source
This commit is contained in:
Oliver Ruebel 2021-05-15 02:00:36 -07:00 committed by Gunther H. Weber
parent ca6f0bba9c
commit bb853befdd
1 changed files with 11 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
vtkm/worklet/contourtree_distributed/HierarchicalAugmenter.h
View File

@ -260,17 +260,21 @@ void HierarchicalAugmenter<FieldType>::Initialize(
vtkm::worklet::contourtree_distributed::hierarchical_augmenter::IsAttachementPointPredicate vtkm::worklet::contourtree_distributed::hierarchical_augmenter::IsAttachementPointPredicate
isAttachementPointPredicate( isAttachementPointPredicate(
this->BaseTree->Superarcs, this->BaseTree->WhichRound, this->BaseTree->NumRounds); this->BaseTree->Superarcs, this->BaseTree->WhichRound, this->BaseTree->NumRounds);
auto tempSupernodeIndex =
vtkm::cont::ArrayHandleIndex(this->BaseTree->Supernodes.GetNumberOfValues());
vtkm::cont::Algorithm::CopyIf( vtkm::cont::Algorithm::CopyIf(
// first we generate a list of all of the supernodes // first a list of all of the supernodes
vtkm::cont::ArrayHandleIndex(this->BaseTree->Supernodes.GetNumberOfValues()), tempSupernodeIndex,
// then our stencil identifies all attachment points // then our stencil
tempSupernodeIndex,
// And the CopyIf compress the supernodes array to eliminate the non-attachement points and
// save to this->AttachmentIds
this->AttachmentIds,
// then our predicate identifies all attachment points
// i.e., an attachment point is defined by having no superarc (NO_SUCH_ELEMENT) and not // i.e., an attachment point is defined by having no superarc (NO_SUCH_ELEMENT) and not
// being in the final round (where this indicates the global root) defined by the condition // being in the final round (where this indicates the global root) defined by the condition
// if (noSuchElement(baseTree->superarcs[supernode]) && (baseTree->whichRound[supernode] < baseTree->nRounds)) // if (noSuchElement(baseTree->superarcs[supernode]) && (baseTree->whichRound[supernode] < baseTree->nRounds))
isAttachementPointPredicate, isAttachementPointPredicate);
// And the CopyIf compress the supernodes array to eliminate the non-attachement points and
// save to this->AttachmentIds
this->AttachmentIds);
} }
// we now resize the working arrays // we now resize the working arrays