Is there a reason to prefer `quat2rotmat(rpy2quat(rpy))` to `rpy2rotmat(rpy)`? If so should we change `rpy2rotmat` to use this two stage approach under the hood?