function [b1, b2, U, V]  = piro_band_lapack(A)
% Reduce a band matrix to bidiagonal form using PIRO_BAND's lapack style
% interface.
% Usage :
% [b1, b2, U, V] = piro_band_lapack(A)
% [b1, b2, U] = piro_band_lapack(A)
% [b1, b2] = piro_band_lapack(A)
%
% or 
% [b1, b2, U] = piro_band_lapack(A, sym, uplo)
% [b1, b2] = piro_band_lapack(A, sym, uplo)
%
% sym, blocksize and uplo are optional arguments. sym=1 for symmetric matrices,
% uplo is present only when the only the lower triangular part of the input 
% symmetric matrix is stored. Upper triangular storage is assumed if uplo is
% not passed.
%

%   Copyright 2009, Sivasankaran Rajamanickam, Timothy A. Davis
%   http://www.cise.ufl.edu/research/sparse
    error('piro_band_lapack is missing\n') ;