Using profile on the code Alex provided, it seems that the slow bits are - [ ] logistic_bisection - specifically/uniroot - [ ] rbind