Thrust快速入门教程（三）——迭代器与静态调度

在这节中我们曾使用了这样的表达式，H.begin() 、H.end()、D.begin() + 7。begin()与end()的返回值在C++中被称为迭代器。vector的迭代器类似于数组的指针，用于指向数组的某个元素。H.begin()是指向H容器中数组第一个元素的迭代器。类似，H.end()指向H容器中的最后一个元素。

虽然说迭代器类似于指针，但它有着更丰富的作用。可以注意到在使用thrust::fill的时候我们并不需要指明这是device_vector的迭代器。这些信息包含在了D.begin()的返回值的迭代器类型中，其类型不同于H.begin()的返回值。当Thrust中的函数调用时，将根据迭代器的类型选择使用主机端还是设备端的算法实现。因为主机/设备调度是在编译时解析，所以这一过程被称为被称为静态调度。请注意，这意味着运行时没有调度进程的开销。

你可能想知道当“原始”指针作为Thrust函数的参数会如何。和STL一样，Thrust允许这种用法，它会调度主机端算法。如果传入的指针是指向设备端内存的指针，那么在调用函数之前需要用thrust::device_ptr封装。例如：

size_t N = 10; // raw pointer to device memory int * raw_ptr ; cudaMalloc (( void **) & raw_ptr , N * sizeof ( int )); // wrap raw pointer with a device_ptr thrust :: device_ptr dev_ptr ( raw_ptr ); // use device_ptr in thrust algorithms thrust :: fill ( dev_ptr , dev_ptr + N, ( int ) 0); 

如需从 device_ptr中提取“原始”指针需要使用raw_pointer_cast，用法如下：

size_t N = 10; // create a device_ptr thrust :: device_ptr dev_ptr = thrust :: device_malloc (N); // extract raw pointer from device_ptr int * raw_ptr = thrust :: raw_pointer_cast ( dev_ptr ); 

迭代器另一个区别于指针的地方在于它可以遍历各种数据结构。例如，STL提供了链表容器（std::list），提供双向的（但不是随机访问）的迭代器。虽然Thrust不提供这类容器的设备端实现，但是与它们兼容。 

# include # include # include # include int main ( void ) { // create an STL list with 4 values std :: list stl_list ; stl_list . push_back (10) ; stl_list . push_back (20) ; stl_list . push_back (30) ; stl_list . push_back (40) ; // initialize a device_vector with the list thrust :: device_vector D( stl_list . begin () , stl_list . end ()); // copy a device_vector into an STL vector std :: vector stl_vector (D. size ()); thrust :: copy (D. begin () , D. end () , stl_vector . begin ()); return 0; } 

 

 

备注：到目前为止，我们所讨论的是十分有用但相当基本的常用迭代器。除了这些常用迭代器，Thrust也提供了counting_iterator和zip_iterator这类特殊迭代器 。虽然他们看起来与常用迭代器一样，但是特殊迭代器能够提供更令人兴奋的特性。我们将在后面的教程讨论这个问题。