trmm.hpp
Go to the documentation of this file.
1 /*!
2  *
3  *
4  * \brief -
5  *
6  * \author O. Krause
7  * \date 2010
8  *
9  *
10  * \par Copyright 1995-2015 Shark Development Team
11  *
12  * <BR><HR>
13  * This file is part of Shark.
14  * <http://image.diku.dk/shark/>
15  *
16  * Shark is free software: you can redistribute it and/or modify
17  * it under the terms of the GNU Lesser General Public License as published
18  * by the Free Software Foundation, either version 3 of the License, or
19  * (at your option) any later version.
20  *
21  * Shark is distributed in the hope that it will be useful,
22  * but WITHOUT ANY WARRANTY; without even the implied warranty of
23  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
24  * GNU Lesser General Public License for more details.
25  *
26  * You should have received a copy of the GNU Lesser General Public License
27  * along with Shark. If not, see <http://www.gnu.org/licenses/>.
28  *
29  */
30 
31 #ifndef SHARK_LINALG_BLAS_KERNELS_DEFAULT_TRMM_HPP
32 #define SHARK_LINALG_BLAS_KERNELS_DEFAULT_TRMM_HPP
33 
35 #include "../../matrix_proxy.hpp"
36 
37 namespace shark { namespace blas { namespace bindings {
38 
39 template <bool Upper,bool Unit,typename TriangularA, typename MatB>
40 void trmm(
41  matrix_expression<TriangularA> const &A,
42  matrix_expression<MatB>& B,
43  boost::mpl::false_
44 ) {
45  SIZE_CHECK(A().size1() == A().size2());
46  SIZE_CHECK(A().size1() == B().size1());
47 
48  std::size_t numCols=B().size2();
49 
50  for(std::size_t i = 0; i != numCols; ++i){
51  matrix_column<MatB> col = column(B,i);
52  kernels::trmv<Upper,Unit>(A,col);
53  }
54 }
55 
56 }}}
57 
58 #endif