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; }
/// Copy construction matrix_view(matrix_view const& X) : IMPL_TYPE(X.indexmap(), X.storage()) {}
/// Copy construction matrix_view( matrix_view const & X): impl_type(X.indexmap(),X.storage()) {}