Remove generic type from precision for simpler use

Author: bytesnake

src/lobpcg/eig.rs

@@ -20,14 +20,14 @@ pub struct TruncatedEig<A: Scalar> {
     problem: Array2<A>,
     pub constraints: Option<Array2<A>>,
     preconditioner: Option<Array2<A>>,
-    precision: A::Real,
+    precision: f32,
     maxiter: usize,
 }
 
 impl<A: Float + Scalar + ScalarOperand + Lapack + PartialOrd + Default> TruncatedEig<A> {
     pub fn new(problem: Array2<A>, order: Order) -> TruncatedEig<A> {
         TruncatedEig {
-            precision: NumCast::from(1e-5).unwrap(),
+            precision: 1e-5,
             maxiter: problem.len_of(Axis(0)) * 2,
             preconditioner: None,
             constraints: None,
@@ -36,7 +36,7 @@ impl<A: Float + Scalar + ScalarOperand + Lapack + PartialOrd + Default> Truncate
         }
     }
 
-    pub fn precision(mut self, precision: A::Real) -> Self {
+    pub fn precision(mut self, precision: f32) -> Self {
         self.precision = precision;
 
         self

src/lobpcg/lobpcg.rs

@@ -146,7 +146,7 @@ pub fn lobpcg<
     mut x: Array2<A>,
     m: G,
     y: Option<Array2<A>>,
-    tol: A::Real,
+    tol: f32,
     maxiter: usize,
     order: Order,
 ) -> LobpcgResult<A> {
@@ -166,6 +166,7 @@ pub fn lobpcg<
 
     // cap the number of iteration
     let mut iter = usize::min(n * 10, maxiter);
+    let tol = NumCast::from(tol).unwrap();
 
     // calculate cholesky factorization of YY' and apply constraints to initial guess
     let cholesky_yy = y.as_ref().map(|y| {

src/lobpcg/svd.rs

@@ -94,21 +94,21 @@ impl<A: Float + PartialOrd + DivAssign<A> + 'static + MagnitudeCorrection> Trunc
 pub struct TruncatedSvd<A: Scalar> {
     order: Order,
     problem: Array2<A>,
-    precision: A::Real,
+    precision: f32,
     maxiter: usize,
 }
 
 impl<A: Float + Scalar + ScalarOperand + Lapack + PartialOrd + Default> TruncatedSvd<A> {
     pub fn new(problem: Array2<A>, order: Order) -> TruncatedSvd<A> {
         TruncatedSvd {
-            precision: NumCast::from(1e-5).unwrap(),
+            precision: 1e-5,
             maxiter: problem.len_of(Axis(0)) * 2,
             order,
             problem,
         }
     }
 
-    pub fn precision(mut self, precision: A::Real) -> Self {
+    pub fn precision(mut self, precision: f32) -> Self {
         self.precision = precision;
 
         self