LIVE / thrust /testing /cuda /adjacent_difference.cu
Xu Ma
update
1c3c0d9
raw
history blame
4.71 kB
#include <unittest/unittest.h>
#include <thrust/adjacent_difference.h>
#include <thrust/execution_policy.h>
#include <thrust/device_malloc.h>
#include <thrust/device_free.h>
template<typename ExecutionPolicy, typename Iterator1, typename Iterator2>
__global__ void adjacent_difference_kernel(ExecutionPolicy exec, Iterator1 first, Iterator1 last, Iterator2 result)
{
thrust::adjacent_difference(exec, first, last, result);
}
template<typename ExecutionPolicy, typename Iterator1, typename Iterator2, typename BinaryFunction>
__global__ void adjacent_difference_kernel(ExecutionPolicy exec, Iterator1 first, Iterator1 last, Iterator2 result, BinaryFunction f)
{
thrust::adjacent_difference(exec, first, last, result, f);
}
template<typename T, typename ExecutionPolicy>
void TestAdjacentDifferenceDevice(ExecutionPolicy exec, const size_t n)
{
thrust::host_vector<T> h_input = unittest::random_samples<T>(n);
thrust::device_vector<T> d_input = h_input;
thrust::host_vector<T> h_output(n);
thrust::device_vector<T> d_output(n);
thrust::adjacent_difference(h_input.begin(), h_input.end(), h_output.begin());
adjacent_difference_kernel<<<1,1>>>(exec, d_input.begin(), d_input.end(), d_output.begin());
{
cudaError_t const err = cudaDeviceSynchronize();
ASSERT_EQUAL(cudaSuccess, err);
}
ASSERT_EQUAL(h_output, d_output);
thrust::adjacent_difference(h_input.begin(), h_input.end(), h_output.begin(), thrust::plus<T>());
adjacent_difference_kernel<<<1,1>>>(exec, d_input.begin(), d_input.end(), d_output.begin(), thrust::plus<T>());
{
cudaError_t const err = cudaDeviceSynchronize();
ASSERT_EQUAL(cudaSuccess, err);
}
ASSERT_EQUAL(h_output, d_output);
// in-place operation
thrust::adjacent_difference(h_input.begin(), h_input.end(), h_input.begin(), thrust::plus<T>());
adjacent_difference_kernel<<<1,1>>>(exec, d_input.begin(), d_input.end(), d_input.begin(), thrust::plus<T>());
{
cudaError_t const err = cudaDeviceSynchronize();
ASSERT_EQUAL(cudaSuccess, err);
}
ASSERT_EQUAL(h_input, h_output); //computed previously
ASSERT_EQUAL(d_input, d_output); //computed previously
}
template<typename T>
void TestAdjacentDifferenceDeviceSeq(const size_t n)
{
TestAdjacentDifferenceDevice<T>(thrust::seq, n);
}
DECLARE_VARIABLE_UNITTEST(TestAdjacentDifferenceDeviceSeq);
template<typename T>
void TestAdjacentDifferenceDeviceDevice(const size_t n)
{
TestAdjacentDifferenceDevice<T>(thrust::device, n);
}
DECLARE_VARIABLE_UNITTEST(TestAdjacentDifferenceDeviceDevice);
void TestAdjacentDifferenceCudaStreams()
{
cudaStream_t s;
cudaStreamCreate(&s);
thrust::device_vector<int> input(3);
thrust::device_vector<int> output(3);
input[0] = 1; input[1] = 4; input[2] = 6;
thrust::adjacent_difference(thrust::cuda::par.on(s), input.begin(), input.end(), output.begin());
cudaStreamSynchronize(s);
ASSERT_EQUAL(output[0], 1);
ASSERT_EQUAL(output[1], 3);
ASSERT_EQUAL(output[2], 2);
cudaStreamDestroy(s);
}
DECLARE_UNITTEST(TestAdjacentDifferenceCudaStreams);
struct detect_wrong_difference
{
bool * flag;
__host__ __device__ detect_wrong_difference operator++() const { return *this; }
__host__ __device__ detect_wrong_difference operator*() const { return *this; }
template<typename Difference>
__host__ __device__ detect_wrong_difference operator+(Difference) const { return *this; }
template<typename Index>
__host__ __device__ detect_wrong_difference operator[](Index) const { return *this; }
__device__
void operator=(long long difference) const
{
if (difference != 1)
{
*flag = false;
}
}
};
void TestAdjacentDifferenceWithBigIndexesHelper(int magnitude)
{
thrust::counting_iterator<long long> begin(1);
thrust::counting_iterator<long long> end = begin + (1ll << magnitude);
ASSERT_EQUAL(thrust::distance(begin, end), 1ll << magnitude);
thrust::device_ptr<bool> all_differences_correct = thrust::device_malloc<bool>(1);
*all_differences_correct = true;
detect_wrong_difference out = { thrust::raw_pointer_cast(all_differences_correct) };
thrust::adjacent_difference(thrust::device, begin, end, out);
bool all_differences_correct_h = *all_differences_correct;
thrust::device_free(all_differences_correct);
ASSERT_EQUAL(all_differences_correct_h, true);
}
void TestAdjacentDifferenceWithBigIndexes()
{
TestAdjacentDifferenceWithBigIndexesHelper(30);
TestAdjacentDifferenceWithBigIndexesHelper(31);
TestAdjacentDifferenceWithBigIndexesHelper(32);
TestAdjacentDifferenceWithBigIndexesHelper(33);
}
DECLARE_UNITTEST(TestAdjacentDifferenceWithBigIndexes);