/*! convert to _dgematrix */ inline _dgematrix dgbmatrix::to_dgematrix() const {VERBOSE_REPORT; dgematrix newmat( dgematrix(m,n).zero() ); for(long i=0; i<m; i++){ for(long j=std::max(long(0),i-kl); j<std::min(n,i+ku+1); j++){ newmat(i,j) =(*this)(i,j); } } return _(newmat); }
int load_dgematrix(dgematrix &A, istream &s){ int i,j,M,N;; s>>M; s>>N; A=dgematrix(M,N); for (i=0; i<M; i++) for(j=0; j<N; j++) s>>A(i,j); return 0; }