Once we have a check for performance regressions (#98) it would be good to - [ ] Unify naming of mask, occupancy - [ ] Replace `(eigvects * occupancy) @ eigvects.T` with `eigvects @ diag(occupancy) @ eigvects.T` Assuming no regressions
Once we have a check for performance regressions (#98) it would be good to
(eigvects * occupancy) @ eigvects.Twitheigvects @ diag(occupancy) @ eigvects.TAssuming no regressions