Spaces:
Runtime error
Runtime error
////////////////////// | |
// Vector Functions // | |
////////////////////// | |
// convert xxx_vector<T1> to xxx_vector<T2> | |
template <class ExampleVector, typename NewType> | |
struct vector_like | |
{ | |
typedef typename ExampleVector::allocator_type alloc; | |
typedef typename thrust::detail::allocator_traits<alloc> alloc_traits; | |
typedef typename alloc_traits::template rebind_alloc<NewType> new_alloc; | |
typedef thrust::detail::vector_base<NewType, new_alloc> type; | |
}; | |
template <class Vector> | |
void TestVectorLowerBoundDescendingSimple(void) | |
{ | |
typedef typename Vector::value_type T; | |
Vector vec(5); | |
vec[0] = 8; | |
vec[1] = 7; | |
vec[2] = 5; | |
vec[3] = 2; | |
vec[4] = 0; | |
Vector input(10); | |
thrust::sequence(input.begin(), input.end()); | |
typedef typename Vector::difference_type int_type; | |
typedef typename vector_like<Vector, int_type>::type IntVector; | |
// test with integral output type | |
IntVector integral_output(10); | |
typename IntVector::iterator output_end = thrust::lower_bound(vec.begin(), vec.end(), input.begin(), input.end(), integral_output.begin(), thrust::greater<T>()); | |
ASSERT_EQUAL_QUIET(integral_output.end(), output_end); | |
ASSERT_EQUAL(4, integral_output[0]); | |
ASSERT_EQUAL(4, integral_output[1]); | |
ASSERT_EQUAL(3, integral_output[2]); | |
ASSERT_EQUAL(3, integral_output[3]); | |
ASSERT_EQUAL(3, integral_output[4]); | |
ASSERT_EQUAL(2, integral_output[5]); | |
ASSERT_EQUAL(2, integral_output[6]); | |
ASSERT_EQUAL(1, integral_output[7]); | |
ASSERT_EQUAL(0, integral_output[8]); | |
ASSERT_EQUAL(0, integral_output[9]); | |
} | |
DECLARE_VECTOR_UNITTEST(TestVectorLowerBoundDescendingSimple); | |
template <class Vector> | |
void TestVectorUpperBoundDescendingSimple(void) | |
{ | |
Vector vec(5); | |
vec[0] = 8; | |
vec[1] = 7; | |
vec[2] = 5; | |
vec[3] = 2; | |
vec[4] = 0; | |
Vector input(10); | |
thrust::sequence(input.begin(), input.end()); | |
typedef typename Vector::difference_type int_type; | |
typedef typename Vector::value_type T; | |
typedef typename vector_like<Vector, int_type>::type IntVector; | |
// test with integral output type | |
IntVector integral_output(10); | |
typename IntVector::iterator output_end = thrust::upper_bound(vec.begin(), vec.end(), input.begin(), input.end(), integral_output.begin(), thrust::greater<T>()); | |
ASSERT_EQUAL_QUIET(output_end, integral_output.end()); | |
ASSERT_EQUAL(5, integral_output[0]); | |
ASSERT_EQUAL(4, integral_output[1]); | |
ASSERT_EQUAL(4, integral_output[2]); | |
ASSERT_EQUAL(3, integral_output[3]); | |
ASSERT_EQUAL(3, integral_output[4]); | |
ASSERT_EQUAL(3, integral_output[5]); | |
ASSERT_EQUAL(2, integral_output[6]); | |
ASSERT_EQUAL(2, integral_output[7]); | |
ASSERT_EQUAL(1, integral_output[8]); | |
ASSERT_EQUAL(0, integral_output[9]); | |
} | |
DECLARE_VECTOR_UNITTEST(TestVectorUpperBoundDescendingSimple); | |
template <class Vector> | |
void TestVectorBinarySearchDescendingSimple(void) | |
{ | |
Vector vec(5); | |
vec[0] = 8; | |
vec[1] = 7; | |
vec[2] = 5; | |
vec[3] = 2; | |
vec[4] = 0; | |
Vector input(10); | |
thrust::sequence(input.begin(), input.end()); | |
typedef typename vector_like<Vector, bool>::type BoolVector; | |
typedef typename Vector::difference_type int_type; | |
typedef typename Vector::value_type T; | |
typedef typename vector_like<Vector, int_type>::type IntVector; | |
// test with boolean output type | |
BoolVector bool_output(10); | |
typename BoolVector::iterator bool_output_end = thrust::binary_search(vec.begin(), vec.end(), input.begin(), input.end(), bool_output.begin(), thrust::greater<T>()); | |
ASSERT_EQUAL_QUIET(bool_output_end, bool_output.end()); | |
ASSERT_EQUAL(true, bool_output[0]); | |
ASSERT_EQUAL(false, bool_output[1]); | |
ASSERT_EQUAL(true, bool_output[2]); | |
ASSERT_EQUAL(false, bool_output[3]); | |
ASSERT_EQUAL(false, bool_output[4]); | |
ASSERT_EQUAL(true, bool_output[5]); | |
ASSERT_EQUAL(false, bool_output[6]); | |
ASSERT_EQUAL(true, bool_output[7]); | |
ASSERT_EQUAL(true, bool_output[8]); | |
ASSERT_EQUAL(false, bool_output[9]); | |
// test with integral output type | |
IntVector integral_output(10, 2); | |
typename IntVector::iterator int_output_end = thrust::binary_search(vec.begin(), vec.end(), input.begin(), input.end(), integral_output.begin(), thrust::greater<T>()); | |
ASSERT_EQUAL_QUIET(int_output_end, integral_output.end()); | |
ASSERT_EQUAL(1, integral_output[0]); | |
ASSERT_EQUAL(0, integral_output[1]); | |
ASSERT_EQUAL(1, integral_output[2]); | |
ASSERT_EQUAL(0, integral_output[3]); | |
ASSERT_EQUAL(0, integral_output[4]); | |
ASSERT_EQUAL(1, integral_output[5]); | |
ASSERT_EQUAL(0, integral_output[6]); | |
ASSERT_EQUAL(1, integral_output[7]); | |
ASSERT_EQUAL(1, integral_output[8]); | |
ASSERT_EQUAL(0, integral_output[9]); | |
} | |
DECLARE_VECTOR_UNITTEST(TestVectorBinarySearchDescendingSimple); | |
template <typename T> | |
struct TestVectorLowerBoundDescending | |
{ | |
void operator()(const size_t n) | |
{ | |
thrust::host_vector<T> h_vec = unittest::random_integers<T>(n); thrust::sort(h_vec.begin(), h_vec.end(), thrust::greater<T>()); | |
thrust::device_vector<T> d_vec = h_vec; | |
thrust::host_vector<T> h_input = unittest::random_integers<T>(2*n); | |
thrust::device_vector<T> d_input = h_input; | |
typedef typename thrust::host_vector<T>::difference_type int_type; | |
thrust::host_vector<int_type> h_output(2*n); | |
thrust::device_vector<int_type> d_output(2*n); | |
thrust::lower_bound(h_vec.begin(), h_vec.end(), h_input.begin(), h_input.end(), h_output.begin(), thrust::greater<T>()); | |
thrust::lower_bound(d_vec.begin(), d_vec.end(), d_input.begin(), d_input.end(), d_output.begin(), thrust::greater<T>()); | |
ASSERT_EQUAL(h_output, d_output); | |
} | |
}; | |
VariableUnitTest<TestVectorLowerBoundDescending, SignedIntegralTypes> TestVectorLowerBoundDescendingInstance; | |
template <typename T> | |
struct TestVectorUpperBoundDescending | |
{ | |
void operator()(const size_t n) | |
{ | |
thrust::host_vector<T> h_vec = unittest::random_integers<T>(n); thrust::sort(h_vec.begin(), h_vec.end(), thrust::greater<T>()); | |
thrust::device_vector<T> d_vec = h_vec; | |
thrust::host_vector<T> h_input = unittest::random_integers<T>(2*n); | |
thrust::device_vector<T> d_input = h_input; | |
typedef typename thrust::host_vector<T>::difference_type int_type; | |
thrust::host_vector<int_type> h_output(2*n); | |
thrust::device_vector<int_type> d_output(2*n); | |
thrust::upper_bound(h_vec.begin(), h_vec.end(), h_input.begin(), h_input.end(), h_output.begin(), thrust::greater<T>()); | |
thrust::upper_bound(d_vec.begin(), d_vec.end(), d_input.begin(), d_input.end(), d_output.begin(), thrust::greater<T>()); | |
ASSERT_EQUAL(h_output, d_output); | |
} | |
}; | |
VariableUnitTest<TestVectorUpperBoundDescending, SignedIntegralTypes> TestVectorUpperBoundDescendingInstance; | |
template <typename T> | |
struct TestVectorBinarySearchDescending | |
{ | |
void operator()(const size_t n) | |
{ | |
thrust::host_vector<T> h_vec = unittest::random_integers<T>(n); thrust::sort(h_vec.begin(), h_vec.end(), thrust::greater<T>()); | |
thrust::device_vector<T> d_vec = h_vec; | |
thrust::host_vector<T> h_input = unittest::random_integers<T>(2*n); | |
thrust::device_vector<T> d_input = h_input; | |
typedef typename thrust::host_vector<T>::difference_type int_type; | |
thrust::host_vector<int_type> h_output(2*n); | |
thrust::device_vector<int_type> d_output(2*n); | |
thrust::binary_search(h_vec.begin(), h_vec.end(), h_input.begin(), h_input.end(), h_output.begin(), thrust::greater<T>()); | |
thrust::binary_search(d_vec.begin(), d_vec.end(), d_input.begin(), d_input.end(), d_output.begin(), thrust::greater<T>()); | |
ASSERT_EQUAL(h_output, d_output); | |
} | |
}; | |
VariableUnitTest<TestVectorBinarySearchDescending, SignedIntegralTypes> TestVectorBinarySearchDescendingInstance; | |