2016-08-05 20:05:32 +00:00
|
|
|
//============================================================================
|
|
|
|
// 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.
|
|
|
|
//
|
2017-09-20 21:33:44 +00:00
|
|
|
// Copyright 2016 National Technology & Engineering Solutions of Sandia, LLC (NTESS).
|
2016-08-05 20:05:32 +00:00
|
|
|
// Copyright 2016 UT-Battelle, LLC.
|
|
|
|
// Copyright 2016 Los Alamos National Security.
|
|
|
|
//
|
2017-09-20 21:33:44 +00:00
|
|
|
// Under the terms of Contract DE-NA0003525 with NTESS,
|
2016-08-05 20:05:32 +00:00
|
|
|
// 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.
|
|
|
|
//============================================================================
|
|
|
|
#ifndef vtk_m_cont_ArrayRangeCompute_h
|
|
|
|
#define vtk_m_cont_ArrayRangeCompute_h
|
|
|
|
|
|
|
|
#include <vtkm/Range.h>
|
|
|
|
|
|
|
|
#include <vtkm/cont/ArrayHandle.h>
|
2017-02-24 00:12:14 +00:00
|
|
|
#include <vtkm/cont/ArrayHandleCartesianProduct.h>
|
|
|
|
#include <vtkm/cont/ArrayHandleCompositeVector.h>
|
2016-08-05 20:05:32 +00:00
|
|
|
#include <vtkm/cont/ArrayHandleUniformPointCoordinates.h>
|
2017-12-12 19:26:23 +00:00
|
|
|
#include <vtkm/cont/ArrayHandleVirtualCoordinates.h>
|
2017-02-24 00:12:14 +00:00
|
|
|
#include <vtkm/cont/RuntimeDeviceTracker.h>
|
2016-08-05 20:05:32 +00:00
|
|
|
|
2017-05-18 14:29:41 +00:00
|
|
|
namespace vtkm
|
|
|
|
{
|
|
|
|
namespace cont
|
|
|
|
{
|
2016-08-05 20:05:32 +00:00
|
|
|
|
|
|
|
/// \brief Compute the range of the data in an array handle.
|
|
|
|
///
|
|
|
|
/// Given an \c ArrayHandle, this function computes the range (min and max) of
|
|
|
|
/// the values in the array. For arrays containing Vec values, the range is
|
|
|
|
/// computed for each component.
|
|
|
|
///
|
2017-09-25 22:28:44 +00:00
|
|
|
/// This method optionally takes a \c RuntimeDeviceTracker to control which
|
|
|
|
/// devices to try.
|
2016-08-05 20:05:32 +00:00
|
|
|
///
|
|
|
|
/// The result is returned in an \c ArrayHandle of \c Range objects. There is
|
|
|
|
/// one value in the returned array for every component of the input's value
|
|
|
|
/// type.
|
|
|
|
///
|
2017-05-18 14:29:41 +00:00
|
|
|
template <typename ArrayHandleType>
|
|
|
|
VTKM_CONT vtkm::cont::ArrayHandle<vtkm::Range> ArrayRangeCompute(
|
|
|
|
const ArrayHandleType& input,
|
|
|
|
vtkm::cont::RuntimeDeviceTracker tracker = vtkm::cont::GetGlobalRuntimeDeviceTracker());
|
2017-02-24 00:12:14 +00:00
|
|
|
|
2018-11-08 15:16:59 +00:00
|
|
|
|
2017-02-24 00:12:14 +00:00
|
|
|
// Precompiled versions of ArrayRangeCompute
|
2017-05-18 14:29:41 +00:00
|
|
|
#define VTKM_ARRAY_RANGE_COMPUTE_EXPORT_T(T, Storage) \
|
|
|
|
VTKM_CONT_EXPORT \
|
|
|
|
VTKM_CONT \
|
|
|
|
vtkm::cont::ArrayHandle<vtkm::Range> ArrayRangeCompute( \
|
|
|
|
const vtkm::cont::ArrayHandle<T, Storage>& input, \
|
|
|
|
vtkm::cont::RuntimeDeviceTracker tracker = vtkm::cont::GetGlobalRuntimeDeviceTracker())
|
|
|
|
#define VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC(T, N, Storage) \
|
|
|
|
VTKM_CONT_EXPORT \
|
|
|
|
VTKM_CONT \
|
|
|
|
vtkm::cont::ArrayHandle<vtkm::Range> ArrayRangeCompute( \
|
|
|
|
const vtkm::cont::ArrayHandle<vtkm::Vec<T, N>, Storage>& input, \
|
|
|
|
vtkm::cont::RuntimeDeviceTracker tracker = vtkm::cont::GetGlobalRuntimeDeviceTracker())
|
2017-02-24 00:12:14 +00:00
|
|
|
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_T(char, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_T(vtkm::Int8, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_T(vtkm::UInt8, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_T(vtkm::Int16, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_T(vtkm::UInt16, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_T(vtkm::Int32, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_T(vtkm::UInt32, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_T(vtkm::Int64, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_T(vtkm::UInt64, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_T(vtkm::Float32, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_T(vtkm::Float64, vtkm::cont::StorageTagBasic);
|
|
|
|
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC(vtkm::Int32, 2, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC(vtkm::Int64, 2, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC(vtkm::Float32, 2, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC(vtkm::Float64, 2, vtkm::cont::StorageTagBasic);
|
|
|
|
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC(vtkm::Int32, 3, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC(vtkm::Int64, 3, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC(vtkm::Float32, 3, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC(vtkm::Float64, 3, vtkm::cont::StorageTagBasic);
|
|
|
|
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC(char, 4, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC(vtkm::Int8, 4, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC(vtkm::UInt8, 4, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC(vtkm::Float32, 4, vtkm::cont::StorageTagBasic);
|
|
|
|
VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC(vtkm::Float64, 4, vtkm::cont::StorageTagBasic);
|
|
|
|
|
2017-12-12 19:26:23 +00:00
|
|
|
|
2017-02-24 00:12:14 +00:00
|
|
|
#undef VTKM_ARRAY_RANGE_COMPUTE_EXPORT_T
|
|
|
|
#undef VTKM_ARRAY_RANGE_COMPUTE_EXPORT_VEC
|
|
|
|
|
2018-11-08 15:16:59 +00:00
|
|
|
// Implementation of uniform point coordinates
|
|
|
|
VTKM_CONT_EXPORT
|
|
|
|
VTKM_CONT
|
|
|
|
vtkm::cont::ArrayHandle<vtkm::Range> ArrayRangeCompute(
|
|
|
|
const vtkm::cont::ArrayHandle<vtkm::Vec<vtkm::FloatDefault, 3>,
|
|
|
|
vtkm::cont::ArrayHandleUniformPointCoordinates::StorageTag>& array,
|
|
|
|
vtkm::cont::RuntimeDeviceTracker = vtkm::cont::GetGlobalRuntimeDeviceTracker());
|
|
|
|
|
|
|
|
|
2017-02-24 00:12:14 +00:00
|
|
|
// Implementation of composite vectors
|
|
|
|
VTKM_CONT_EXPORT
|
|
|
|
VTKM_CONT
|
2017-05-18 14:29:41 +00:00
|
|
|
vtkm::cont::ArrayHandle<vtkm::Range> ArrayRangeCompute(
|
2018-05-10 17:06:08 +00:00
|
|
|
const vtkm::cont::ArrayHandle<vtkm::Vec<vtkm::Float32, 3>,
|
|
|
|
typename vtkm::cont::ArrayHandleCompositeVector<
|
|
|
|
vtkm::cont::ArrayHandle<vtkm::Float32>,
|
2017-05-26 17:53:28 +00:00
|
|
|
vtkm::cont::ArrayHandle<vtkm::Float32>,
|
2018-05-10 17:06:08 +00:00
|
|
|
vtkm::cont::ArrayHandle<vtkm::Float32>>::StorageTag>& input,
|
2017-05-18 14:29:41 +00:00
|
|
|
vtkm::cont::RuntimeDeviceTracker tracker = vtkm::cont::GetGlobalRuntimeDeviceTracker());
|
2017-02-24 00:12:14 +00:00
|
|
|
|
|
|
|
VTKM_CONT_EXPORT
|
|
|
|
VTKM_CONT
|
2017-05-18 14:29:41 +00:00
|
|
|
vtkm::cont::ArrayHandle<vtkm::Range> ArrayRangeCompute(
|
2018-05-10 17:06:08 +00:00
|
|
|
const vtkm::cont::ArrayHandle<vtkm::Vec<vtkm::Float64, 3>,
|
|
|
|
typename vtkm::cont::ArrayHandleCompositeVector<
|
|
|
|
vtkm::cont::ArrayHandle<vtkm::Float64>,
|
2017-05-26 17:53:28 +00:00
|
|
|
vtkm::cont::ArrayHandle<vtkm::Float64>,
|
2018-05-10 17:06:08 +00:00
|
|
|
vtkm::cont::ArrayHandle<vtkm::Float64>>::StorageTag>& input,
|
2017-05-18 14:29:41 +00:00
|
|
|
vtkm::cont::RuntimeDeviceTracker tracker = vtkm::cont::GetGlobalRuntimeDeviceTracker());
|
2017-02-24 00:12:14 +00:00
|
|
|
|
|
|
|
// Implementation of cartesian products
|
2017-05-18 14:29:41 +00:00
|
|
|
template <typename T, typename ArrayType1, typename ArrayType2, typename ArrayType3>
|
|
|
|
VTKM_CONT inline vtkm::cont::ArrayHandle<vtkm::Range> ArrayRangeCompute(
|
|
|
|
const vtkm::cont::ArrayHandle<
|
2017-05-26 17:53:28 +00:00
|
|
|
T,
|
|
|
|
vtkm::cont::internal::StorageTagCartesianProduct<ArrayType1, ArrayType2, ArrayType3>>& input,
|
2017-05-18 14:29:41 +00:00
|
|
|
vtkm::cont::RuntimeDeviceTracker tracker = vtkm::cont::GetGlobalRuntimeDeviceTracker())
|
2016-08-05 20:05:32 +00:00
|
|
|
{
|
2017-02-24 00:12:14 +00:00
|
|
|
vtkm::cont::ArrayHandle<vtkm::Range> result;
|
|
|
|
result.Allocate(3);
|
2016-08-05 20:05:32 +00:00
|
|
|
|
2017-02-24 00:12:14 +00:00
|
|
|
vtkm::cont::ArrayHandle<vtkm::Range> componentRangeArray;
|
|
|
|
vtkm::Range componentRange;
|
2016-08-05 20:05:32 +00:00
|
|
|
|
2017-02-24 00:12:14 +00:00
|
|
|
ArrayType1 firstArray = input.GetStorage().GetFirstArray();
|
|
|
|
componentRangeArray = vtkm::cont::ArrayRangeCompute(firstArray, tracker);
|
|
|
|
componentRange = componentRangeArray.GetPortalConstControl().Get(0);
|
|
|
|
result.GetPortalControl().Set(0, componentRange);
|
2016-08-05 20:05:32 +00:00
|
|
|
|
2017-02-24 00:12:14 +00:00
|
|
|
ArrayType2 secondArray = input.GetStorage().GetSecondArray();
|
|
|
|
componentRangeArray = vtkm::cont::ArrayRangeCompute(secondArray, tracker);
|
|
|
|
componentRange = componentRangeArray.GetPortalConstControl().Get(0);
|
|
|
|
result.GetPortalControl().Set(1, componentRange);
|
2016-08-05 20:05:32 +00:00
|
|
|
|
2017-02-24 00:12:14 +00:00
|
|
|
ArrayType3 thirdArray = input.GetStorage().GetThirdArray();
|
|
|
|
componentRangeArray = vtkm::cont::ArrayRangeCompute(thirdArray, tracker);
|
|
|
|
componentRange = componentRangeArray.GetPortalConstControl().Get(0);
|
|
|
|
result.GetPortalControl().Set(2, componentRange);
|
2016-08-05 20:05:32 +00:00
|
|
|
|
2017-02-24 00:12:14 +00:00
|
|
|
return result;
|
2016-08-05 20:05:32 +00:00
|
|
|
}
|
2018-02-14 15:46:12 +00:00
|
|
|
|
|
|
|
VTKM_CONT_EXPORT void ThrowArrayRangeComputeFailed();
|
2016-08-05 20:05:32 +00:00
|
|
|
}
|
|
|
|
} // namespace vtkm::cont
|
|
|
|
|
|
|
|
#endif //vtk_m_cont_ArrayRangeCompute_h
|