2017-02-24 00:12:14 +00:00
|
|
|
//============================================================================
|
|
|
|
// Copyright (c) Kitware, Inc.
|
|
|
|
// All rights reserved.
|
|
|
|
// See LICENSE.txt for details.
|
2019-04-15 23:24:21 +00:00
|
|
|
//
|
2017-02-24 00:12:14 +00:00
|
|
|
// 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.
|
|
|
|
//============================================================================
|
|
|
|
#ifndef vtk_m_cont_ArrayRangeCompute_hxx
|
|
|
|
#define vtk_m_cont_ArrayRangeCompute_hxx
|
|
|
|
|
|
|
|
#include <vtkm/cont/ArrayRangeCompute.h>
|
|
|
|
|
|
|
|
#include <vtkm/BinaryOperators.h>
|
|
|
|
#include <vtkm/VecTraits.h>
|
|
|
|
|
2018-02-14 15:46:12 +00:00
|
|
|
#include <vtkm/cont/Algorithm.h>
|
|
|
|
|
|
|
|
#include <limits>
|
2017-02-24 00:12:14 +00:00
|
|
|
|
2017-05-18 14:29:41 +00:00
|
|
|
namespace vtkm
|
|
|
|
{
|
|
|
|
namespace cont
|
|
|
|
{
|
2017-02-24 00:12:14 +00:00
|
|
|
|
2017-05-18 14:29:41 +00:00
|
|
|
namespace detail
|
|
|
|
{
|
2017-02-24 00:12:14 +00:00
|
|
|
|
|
|
|
struct ArrayRangeComputeFunctor
|
|
|
|
{
|
2018-02-14 15:46:12 +00:00
|
|
|
template <typename Device, typename T, typename S>
|
2018-02-14 14:25:02 +00:00
|
|
|
VTKM_CONT bool operator()(Device,
|
2018-02-14 15:46:12 +00:00
|
|
|
const vtkm::cont::ArrayHandle<T, S>& handle,
|
|
|
|
const vtkm::Vec<T, 2>& initialValue,
|
|
|
|
vtkm::Vec<T, 2>& result) const
|
2017-02-24 00:12:14 +00:00
|
|
|
{
|
|
|
|
VTKM_IS_DEVICE_ADAPTER_TAG(Device);
|
2017-08-16 15:34:21 +00:00
|
|
|
using Algorithm = vtkm::cont::DeviceAdapterAlgorithm<Device>;
|
2018-02-14 15:46:12 +00:00
|
|
|
result = Algorithm::Reduce(handle, initialValue, vtkm::MinAndMax<T>());
|
2017-02-24 00:12:14 +00:00
|
|
|
return true;
|
|
|
|
}
|
|
|
|
};
|
|
|
|
|
2018-02-14 15:46:12 +00:00
|
|
|
template <typename T, typename S>
|
2017-05-18 14:29:41 +00:00
|
|
|
inline vtkm::cont::ArrayHandle<vtkm::Range> ArrayRangeComputeImpl(
|
2018-02-14 15:46:12 +00:00
|
|
|
const vtkm::cont::ArrayHandle<T, S>& input,
|
2019-04-01 21:25:44 +00:00
|
|
|
vtkm::cont::DeviceAdapterId device)
|
2017-02-24 00:12:14 +00:00
|
|
|
{
|
2018-02-14 15:46:12 +00:00
|
|
|
using VecTraits = vtkm::VecTraits<T>;
|
|
|
|
using CT = typename VecTraits::ComponentType;
|
|
|
|
//We want to minimize the amount of code that we do in try execute as
|
|
|
|
//it is repeated for each
|
2018-02-14 14:25:02 +00:00
|
|
|
vtkm::cont::ArrayHandle<vtkm::Range> range;
|
2018-02-14 15:46:12 +00:00
|
|
|
range.Allocate(VecTraits::NUM_COMPONENTS);
|
2017-02-24 00:12:14 +00:00
|
|
|
|
2018-02-14 15:46:12 +00:00
|
|
|
if (input.GetNumberOfValues() < 1)
|
2017-02-24 00:12:14 +00:00
|
|
|
{
|
2018-02-14 15:46:12 +00:00
|
|
|
auto portal = range.GetPortalControl();
|
|
|
|
for (vtkm::IdComponent i = 0; i < VecTraits::NUM_COMPONENTS; ++i)
|
|
|
|
{
|
|
|
|
portal.Set(i, vtkm::Range());
|
|
|
|
}
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
//We used the limits, so that we don't need to sync the array handle
|
|
|
|
//
|
|
|
|
vtkm::Vec<T, 2> result;
|
|
|
|
vtkm::Vec<T, 2> initial;
|
|
|
|
initial[0] = T(std::numeric_limits<CT>::max());
|
|
|
|
initial[1] = T(std::numeric_limits<CT>::lowest());
|
|
|
|
|
2019-04-01 21:25:44 +00:00
|
|
|
const bool rangeComputed = vtkm::cont::TryExecuteOnDevice(
|
|
|
|
device, detail::ArrayRangeComputeFunctor{}, input, initial, result);
|
2018-02-14 15:46:12 +00:00
|
|
|
if (!rangeComputed)
|
|
|
|
{
|
|
|
|
ThrowArrayRangeComputeFailed();
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
auto portal = range.GetPortalControl();
|
|
|
|
for (vtkm::IdComponent i = 0; i < VecTraits::NUM_COMPONENTS; ++i)
|
|
|
|
{
|
|
|
|
portal.Set(i,
|
|
|
|
vtkm::Range(VecTraits::GetComponent(result[0], i),
|
|
|
|
VecTraits::GetComponent(result[1], i)));
|
|
|
|
}
|
|
|
|
}
|
2017-02-24 00:12:14 +00:00
|
|
|
}
|
2018-02-14 14:25:02 +00:00
|
|
|
return range;
|
2017-02-24 00:12:14 +00:00
|
|
|
}
|
|
|
|
|
|
|
|
} // namespace detail
|
|
|
|
|
2018-12-11 19:12:47 +00:00
|
|
|
|
|
|
|
VTKM_CONT
|
|
|
|
inline vtkm::cont::ArrayHandle<vtkm::Range> ArrayRangeCompute(
|
2019-07-31 16:20:38 +00:00
|
|
|
const vtkm::cont::ArrayHandleVirtual<vtkm::Vec3f>& input,
|
2019-04-01 21:25:44 +00:00
|
|
|
vtkm::cont::DeviceAdapterId device)
|
2018-12-11 19:12:47 +00:00
|
|
|
{
|
|
|
|
using UniformHandleType = ArrayHandleUniformPointCoordinates;
|
|
|
|
using RectilinearHandleType =
|
|
|
|
vtkm::cont::ArrayHandleCartesianProduct<vtkm::cont::ArrayHandle<vtkm::FloatDefault>,
|
|
|
|
vtkm::cont::ArrayHandle<vtkm::FloatDefault>,
|
|
|
|
vtkm::cont::ArrayHandle<vtkm::FloatDefault>>;
|
|
|
|
|
|
|
|
if (input.IsType<UniformHandleType>())
|
|
|
|
{
|
|
|
|
using T = typename UniformHandleType::ValueType;
|
|
|
|
using S = typename UniformHandleType::StorageTag;
|
2019-02-10 08:46:43 +00:00
|
|
|
const vtkm::cont::internal::detail::StorageVirtual* storage =
|
|
|
|
input.GetStorage().GetStorageVirtual();
|
|
|
|
const auto* castStorage =
|
|
|
|
storage->Cast<vtkm::cont::internal::detail::StorageVirtualImpl<T, S>>();
|
2018-12-11 19:12:47 +00:00
|
|
|
|
2019-04-01 21:25:44 +00:00
|
|
|
return ArrayRangeCompute(castStorage->GetHandle(), device);
|
2018-12-11 19:12:47 +00:00
|
|
|
}
|
|
|
|
else if (input.IsType<RectilinearHandleType>())
|
|
|
|
{
|
|
|
|
using T = typename RectilinearHandleType::ValueType;
|
|
|
|
using S = typename RectilinearHandleType::StorageTag;
|
2019-02-10 08:46:43 +00:00
|
|
|
const vtkm::cont::internal::detail::StorageVirtual* storage =
|
|
|
|
input.GetStorage().GetStorageVirtual();
|
|
|
|
const auto* castStorage =
|
|
|
|
storage->Cast<vtkm::cont::internal::detail::StorageVirtualImpl<T, S>>();
|
2018-12-11 19:12:47 +00:00
|
|
|
|
2019-04-01 21:25:44 +00:00
|
|
|
return ArrayRangeCompute(castStorage->GetHandle(), device);
|
2018-12-11 19:12:47 +00:00
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
2019-04-01 21:25:44 +00:00
|
|
|
return detail::ArrayRangeComputeImpl(input, device);
|
2018-12-11 19:12:47 +00:00
|
|
|
}
|
|
|
|
}
|
|
|
|
|
2017-05-18 14:29:41 +00:00
|
|
|
template <typename ArrayHandleType>
|
2019-04-01 21:25:44 +00:00
|
|
|
inline vtkm::cont::ArrayHandle<vtkm::Range> ArrayRangeCompute(const ArrayHandleType& input,
|
|
|
|
vtkm::cont::DeviceAdapterId device)
|
2017-02-24 00:12:14 +00:00
|
|
|
{
|
|
|
|
VTKM_IS_ARRAY_HANDLE(ArrayHandleType);
|
2019-04-01 21:25:44 +00:00
|
|
|
return detail::ArrayRangeComputeImpl(input, device);
|
2017-02-24 00:12:14 +00:00
|
|
|
}
|
|
|
|
}
|
|
|
|
} // namespace vtkm::cont
|
|
|
|
|
|
|
|
#endif //vtk_m_cont_ArrayRangeCompute_hxx
|