the symplectic diagonal
  matrix D1 = diag(D_A,D_A**(-1)).  Next, the off-diagonal blocks of
  the resulting Hamiltonian matrix are equilibrated in the 1-norm
  using the symplectic diagonal matrix D2 of the form

              ( I/rho    0   )
         D2 = (              )
              (   0    rho*I )

  where rho is a real scalar. Thus, in (2), D = D1*D2.

  2. Norm scaling (JOBSCL = '1' or 'O'):

  The norm of the matrices A and G of (1) is reduced by setting
  A := A/tau  and  G := G/(tau**2) where t