This is largely taken care of in #25, but not all: - [ ] scafacos - [ ] direct sum - [ ] Use `neighbor_parallel_for`