void test_with (MP &m1, MP &m2, MP &m3) const { { value_type t; // Default Construct default_construct<MP>::test (); // Copy and swap initialize_matrix (m1); initialize_matrix (m2); m1 = m2; std::cout << "m1 = m2 = " << m1 << std::endl; m1.assign_temporary (m2); std::cout << "m1.assign_temporary (m2) = " << m1 << std::endl; m1.swap (m2); std::cout << "m1.swap (m2) = " << m1 << " " << m2 << std::endl; // Unary matrix operations resulting in a matrix initialize_matrix (m1); m2 = - m1; std::cout << "- m1 = " << m2 << std::endl; m2 = ublas::conj (m1); std::cout << "conj (m1) = " << m2 << std::endl; // Binary matrix operations resulting in a matrix initialize_matrix (m1); initialize_matrix (m2); m3 = m1 + m2; std::cout << "m1 + m2 = " << m3 << std::endl; m3 = m1 - m2; std::cout << "m1 - m2 = " << m3 << std::endl; // Scaling a matrix t = N; initialize_matrix (m1); m2 = value_type (1.) * m1; std::cout << "1. * m1 = " << m2 << std::endl; m2 = t * m1; std::cout << "N * m1 = " << m2 << std::endl; initialize_matrix (m1); m2 = m1 * value_type (1.); std::cout << "m1 * 1. = " << m2 << std::endl; m2 = m1 * t; std::cout << "m1 * N = " << m2 << std::endl; // Some assignments initialize_matrix (m1); initialize_matrix (m2); m2 += m1; std::cout << "m2 += m1 = " << m2 << std::endl; m2 -= m1; std::cout << "m2 -= m1 = " << m2 << std::endl; m2 = m2 + m1; std::cout << "m2 = m2 + m1 = " << m2 << std::endl; m2 = m2 - m1; std::cout << "m2 = m2 - m1 = " << m2 << std::endl; m1 *= value_type (1.); std::cout << "m1 *= 1. = " << m1 << std::endl; m1 *= t; std::cout << "m1 *= N = " << m1 << std::endl; // Transpose initialize_matrix (m1); // Transpose of a triangular isn't triangular of the same kind std::cout << "trans (m1) = " << ublas::trans (m1) << std::endl; // Hermitian initialize_matrix (m1); // Hermitian of a triangular isn't hermitian of the same kind std::cout << "herm (m1) = " << ublas::herm (m1) << std::endl; // Matrix multiplication initialize_matrix (m1); initialize_matrix (m2); m3 = ublas::prod (m1, m2); std::cout << "prod (m1, m2) = " << m3 << std::endl; } }
void test_expression_with (MP &m1, MP &m2, MP &m3) const { value_type t; // Default Construct default_construct<MP>::test (); // Copy and swap initialize_matrix (m1); initialize_matrix (m2); m1 = m2; std::cout << "m1 = m2 = " << m1 << std::endl; m1.assign_temporary (m2); std::cout << "m1.assign_temporary (m2) = " << m1 << std::endl; m1.swap (m2); std::cout << "m1.swap (m2) = " << m1 << " " << m2 << std::endl; // Zero assignment m1 = ublas::zero_matrix<> (m1.size1 (), m1.size2 ()); std::cout << "m1.zero_matrix = " << m1 << std::endl; m1 = m2; #ifndef BOOST_NO_FUNCTION_TEMPLATE_ORDERING // Project range and slice initialize_matrix (m1); initialize_matrix (m2); project (m1, ublas::range(0,1),ublas::range(0,1)) = project (m2, ublas::range(0,1),ublas::range(0,1)); project (m1, ublas::range(0,1),ublas::range(0,1)) = project (m2, ublas::slice(0,1,1),ublas::slice(0,1,1)); project (m1, ublas::slice(2,-1,2),ublas::slice(2,-1,2)) = project (m2, ublas::slice(0,1,2),ublas::slice(0,1,2)); project (m1, ublas::slice(2,-1,2),ublas::slice(2,-1,2)) = project (m2, ublas::range(0,2),ublas::range(0,2)); std::cout << "m1 = range/slice " << m1 << std::endl; #endif // Unary matrix operations resulting in a matrix initialize_matrix (m1); m2 = - m1; std::cout << "- m1 = " << m2 << std::endl; m2 = ublas::conj (m1); std::cout << "conj (m1) = " << m2 << std::endl; // Binary matrix operations resulting in a matrix initialize_matrix (m1); initialize_matrix (m2); m3 = m1 + m2; std::cout << "m1 + m2 = " << m3 << std::endl; m3 = m1 - m2; std::cout << "m1 - m2 = " << m3 << std::endl; m3 = ublas::element_prod (m1, m2); std::cout << "element_prod (m1, m2) = " << m3 << std::endl; // Scaling a matrix t = N; initialize_matrix (m1); m2 = value_type (1.) * m1; std::cout << "1. * m1 = " << m2 << std::endl; m2 = t * m1; std::cout << "N * m1 = " << m2 << std::endl; initialize_matrix (m1); m2 = m1 * value_type (1.); std::cout << "m1 * 1. = " << m2 << std::endl; m2 = m1 * t; std::cout << "m1 * N = " << m2 << std::endl; m2 = m1 / value_type (2.); std::cout << "m1 / 2. = " << m2 << std::endl; m2 = m1 / t; std::cout << "m1 / N = " << m2 << std::endl; // Some assignments initialize_matrix (m1); initialize_matrix (m2); m2 += m1; std::cout << "m2 += m1 = " << m2 << std::endl; m2 -= m1; std::cout << "m2 -= m1 = " << m2 << std::endl; m2 = m2 + m1; std::cout << "m2 = m2 + m1 = " << m2 << std::endl; m2 = m2 - m1; std::cout << "m2 = m2 - m1 = " << m2 << std::endl; m1 *= value_type (1.); std::cout << "m1 *= 1. = " << m1 << std::endl; m1 *= t; std::cout << "m1 *= N = " << m1 << std::endl; // Transpose initialize_matrix (m1); m2 = ublas::trans (m1); std::cout << "trans (m1) = " << m2 << std::endl; // Hermitean initialize_matrix (m1); m2 = ublas::herm (m1); std::cout << "herm (m1) = " << m2 << std::endl; // Matrix multiplication initialize_matrix (m1); initialize_matrix (m2); m3 = ublas::prod (m1, m2); std::cout << "prod (m1, m2) = " << m3 << std::endl; }
void operator () (MP &m1, MP &m2, MP &m3) const { try { value_type t; // Copy and swap initialize_matrix (m1); initialize_matrix (m2); m1 = m2; std::cout << "m1 = m2 = " << m1 << std::endl; m1.assign_temporary (m2); std::cout << "m1.assign_temporary (m2) = " << m1 << std::endl; m1.swap (m2); std::cout << "m1.swap (m2) = " << m1 << " " << m2 << std::endl; // Unary matrix operations resulting in a matrix initialize_matrix (m1); m2 = - m1; std::cout << "- m1 = " << m2 << std::endl; m2 = ublas::conj (m1); std::cout << "conj (m1) = " << m2 << std::endl; // Binary matrix operations resulting in a matrix initialize_matrix (m1); initialize_matrix (m2); m3 = m1 + m2; std::cout << "m1 + m2 = " << m3 << std::endl; m3 = m1 - m2; std::cout << "m1 - m2 = " << m3 << std::endl; // Scaling a matrix t = N; initialize_matrix (m1); m2 = value_type (1.) * m1; std::cout << "1. * m1 = " << m2 << std::endl; m2 = t * m1; std::cout << "N * m1 = " << m2 << std::endl; initialize_matrix (m1); m2 = m1 * value_type (1.); std::cout << "m1 * 1. = " << m2 << std::endl; m2 = m1 * t; std::cout << "m1 * N = " << m2 << std::endl; // Some assignments initialize_matrix (m1); initialize_matrix (m2); #ifdef BOOST_UBLAS_USE_ET m2 += m1; std::cout << "m2 += m1 = " << m2 << std::endl; m2 -= m1; std::cout << "m2 -= m1 = " << m2 << std::endl; #else m2 = m2 + m1; std::cout << "m2 += m1 = " << m2 << std::endl; m2 = m2 - m1; std::cout << "m2 -= m1 = " << m2 << std::endl; #endif m1 *= value_type (1.); std::cout << "m1 *= 1. = " << m1 << std::endl; m1 *= t; std::cout << "m1 *= N = " << m1 << std::endl; // Transpose initialize_matrix (m1); m2 = ublas::trans (m1); std::cout << "trans (m1) = " << m2 << std::endl; // Hermitean initialize_matrix (m1); m2 = ublas::herm (m1); std::cout << "herm (m1) = " << m2 << std::endl; // Matrix multiplication initialize_matrix (m1); initialize_matrix (m2); m3 = ublas::prod (m1, m2); std::cout << "prod (m1, m2) = " << m3 << std::endl; } catch (std::exception &e) { std::cout << e.what () << std::endl; } catch (...) { std::cout << "unknown exception" << std::endl; } }
void test_with (MP &m1, MP &m2, MP &m3) const { { value_type t; // Copy and swap initialize_matrix (m1); initialize_matrix (m2); m1 = m2; std::cout << "m1 = m2 = " << m1 << std::endl; m1.assign_temporary (m2); std::cout << "m1.assign_temporary (m2) = " << m1 << std::endl; m1.swap (m2); std::cout << "m1.swap (m2) = " << m1 << " " << m2 << std::endl; // Zero assignment m1 = ublas::zero_matrix<value_type> (m1.size1 (), m1.size2 ()); std::cout << "m1.zero_matrix = " << m1 << std::endl; m1 = m2; // Unary matrix operations resulting in a matrix initialize_matrix (m1); m2 = - m1; std::cout << "- m1 = " << m2 << std::endl; m2 = ublas::conj (m1); std::cout << "conj (m1) = " << m2 << std::endl; // Binary matrix operations resulting in a matrix initialize_matrix (m1); initialize_matrix (m2); m3 = m1 + m2; std::cout << "m1 + m2 = " << m3 << std::endl; m3 = m1 - m2; std::cout << "m1 - m2 = " << m3 << std::endl; // Scaling a matrix t = N; initialize_matrix (m1); m2 = value_type (1.) * m1; std::cout << "1. * m1 = " << m2 << std::endl; m2 = t * m1; std::cout << "N * m1 = " << m2 << std::endl; initialize_matrix (m1); m2 = m1 * value_type (1.); std::cout << "m1 * 1. = " << m2 << std::endl; m2 = m1 * t; std::cout << "m1 * N = " << m2 << std::endl; // Some assignments initialize_matrix (m1); initialize_matrix (m2); m2 += m1; std::cout << "m2 += m1 = " << m2 << std::endl; m2 -= m1; std::cout << "m2 -= m1 = " << m2 << std::endl; m2 = m2 + m1; std::cout << "m2 = m2 + m1 = " << m2 << std::endl; m2 = m2 - m1; std::cout << "m2 = m1 - m1 = " << m2 << std::endl; m1 *= value_type (1.); std::cout << "m1 *= 1. = " << m1 << std::endl; m1 *= t; std::cout << "m1 *= N = " << m1 << std::endl; // Transpose initialize_matrix (m1); m2 = ublas::trans (m1); std::cout << "trans (m1) = " << m2 << std::endl; // Hermitean initialize_matrix (m1); m2 = ublas::herm (m1); std::cout << "herm (m1) = " << m2 << std::endl; // Matrix multiplication initialize_matrix (m1); initialize_matrix (m2); m3 = ublas::prod (m1, m2); std::cout << "prod (m1, m2) = " << m3 << std::endl; } }