thrust

I have the following (compilable and executable) code using CUDA Thrust to perform reductions of float2 arrays. It works correctly using namespace std; // includes, system #include <stdlib.h> #include <stdio.h> #include <string.h> #include <math.h> #include <conio.h> #include <typeinfo> #include <iostream> // includes CUDA #include <cuda.h> #include <cuda_runtime.h> // includes Thrust #include <thrust/host_vector.h> #include <thrust/device_vector.h> #include <thrust/reduce.h> // float2 + struct struct add_float2 { __device__ float2 operator()(const float2& a, const float2& b) const { float2 r;

I am using visual studio 2010, parallel nsight 2.2 and cuda 4.2 for learning. My system is Windows 8 pro x64. I opened the radix sort project which included by cuda computing SDK in VS, and compiled it with no error. The sort code uses thrust library: if(keysOnly) thrust::sort(d_keys.begin(), d_keys.end()); else thrust::sort_by_key(d_keys.begin(), d_keys.end(), d_values.begin()); I want to know how thrust dispatch the sort function to cuda kernels, so I tried to add breakpoints in front of lines above and compiled the project in debug mode. But when I use parallel nsight for cuda debugging,

I rendered a scene with opengl (I can also render it to a texture) I want to use CUDA / Thrust to sort this rendered image How do I link the texture I made from : cudaGraphicsGLRegisterImage to be used via thrust? maybe something like this ? how to calculate an average from a int2 array using Thrust Robert Crovella I'm not sure it makes sense to try and use textures directly with thrust. However using an ordinary GL pixel buffer can be made to work directly with thrust. The following example creates an openGL pixel buffer with a particular green/black pattern, and then displays it. When you

How can I impliment this function in Thrust? for (i=0;i<n;i++) if (i==pos) h1[i]=1/h1[i]; else h1[i]=-h1[i]/value; In CUDA I did it like: __global__ void inverse_1(double* h1, double value, int pos, int N) { int i = blockDim.x * blockIdx.x + threadIdx.x; if (i < N){ if (i == pos) h1[i] = 1 / h1[i]; else h1[i] = -h1[i] / value; } } Thanks! You need to create a binary functor to apply the operation, then use a counting iterator as the second input. You can pass pos and value into the functor's constructor. It'd look something like: struct inv1_functor { const int pos; const double value; inv1

I'm trying to introduce some CUB into my "old" Thrust code, and so have started with a small example to compare thrust::reduce_by_key with cub::DeviceReduce::ReduceByKey , both applied to thrust::device_vectors . The thrust part of the code is fine, but the CUB part, which naively uses raw pointers obtained via thrust::raw_pointer_cast, crashes after the CUB calls. I put in a cudaDeviceSynchronize() to try to solve this problem, but it didn't help. The CUB part of the code was cribbed from the CUB web pages. On OSX the runtime error is: libc++abi.dylib: terminate called throwing an exception

I'm using the library Thrust to get intersection set of a two larger sets of integers. In test with 2 small inputs i got correct results, but when i use two sets with 10^8 and 65535*1024 elements i got a empty set. Who's can explain this problem? Changing the two first variables to smaller values the thrust returns a expected intersection set. My code is following. #include <thrust/set_operations.h> #include <thrust/device_vector.h> #include <thrust/device_ptr.h> #include <iostream> #include <stdio.h> int main() { int sizeArrayLonger = 100*1000*1000; int sizeArraySmaller = 65535*1024; int