| // |
| // Copyright (c) 2000-2002 |
| // Joerg Walter, Mathias Koch |
| // |
| // Distributed under the Boost Software License, Version 1.0. (See |
| // accompanying file LICENSE_1_0.txt or copy at |
| // http://www.boost.org/LICENSE_1_0.txt) |
| // |
| // The authors gratefully acknowledge the support of |
| // GeNeSys mbH & Co. KG in producing this work. |
| // |
| |
| #include <boost/numeric/ublas/matrix.hpp> |
| #include <boost/numeric/ublas/io.hpp> |
| |
| int main () { |
| using namespace boost::numeric::ublas; |
| matrix<double> m1 (3, 3), m2 (3, 3); |
| for (unsigned i = 0; i < (std::min) (m1.size1 (), m2.size1 ()); ++ i) |
| for (unsigned j = 0; j < (std::min) (m1.size2 (), m2.size2 ()); ++ j) |
| m1 (i, j) = m2 (i, j) = 3 * i + j; |
| |
| std::cout << prod (m1, m2) << std::endl; |
| } |
| |