add outerproduct coeff(int,int) method.
This is needed to make this expression work:
(vec1*vec2.transpose())*vec3
Gael, no objection? Seems to make sense as that's fast.

Eigen/src/Core/Product.h

diff --git a/Eigen/src/Core/Product.h b/Eigen/src/Core/Product.h
index 287fc9c..3c5c56f 100644
--- a/Eigen/src/Core/Product.h
+++ b/Eigen/src/Core/Product.h
@@ -212,6 +212,11 @@
       ei_outer_product_selector<(int(Dest::Flags)&RowMajorBit) ? RowMajor : ColMajor>::run(*this, dest, alpha);
     }
 
+    Scalar coeff(int row, int col) const
+    {
+      return this->lhs().coeff(row) * this->rhs().coeff(col);
+    }
+
   private:
     GeneralProduct& operator=(const GeneralProduct&);
 };