function B=storeband(A)
% Usage B=storeband(A).
% Stores the n-by-n banded upper triangular matrix A in with bandwidth b, in B 
% such that A(i,j)=B((b+1)-(j-i),j) where B is b+1-by-n matrix. A can be sparse
% or full.
%

%   Copyright 2009, Sivasankaran Rajamanickam, Timothy A. Davis
%   http://www.cise.ufl.edu/research/sparse

    error('storeband is missing\n') ;