c++ - std::reduce 似乎将结果转换为整数-6ren

c++ - std::reduce 似乎将结果转换为整数

转载作者：塔克拉玛干更新时间：2023-11-02 23:54:53

28

4

<分区>

我正在比较这两个函数:

double polynomials(const vector<double> & coeffs, double x) {
    double sum = 0.0;
    double factor = 1.0;
    for (double coeff: coeffs) {
        sum += coeff * factor;
        factor *= x;
    }
    return sum;
}

和

double algorithm_polynomials(const vector<double> & coeffs, double x) {
    return reduce(execution::seq, cbegin(coeffs), end(coeffs), 0, [factor=1.0, x](double sum, double coeff) mutable {
        double curr_sum = sum + coeff * factor;
        factor *= x;
        return curr_sum;
    });
}

对于整数系数和整数x，这两个函数的结果是等价的。但对于实际值，结果有所不同。第一个函数给出的是正确的，而第二个函数给出的总是整数(不一定接近真实答案)。

我的猜测是每个中间值在求和之前都被转换为整数，但我不明白为什么会发生这种情况，或者我该如何解决它。

要运行此代码，请使用

#include <numeric>
#include <vector>
#include <algorithm>
#include <execution>

using namespace std;

需要带有 -std=c++17 标志的 gcc 9 或更高版本以及从 2019 年到 tbb 的链接。

编辑:感谢大家指出我忽略的初始值类型。

我感到有点惭愧，既因为没有看到它，也因为 reduce 函数使用不当(由于可变的 lambda 而无法将其并行化)。

为了救赎自己，我在这里添加了一种使其并行的方法...

class counter: public std::iterator<
        std::random_access_iterator_tag,   // iterator_category
        size_t,                      // value_type
        size_t,                      // difference_type
        const size_t*,               // pointer
        size_t                       // reference
>{
    size_t num = 0;
public:
    explicit counter(size_t _num) : num(_num) {}
    counter& operator++() {num += 1; return *this;}
    counter operator++(int) {counter retval = *this; ++(*this); return retval;}
    bool operator==(counter other) const {return num == other.num;}
    bool operator!=(counter other) const {return !(*this == other);}
    counter& operator+=(size_t i) { num += i; return *this; }
    counter& operator-=(size_t i) { num -= i; return *this; }
    counter operator +(counter &other) const { return counter(num + other.num);}
    counter operator -(counter &other) const { return counter(num - other.num); }
    counter operator +(size_t i) const { return counter(num + i); }
    counter operator -(size_t i) const {return counter(num - i); }
    reference operator*() const {return num;}
};


double better_algorithm_polinomials(const vector<double> & coeffs, double x) {
    //this has the advantage of being easily parallelized
    return transform_reduce(execution::par, cbegin(coeffs), end(coeffs), counter(0), 0.0, plus{}, [x](double coeff, size_t index) { return coeff * pow<double>(x, index); });
}

28

4

0

文章推荐： linux - 连字符在 Linux 命令选项中的用法

文章推荐： linux - Linux内核中的原子操作

文章推荐： c++ - 如何在命令行界面中处理 ctrl-break 信号

c++ - 为什么 `std::common_type_t` 等于 `std::ostream` 而不是 `std::ostream &` ？
我正在开发一个小型图书馆，我需要做的一件事是让访问者访问一些数据并返回结果。在一些较旧的 C++ 代码中，访问者需要声明一个 typedef return_type .例如，boost::stati
c++ - std::map 麻烦
我正在尝试使用std:map类型的键和值制作std::any Visual Studio 2017 std::map m("lastname", "Ivanov"); std::cout (m["la
C++ std::map> 。如何循环设定值？
我已经在 C++ 的 map 中声明了一个集合为 std::map> .如何循环访问或打印设定值？最佳答案如果你知道如何迭代 std::map或 std::set单独地，您应该可以毫无问题地组合迭
C++ 循环 std::vector>
如何循环？我已经试过了: //----- code std::vector >::iterator it; for ( it = users.begin(); it != users.end();
c++ - std::unique_lock 还是 std::lock_guard？
我有两个用例。 A.我想同步访问两个线程的队列。 B.我想同步两个线程对队列的访问并使用条件变量，因为其中一个线程将等待另一个线程将内容存储到队列中。对于用例 A，我看到了使用 std::lock_
c++ - std::trivially_copyable_v 和 std::is_pod_v 之间有什么区别(std::is_standard_layout && std::is_trivial_v)
我正在查看这两种类型特征的文档，但不确定有什么区别。我不是语言律师，但据我所知，它们都适用于“memcpy-able”类型。它们可以互换使用吗？最佳答案不，这些术语不能互换使用。这两个术语都表示
c++ - 为什么我可以有一个 std::vector 而不是 std::vector？
我有以下测试代码，其中有一个参数 fS，它是 ofstream 的容器: #include #include #include #include int
c++ - std::unordered_map
这是这个问题的延续 c++ function ptr in unorderer_map, compile time error 我试图使用 std::function 而不是函数指针，并且只有当函数是

c++ - 将 std::any_of、std::all_of、std::none_of 等与 std::map 一起使用
std::unordered_map str_bool_map = { {"a", true}, {"b", false}, {"c", true} }; 我们可以在此映射上使
c++ - 使用 std::find 检查 std::vector> 中的项目
我有以下对象 std::vector> vectorList; 然后我添加到这个使用 std::vector vec_tmp; vec_tmp.push_back(strDRG); vec_tmp.p
c++ - 为什么 std::initializer_list 不支持 std::get<>、std::tuple_size 和 std::tuple_element
为什么 std::initializer_list不支持std::get<> , std::tuple_size和 std::tuple_element ？在constexpr中用得很多现在的表达式，
c++ - std::tuple 和 std::tuple 是否被 std::vector 视为同一类型？
我有一个像这样定义的变量 auto drum = std::make_tuple ( std::make_tuple ( 0.3f , Ex
c++ :将 std::map 转换为 std::map
假设我有一个私有(private)std::map在我的类(class)里std::map 。我怎样才能将其转换为std::map返回给用户？我想要下面的原型(prototype) const std
c++ :将 std::map 转换为 std::map
假设我有一个私有(private)std::map在我的类(class)里std::map 。我怎样才能将其转换为std::map返回给用户？我想要下面的原型(prototype) const std
c++ - 在带有 std::ref 的 std::thread 中使用地址清理调用 std::invoke(std::forward(...)) 时的奇怪行为
问题我正在尝试将 lambda 闭包传递给 std::thread，它使用任意封闭参数调用任意封闭函数。 template std::thread timed_thread(Function&& f
c++ - 具有模板模板参数的模板定义，可以专门化为类，例如，std::vector 或 std::map
我想创建一个模板类，可以容纳容器和容器的任意组合。例如，std::vector或 std::map ，例如。我尝试了很多组合，但我必须承认模板的复杂性让我不知所措。我编译的关闭是这样的: templ
c++ - 将 std::vector> 分配给另一个 std::vector>
我有一个 std::vector>我将其分配给相同类型的第二个 vector 。我收到这个编译器错误: /opt/gcc-8.2.0/include/c++/8.2.0/bits/stl_algob
c++ - 将 std::vector> 移动到 std::vector>
有时候，我们有一个工厂可以生成一个 std::unique_ptr vector ，后来我们想在类/线程/你命名的之间共享这些指针。因此，最好改用 std::shared_ptr 。当然有一种方法可以
c++ - 为什么 std::sort 假定 std::vector< std::vector> 默认为 std::vector，从而产生错误的结果？
这个问题在这里已经有了答案: Sorting a vector of custom objects (14 个答案) 关闭 6 年前。我创建了一个 vector vector ，我想根据我定义的参
c++ - 将 std::vector> 转换为 std::vector>
我有三个类(class)成员: public: std::vector > getObjects(); std::vector > getObjects() const; privat

首页

博学

6Ren·AI

商城

c++ - std::reduce 似乎将结果转换为整数