Eig computational routines, packed

Modules
	{hp,sp}trd: reduction to tridiagonal
	{up,op}gtr: generate Q from hetrd
	{up,op}mtr: multiply by Q from hptrd

Detailed Description