Merge pull request #151 from FabienPean/feature/array_indexing

Allow indexing by operator[] when dimension is 1

Author: romeric

Fastor/tensor/ScalarIndexing.h

@@ -9,6 +9,11 @@ template<typename... Args, typename std::enable_if<sizeof...(Args)==dimension_t:
 FASTOR_INLINE T& operator()(Args ... args) {
     return _data[get_flat_index(args...)];
 }
+template<typename Arg, typename std::enable_if<1==dimension_t::value &&
+                                is_arithmetic_pack<Arg>::value,bool>::type =0>
+FASTOR_INLINE T& operator[](Arg arg) {
+    return _data[get_flat_index(arg)];
+}
 //----------------------------------------------------------------------------------------------------------//
 
 #endif // SCALAR_INDEXING_NONCONST_H
@@ -24,7 +29,11 @@ template<typename... Args, typename std::enable_if<sizeof...(Args)==dimension_t:
 constexpr FASTOR_INLINE const T& operator()(Args ... args) const {
     return _data[get_flat_index(args...)];
 }
-
+template<typename Arg, typename std::enable_if<1==dimension_t::value &&
+                                is_arithmetic_pack<Arg>::value,bool>::type =0>
+constexpr FASTOR_INLINE const T& operator[](Arg arg) const {
+    return _data[get_flat_index(arg)];
+}
 //----------------------------------------------------------------------------------------------------------//
 
 #endif // SCALAR_INDEXING_CONST_H