  CS      (output) DOUBLE PRECISION array, dimension (2 * min(M,N))
          On exit, the first 2*min(M,N) elements of this array
          contain the cosines and sines of the symplectic Givens
          rotations used to compute the symplectic QR decomposition.

  TAU     (output) DOUBLE PRECISION array, dimension (min(M,N))
          On exit, the first min(M,N) elements of this array
          contain the scalar factors of some of the elementary
          reflectors.

</PRE>
<B>Workspace</B>
<PRE>
  DWO