--- a/src/loc.rs Sat Dec 21 14:27:14 2024 -0500 +++ b/src/loc.rs Fri Dec 20 16:14:17 2024 -0500 @@ -709,9 +709,12 @@ } } -impl<F : Num, const N : usize> Linear<Loc<F, N>> for Loc<F, N> { } +impl<F : Num + AXPY, const N : usize> Linear<Loc<F, N>> for Loc<F, N> { + type LinCodomain = F; +} -impl<F : Num, const N : usize> AXPY<F, Loc<F, N>> for Loc<F, N> { +impl<F : Num + AXPY, const N : usize> AXPY for Loc<F, N> { + type Field = F; #[inline] fn axpy<I : Instance<Loc<F, N>>>(&mut self, α : F, x : I, β : F) { x.eval(|x̃| {