hemm: Hermitian matrix multiply

\(C = \alpha A B + \beta C\) or \(C = \alpha B A + \beta C\) where \(A\) is Hermitian More...

Functions
template<typename scalar_t >
void	slate::hemm (Side side, scalar_t alpha, HermitianMatrix< scalar_t > &A, Matrix< scalar_t > &B, scalar_t beta, Matrix< scalar_t > &C, Options const &opts)
	Distributed parallel Hermitian matrix-matrix multiplication.
template<typename scalar_t >
void	slate::hemmA (Side side, scalar_t alpha, HermitianMatrix< scalar_t > &A, Matrix< scalar_t > &B, scalar_t beta, Matrix< scalar_t > &C, Options const &opts)
	Distributed parallel Hermitian matrix-matrix multiplication.
template<typename scalar_t >
void	slate::hemmC (Side side, scalar_t alpha, HermitianMatrix< scalar_t > &A, Matrix< scalar_t > &B, scalar_t beta, Matrix< scalar_t > &C, Options const &opts)
	Distributed parallel Hermitian matrix-matrix multiplication.

Detailed Description

\(C = \alpha A B + \beta C\) or \(C = \alpha B A + \beta C\) where \(A\) is Hermitian

Function Documentation

hemm()

template<typename scalar_t >

void slate::hemm	(	Side	side,
		scalar_t	alpha,
		HermitianMatrix< scalar_t > &	A,
		Matrix< scalar_t > &	B,
		scalar_t	beta,
		Matrix< scalar_t > &	C,
		Options const &	opts
	)

Distributed parallel Hermitian matrix-matrix multiplication.

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 Hermitian matrix and B and C are m-by-n matrices.

Complexity (in real): \(2 m^{2} n\) flops.

Template Parameters

scalar_t

One of float, double, std::complex<float>, std::complex<double>.

Parameters

[in]	side	Whether the Hermitian matrix A appears on the left or right: Side::Left: \(C = \alpha A B + \beta C\) Side::Right: \(C = \alpha B A + \beta C\)
[in]	alpha	The scalar alpha.
[in]	A	If side = left, the m-by-m Hermitian matrix A; if side = right, the n-by-n Hermitian matrix A.
[in]	B	The m-by-n matrix B.
[in]	beta	The scalar beta.
[in,out]	C	On entry, the m-by-n matrix C. On exit, overwritten by the result \(\alpha A B + \beta C\) or \(\alpha B A + \beta C\).
[in]	opts	Additional options, as map of name = value pairs. Possible options: Option::Lookahead: Number of blocks to overlap communication and computation. lookahead >= 0. Default 1. Option::MethodHemm: Select the right routine to call. Possible values: Auto: let the routine decides [default] hemmA: select hemmA routine hemmC: select hemmC routine Option::Target: Implementation to target. Possible values: HostTask: OpenMP tasks on CPU host [default]. HostNest: nested OpenMP parallel for loop on CPU host. HostBatch: batched BLAS on CPU host. Devices: batched BLAS on GPU device.

hemmA()

template<typename scalar_t >

void slate::hemmA	(	Side	side,
		scalar_t	alpha,
		HermitianMatrix< scalar_t > &	A,
		Matrix< scalar_t > &	B,
		scalar_t	beta,
		Matrix< scalar_t > &	C,
		Options const &	opts
	)

Distributed parallel Hermitian matrix-matrix multiplication.

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 Hermitian matrix and B and C are m-by-n matrices.

Template Parameters

scalar_t

One of float, double, std::complex<float>, std::complex<double>.

Parameters

[in]	side	Whether the Hermitian matrix A appears on the left or right: Side::Left: \(C = \alpha A B + \beta C\) Side::Right: \(C = \alpha B A + \beta C\)
[in]	alpha	The scalar alpha.
[in]	A	If side = left, the m-by-m Hermitian matrix A; if side = right, the n-by-n Hermitian matrix A.
[in]	B	The m-by-n matrix B.
[in]	beta	The scalar beta.
[in,out]	C	On entry, the m-by-n matrix C. On exit, overwritten by the result \(\alpha A B + \beta C\) or \(\alpha B A + \beta C\).
[in]	opts	Additional options, as map of name = value pairs. Possible options: Option::Lookahead: Number of blocks to overlap communication and computation. lookahead >= 0. Default 1. Option::Target: Implementation to target. Possible values: HostTask: OpenMP tasks on CPU host [default]. HostNest: nested OpenMP parallel for loop on CPU host. HostBatch: batched BLAS on CPU host. Devices: batched BLAS on GPU device.

hemmC()

template<typename scalar_t >

void slate::hemmC	(	Side	side,
		scalar_t	alpha,
		HermitianMatrix< scalar_t > &	A,
		Matrix< scalar_t > &	B,
		scalar_t	beta,
		Matrix< scalar_t > &	C,
		Options const &	opts
	)

Distributed parallel Hermitian matrix-matrix multiplication.

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 Hermitian matrix and B and C are m-by-n matrices.

Template Parameters

scalar_t

One of float, double, std::complex<float>, std::complex<double>.

Parameters

[in]	side	Whether the Hermitian matrix A appears on the left or right: Side::Left: \(C = \alpha A B + \beta C\) Side::Right: \(C = \alpha B A + \beta C\)
[in]	alpha	The scalar alpha.
[in]	A	If side = left, the m-by-m Hermitian matrix A; if side = right, the n-by-n Hermitian matrix A.
[in]	B	The m-by-n matrix B.
[in]	beta	The scalar beta.
[in,out]	C	On entry, the m-by-n matrix C. On exit, overwritten by the result \(\alpha A B + \beta C\) or \(\alpha B A + \beta C\).
[in]	opts	Additional options, as map of name = value pairs. Possible options: Option::Lookahead: Number of blocks to overlap communication and computation. lookahead >= 0. Default 1. Option::Target: Implementation to target. Possible values: HostTask: OpenMP tasks on CPU host [default]. HostNest: nested OpenMP parallel for loop on CPU host. HostBatch: batched BLAS on CPU host. Devices: batched BLAS on GPU device.