スラストデバイスベクトルイテレータの位置

Question

イテレータのインデックス/位置を検索しようとしていますが、私はthrust :: distance（）を使用します。しかし、それは奇妙な値を返します。 ベクトルサイズは10です。このメソッドを使用すると、 "131"という値が返されます。 ここでは完全に動作する例です。

#include <thrust/device_vector.h> 
#include <thrust/host_vector.h> 
#include <thrust/reduce.h> 
#include <thrust/extrema.h> 
#include <iostream> 
#include <iomanip> 
#include <thrust/sort.h> 
#include <thrust/copy.h> 
#include <thrust/random.h> 
#include <thrust/unique.h> 
#include <thrust/reduce.h> 
#include <thrust/iterator/constant_iterator.h> 

using namespace std; 
template <typename Vector> 
void print_vector(const std::string& name, const Vector& v) 
{ 
    typedef typename Vector::value_type T; 
    std::cout << " " << std::setw(20) << name << " "; 
    thrust::copy(v.begin(), v.end(), std::ostream_iterator<T>(std::cout, "")); 
    std::cout << std::endl; 
} 

int main() 
{ 
thrust::device_vector<int> x; 
x.push_back(1); 
x.push_back(10); 
x.push_back(1); 
x.push_back(11); 
x.push_back(1); 
x.push_back(11); 
thrust::device_vector<int> y(10); 

print_vector("Original",x); 

thrust::sort(x.begin(),x.end()); 
print_vector("sort",x); 

thrust::device_vector<int>::iterator it=thrust::unique(x.begin(),x.end()); 
std::cout<<*it<<std::endl; 

//int newsize=it-y.begin(); 
int newsize=thrust::distance(y.begin(),it); 
cout<<"nsz:"<<newsize<<endl; 

return 0; 
} 

Answer 1

イテレータitはベクトルxに関して確立されています

thrust::device_vector<int>::iterator it=thrust::unique(x.begin(),x.end()); 
                ^  ^

しかし、あなたはベクトルyの初めに、このイテレータからの距離を求めている：

int newsize=thrust::distance(y.begin(),it); 
          ^

それは意味をなさない。 itとベクトルyの間に定義された関係はありません。

あなたの代わりにベクトルxの先頭までの距離を求める場合は、より賢明な結果を得るでしょう：

$ cat t1244.cu 
#include <thrust/device_vector.h> 
#include <thrust/host_vector.h> 
#include <thrust/reduce.h> 
#include <thrust/extrema.h> 
#include <iostream> 
#include <iomanip> 
#include <thrust/sort.h> 
#include <thrust/copy.h> 
#include <thrust/random.h> 
#include <thrust/unique.h> 
#include <thrust/reduce.h> 
#include <thrust/iterator/constant_iterator.h> 

using namespace std; 
template <typename Vector> 
void print_vector(const std::string& name, const Vector& v) 
{ 
    typedef typename Vector::value_type T; 
    std::cout << " " << std::setw(20) << name << " "; 
    thrust::copy(v.begin(), v.end(), std::ostream_iterator<T>(std::cout, " ")); 
    std::cout << std::endl; 
} 

int main() 
{ 
thrust::device_vector<int> x; 
x.push_back(1); 
x.push_back(10); 
x.push_back(1); 
x.push_back(11); 
x.push_back(1); 
x.push_back(11); 
thrust::device_vector<int> y(10); 

print_vector("Original",x); 

thrust::sort(x.begin(),x.end()); 
print_vector("sort",x); 

thrust::device_vector<int>::iterator it=thrust::unique(x.begin(),x.end()); 
std::cout<<*it<<std::endl; 

//int newsize=it-y.begin(); 
int newsize=thrust::distance(x.begin(),it); 
cout<<"nsz:"<<newsize<<endl; 

return 0; 
} 
$ nvcc -o t1244 t1244.cu 
$ ./t1244 
       Original 1 10 1 11 1 11 
        sort 1 1 1 10 11 11 
10 
nsz:3 
$