wal*_*mat 11 c++ stl vector c++11
我有一个期望std::vector
(调用它A
)的算法.不过,我已经B
有N + 2
条目和我基本上要为通过B.data() + 2
,因此该算法得到最后N
的条目B
.如果A
被修改,那么也是如此B
.
当使用double*
指针时,我非常清楚我应该如何制作它,但这对std::vector
s 也是可能的吗?我的意思是,矢量的好处是它为我处理内存,我现在想要的是禁止它(如果B
或被A
破坏,它们应该保持指向的数据不变).
像这样:
std::vector< double > B({1,2,3,4,5});
std::vector< double > A(B.data() + 2, B.size() - 2);
// A and B share data now. A is 3, 4, 5
Run Code Online (Sandbox Code Playgroud)
我知道通过采用一对迭代器可以更好地为此目的设计算法,但这不在我手中.
UPDATE
(在评论中,有人希望看到签名,这里是)
nlopt::result nlopt::opt::optimize(std::vector<double> &x, double &opt_f);
Run Code Online (Sandbox Code Playgroud)
但是,我的初衷是非常聪明,让算法在我的向量中优化direclty B
,所以我最终得到的是这样的:
std::vector< double > B(N + 2);
// do something with it here, like put initial values
std::vector< double > A( B.begin() + 2, B.end() );
optimize( A );
std::copy(A.begin(), A.end(), B.begin() + 2);
Run Code Online (Sandbox Code Playgroud)
我真的不担心这个解决方法,我还在文档中读到,nlopt
如果使用C++接口而不是C接口,那么向量内部副本无论多少次.
但同样,这个例子真的让我开始关注算法界面设计,并且值得投入一些时间浏览一下像升级库range
等等.
实际上,我在项目中遇到了类似的问题,这是我提出的解决方案:
#include <iterator>
#include <type_traits>
/** \brief Provides an interface to random accessible ranges
* \tparam RAIterator must be a random access iterator
*
* This class doesn't allocate any memory by itself and
* provides only an interface into the data of another
* container.
* \attention Keep in mind that although all methods are
* \c const, the returned RAIterator might not be a \c const_iterator
* at all. It is your responsibility to make sure that
* you don't invalidate the given range while working on it.
**/
template <class RAIterator>
class Slice
{
public:
//! Typedef for convenience when working with standard algorithms
typedef RAIterator iterator;
//! Alias to the iterator's reference type
typedef typename std::iterator_traits<RAIterator>::reference reference;
//! Creates the slice.
//! \param first, last a valid range
Slice(RAIterator first, RAIterator last) : first(first), last(last){}
//! Creates the slice.
//! \param first iterator to the first element
//! \param length of the range [first, first+length)
//! \remark if length is negative, an empty slice is assumed.
Slice(RAIterator first, int length) :
first(first),
last( length > 0 ? first + length : first)
{}
//! The default constructor is deleted, as it would resemble an empty slice
Slice() = delete;
///@{
//! \brief Defaulted construcors.
Slice(const Slice&) = default;
Slice(Slice&&) = default;
Slice& operator=(const Slice&)= default;
Slice& operator=(Slice&&)= default;
/**@}*/
//! Returns an iterator to the begin of the range
RAIterator begin() const{ return first; }
//! Returns an iterator to the end of the range
RAIterator end() const{ return last; }
//! Returns the size of the slice.
typename std::iterator_traits<RAIterator>::difference_type size() const{
return std::distance(first, last);
}
//! Provides random access for the values interfaced by Slice
reference operator[](size_t index) const { return first[index]; }
private:
RAIterator first; //!< begin of the range
RAIterator last; //!< end of the range
};
/** \brief Creates a slice from the given range
* \tparam RAIterator should be an random access iterator
* \returns a slice [first,last)
* \param first, last is the range
**/
template <class RAIterator>
Slice<RAIterator> make_slice(RAIterator first, RAIterator last){
return Slice<RAIterator>(first, last);
}
Run Code Online (Sandbox Code Playgroud)
现在您可以使用Slice
示例中的just:
std::vector< double > B({1,2,3,4,5});
Slice A(B.begin() + 2, B.size() - 2);
A[0] = 5;
// B == {1,2,5,4,5}
Run Code Online (Sandbox Code Playgroud)
编辑:如果你想要一个更成熟的切片,boost::adaptors::slice
请尽可能使用.
归档时间: |
|
查看次数: |
953 次 |
最近记录: |