Rewrite opnorm by LAPACK

Author: termoshtt

lax/src/opnorm.rs

@@ -2,7 +2,7 @@
 
 use crate::layout::MatrixLayout;
 use cauchy::*;
-use lapacke::Layout::ColumnMajor as cm;
+use num_traits::Zero;
 
 pub use super::NormType;
 
@@ -14,18 +14,24 @@ macro_rules! impl_opnorm {
     ($scalar:ty, $lange:path) => {
         impl OperatorNorm_ for $scalar {
             unsafe fn opnorm(t: NormType, l: MatrixLayout, a: &[Self]) -> Self::Real {
-                match l {
-                    MatrixLayout::F { col, lda } => $lange(cm, t as u8, lda, col, a, lda),
-                    MatrixLayout::C { row, lda } => {
-                        $lange(cm, t.transpose() as u8, lda, row, a, lda)
-                    }
-                }
+                let m = l.lda();
+                let n = l.len();
+                let t = match l {
+                    MatrixLayout::F { .. } => t,
+                    MatrixLayout::C { .. } => t.transpose(),
+                };
+                let mut work = if matches!(t, NormType::Infinity) {
+                    vec![Self::Real::zero(); m as usize]
+                } else {
+                    Vec::new()
+                };
+                $lange(t as u8, m, n, a, m, &mut work)
             }
         }
     };
 } // impl_opnorm!
 
-impl_opnorm!(f64, lapacke::dlange);
-impl_opnorm!(f32, lapacke::slange);
-impl_opnorm!(c64, lapacke::zlange);
-impl_opnorm!(c32, lapacke::clange);
+impl_opnorm!(f64, lapack::dlange);
+impl_opnorm!(f32, lapack::slange);
+impl_opnorm!(c64, lapack::zlange);
+impl_opnorm!(c32, lapack::clange);