Added routines to set dense submatrices of a sparse matrix
CS and ellpack matrices have the routine for setting a dense submatrix overridden for extra performance
Started writing test for setting multiple matrix entries
Added more to test for multiple matrix entries
All matrices pass elementary test for setting multiple entries
More informative error message in test for setting multiple matrix entries
Setting multiple matrix entries with realloction works