Follow up to #102 and #113 (comment).
We have some outstanding issues with handling optional return of the truncation error for svd_trunc and friends.
These are:
- Should
epsilon always be passed by reference, as an array or Ref? This would allow us to pass it as a reference to truncation_error!
- How to handle the
compute_error keyword argument in svd_trunc (without a !) for algorithm selection?