아이겐 : 나는 <code>Eigen</code>에서 작성하려는이 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::MatrixXcf에 std::sqrt(v)을 적용 할 수 없습니다) 분명히 잘못된 것입니다,하지만 난 그것을 작성하는 방법을 알아낼하지 않았다

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