--- a/src/loc.rs Mon Dec 30 15:46:28 2024 -0500 +++ b/src/loc.rs Tue Dec 31 10:51:32 2024 -0500 @@ -709,6 +709,12 @@ type Owned = Self; #[inline] + fn add_mul<I : Instance<Loc<F, N>>>(mut self, α : F, x : I, β : F) -> Self { + self.axpy(α, x, β); + self + } + + #[inline] fn axpy<I : Instance<Loc<F, N>>>(&mut self, α : F, x : I, β : F) { x.eval(|x̃| { if β == F::ZERO {