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

Correction for not calculating absolurte value

This commit is contained in:
metric-space 2022-01-21 06:41:06 -05:00
parent 6a28306074
commit 7e9345d91e
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
nalgebra-lapack/src/qz.rs
View File

@ -167,6 +167,8 @@ where
(self.vsl, self.s, self.t, self.vsr) (self.vsl, self.s, self.t, self.vsr)
} }
/// computes the generalized eigenvalues /// computes the generalized eigenvalues
#[must_use] #[must_use]
pub fn eigenvalues(&self) -> OVector<Complex<T>, D> pub fn eigenvalues(&self) -> OVector<Complex<T>, D>
@ -176,20 +178,20 @@ where
let mut out = Matrix::zeros_generic(self.t.shape_generic().0, Const::<1>); let mut out = Matrix::zeros_generic(self.t.shape_generic().0, Const::<1>);
for i in 0..out.len() { for i in 0..out.len() {
out[i] = if self.beta[i].clone() < T::RealField::default_epsilon() { out[i] = if self.beta[i].clone().abs() < T::RealField::default_epsilon() {
Complex::zero() Complex::zero()
} else { } else {
let mut cr = self.alphar[i].clone(); let mut cr = self.alphar[i].clone();
let mut ci = self.alphai[i].clone(); let mut ci = self.alphai[i].clone();
let b = self.beta[i].clone(); let b = self.beta[i].clone();
if cr < T::RealField::default_epsilon() { if cr.clone().abs() < T::RealField::default_epsilon() {
cr = T::RealField::zero() cr = T::RealField::zero()
} else { } else {
cr = cr / b.clone() cr = cr / b.clone()
}; };
if ci < T::RealField::default_epsilon() { if ci.clone().abs() < T::RealField::default_epsilon() {
ci = T::RealField::zero() ci = T::RealField::zero()
} else { } else {
ci = ci / b ci = ci / b