There are two prototypes of symm
            available, please see below. 
symm( const Side side, const Scalar >, const MatrixA& a, const MatrixB& b, const Scalar >, MatrixC& c );
symm( const Side side, const Scalar alpha, const MatrixA& a, const MatrixB& b, const Scalar beta, MatrixC& c );
            symm (short for symmetric
            matrix-matrix operation) provides a C++ interface to BLAS routines SSYMM,
            DSYMM, CSYMM, and ZSYMM. symm
            performs one of the matrix-matrix operations
          
C := alpha*A*B + beta*C,
or
C := alpha*B*A + beta*C,
where alpha and beta are scalars, A is a symmetric matrix and B and C are m by n matrices.
            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.90. Dispatching of symm.
| Value type of MatrixA | BLAS routine | CBLAS routine | CUBLAS routine | 
|---|---|---|---|
| 
                       | SSYMM | cblas_ssymm | cublasSsymm | 
| 
                       | DSYMM | cblas_dsymm | cublasDsymm | 
| 
                       | CSYMM | cblas_csymm | cublasCsymm | 
| 
                       | ZSYMM | cblas_zsymm | Unavailable | 
            The original routines SSYMM, DSYMM, CSYMM, and ZSYMM have twelve arguments,
            whereas symm requires
            six arguments.
          
            Defined in header boost/numeric/bindings/blas/level3/symm.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/symm.hpp> using namespace boost::numeric::bindings; blas::symm( x, y, z );
this will output
[5] 0 1 2 3 4 5