implicit conversion to scalar for inner product

1 files changed, 7 insertions, 2 deletions

diff --git a/test/product.h b/test/product.h
index 277b73c45..71dc4bde2 100644
--- a/test/product.h
+++ b/test/product.h
@@ -71,8 +71,9 @@ template<typename MatrixType> void product(const MatrixType& m)
 
   Scalar s1 = ei_random<Scalar>();
 
-  int r = ei_random<int>(0, rows-1),
-      c = ei_random<int>(0, cols-1);
+  int r  = ei_random<int>(0, rows-1),
+      c  = ei_random<int>(0, cols-1),
+      c2 = ei_random<int>(0, cols-1);
 
   // begin testing Product.h: only associativity for now
   // (we use Transpose.h but this doesn't count as a test for it)
@@ -150,4 +151,8 @@ template<typename MatrixType> void product(const MatrixType& m)
   {
     VERIFY(areNotApprox(res2,square2 + m2.transpose() * m1));
   }
+
+  // inner product
+  Scalar x = square2.row(c) * square2.col(c2);
+  VERIFY_IS_APPROX(x, square2.row(c).transpose().cwiseProduct(square2.col(c2)).sum());
 }