As we are doing PR's we should add notes to this issue regarding parts of the code we may later optimize through numba or other means. * [ ] @gbrener and I discussed this for PR #39