ion and
<LI> orthogonalization.
</UL>
To this list our quadratically convergent algorithm adds
<UL>
<LI> matrix exponentiation.
</UL>
With <A HREF="../peigs/peigs.html">PeIGS</A> and <A
HREF="http://www.netlib.org/scalapack/">ScaLAPACK</A> we now have
quite efficient matrix diagonalization and orthogonalization
routines, but these operations still take much longer than a matrix
multiplication.  This is because the latter both parallelizes more
readily and also drives a single processor at peak speed.

 <P> 