I am trying to do a simple matrix multiplication of 2 complex matrices. C = A*B
The first matrix has a complex value of 1.0+0i in rows 1 and 10 and zeros in all other rows.
The second matrix has values in all rows for column 1.
The values in row1 are -6.6617e-9+i1.92738e-9
and in row 10 I have -6.6617e-9-i1.92738e-10

These matrices are 20 by 20.

The problem I have is that in Matlab, and by hand calculation, the result for C row 1 column 1 should be -1.3324e-8+i0.0
But from the eigen package I get: -2.49229e-21+i7.206811e-23

I've tried adding the noalias (eg: C.noalias() = A*B;), but it made no difference.
My matrices were declared using the MatrixXcd. Am I missing something else?
Does anyone have any ideas on what I can try?

Thanks,
Steve

Ok, never mind. I found my problem. It wasn't an EIGEN problem. It was a problem with my data.