added functions to return random orthogonal matrices in double double and
quad double precision arithmetic