ontinuous-time Lyapunov equation to the canonical form
         H      H              H                     2      H
    op(S) *op(V) *op(V) + op(V) *op(V)*op(S) = -scale *op(F) *op(F),

  or the discrete-time Lyapunov equation to the canonical form
         H      H                    H               2      H
    op(S) *op(V) *op(V)*op(S) - op(V) *op(V) = -scale *op(F) *op(F),

  where V and F are upper triangular, and

     F = R,  M &gt;= N,   F = ( R  Z ),  M &lt; N,  if TRANS = 'N';
                   