- [ ] Inject factory function from purify into [sopt](https://github.com/astro-informatics/sopt/blob/937bfbfc9df4dc8939369391183b444e4b550deb/cpp/sopt/forward_backward.h#L266) to construct measurement operator each step of iteration - [ ] Investigate how to construct partial measurement operators - Random sampling of data is done at the beginning of the iteration loop. - Sub-operators should be implemented (by Luke) somewhere in the primal dual algorithm (maybe in `fly_operators.h`)