Added new matrix type 'DiagMatrix' for diagonal matrices. Various minor other changes in BLAS function for band and packed storage.