- [ ] Use an orthogonality convergence criterion in pnls! - [ ] Store nm.u after pnls! at initial beta and theta and use it for starting nm.u in later pnls! - [ ] Add an external constructor taking a NonlinearLS object.