Fix bug in outerproducts for tensormaps

Author: romeric

Fastor/tensor_algebra/outerproduct.h

@@ -97,7 +97,8 @@ dyadic(const AbstractTensor<Derived0,DIM0> &a, const AbstractTensor<Derived1,DIM
 //     return a;
 // }
 
-template<typename AbstractTensorType0, typename AbstractTensorType1, typename ... AbstractTensorTypes>
+template<typename AbstractTensorType0, typename AbstractTensorType1, typename ... AbstractTensorTypes,
+    enable_if_t_<is_greater_equal_v_<sizeof...(AbstractTensorTypes),1>,bool> >
 FASTOR_INLINE
 auto
 outer(const AbstractTensorType0& a, const AbstractTensorType1& b, const AbstractTensorTypes& ... rest)
@@ -106,7 +107,8 @@ outer(const AbstractTensorType0& a, const AbstractTensorType1& b, const Abstract
     return outer(res, rest...);
 }
 
-template<typename AbstractTensorType0, typename AbstractTensorType1, typename ... AbstractTensorTypes>
+template<typename AbstractTensorType0, typename AbstractTensorType1, typename ... AbstractTensorTypes,
+    enable_if_t_<is_greater_equal_v_<sizeof...(AbstractTensorTypes),1>,bool> >
 FASTOR_INLINE
 auto
 dyadic(const AbstractTensorType0& a, const AbstractTensorType1& b, const AbstractTensorTypes& ... rest)

tests/test_tensormap/test_tensormap.cpp

@@ -88,6 +88,21 @@ void run() {
         FASTOR_EXIT_ASSERT(std::abs(res1.sum() + 28 ) < Tol);
     }
 
+    // Bug 132 - add 1D tensormaps, outer tensormaps
+    {
+        std::array<double,3> a = {1,1,1};
+        TensorMap<double,3> at(a.data());
+
+        std::array<double,3> b = {2,2,2};
+        TensorMap<double,3> bt(b.data());
+
+        Tensor<double,3> ct = bt + at;
+        FASTOR_EXIT_ASSERT(std::abs(ct.sum() - 9 ) < Tol);
+
+        Tensor<double,3,3> dt = outer(at, bt);
+        FASTOR_EXIT_ASSERT(std::abs(dt.sum() - 18 ) < Tol);
+    }
+
     print(FGRN(BOLD("All tests passed successfully")));
 }