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); }