固有：私は<code>Eigen</code>で書きたい、このMathWorks社のMATLABコードを持って

Question

複雑なマトリックスの選択機能：

[V_K,D_K] = eig(K); 
d_k = diag(D_K); 
ind_k = find(d_k > 1e-8); 
d_k(ind_k) = d_k(ind_k).^(-1/2); 
K_half = V_K*diag(d_k)*V_K'; 

これまでのところ、私が書いたコードは（最後の部分は、this例以下基づいています）されます：

Eigen::EigenSolver<Eigen::MatrixXf> es (K,true); 
Eigen::MatrixXcf v = es.eigenvalues(); 
v = (v.array()).select(1/std::sqrt(v),v); 

これは（最初にすべての、あなたはEigen::MatrixXcfにstd::sqrt(v)を適用することはできません）、明らかに間違っているが、私はそれを記述する方法を見つけ出すませんでした。手伝って頂けますか？

Answer 1

あなたはにリンクされている例のように、あなたは> 1e-8に.array()（.real()？）を比較すると、行列の値を操作するEigens機能（またはunaryExpr）を使用する必要があります。本当にありがとうございましたを

v = (v.array().real() 
//   ^^^ 
//  can't compare complex numbers 

    > 1e-8) 
// ^^^ 
// Actually compare the number 
    .select(v.cwiseSqrt().cwiseInverse(), v); 
//   ^^^ 
//   Use Eigens functions or alternatively unaryExpr 

// In one line: 
v = (v.array().real() > 1e-8).select(v.cwiseSqrt().cwiseInverse(), v);