Bug fix in using reconstruct_eigenvectors

termoshtt · termoshtt · commit a94bd45e0e33 · 2022-09-24T15:57:16.000+09:00
diff --git a/lax/src/eig.rs b/lax/src/eig.rs
@@ -357,11 +357,11 @@ macro_rules! impl_eig_work_r {
 
                 if let Some(v) = self.vr_l.as_ref() {
                     let v = unsafe { v.slice_assume_init_ref() };
-                    reconstruct_eigenvectors(false, eigs_im, v, self.vc_l.as_mut().unwrap());
+                    reconstruct_eigenvectors(true, eigs_im, v, self.vc_l.as_mut().unwrap());
                 }
                 if let Some(v) = self.vr_r.as_ref() {
                     let v = unsafe { v.slice_assume_init_ref() };
-                    reconstruct_eigenvectors(false, eigs_im, v, self.vc_l.as_mut().unwrap());
+                    reconstruct_eigenvectors(false, eigs_im, v, self.vc_r.as_mut().unwrap());
                 }
 
                 Ok(EigRef {
@@ -414,11 +414,11 @@ macro_rules! impl_eig_work_r {
 
                 if let Some(v) = self.vr_l.as_ref() {
                     let v = unsafe { v.slice_assume_init_ref() };
-                    reconstruct_eigenvectors(false, eigs_im, v, self.vc_l.as_mut().unwrap());
+                    reconstruct_eigenvectors(true, eigs_im, v, self.vc_l.as_mut().unwrap());
                 }
                 if let Some(v) = self.vr_r.as_ref() {
                     let v = unsafe { v.slice_assume_init_ref() };
-                    reconstruct_eigenvectors(false, eigs_im, v, self.vc_l.as_mut().unwrap());
+                    reconstruct_eigenvectors(false, eigs_im, v, self.vc_r.as_mut().unwrap());
                 }
 
                 Ok(Eig {

Original file line number	Diff line number	Diff line change
`@@ -357,11 +357,11 @@ macro_rules! impl_eig_work_r {`
`357`	`357`
`358`	`358`	`if let Some(v) = self.vr_l.as_ref() {`
`359`	`359`	`let v = unsafe { v.slice_assume_init_ref() };`
`360`		`- reconstruct_eigenvectors(false, eigs_im, v, self.vc_l.as_mut().unwrap());`
	`360`	`+ reconstruct_eigenvectors(true, eigs_im, v, self.vc_l.as_mut().unwrap());`
`361`	`361`	`}`
`362`	`362`	`if let Some(v) = self.vr_r.as_ref() {`
`363`	`363`	`let v = unsafe { v.slice_assume_init_ref() };`
`364`		`- reconstruct_eigenvectors(false, eigs_im, v, self.vc_l.as_mut().unwrap());`
	`364`	`+ reconstruct_eigenvectors(false, eigs_im, v, self.vc_r.as_mut().unwrap());`
`365`	`365`	`}`
`366`	`366`
`367`	`367`	`Ok(EigRef {`
`@@ -414,11 +414,11 @@ macro_rules! impl_eig_work_r {`
`414`	`414`
`415`	`415`	`if let Some(v) = self.vr_l.as_ref() {`
`416`	`416`	`let v = unsafe { v.slice_assume_init_ref() };`
`417`		`- reconstruct_eigenvectors(false, eigs_im, v, self.vc_l.as_mut().unwrap());`
	`417`	`+ reconstruct_eigenvectors(true, eigs_im, v, self.vc_l.as_mut().unwrap());`
`418`	`418`	`}`
`419`	`419`	`if let Some(v) = self.vr_r.as_ref() {`
`420`	`420`	`let v = unsafe { v.slice_assume_init_ref() };`
`421`		`- reconstruct_eigenvectors(false, eigs_im, v, self.vc_l.as_mut().unwrap());`
	`421`	`+ reconstruct_eigenvectors(false, eigs_im, v, self.vc_r.as_mut().unwrap());`
`422`	`422`	`}`
`423`	`423`
`424`	`424`	`Ok(Eig {`