Lev

Does anyone understand why does the very simple example below fail ? #define VIENNACL_WITH_UBLAS #include "viennacl/compressed_matrix.hpp" #include "viennacl/linalg/power_iter.hpp" using namespace boost::numeric; int main(){ double lambda; ublas::compressed_matrix<double> M(10, 10), P(10,10); viennacl::compressed_matrix<double> C(10, 10); viennacl::linalg::power_iter_tag ptag(1e-6); M(2,2) = M(2,4) = 1; viennacl::copy(M,C); lambda = viennacl::linalg::eig(C, ptag); P = M / 2; viennacl::copy(P,C);...