From e48a3fcffbda34566f33d4a40a9c1b66fe050b10 Mon Sep 17 00:00:00 2001
From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Mon, 8 Jul 2013 18:13:17 +0000
Subject: [PATCH] Add accessor functions that take a whole set of indices all
 at once.

git-svn-id: https://svn.dealii.org/trunk@29954 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/include/deal.II/base/tensor.h      | 32 ++++++++++++++++++++++
 deal.II/include/deal.II/base/tensor_base.h | 26 ++++++++++++++++++
 2 files changed, 58 insertions(+)

diff --git a/deal.II/include/deal.II/base/tensor.h b/deal.II/include/deal.II/base/tensor.h
index bb5a41193f..bb1e0ee704 100644
--- a/deal.II/include/deal.II/base/tensor.h
+++ b/deal.II/include/deal.II/base/tensor.h
@@ -143,6 +143,16 @@ public:
    */
   const Tensor<rank_-1,dim,Number> &operator [] (const unsigned int i) const;
 
+  /**
+   * Read access using TableIndices <tt>indices</tt>
+   */
+  Number operator [](const TableIndices<rank_> & indices) const; 
+  
+  /**
+   * Read and write access using TableIndices <tt>indices</tt>
+   */
+  Number &operator [](const TableIndices<rank_> & indices); 
+  
   /**
    *  Assignment operator.
    */
@@ -405,7 +415,29 @@ Tensor<rank_,dim,Number>::operator[] (const unsigned int i) const
   return subtensor[i];
 }
 
+template <int rank_, int dim, typename Number>
+inline
+Number 
+Tensor<rank_,dim,Number>::operator[] (const TableIndices<rank_> & indices) const
+{
+    const unsigned int outer_ind = indices[rank_-1];
+	TableIndices<rank_-1> indices1;
+	for (unsigned int i = 0; i < rank_-1;i++)
+       indices1[i] = indices[i];
+	return (subtensor[outer_ind])[indices1];	
+}
 
+template <int rank_, int dim, typename Number>
+inline
+Number & 
+Tensor<rank_,dim,Number>::operator[] (const TableIndices<rank_> & indices)
+{
+    const unsigned int outer_ind = indices[rank_-1];
+	TableIndices<rank_-1> indices1;
+	for (unsigned int i = 0; i < rank_-1;i++)
+       indices1[i] = indices[i];
+	return (subtensor[outer_ind])[indices1];	
+}
 
 template <int rank_, int dim, typename Number>
 inline
diff --git a/deal.II/include/deal.II/base/tensor_base.h b/deal.II/include/deal.II/base/tensor_base.h
index 6a56daa4dd..fe44dc18b0 100644
--- a/deal.II/include/deal.II/base/tensor_base.h
+++ b/deal.II/include/deal.II/base/tensor_base.h
@@ -434,6 +434,16 @@ public:
    * backcompatibility.
    */
   Number &operator [] (const unsigned int index);
+  
+  /**
+   * Read access using TableIndices <tt>indices</tt>
+   */
+  Number operator [](const TableIndices<1> & indices) const; 
+  
+  /**
+   * Read and write access using TableIndices <tt>indices</tt>
+   */
+  Number &operator [](const TableIndices<1> & indices); 
 
   /**
    * Assignment operator.
@@ -982,6 +992,22 @@ Number &Tensor<1,dim,Number>::operator [] (const unsigned int index)
   return values[index];
 }
 
+template <int dim, typename Number>
+inline
+Number Tensor<1,dim,Number>::operator [] (const TableIndices<1> & indices) const
+{
+  Assert (indices[0]<dim, ExcIndexRange (indices[0], 0, dim));
+  return values[indices[0]];
+}
+
+template <int dim, typename Number>
+inline
+Number &Tensor<1,dim,Number>::operator [] (const TableIndices<1> & indices) 
+{
+  Assert (indices[0]<dim, ExcIndexRange (indices[0], 0, dim));
+  return values[indices[0]];
+}
+
 
 
 template <>
-- 
2.39.5