A computationally efficient method for solving SUR models with orthogonal regressors