예제 #1
0
 eigenelements_worker_base ( matrix_view <T,Opt> the_matrix) : mat(the_matrix) { 
  if (mat.is_empty())   TRIQS_RUNTIME_ERROR<<"eigenelements_worker : the matrix is empty : matrix =  "<<mat<<"  ";
  if (!mat.is_square())   TRIQS_RUNTIME_ERROR<<"eigenelements_worker : the matrix "<<mat<<" is not square ";
  if (!mat.indexmap().is_contiguous())   TRIQS_RUNTIME_ERROR<<"eigenelements_worker : the matrix "<<mat<<" is not contiguous in memory";
  dim = mat.dim0();
  ev.resize(dim);
  lwork = 64*dim;
  work.resize(lwork);
  uplo='U';compz='N' ;
  has_run = false;
 }
예제 #2
0
파일: matrix.hpp 프로젝트: Titan-C/triqs
 /// Copy construction
 matrix_view(matrix_view const& X) : IMPL_TYPE(X.indexmap(), X.storage()) {}
예제 #3
0
파일: matrix.hpp 프로젝트: pradhanphy/TRIQS
 /// Copy construction
 matrix_view( matrix_view const & X): impl_type(X.indexmap(),X.storage()) {}