int load_dcovector(dcovector &X, istream &s){
int i , N;
s>>N;
X = dcovector(N);
for(i=0;i<N;i++)
s>>X(i);
return 0;
}
/*! get column of the matrix */
inline _dcovector dgbmatrix::col(const long& _n) const
{VERBOSE_REPORT;
#ifdef CPPL_DEBUG
if( _n<0 || _n>n ){
ERROR_REPORT;
std::cerr << "Input row number must be between 0 and " << n << "." << std::endl
<< "Your input was " << _n << "." << std::endl;
exit(1);
}
#endif//CPPL_DEBUG
dcovector v( dcovector(m).zero() );
for(long i=std::max(long(0),_n-ku); i<std::min(m,_n+kl+1); i++){ v(i)=(*this)(i,_n); }
return _(v);
}
