- [ ] numerical improvements (less sqrt, no branching etc.) - [ ] keeping RegionMatrix alive and avoid reallocating - [ ] pre-scale weight