M-Labs
/
nalgebra
3
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Added comment on logic

This commit is contained in:
metric-space 2022-02-04 00:09:29 -05:00
parent 55fdd84e1d
commit 4362f0004c
1 changed files with 5 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
nalgebra-lapack/src/generalized_eigenvalues.rs
View File

@ -194,6 +194,11 @@ where
/// ///
/// u(j)**H * A = lambda(j) * u(j)**H * B . /// u(j)**H * A = lambda(j) * u(j)**H * B .
/// where u(j)**H is the conjugate-transpose of u(j). /// where u(j)**H is the conjugate-transpose of u(j).
///
/// What is going on below?
/// If the j-th and (j+1)-th eigenvalues form a complex conjugate pair,
/// then v_l(j) = VSL(:,j)+i*VSL(:,j+1) and v_l(j+1) = VSL(:,j)-i*VSL(:,j+1).
/// and then v_r(j) = VSR(:,j)+i*VSR(:,j+1) and v_r(j+1) = VSR(:,j)-i*VSR(:,j+1).
pub fn eigenvectors(self) -> (OMatrix<Complex<T>, D, D>, OMatrix<Complex<T>, D, D>) pub fn eigenvectors(self) -> (OMatrix<Complex<T>, D, D>, OMatrix<Complex<T>, D, D>)
where where
DefaultAllocator: DefaultAllocator: