Question

我有std::vector<double>，我需要插入其值。例如，只有1个中间值，并给出一个用

填充的向量

1 / 2 / 3 / 4

我想访问以下值

1 / 1.5 / 2 / 2.5 / 3 / 3.5 / 4

当然我不需要存储这个中间值（简单的线性插值，我不必经常阅读它们），所以我写了这个简单的类：

typedef std::vector<double> DVector;
class InterpolatedVector {
public:
    InterpolatedVector(const DVector& v,int steps) : v(v),steps(steps){}
    double at(int i){
        int j = i%steps;
        int k = (int)i/steps;
        if (i==0){return v[0];}
        else if (i==v.size()*steps){return v.back();}
        else {return  ((steps-j)*v[k] + j*v[k+1])/steps;}
    }
    int size(){return steps*(v.size()-1) + 1;}
private:
    DVector v;
    int steps;
};

它工作正常，我得到（几乎）我想要的东西。但是，这个“容器”我不能使用std :: algorithms而我没有迭代器。（当然，我不能写中间值，但至少在阅读时，我想使用算法。）我应该提到我仍然对迭代器等缺乏一点了解。

我应该如何实现这个“InterpolatedVector”，以便我可以做

之类的事情

std::accumulate(/* passing Interpolated iterators? */ );

吗

Answer 1

鉴于您已经拥有处理索引本身的代码，将其包装为迭代器非常简单。如果你原谅我，我也会更加通用。

#include <vector>
#include <iterator>

template <class T>
class InterpolatedVector {
    typedef std::vector<T> DVector;
public:
    InterpolatedVector(const DVector& v,int steps) : v(v),steps(steps){}
    T at(int i){
        int j = i%steps;
        int k = (int)i/steps;
        if (i==0){return v[0];}
        else if (i==v.size()*steps){return v.back();}
        else {return  ((steps-j)*v[k] + j*v[k+1])/steps;}
    }
    int size(){return steps*(v.size()-1) + 1;}

    class iterator : public std::iterator < std::random_access_iterator_tag, T > {
        InterpolatedVector *vec;
        int index;
    public:
        iterator(InterpolatedVector &d, int index) : vec(&d), index(index) {}
        iterator &operator++() { ++index; return *this; }
        iterator operator++(int) { iterator tmp{ *vec, index }; ++index; return tmp; }
        iterator operator+(int off) { return iterator(*vec, index + off); }
        iterator operator-(int off) { return iterator(*vec, index - off); }
        value_type operator*() { return (*vec).at(index);   }
        bool operator!=(iterator const &other) { return index != other.index; }
        bool operator<(iterator const &other) { return index < other.index; }
    };

    iterator begin() { return iterator(*this, 0); }
    iterator end() { return iterator(*this, size()); }
private:

    DVector v;
    int steps;
};

...以及一些快速的演示代码来测试它：

#include <iostream>

int main() {
    std::vector<double> d{ 1, 2, 3, 4 };
    InterpolatedVector<double> id(d, 2);

    std::copy(id.begin(), id.end(), std::ostream_iterator<double>(std::cout, "\t"));
    std::cout << "\n";

    std::vector<int> i{ 0, 5 };
    InterpolatedVector<int> ii(i, 5);

    std::copy(ii.begin(), ii.end(), std::ostream_iterator<int>(std::cout, "\t"));
}

输出：

1       1.5     2       2.5     3       3.5     4
0       1       2       3       4       5

当然，有些算法仍然无法通过这种＆＃34;收集＆＃34;来做很多事情。尝试将插值集合提供给std::sort并不是很有意义（您显然需要对底层容器进行排序）。只要算法只需要读取数据，这应该没问题。

Answer 2

如果要将其与算法一起使用，请使用中间值填充向量。这会简单得多。

如何实现“InterpolatedVector”？

2 个答案: