function [b1, b2, U, V, fl]  = piro_band_flops(A, sym, blks, wantflops)
% Usage [b1, b2]  = piro_band_flops(A, sym, blks, wantflops)
% Band reduction of banded matrix to upper bidiagonal form.
% Uses blocked and interleaved Given's rotations.
% piro_bandmex should be compiled with -DBENCHMARK for this function. This 
% function is for internal benchmarking purposes. Use piro_band.m for normal
% usage. 
% Usage :
% [b1, b2, U, V, fl] = piro_band_flops(A, sym, blksize, 3.14)
% [b1, b2, U, fl] = piro_band_flops(A, sym, blksize, 3.14)
% [b1, b2, fl] = piro_band_flops(A, sym, blksize, 3.14)
% When the 4th input argument is present the value of it doesn't matter.
% Without the 4th input argument the usage will remain the same as piro_band.
%
% When compiled with -DBENCHMARK for this function and when there are four
% arguments in the input, the last output argument of piro_band_flops will be
% the flops for the reduction.
%

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

    error ('piro_band is missing \n')