Hi,

I'd like to compute the inner product of two Mapped vectors : one vector is made of T [T=double], the other is made of std::complex<T>, but Eigen does not like it. It cannot cast the result to a std::complex<T>.

Is it something planned ?

Cheers
Romain

Works for me:

Code: Select all

#include <Eigen/Dense>
using namespace Eigen;
int main()
{
  double da[10];
  std::complex<double> db[10];
  Map<VectorXd> a(da,10);
  a.setRandom();
  Map<VectorXcd> b(db,10);
  b.setRandom();
  std::complex<double> res1 = a.transpose() * b;
  std::complex<double> res2 = a.dot(b);
}


And it works! Thanks Gaël.
(I was using the res = X*B notation, but with 'dot', it works perfectly)

By the way, in 3.1.0 alpha2, the speed improvement for product of small dynamic matrices is *impressive*.
For my problems, I get a 30% improvement over 3.0.5, which was already fast.
Many thanks

Romain