21     void test_with (MP &m1, MP &m2, MP &m3) const {  in test_with()
30             initialize_matrix (m2);  in test_with()
31             m1 = m2;  in test_with()
33             m1.assign_temporary (m2);  in test_with()
35             m1.swap (m2);  in test_with()
36             std::cout << "m1.swap (m2) = " << m1 << " " << m2 << std::endl;  in test_with()
41             m1 = m2;  in test_with()
45             m2 = - m1;  in test_with()
46             std::cout << "- m1 = " << m2 << std::endl;  in test_with()
47             m2 = ublas::conj (m1);  in test_with()
48             std::cout << "conj (m1) = " << m2 << std::endl;  in test_with()
52             initialize_matrix (m2);  in test_with()
53             m3 = m1 + m2;  in test_with()
55             m3 = m1 - m2;  in test_with()
61             m2 = value_type (1.) * m1;  in test_with()
62             std::cout << "1. * m1 = " << m2 << std::endl;  in test_with()
63             m2 = t * m1;  in test_with()
64             std::cout << "N * m1 = " << m2 << std::endl;  in test_with()
66             m2 = m1 * value_type (1.);  in test_with()
67             std::cout << "m1 * 1. = " << m2 << std::endl;  in test_with()
68             m2 = m1 * t;  in test_with()
69             std::cout << "m1 * N = " << m2 << std::endl;  in test_with()
73             initialize_matrix (m2);  in test_with()
74             m2 += m1;  in test_with()
75             std::cout << "m2 += m1 = " << m2 << std::endl;  in test_with()
76             m2 -= m1;  in test_with()
77             std::cout << "m2 -= m1 = " << m2 << std::endl;  in test_with()
78             m2 = m2 + m1;  in test_with()
79             std::cout << "m2 = m2 + m1 = " << m2 << std::endl;  in test_with()
80             m2 = m2 - m1;  in test_with()
81             std::cout << "m2 = m1 - m1 = " << m2 << std::endl;  in test_with()
89             m2 = ublas::trans (m1);  in test_with()
90             std::cout << "trans (m1) = " << m2 << std::endl;  in test_with()
94             m2 = ublas::herm (m1);  in test_with()
95             std::cout << "herm (m1) = " << m2 << std::endl;  in test_with()
99             initialize_matrix (m2);  in test_with()
100             m3 = ublas::prod (m1, m2);  in test_with()
106             M m1 (N, N), m2 (N, N), m3 (N, N);  in operator ()()  local
107             test_with (m1, m2, m3);  in operator ()()
111                                    mr2 (m2, ublas::range (0, N), ublas::range (0, N)),  in operator ()()
118                                    ms2 (m2, ublas::slice (0, 1, N), ublas::slice (0, 1, N)),  in operator ()()
126             M m1 (N, N), m2 (N, N), m3 (N, N);  in operator ()()  local
127             ublas::symmetric_adaptor<M> sam1 (m1), sam2 (m2), sam3 (m3);  in operator ()()