次のコード(crtp.ccに格納されている)でコンパイラに依存する問題が発生しました:
#include <vector>
#include <cassert>
#include <iostream>
template < class Derived >
class AlgebraicVectorExpression {
public:
typedef std::vector<double>::size_type SizeType;
typedef std::vector<double>::value_type ValueType;
typedef std::vector<double>::reference ReferenceType;
SizeType size() const {
return static_cast<const Derived&>(*this).size();
}
ValueType operator[](SizeType ii) const {
return static_cast<const Derived&>(*this)[ii];
}
operator Derived&() {
return static_cast<Derived&>(*this);
}
operator const Derived&() const {
return static_cast< const Derived& >(*this);
}
};
template< class T1, class T2>
class AlgebraicVectorSum : public AlgebraicVectorExpression< AlgebraicVectorSum<T1,T2> > {
const T1 & a_;
const T2 & b_;
typedef typename AlgebraicVectorExpression< AlgebraicVectorSum<T1,T2> >::SizeType SizeType;
typedef typename AlgebraicVectorExpression< AlgebraicVectorSum<T1,T2> >::ValueType ValueType;
public:
AlgebraicVectorSum(const AlgebraicVectorExpression<T1>& a, const AlgebraicVectorExpression<T1>& b) :
a_(a), b_(b) {
assert(a_.size() == b_.size());
}
SizeType size() const {
return a_.size();
}
ValueType operator[](SizeType ii) const {
return (a_[ii] + b_[ii]);
}
};
template< class T1, class T2>
const AlgebraicVectorSum<T1,T2>
operator+(const AlgebraicVectorExpression<T1>& a, const AlgebraicVectorExpression<T2>& b) {
return AlgebraicVectorSum<T1,T2>(a,b);
}
class AlgebraicVector : public AlgebraicVectorExpression<AlgebraicVector>{
std::vector<double> data_;
public:
SizeType size() const {
return data_.size();
}
ValueType operator[](SizeType ii) const {
return data_[ii];
}
ValueType& operator[](SizeType ii) {
return data_[ii];
}
AlgebraicVector(SizeType n) : data_(n,0.0) {
};
template< class T>
AlgebraicVector(const AlgebraicVectorExpression<T>& vec) {
const T& v = vec;
data_.resize(v.size());
for( SizeType idx = 0; idx != v.size(); ++idx) {
data_[idx] = v[idx];
}
}
};
int main() {
AlgebraicVector x(10);
AlgebraicVector y(10);
for (int ii = 0; ii != 10; ++ii)
x[ii] = y[ii] = ii;
AlgebraicVector z(10);
z = x + y;
for(int ii = 0; ii != 10; ++ii)
std::cout << z[ii] << std::endl;
return 0;
}
実際、私がそれをコンパイルするとき:
$ g++ --version
g++ (Ubuntu 4.4.3-4ubuntu5) 4.4.3
Copyright (C) 2009 Free Software Foundation, Inc.
This is free software; see the source for copying conditions. There is NO
warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
$ g++ -O0 -g crtp.cc
私が得る:
$ ./a.out
0
2
4
6
8
10
12
14
16
18
これは予想される動作です。icpcを使用する場合:
$ icpc --version
icpc (ICC) 12.1.0 20110811
Copyright (C) 1985-2011 Intel Corporation. All rights reserved.
$ icpc -g -O0 crtp.cc
代わりにを取得しSegmentation fault
ます。ランニング
valgrind --tool=memcheck ./a.out
ソースの29行目を指している
AlgebraicVectorExpression<AlgebraicVector>::operator AlgebraicVector const&() const (crtp.cc:29)
私はC++にまったく慣れておらず、結果が得られずにバグを探すのにかなりの時間を費やしたので、この問題が私が導入したエラーによるものかどうかを理解するために、より経験豊富な人の意見を聞きたいと思います。コンパイラのバグに。
編集:Mike Seymourからの回答の後、コードをそのまま変更しました。これでコンパイラの警告は表示されませんが、以前と同じ動作が得られます(同じvalgrind応答で)。誰かがIntelでコンパイルしようとしましたか?
編集:ウィキペディアの式テンプレートページでコードをコンパイルしようとしました。提供した例とまったく同じ動作が得られました。
編集:私は問題をさらに調査しました、そしてインテルと一緒にコンパイルすることicpc
はオペレーターのようです
operator const Derived&() const {
return static_cast< const Derived& >(*this);
}
再帰的に自分自身を呼び出します。私が見つけた回避策の1つは、この演算子をメソッドに置き換えることです。
const Derived& get_ref() const {
return static_cast< const Derived& >(*this);
}
それに応じて、さまざまなクラスのコンストラクターを変更します。誰かがこれらの2つの振る舞いのどちらが正しいのか、それを説明するための標準を指している可能性があるかを知ることができますか?