31 #ifndef SHARK_LINALG_BLAS_KERNELS_DEFAULT_TRMM_HPP 32 #define SHARK_LINALG_BLAS_KERNELS_DEFAULT_TRMM_HPP 35 #include "../../matrix_proxy.hpp" 37 namespace shark {
namespace blas {
namespace bindings {
39 template <
bool Upper,
bool Unit,
typename TriangularA,
typename MatB>
41 matrix_expression<TriangularA>
const &A,
42 matrix_expression<MatB>& B,
48 std::size_t numCols=B().size2();
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);