--- a/src/linops.rs	Mon Dec 30 15:46:28 2024 -0500
+++ b/src/linops.rs	Tue Dec 31 10:51:32 2024 -0500
@@ -24,6 +24,9 @@
 {
     type Owned : AXPY<F, X>;
 
+    /// Computes  `βy + αx`, where `y` is `Self`.
+    fn add_mul<I : Instance<X>>(self, α : F, x : I, β : F) -> Self::Owned ;
+
     /// Computes  `y = βy + αx`, where `y` is `Self`.
     fn axpy<I : Instance<X>>(&mut self, α : F, x : I, β : F);
 
@@ -47,6 +50,7 @@
 /// Efficient in-place application for [`Linear`] operators.
 #[replace_float_literals(F::cast_from(literal))]
 pub trait GEMV<F : Num, X : Space, Y = <Self as Mapping<X>>::Codomain> : Linear<X> {
+    
     /// Computes  `y = αAx + βy`, where `A` is `Self`.
     fn gemv<I : Instance<X>>(&self, y : &mut Y, α : F, x : I, β : F);