00001 //============================================================================= 00002 /*! convert to _zgematrix */ 00003 inline _zgematrix zgbmatrix::to_zgematrix() const 00004 {VERBOSE_REPORT; 00005 zgematrix newmat( zgematrix(m,n).zero() ); 00006 for(long i=0; i<m; i++){ 00007 for(long j=std::max(long(0),i-kl); j<std::min(n,i+ku+1); j++){ 00008 newmat(i,j) =(*this)(i,j); 00009 } 00010 } 00011 return _(newmat); 00012 }