cds_mat.m
View Options

			function [mata, diag] = cds_mat(file, dset)
			n=double(hdf5read(file,dset, 'RANK'));
			dists=double(hdf5read(file, strcat(dset,'/dists')));
			val=hdf5read(file, strcat(dset,'/vals'));

			%% Shift the off-diagonals %%
			for k=1:length(dists)
			d=dists(k);
			if d < 0
			val(1:n+d,k) = val(1-d:n,k);
			elseif d > 0
			val(n:-1:d+1,k) = val(n-d:-1:1,k);
			end
			end
			mata = spdiags(val, dists, n,n);
			if nargout == 2
			idiag = find(dists==0);
			diag = val(:,idiag);
			end

R8158/matlab/cds_mat.mcb8c37ade3f4master