ute || Q^H Q - I ||_F for a complex matrix of the form

                    [  op( Q1 )  op( Q2 ) ]
               Q =  [                     ],
                    [ -op( Q2 )  op( Q1 ) ]

  where Q1 and Q2 are N-by-N matrices. This residual can be used to
  test wether Q is numerically a unitary symplectic matrix.

</PRE>
<A name="Specification"><B><FONT SIZE="+1">Specification</FONT></B></A>
<PRE>
      DOUBLE PRECISION FUNCTION MA02JZ( LTRAN1, LTRAN2, N, Q1, LDQ1, Q2,
     $                             