There are two prototypes of trsm
            available, please see below. 
trsm( const Side side, const Scalar >, const MatrixA& a, MatrixB& b );
trsm( const Side side, const Scalar alpha, const MatrixA& a, MatrixB& b );
            trsm (short for TODO)
            provides a C++ interface to BLAS routines STRSM, DTRSM, CTRSM, and ZTRSM.
            trsm solves one of the
            matrix equations
          
op( A )*X = alpha*B, or X*op( A ) = alpha*B,
where alpha is a scalar, X and B are m by n matrices, A is a unit, or non-unit, upper or lower triangular matrix and op( A ) is one of
op( A ) = A or op( A ) = A' or op( A ) = conjg( A' ).
The matrix X is overwritten on B.
            The selection of the BLAS routine is done during compile-time, and is
            determined by the type of values contained in type MatrixA.
            The type of values is obtained through the value_type
            meta-function typename value_type<MatrixA>::type. Table X below illustrates to
            which specific routine this dispatching will take place.
          
Table 1.98. Dispatching of trsm.
| Value type of MatrixA | BLAS routine | CBLAS routine | CUBLAS routine | 
|---|---|---|---|
| 
                       | STRSM | cblas_strsm | cublasStrsm | 
| 
                       | DTRSM | cblas_dtrsm | cublasDtrsm | 
| 
                       | CTRSM | cblas_ctrsm | cublasCtrsm | 
| 
                       | ZTRSM | cblas_ztrsm | cublasZtrsm | 
            The original routines STRSM, DTRSM, CTRSM, and ZTRSM have eleven arguments,
            whereas trsm requires
            four arguments.
          
            Defined in header boost/numeric/bindings/blas/level3/trsm.hpp.
          
Parameters
The definition of term 1
The definition of term 2
The definition of term 3.
Definitions may contain paragraphs.
#include <boost/numeric/bindings/blas/level3/trsm.hpp> using namespace boost::numeric::bindings; blas::trsm( x, y, z );
this will output
[5] 0 1 2 3 4 5