vtk-m2/vtkm/cont/Field.cxx
Kenneth Moreland 7d681fb585 Deprecate templated versions of Field::GetRange
Instead, always use precompiled versions of range computing. This means
you won't be able to specify the type. Currently, types are limited to
scalars vecs up to size 4.

The main motivation for this change is to allow you to include Field.h
with a non-device compiler. This is an important feature for our
customers.

I plan in the future to implement a mechanism to pull out a component of
most ArrayHandle's as a single array. This would enable us to support a
precompiled version that can compute the range of arbitrarily sized
Vecs.
2020-11-09 12:28:29 -07:00

156 lines
3.4 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.
//============================================================================
#include <vtkm/cont/Field.h>
#include <vtkm/TypeList.h>
namespace
{
using ComputeRangeTypes = vtkm::TypeListAll;
struct ComputeRange
{
template <typename ArrayHandleType>
void operator()(const ArrayHandleType& input, vtkm::cont::ArrayHandle<vtkm::Range>& range) const
{
range = vtkm::cont::ArrayRangeCompute(input);
}
};
} // anonymous namespace
namespace vtkm
{
namespace cont
{
/// constructors for points / whole mesh
VTKM_CONT
Field::Field(std::string name, Association association, const vtkm::cont::VariantArrayHandle& data)
: Name(name)
, FieldAssociation(association)
, Data(data)
, Range()
, ModifiedFlag(true)
{
}
VTKM_CONT
Field::Field(const vtkm::cont::Field& src)
: Name(src.Name)
, FieldAssociation(src.FieldAssociation)
, Data(src.Data)
, Range(src.Range)
, ModifiedFlag(src.ModifiedFlag)
{
}
VTKM_CONT
Field::Field(vtkm::cont::Field&& src) noexcept
: Name(std::move(src.Name))
, FieldAssociation(std::move(src.FieldAssociation))
, Data(std::move(src.Data))
, Range(std::move(src.Range))
, ModifiedFlag(std::move(src.ModifiedFlag))
{
}
VTKM_CONT
Field& Field::operator=(const vtkm::cont::Field& src)
{
this->Name = src.Name;
this->FieldAssociation = src.FieldAssociation;
this->Data = src.Data;
this->Range = src.Range;
this->ModifiedFlag = src.ModifiedFlag;
return *this;
}
VTKM_CONT
Field& Field::operator=(vtkm::cont::Field&& src) noexcept
{
this->Name = std::move(src.Name);
this->FieldAssociation = std::move(src.FieldAssociation);
this->Data = std::move(src.Data);
this->Range = std::move(src.Range);
this->ModifiedFlag = std::move(src.ModifiedFlag);
return *this;
}
VTKM_CONT
void Field::PrintSummary(std::ostream& out) const
{
out << " " << this->Name;
out << " assoc= ";
switch (this->GetAssociation())
{
case Association::ANY:
out << "Any ";
break;
case Association::WHOLE_MESH:
out << "Mesh ";
break;
case Association::POINTS:
out << "Points ";
break;
case Association::CELL_SET:
out << "Cells ";
break;
}
this->Data.PrintSummary(out);
}
VTKM_CONT
Field::~Field() {}
VTKM_CONT
const vtkm::cont::VariantArrayHandle& Field::GetData() const
{
return this->Data;
}
VTKM_CONT
vtkm::cont::VariantArrayHandle& Field::GetData()
{
this->ModifiedFlag = true;
return this->Data;
}
VTKM_CONT const vtkm::cont::ArrayHandle<vtkm::Range>& Field::GetRange() const
{
VTKM_LOG_SCOPE(vtkm::cont::LogLevel::Perf, "Field::GetRange");
if (this->ModifiedFlag)
{
vtkm::cont::CastAndCall(
this->Data.ResetTypes(ComputeRangeTypes{}), ComputeRange{}, this->Range);
this->ModifiedFlag = false;
}
return this->Range;
}
VTKM_CONT void Field::GetRange(vtkm::Range* range) const
{
this->GetRange();
const vtkm::Id length = this->Range.GetNumberOfValues();
auto portal = this->Range.ReadPortal();
for (vtkm::Id i = 0; i < length; ++i)
{
range[i] = portal.Get(i);
}
}
}
} // namespace vtkm::cont