diff -r 99ad55974e62 -r 981069ef919b src/nalgebra_support.rs
--- a/src/nalgebra_support.rs	Mon Dec 30 15:46:28 2024 -0500
+++ b/src/nalgebra_support.rs	Tue Dec 31 10:51:32 2024 -0500
@@ -91,6 +91,13 @@
     type Owned = OVector<E, M>;
 
     #[inline]
+    fn add_mul<I : Instance<Vector<E,M,SV1>>>(self, α : E, x : I, β : E) -> Self::Owned {
+        let mut owned = self.into_owned();
+        x.eval(|x̃| Matrix::axpy(&mut owned, α, x̃, β));
+        owned
+    }
+
+    #[inline]
     fn axpy<I : Instance<Vector<E,M,SV1>>>(&mut self, α : E, x : I, β : E) {
         x.eval(|x̃| Matrix::axpy(self, α, x̃, β))
     }