]>
git.proxmox.com Git - ceph.git/blob - ceph/src/boost/libs/numeric/ublas/doc/samples/matrix_vector_solve.cpp
2 // Copyright (c) 2000-2002
3 // Joerg Walter, Mathias Koch
5 // Distributed under the Boost Software License, Version 1.0. (See
6 // accompanying file LICENSE_1_0.txt or copy at
7 // http://www.boost.org/LICENSE_1_0.txt)
9 // The authors gratefully acknowledge the support of
10 // GeNeSys mbH & Co. KG in producing this work.
13 #include <boost/numeric/ublas/triangular.hpp>
14 #include <boost/numeric/ublas/io.hpp>
17 using namespace boost::numeric::ublas
;
18 matrix
<double> m (3, 3);
20 for (unsigned i
= 0; i
< (std::min
) (m
.size1 (), v
.size ()); ++ i
) {
21 for (unsigned j
= 0; j
<= i
; ++ j
)
22 m (i
, j
) = 3 * i
+ j
+ 1;
26 std::cout
<< solve (m
, v
, lower_tag ()) << std::endl
;
27 std::cout
<< solve (v
, m
, lower_tag ()) << std::endl
;