tion
<UL>
  <LI> one general real symmetric diagonalization is performed (to
       generate orthonormal, energy weighted starting guess orbitals),
  <LI> one or two real symmetric diagonalizations (to make the
       orbital-Hessian more diagonally dominant, and to
       canonicalize the final orbitals), 
  <LI> one or more matrix exponentials are performed each
       macro-iteration, and
  <LI> no orthogonalizations are performed.
</UL>
As a result, small and in-core calculations are up to 30% improved 