Tensor::operator[] should always return by reference.

One of the access operators returned by value, instead of by reference.

diff --git a/doc/news/changes.h b/doc/news/changes.h
index 0e879fb5cff7958c81e35a39d5e9b88cb6280530..372805040d5800b883f73f2e447c59a152046077 100644 (file)
--- a/doc/news/changes.h
+++ b/doc/news/changes.h
@@ -198,6 +198,12 @@ inconvenience this causes.
 <a name="general"></a>
 <h3>General</h3>
 <ol>
+  <li> Fixed: Tensor::operator[] that takes TableIndices as a parameter no
+  longer returns by value, but rather by reference.
+  <br>
+  (Jean-Paul Pelteret, 2016/01/08)
+  </li>
+
   <li> New: constrained_linear_operator() and constrained_right_hand_side()
   provide a generic mechanism of applying constraints to a LinearOperator.
   A detailed explanation with example code is given in the @ref constraints

diff --git a/include/deal.II/base/tensor.h b/include/deal.II/base/tensor.h
index 125f964860b1d30f5be26a5d695dc8abc68b8a4c..7aa7c52fc50a4f2e829b7a10218464cd0528a676 100644 (file)
--- a/include/deal.II/base/tensor.h
+++ b/include/deal.II/base/tensor.h
@@ -415,7 +415,7 @@ public:
   /**
    * Read access using TableIndices <tt>indices</tt>
    */
-  Number operator [] (const TableIndices<rank_> &indices) const;
+  const Number &operator [] (const TableIndices<rank_> &indices) const;
 
   /**
    * Read and write access using TableIndices <tt>indices</tt>
@@ -887,7 +887,7 @@ Tensor<rank_,dim,Number>::operator[] (const unsigned int i) const
 
 template <int rank_, int dim, typename Number>
 inline
-Number
+const Number &
 Tensor<rank_,dim,Number>::operator[] (const TableIndices<rank_> &indices) const
 {
   Assert(dim != 0, ExcMessage("Cannot access an object of type Tensor<rank_,0,Number>"));