使用特征：在向量和矩阵之间执行元素乘法的错误

Question

我正在尝试使用矩阵对行矢量进行元素乘法。在MATLAB中，这将简单地由“点”运算符或：

  deriv  = 1i * k .* fk;

其中k是行矢量，而fk是矩阵。 现在，在C ++中，我有此代码：

static const int nx = 10;
static const int ny = 10; 
static const int nyk = ny/2 + 1;
static const int nxk = nx/2 + 1;
static const int ncomp = 2;

Matrix <double, 1, nx> eK; 
eK.setZero();

for(int i = 0; i < nx; i++){
            eK[i] = //some expression
    }
fftw_complex *UOut; 
    UOut= (fftw_complex*) fftw_malloc((((nx)*(ny+1))*nyk)* sizeof(fftw_complex));

for (int i = 0; i < nx; i++){
        for (int j = 0; j < ny+1; j++){ 
            for (int k = 0; k < ncomp; k++){
                UOut[i*(ny+1)+j][k] = //FFT of some expression 
            }
        }
    }

Eigen::Map<Eigen::MatrixXcd, Eigen::Unaligned> U(reinterpret_cast<std::complex<double>*>(UOut),(ny+1),nx); 

现在，我正在尝试将ek的产品摄取，该产品是1 x 10的行矢量和11 x的矩阵 u 10。我尝试了几件事，似乎没有任何事情。

U = 1i * eKX.array() * euhX.array() ; //ERROR

static assertion failed: YOU_MIXED_MATRICES_OF_DIFFERENT_SIZES

( \
      |      ~~~
  176 |         (int(Eigen::internal::size_of_xpr_at_compile_time<TYPE0>::ret)==0 && int(Eigen::internal::size_of_xpr_at_compile_time<TYPE1>::ret)==0) \
      |         ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  177 |     || (\
      |     ^~~~~
  178 |           (int(TYPE0::RowsAtCompileTime)==Eigen::Dynamic \

Answer 1

EIGEN不会像Matlab或Numpy那样广播，例如使用matrix.Array（）。rowwise（） * vector.aray（）

除非您明确要求它，否则 将矢量解释为对角线矩阵。

Eigen::VectorXd eK = ...;
Eigen::Map<Eigen::MatrixXcd, Eigen::Unaligned> U = ...;

Eigen::MatrixXcd result = U * (eK * 1i).asDiagonal();