Add EigWork<T>

Author: termoshtt

lax/src/eig.rs

@@ -32,6 +32,31 @@ pub trait Eig_: Scalar {
     ) -> Result<(Vec<Self::Complex>, Vec<Self::Complex>)>;
 }
 
+/// Working memory for [Eig_]
+#[derive(Debug, Clone)]
+pub struct EigWork<T: Scalar> {
+    pub n: i32,
+    pub jobvr: JobEv,
+    pub jobvl: JobEv,
+
+    /// Eigenvalues used in complex routines
+    pub eigs: Option<Vec<MaybeUninit<T>>>,
+    /// Real part of eigenvalues used in real routines
+    pub eigs_re: Option<Vec<MaybeUninit<T>>>,
+    /// Imaginary part of eigenvalues used in real routines
+    pub eigs_im: Option<Vec<MaybeUninit<T>>>,
+
+    /// Left eigenvectors
+    pub vl: Option<Vec<MaybeUninit<T>>>,
+    /// Right eigenvectors
+    pub vr: Option<Vec<MaybeUninit<T>>>,
+
+    /// Working memory
+    pub work: Vec<MaybeUninit<T>>,
+    /// Working memory with `T::Real`
+    pub rwork: Option<Vec<MaybeUninit<T::Real>>>,
+}
+
 macro_rules! impl_eig_complex {
     ($scalar:ty, $ev:path) => {
         impl Eig_ for $scalar {