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.
|
|
|
|
//============================================================================
|
|
|
|
|
2020-10-20 19:22:43 +00:00
|
|
|
#include <vtkm/cont/ArrayRangeComputeTemplate.h>
|
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
|
|
|
|
2018-02-14 15:46:12 +00:00
|
|
|
void ThrowArrayRangeComputeFailed()
|
|
|
|
{
|
|
|
|
throw vtkm::cont::ErrorExecution("Failed to run ArrayRangeComputation on any device.");
|
|
|
|
}
|
|
|
|
|
2020-08-17 14:18:24 +00:00
|
|
|
#define VTKM_ARRAY_RANGE_COMPUTE_IMPL_T(T, Storage) \
|
|
|
|
VTKM_CONT \
|
|
|
|
vtkm::cont::ArrayHandle<vtkm::Range> ArrayRangeCompute( \
|
|
|
|
const vtkm::cont::ArrayHandle<T, Storage>& input, vtkm::cont::DeviceAdapterId device) \
|
|
|
|
{ \
|
|
|
|
return detail::ArrayRangeComputeImpl(input, device); \
|
|
|
|
} \
|
2017-02-24 00:12:14 +00:00
|
|
|
struct SwallowSemicolon
|
2020-08-17 14:18:24 +00:00
|
|
|
#define VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC(T, N, Storage) \
|
|
|
|
VTKM_CONT \
|
|
|
|
vtkm::cont::ArrayHandle<vtkm::Range> ArrayRangeCompute( \
|
|
|
|
const vtkm::cont::ArrayHandle<vtkm::Vec<T, N>, Storage>& input, \
|
|
|
|
vtkm::cont::DeviceAdapterId device) \
|
|
|
|
{ \
|
|
|
|
return detail::ArrayRangeComputeImpl(input, device); \
|
|
|
|
} \
|
2017-02-24 00:12:14 +00:00
|
|
|
struct SwallowSemicolon
|
|
|
|
|
2020-11-11 23:15:14 +00:00
|
|
|
#define VTKM_ARRAY_RANGE_COMPUTE_IMPL_ALL_SCALAR_T(Storage) \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_T(vtkm::Int8, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_T(vtkm::UInt8, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_T(vtkm::Int16, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_T(vtkm::UInt16, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_T(vtkm::Int32, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_T(vtkm::UInt32, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_T(vtkm::Int64, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_T(vtkm::UInt64, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_T(vtkm::Float32, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_T(vtkm::Float64, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_T(char, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_T(signed VTKM_UNUSED_INT_TYPE, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_T(unsigned VTKM_UNUSED_INT_TYPE, Storage)
|
2017-02-24 00:12:14 +00:00
|
|
|
|
2020-11-11 23:15:14 +00:00
|
|
|
#define VTKM_ARRAY_RANGE_COMPUTE_IMPL_ALL_VEC(N, Storage) \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC(vtkm::Int8, N, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC(vtkm::UInt8, N, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC(vtkm::Int16, N, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC(vtkm::UInt16, N, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC(vtkm::Int32, N, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC(vtkm::UInt32, N, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC(vtkm::Int64, N, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC(vtkm::UInt64, N, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC(vtkm::Float32, N, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC(vtkm::Float64, N, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC(char, N, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC(signed VTKM_UNUSED_INT_TYPE, N, Storage); \
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC(unsigned VTKM_UNUSED_INT_TYPE, N, Storage)
|
2017-02-24 00:12:14 +00:00
|
|
|
|
2020-10-20 19:22:43 +00:00
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_ALL_SCALAR_T(vtkm::cont::StorageTagBasic);
|
2017-02-24 00:12:14 +00:00
|
|
|
|
2020-10-20 19:22:43 +00:00
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_ALL_VEC(2, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_ALL_VEC(3, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_ALL_VEC(4, vtkm::cont::StorageTagBasic);
|
2017-02-24 00:12:14 +00:00
|
|
|
|
2020-12-09 22:49:54 +00:00
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_ALL_VEC(2, vtkm::cont::StorageTagSOA);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_ALL_VEC(3, vtkm::cont::StorageTagSOA);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_IMPL_ALL_VEC(4, vtkm::cont::StorageTagSOA);
|
|
|
|
|
2017-02-24 00:12:14 +00:00
|
|
|
#undef VTKM_ARRAY_RANGE_COMPUTE_IMPL_T
|
|
|
|
#undef VTKM_ARRAY_RANGE_COMPUTE_IMPL_VEC
|
2020-10-20 19:22:43 +00:00
|
|
|
#undef VTKM_ARRAY_RANGE_COMPUTE_IMPL_ALL_SCALAR_T
|
|
|
|
#undef VTKM_ARRAY_RANGE_COMPUTE_IMPL_ALL_VEC
|
2017-02-24 00:12:14 +00:00
|
|
|
|
|
|
|
// Special implementation for regular point coordinates, which are easy
|
|
|
|
// to determine.
|
|
|
|
VTKM_CONT
|
2017-05-18 14:29:41 +00:00
|
|
|
vtkm::cont::ArrayHandle<vtkm::Range> ArrayRangeCompute(
|
2019-07-31 16:20:38 +00:00
|
|
|
const vtkm::cont::ArrayHandle<vtkm::Vec3f,
|
2017-05-18 14:29:41 +00:00
|
|
|
vtkm::cont::ArrayHandleUniformPointCoordinates::StorageTag>& array,
|
2019-04-01 21:25:44 +00:00
|
|
|
vtkm::cont::DeviceAdapterId)
|
2017-02-24 00:12:14 +00:00
|
|
|
{
|
2020-01-28 19:14:32 +00:00
|
|
|
vtkm::internal::ArrayPortalUniformPointCoordinates portal = array.ReadPortal();
|
2017-02-24 00:12:14 +00:00
|
|
|
|
|
|
|
// In this portal we know that the min value is the first entry and the
|
|
|
|
// max value is the last entry.
|
2019-07-31 16:20:38 +00:00
|
|
|
vtkm::Vec3f minimum = portal.Get(0);
|
|
|
|
vtkm::Vec3f maximum = portal.Get(portal.GetNumberOfValues() - 1);
|
2017-02-24 00:12:14 +00:00
|
|
|
|
|
|
|
vtkm::cont::ArrayHandle<vtkm::Range> rangeArray;
|
|
|
|
rangeArray.Allocate(3);
|
2020-01-28 19:14:32 +00:00
|
|
|
vtkm::cont::ArrayHandle<vtkm::Range>::WritePortalType outPortal = rangeArray.WritePortal();
|
2017-02-24 00:12:14 +00:00
|
|
|
outPortal.Set(0, vtkm::Range(minimum[0], maximum[0]));
|
|
|
|
outPortal.Set(1, vtkm::Range(minimum[1], maximum[1]));
|
|
|
|
outPortal.Set(2, vtkm::Range(minimum[2], maximum[2]));
|
|
|
|
|
|
|
|
return rangeArray;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
} // namespace vtkm::cont
|