---------------------------------------------------------------------------------
fixed bug in matvec method of csr_mat, occuring if stride of vector y is not one
ll_mat.shift now supports shift of a general matrix by a symmetric matrix
new matvec and matvec:transp object methods for ll_mat type
fixed bug when assigning a symmetric matrix to a block of a general matrix
delete_cols method added for ll_mat objects