From 0c89aceed3e5c7c44d4666154396fe6a1f7ae469 Mon Sep 17 00:00:00 2001 From: metric-space Date: Fri, 4 Feb 2022 00:09:29 -0500 Subject: [PATCH] Added comment on logic --- nalgebra-lapack/src/generalized_eigenvalues.rs | 5 +++++ 1 file changed, 5 insertions(+) diff --git a/nalgebra-lapack/src/generalized_eigenvalues.rs b/nalgebra-lapack/src/generalized_eigenvalues.rs index 5c273e9b..4a2a293f 100644 --- a/nalgebra-lapack/src/generalized_eigenvalues.rs +++ b/nalgebra-lapack/src/generalized_eigenvalues.rs @@ -194,6 +194,11 @@ where /// /// u(j)**H * A = lambda(j) * u(j)**H * B . /// 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, D, D>, OMatrix, D, D>) where DefaultAllocator: