Numerical Methods for computing SVD in the D-orthogonal group