perf: vectorized small lblt prototype

faer-traits/src/lib.rs

@@ -2394,6 +2394,46 @@ impl<T: RealField> core::ops::Mul<&ComplexImpl<T>> for &ComplexImpl<T> {
 	}
 }
 
+impl<T: RealField> core::ops::Neg for ComplexImpl<T> {
+	type Output = ComplexImpl<T>;
+
+	#[inline]
+	fn neg(self) -> Self::Output {
+		use math_utils::*;
+
+		ComplexImpl(neg(&self.0))
+	}
+}
+impl<T: RealField> core::ops::Add<ComplexImpl<T>> for ComplexImpl<T> {
+	type Output = ComplexImpl<T>;
+
+	#[inline]
+	fn add(self, rhs: ComplexImpl<T>) -> Self::Output {
+		use math_utils::*;
+
+		ComplexImpl(add(&self.0, &rhs.0))
+	}
+}
+impl<T: RealField> core::ops::Sub<ComplexImpl<T>> for ComplexImpl<T> {
+	type Output = ComplexImpl<T>;
+
+	#[inline]
+	fn sub(self, rhs: ComplexImpl<T>) -> Self::Output {
+		use math_utils::*;
+
+		ComplexImpl(sub(&self.0, &rhs.0))
+	}
+}
+impl<T: RealField> core::ops::Mul<ComplexImpl<T>> for ComplexImpl<T> {
+	type Output = ComplexImpl<T>;
+
+	#[inline]
+	fn mul(self, rhs: ComplexImpl<T>) -> Self::Output {
+		use math_utils::*;
+
+		ComplexImpl(mul(&self.0, &rhs.0))
+	}
+}
 impl<T> From<Complex<T>> for ComplexImpl<T> {
 	#[inline]
 	fn from(value: Complex<T>) -> Self {

faer/examples/lblt.rs

@@ -4,6 +4,7 @@ use diol::prelude::*;
 use dyn_stack::{MemBuffer, MemStack};
 use faer::diag::Diag;
 use faer::linalg::cholesky::lblt;
+use faer::linalg::cholesky::lblt::factor::PivotingStrategy;
 use faer::prelude::*;
 use faer::reborrow::ReborrowMut;
 use faer::stats::prelude::*;
@@ -88,7 +89,8 @@ fn faer(bencher: Bencher, PlotArg(n): PlotArg) {
 	let perm_inv = &mut *vec![0usize; n];
 
 	let par = Par::Seq;
-	let params = default();
+	let mut params: faer::Spec<faer::linalg::cholesky::lblt::factor::LbltParams, f64> = default();
+	params.pivoting = PivotingStrategy::Partial;
 
 	let scratch = lblt::factor::cholesky_in_place_scratch::<usize, f64>(n, par, params);
 	let mut mem = MemBuffer::new(scratch);
@@ -113,7 +115,7 @@ fn main() -> eyre::Result<()> {
 
 			f
 		},
-		[64, 96, 128, 192, 256, 1024, 2048, 4096, 8192].map(PlotArg),
+		[8, 16, 32, 64, 96, 128, 192, 256, 1024, 2048, 4096, 8192].map(PlotArg),
 	);
 
 	bench.run()?;

faer/src/lib.rs

@@ -166,16 +166,19 @@ extern crate alloc;
 #[cfg(feature = "std")]
 extern crate std;
 
+#[doc(hidden)]
+pub extern crate generativity;
+
 /// see: [`generativity::make_guard`]
 #[macro_export]
 macro_rules! make_guard {
     ($($name:ident),* $(,)?) => {$(
         #[allow(unused_unsafe)]
-        let $name = unsafe { extern crate generativity; ::generativity::Id::new() };
+        let $name = unsafe { $crate::generativity::Id::new() };
         #[allow(unused, unused_unsafe)]
-        let lifetime_brand = unsafe { extern crate generativity; ::generativity::LifetimeBrand::new(&$name) };
+        let lifetime_brand = unsafe {  $crate::generativity::LifetimeBrand::new(&$name) };
         #[allow(unused_unsafe)]
-        let $name = unsafe { extern crate generativity; ::generativity::Guard::new($name) };
+        let $name = unsafe {  $crate::generativity::Guard::new($name) };
     )*};
 }