Пример #1
0
/*! 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);
}
Пример #2
0
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;
}