mirror of
https://gitlab.kitware.com/vtk/vtk-m
synced 2024-10-06 02:18:58 +00:00
180 lines
5.7 KiB
C++
180 lines
5.7 KiB
C++
|
//============================================================================
|
||
|
// Copyright (c) Kitware, Inc.
|
||
|
// All rights reserved.
|
||
|
// See LICENSE.txt for details.
|
||
|
// This software is distributed WITHOUT ANY WARRANTY; without even
|
||
|
// the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
|
||
|
// PURPOSE. See the above copyright notice for more information.
|
||
|
//
|
||
|
// Copyright 2016 Sandia Corporation.
|
||
|
// Copyright 2016 UT-Battelle, LLC.
|
||
|
// Copyright 2016 Los Alamos National Security.
|
||
|
//
|
||
|
// Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
|
||
|
// the U.S. Government retains certain rights in this software.
|
||
|
//
|
||
|
// Under the terms of Contract DE-AC52-06NA25396 with Los Alamos National
|
||
|
// Laboratory (LANL), the U.S. Government retains certain rights in
|
||
|
// this software.
|
||
|
//============================================================================
|
||
|
|
||
|
#include <vtkm/filter/CleanGrid.h>
|
||
|
|
||
|
#include <vtkm/worklet/CellDeepCopy.h>
|
||
|
#include <vtkm/worklet/RemoveUnusedPoints.h>
|
||
|
|
||
|
#include <vector>
|
||
|
|
||
|
namespace vtkm {
|
||
|
namespace filter {
|
||
|
|
||
|
namespace detail {
|
||
|
|
||
|
template<typename Device>
|
||
|
struct CleanCompactPointArrayFunctor
|
||
|
{
|
||
|
vtkm::cont::DataSet &OutDataSet;
|
||
|
std::string Name;
|
||
|
const vtkm::filter::CleanGrid *Self;
|
||
|
|
||
|
CleanCompactPointArrayFunctor(vtkm::cont::DataSet &outDataSet,
|
||
|
const std::string &name,
|
||
|
const vtkm::filter::CleanGrid *self)
|
||
|
: OutDataSet(outDataSet), Name(name), Self(self)
|
||
|
{ }
|
||
|
|
||
|
template<typename ArrayHandleType>
|
||
|
void operator()(const ArrayHandleType &coordSystemArray) const
|
||
|
{
|
||
|
VTKM_IS_ARRAY_HANDLE(ArrayHandleType);
|
||
|
|
||
|
vtkm::cont::ArrayHandle<typename ArrayHandleType::ValueType> outArray =
|
||
|
this->Self->CompactPointArray(coordSystemArray, Device());
|
||
|
this->OutDataSet.AddCoordinateSystem(
|
||
|
vtkm::cont::CoordinateSystem(this->Name, outArray));
|
||
|
}
|
||
|
};
|
||
|
|
||
|
} // namespace detail
|
||
|
|
||
|
VTKM_CONT
|
||
|
CleanGrid::CleanGrid()
|
||
|
: CompactPointFields(true)
|
||
|
{ }
|
||
|
|
||
|
template<typename Policy, typename Device>
|
||
|
VTKM_CONT
|
||
|
vtkm::filter::ResultDataSet
|
||
|
CleanGrid::DoExecute(const vtkm::cont::DataSet &inData,
|
||
|
vtkm::filter::PolicyBase<Policy>,
|
||
|
Device)
|
||
|
{
|
||
|
VTKM_IS_DEVICE_ADAPTER_TAG(Device);
|
||
|
|
||
|
using CellSetType = vtkm::cont::CellSetExplicit<>;
|
||
|
using VecId = std::vector<CellSetType>::size_type;
|
||
|
|
||
|
VecId numCellSets = static_cast<VecId>(inData.GetNumberOfCellSets());
|
||
|
|
||
|
std::vector<CellSetType> outputCellSets(numCellSets);
|
||
|
|
||
|
// Do a deep copy of the cells to new CellSetExplicit structures
|
||
|
for (VecId cellSetIndex = 0; cellSetIndex < numCellSets; cellSetIndex++)
|
||
|
{
|
||
|
vtkm::cont::DynamicCellSet inCellSet =
|
||
|
inData.GetCellSet(static_cast<vtkm::IdComponent>(cellSetIndex));
|
||
|
|
||
|
vtkm::worklet::CellDeepCopy::Run(
|
||
|
vtkm::filter::ApplyPolicy(inCellSet, Policy()),
|
||
|
outputCellSets[cellSetIndex],
|
||
|
Device());
|
||
|
}
|
||
|
|
||
|
// Optionally adjust the cell set indices to remove all unused points
|
||
|
if (this->GetCompactPointFields())
|
||
|
{
|
||
|
this->PointCompactor.FindPointsStart(Device());
|
||
|
for (VecId cellSetIndex = 0; cellSetIndex < numCellSets; cellSetIndex++)
|
||
|
{
|
||
|
this->PointCompactor.FindPoints(outputCellSets[cellSetIndex], Device());
|
||
|
}
|
||
|
this->PointCompactor.FindPointsEnd(Device());
|
||
|
|
||
|
for (VecId cellSetIndex = 0; cellSetIndex < numCellSets; cellSetIndex++)
|
||
|
{
|
||
|
outputCellSets[cellSetIndex] =
|
||
|
this->PointCompactor.MapCellSet(outputCellSets[cellSetIndex], Device());
|
||
|
}
|
||
|
}
|
||
|
|
||
|
// Construct resulting data set with new cell sets
|
||
|
vtkm::cont::DataSet outData;
|
||
|
for (VecId cellSetIndex = 0; cellSetIndex < numCellSets; cellSetIndex++)
|
||
|
{
|
||
|
outData.AddCellSet(outputCellSets[cellSetIndex]);
|
||
|
}
|
||
|
|
||
|
// Pass the coordinate systems
|
||
|
// TODO: This is very awkward. First of all, there is no support for dealing
|
||
|
// with coordinate systems at all. That is fine if you are computing a new
|
||
|
// coordinate system, but a pain if you are deriving the coordinate system
|
||
|
// array. Second, why is it that coordinate systems are automtically mapped
|
||
|
// but other fields are not? Why shouldn't the Execute of a filter also set
|
||
|
// up all the fields of the output data set?
|
||
|
for (vtkm::IdComponent coordSystemIndex = 0;
|
||
|
coordSystemIndex < inData.GetNumberOfCoordinateSystems();
|
||
|
coordSystemIndex++)
|
||
|
{
|
||
|
vtkm::cont::CoordinateSystem coordSystem =
|
||
|
inData.GetCoordinateSystem(coordSystemIndex);
|
||
|
|
||
|
vtkm::filter::ApplyPolicy(coordSystem,
|
||
|
Policy(),
|
||
|
vtkm::filter::FilterTraits<CleanGrid>())
|
||
|
.CastAndCall(detail::CleanCompactPointArrayFunctor<Device>(
|
||
|
outData,coordSystem.GetName(),this));
|
||
|
}
|
||
|
|
||
|
return outData;
|
||
|
}
|
||
|
|
||
|
template<typename ValueType,
|
||
|
typename Storage,
|
||
|
typename Policy,
|
||
|
typename Device>
|
||
|
VTKM_CONT
|
||
|
bool CleanGrid::DoMapField(
|
||
|
vtkm::filter::ResultDataSet& result,
|
||
|
const vtkm::cont::ArrayHandle<ValueType, Storage>& input,
|
||
|
const vtkm::filter::FieldMetadata& fieldMeta,
|
||
|
vtkm::filter::PolicyBase<Policy>,
|
||
|
Device)
|
||
|
{
|
||
|
if (this->GetCompactPointFields() && fieldMeta.IsPointField())
|
||
|
{
|
||
|
vtkm::cont::ArrayHandle<ValueType> compactedArray =
|
||
|
this->CompactPointArray(input, Device());
|
||
|
result.GetDataSet().AddField(fieldMeta.AsField(compactedArray));
|
||
|
}
|
||
|
else
|
||
|
{
|
||
|
result.GetDataSet().AddField(fieldMeta.AsField(input));
|
||
|
}
|
||
|
|
||
|
return true;
|
||
|
}
|
||
|
|
||
|
template<typename ValueType, typename Storage, typename Device>
|
||
|
VTKM_CONT
|
||
|
vtkm::cont::ArrayHandle<ValueType>
|
||
|
CleanGrid::CompactPointArray(
|
||
|
const vtkm::cont::ArrayHandle<ValueType,Storage> &inArray, Device) const
|
||
|
{
|
||
|
VTKM_ASSERT(this->GetCompactPointFields());
|
||
|
|
||
|
return this->PointCompactor.MapPointFieldDeep(inArray, Device());
|
||
|
}
|
||
|
|
||
|
}
|
||
|
}
|