update gemm

Author: sarah

Cargo.toml

@@ -21,7 +21,7 @@ dyn-stack = { version = "0.10.0", default-features = false }
 equator = "0.2.2"
 faer-entity = { version ="0.18.0", default-features = false, path = "./faer-entity" }
 
-gemm = { version = "0.17.1", default-features = false }
+gemm = { version = "0.18.0", default-features = false }
 nano-gemm = { version = "0.1.2", default-features = false }
 num-complex = { version = "0.4.5", default-features = false }
 num-traits = { version = "0.2.18", default-features = false }
@@ -36,7 +36,6 @@ npyz = { version = "0.8", optional = true }
 rand = { version = "0.8.5", default-features = false, optional = true }
 rand_distr = { version = "0.4.3", default-features = false, optional = true }
 libm = "0.2.8"
-likely_stable = "0.1.2"
 
 [features]
 default = ["std", "rayon", "serde", "rand", "npy"]
@@ -95,10 +94,6 @@ harness = false
 name = "meanvar"
 harness = false
 
-[[bench]]
-name = "sparse_solve"
-harness = false
-
 [[bench]]
 name = "bench"
 harness = false

benches/qr.rs

@@ -51,7 +51,10 @@ where
 
 fn main() -> std::io::Result<()> {
     let mut bench = Bench::new(BenchConfig::from_args()?);
-    let args = [1, 2, 4, 8, 16, 24, 32, 48, 64, 128].map(PlotArg);
+    let args = [
+        1, 2, 4, 8, 16, 24, 32, 48, 64, 128, 256, 512, 1024, 2048, 4096,
+    ]
+    .map(PlotArg);
     bench.register_many(list![faer::<f32>, nalgebra::<f32>], args);
     bench.register_many(list![faer::<f64>, nalgebra::<f64>], args);
     bench.run()?;

benches/sparse_solve.rs

@@ -169,9 +169,9 @@ pub fn recommended_blocksize<E: Entity>(nrows: usize, ncols: usize) -> usize {
     } else if prod > 128 * 128 {
         32
     } else if prod > 32 * 32 {
-        16
-    } else if prod > 16 * 16 {
         8
+    } else if prod > 16 * 16 {
+        4
     } else {
         1
     })