1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
|
#include <thrust/functional.h>
#include <thrust/set_operations.h>
#include <thrust/sort.h>
#include <unittest/unittest.h>
template <typename Vector>
void TestSetDifferenceDescendingSimple()
{
typedef typename Vector::value_type T;
typedef typename Vector::iterator Iterator;
Vector a(4), b(5);
// clang-format off
a[0] = 5; a[1] = 4; a[2] = 2; a[3] = 0;
b[0] = 6; b[1] = 4; b[2] = 3; b[3] = 3; b[4] = 0;
// clang-format on
Vector ref(2);
ref[0] = 5;
ref[1] = 2;
Vector result(2);
Iterator end = thrust::set_difference(a.begin(), a.end(), b.begin(), b.end(), result.begin(), thrust::greater<T>());
ASSERT_EQUAL_QUIET(result.end(), end);
ASSERT_EQUAL(ref, result);
}
DECLARE_VECTOR_UNITTEST(TestSetDifferenceDescendingSimple);
template <typename T>
void TestSetDifferenceDescending(const size_t n)
{
thrust::host_vector<T> temp = unittest::random_integers<T>(2 * n);
thrust::host_vector<T> h_a(temp.begin(), temp.begin() + n);
thrust::host_vector<T> h_b(temp.begin() + n, temp.end());
thrust::sort(h_a.begin(), h_a.end(), thrust::greater<T>());
thrust::sort(h_b.begin(), h_b.end(), thrust::greater<T>());
thrust::device_vector<T> d_a = h_a;
thrust::device_vector<T> d_b = h_b;
thrust::host_vector<T> h_result(n);
thrust::device_vector<T> d_result(n);
typename thrust::host_vector<T>::iterator h_end;
typename thrust::device_vector<T>::iterator d_end;
h_end =
thrust::set_difference(h_a.begin(), h_a.end(), h_b.begin(), h_b.end(), h_result.begin(), thrust::greater<T>());
h_result.resize(h_end - h_result.begin());
d_end =
thrust::set_difference(d_a.begin(), d_a.end(), d_b.begin(), d_b.end(), d_result.begin(), thrust::greater<T>());
d_result.resize(d_end - d_result.begin());
ASSERT_EQUAL(h_result, d_result);
}
DECLARE_VARIABLE_UNITTEST(TestSetDifferenceDescending);
|
