我想做的是通过thrust::reduce_by_key通过键获取值的平均值。我首先使用sort_by_key,按reduce_by_key的连续键进行分组就可以了。我使用this来帮助我达到目标。但是,我遇到了许多我无法理解的错误(这也是我第一次使用reduce_by_key),并且我想不出一种更好的方法来实现此目的,而无需使用大量临时分配给(1)来获得按键计算值,然后按键计数和(2)将两者相除得出平均值。

input keys:   1,   1,   1,  2,   3,   5,  5,  2
input values: 120, 477, 42, 106, 143, 53, 83, 24

expected output values: 213, 65, 143, 68


我有以下自定义函子:

struct GetAverage
{
    template<typename Tuple>
    __host__ __device__
    int operator()(const Tuple& t)
    {
        //SumByKey / CountByKey
        return thrust::get<0>(t) / thrust::get<1>(t);
    }
};


函子是从以下代码(位于main()中)调用的

thrust::device_vector<unsigned int> tempKey(8);
thrust::device_vector<unsigned int> tempValue(8);

tempKey[0] = 1;
tempKey[1] = 1;
tempKey[2] = 1;
tempKey[3] = 2;
tempKey[4] = 3;
tempKey[5] = 5;
tempKey[6] = 5;
tempKey[7] = 2;

tempValue[0] = 120;
tempValue[1] = 477;
tempValue[2] = 42;
tempValue[3] = 106;
tempValue[4] = 143;
tempValue[5] = 53;
tempValue[6] = 83;
tempValue[7] = 24;

thrust::sort_by_key(tempKey.begin(), tempKey.end(), tempValue.begin());

thrust::equal_to<int> binary_pred;
thrust::reduce_by_key(
    tempKey.begin(),
    tempKey.end(),
    thrust::make_zip_iterator(
        thrust::make_tuple(
            tempValue.begin(),
            thrust::make_constant_iterator(1)
        )
    ), //values_first; Should go into GetAverage() custom functor as a zipped tuple <tempValue, 1>
    tempKey.begin(), //keys_output; Should be returning the unique keys
    tempValue.begin(), //values_output; Should be returning the average by key
    binary_pred,
    GetAverage()
);


错误示例:
-no instance of function template "GetAverage::operator()" matches the argument list
-no operator "=" matches these operands
-no suitable conversion function from "InputValueType" to "TemporaryType" exists
-no suitable conversion function from "thrust::detail::tuple_of_iterator_references<thrust::device_reference<int>, int, thrust::null_type, thrust::null_type, thrust::null_type, thrust::null_type, thrust::null_type, thrust::null_type, thrust::null_type, thrust::null_type>" to "TemporaryType" exists

有人对如何解决这个问题有任何想法吗?还是链接?我阅读了此处使用的所有内容的文档,并尽力尝试尽可能地理解它,但没有解决方案。谢谢!

更新
请参阅Eric的答案。结合他的说法,这是新的源代码。创建了一个操作来处​​理元组的加号。该代码唯一不做的事情是在reduce_by_key调用之后,应在结果上使用thrust::transform以通过将总和除以计数来获得平均值。

// --- Defining key tuple type
typedef thrust::tuple<int, int> Tuple;

/* PLUS OPERATOR BETWEEN TUPLES */
struct TuplePlus
{
    __host__ __device__
    Tuple operator ()(const Tuple& lhs, const Tuple& rhs)
    {
        return thrust::make_tuple(
            thrust::get<0>(lhs) + thrust::get<0>(rhs),
            thrust::get<1>(lhs) + thrust::get<1>(rhs)
        );
    }
};


main()内部,现在有以下内容。

thrust::equal_to<int> binary_pred;
thrust::reduce_by_key(
    tempKey.begin(),
    tempKey.end(),
    thrust::make_zip_iterator(
        thrust::make_tuple(
            tempValue.begin(),
            thrust::make_constant_iterator(1)
        )
    ), //values_first; goes in as a zipped up tuple <value, 1>
    tempKey.begin(), //keys_output
    thrust::make_zip_iterator(
        thrust::make_tuple(
            tempValue.begin(),
            tempCount.begin()
        )
    ), //values_output; ZipIterator<Sum, Count> by key
    binary_pred,
    TuplePlus()
);

最佳答案

有两个问题。

元组序列减少的结果应该是元组,而不是int。根据文件

https://thrust.github.io/doc/group__reductions.html#ga633d78d4cb2650624ec354c9abd0c97f

最后一个参数binary_op的类型应为


  BinaryFunction是Binary Function的模型,BinaryFunction的result_type可转换为OutputIterator2的value_type。


这意味着您的还原操作应类似于

struct GetSum
{
  template<typename Tuple>
  __host__ __device__
  Tuple operator()(const Tuple& a, construction Tuple& b)
  {
    ...
  }
}


另一方面,在缩减阶段,您只能有效地计算总和,而不能有效地计算平均值。这意味着您的values_output也应该是与values_first相同类型的zip迭代器。


  OutputIterator2是Output Iterator的模型,InputIterator2的value_type可转换为OutputIterator2的value_type。


因此,您需要两个结果数组,一个用于键的总和,另一个用于键的计数。它们应压缩在一起并用作values_output

然后,您需要另一个thrust::transform来计算最终结果-按键求平均值。



您也可以尝试@RobertCrovella提出的方法,该方法使用单个thrust::reduce_by_key来计算平均值。

Output from reduce_by_key() as a function of two reduced vectors

关于c++ - 推力:reduce_by_key将zip_iterator(tuple)传递到自定义仿函数以按键检索平均值,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/37890860/

10-12 15:23