Added read_indexed_matrix routine
Added write_indexed_matrix routine
Added is_hermitian function