commit 6544b49e59f6f43c9183bf0d15d74267f5070137
author Gael Guennebaud <g.gael@free.fr> Mon Jul 20 13:57:55 2015 +0200
committer Gael Guennebaud <g.gael@free.fr> Mon Jul 20 13:57:55 2015 +0200
tree 140eff6e3e17fc7681dfefda759ec0a2f3591836
parent 2d93060291f3373a3325413159062a3b4d1e95e8

Generalize pow(x,e) such that x and e can be a different expression type or a scalar for either x or e. Add x.pow(e) with e an array expression.

test/array.cpp

diff --git a/test/array.cpp b/test/array.cpp
index 8cd4fd8..3d00d1d 100644
--- a/test/array.cpp
+++ b/test/array.cpp
@@ -273,9 +273,18 @@
   VERIFY_IS_APPROX(pow(m1,2), m1.square());
   VERIFY_IS_APPROX(m1.pow(3), m1.cube());
   VERIFY_IS_APPROX(pow(m1,3), m1.cube());
+  VERIFY_IS_APPROX((-m1).pow(3), -m1.cube());
+  VERIFY_IS_APPROX(pow(2*m1,3), 8*m1.cube());
 
   ArrayType exponents = ArrayType::Constant(rows, cols, RealScalar(2));
   VERIFY_IS_APPROX(Eigen::pow(m1,exponents), m1.square());
+  VERIFY_IS_APPROX(m1.pow(exponents), m1.square());
+  VERIFY_IS_APPROX(Eigen::pow(2*m1,exponents), 4*m1.square());
+  VERIFY_IS_APPROX((2*m1).pow(exponents), 4*m1.square());
+  VERIFY_IS_APPROX(Eigen::pow(m1,2*exponents), m1.square().square());
+  VERIFY_IS_APPROX(m1.pow(2*exponents), m1.square().square());
+  VERIFY_IS_APPROX(pow(m1(0,0), exponents), ArrayType::Constant(rows,cols,m1(0,0)*m1(0,0)));
+  
 
   VERIFY_IS_APPROX(m3.pow(RealScalar(0.5)), m3.sqrt());
   VERIFY_IS_APPROX(pow(m3,RealScalar(0.5)), m3.sqrt());