c-CRTP：表达式模板与编译器有关的问题

我发生了以下代码(存储在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

我获得了细分错误.运行

valgrind --tool=memcheck ./a.out

指向源代码中的第29行

AlgebraicVectorExpression<AlgebraicVector>::operator AlgebraicVector const&() const (crtp.cc:29)

由于我是C语言的新手,我花了相当长的时间查找没有任何结果的错误,因此我想请一个更有经验的人发表意见,以了解此问题是否是由于我引入的错误(如我预期)或到编译器错误.

编辑：
根据Mike Seymour的回答,我现在更改了代码.现在我没有得到编译器警告,但是我仍然得到与以前相同的行为(具有相同的valgrind响应).是否有人尝试与Intel进行编译？

编辑：
我试图在Wikipedia的Expression Templates页中编译代码.我获得的行为与我提供的示例相同.

编辑：
我已经对问题进行了进一步调查,并且似乎可以使用Intel icpc进行编译

operator const Derived&() const {
    return static_cast< const Derived& >(*this);
  }

递归调用自己.我发现的一种解决方法是将该操作符替换为一种方法：

const Derived& get_ref() const {
    return static_cast< const Derived& >(*this);
  }

并相应地修改各个类的构造函数.有人能指出这两种行为中的哪一种是正确的吗？

解决方法:

您应该始终启用编译器警告.他们通常可以发现细微的问题.在这种情况下：

g++ -Wall -Wextra test.cpp
test.cpp: In member function ‘const typename AlgebraicVectorExpression<AlgebraicVectorSum<T1, T2> >::ValueType& AlgebraicVectorSum<T1, T2>::operator[](typename AlgebraicVectorExpression<AlgebraicVectorSum<T1, T2> >::SizeType) const [with T1 = AlgebraicVector, T2 = AlgebraicVector]’:
test.cpp:90:   instantiated from ‘AlgebraicVector::AlgebraicVector(const AlgebraicVectorExpression<T1>&) [with T = AlgebraicVectorSum<AlgebraicVector, AlgebraicVector>]’
test.cpp:103:   instantiated from here
test.cpp:52: warning: returning reference to temporary

这告诉您问题：

const ValueType& operator[](SizeType ii) const {
    return (a_[ii] + b_[ii]);
}

表达式的结果是一个临时的,在该行的末尾被销毁,因此该函数将一个悬空的引用返回到不存在的对象.该运算符将不得不按值返回,并且您不应该实现非常量重载,因为没有要修改的值.

标签：crtp,c,expression-templates
来源： https://codeday.me/bug/20191009/1881758.html