   _ _ _         _
      X = op(T)'*X*op(T) + op(T)'*X*G*X*op(T) + Q,               (3)
        _           _               _
  where X = U'*X*U, Q = U'*Q*U, and G = U'*G*U, with U the
  orthogonal matrix reducing Ac to a real Schur form, T = U'*Ac*U.

  The routine estimates the quantities

  sepd(op(Ac),op(Ac)') = 1 / norm(inv(Omega)),

  norm(Theta) and norm(Pi) using 1-norm condition estimator.

  The forward error bound is estimated using a practical error bound
  similar to the one proposed in [2].

<