Developer Reference

  • 098
  • 03/30/2020
  • Public Content
Contents

p?trmm

Computes a scalar-matrix-matrix product (one matrix operand is triangular) for distributed matrices.

Syntax

call pstrmm
(
side
,
uplo
,
transa
,
diag
,
m
,
n
,
alpha
,
a
,
ia
,
ja
,
desca
,
b
,
ib
,
jb
,
descb
)
call pdtrmm
(
side
,
uplo
,
transa
,
diag
,
m
,
n
,
alpha
,
a
,
ia
,
ja
,
desca
,
b
,
ib
,
jb
,
descb
)
call pctrmm
(
side
,
uplo
,
transa
,
diag
,
m
,
n
,
alpha
,
a
,
ia
,
ja
,
desca
,
b
,
ib
,
jb
,
descb
)
call pztrmm
(
side
,
uplo
,
transa
,
diag
,
m
,
n
,
alpha
,
a
,
ia
,
ja
,
desca
,
b
,
ib
,
jb
,
descb
)
Include Files
  • mkl_pblas.h
Description
The
p?trmm
routines perform a matrix-matrix operation using triangular matrices. The operation is defined as
sub(
B
) :=
alpha
*op(sub(
A
))*sub(
B
)
or
sub(
B
) :=
alpha
*sub(
B
)*op(sub(
A
))
where:
alpha
is a scalar,
sub(
B
)
is an
m
-by-
n
distributed matrix,
sub(
B
)=
B
(
ib
:
ib
+
m
-1,
jb
:
jb
+
n
-1)
.
A
is a unit, or non-unit, upper or lower triangular distributed matrix,
sub(
A
)=
A
(
ia
:
ia
+
m
-1,
ja
:
ja
+
m
-1)
, if
side
= '
L
'
or
'
l
'
, and
sub(
A
)=
A
(
ia
:
ia
+
n
-1,
ja
:
ja
+
n
-1)
, if
side
= '
R
'
or
'
r
'
.
op(sub(
A
))
is one of
op(sub(
A
)) = sub(
A
)
, or
op(sub(
A
)) = sub(
A
)'
, or
op(sub(
A
)) = conjg(sub(
A
)')
.
Input Parameters
side
(global)
CHARACTER*1
.
Specifies whether
op(sub(
A
))
appears on the left or right of
sub(
B
)
in the operation:
if
side
=
'L'
or
'l'
, then
sub(
B
) :=
alpha
*op(sub(
A
))*sub(
B
)
;
if
side
=
'R'
or
'r'
, then
sub(
B
) :=
alpha
*sub(
B
)*op(sub(
A
))
.
uplo
(global)
CHARACTER*1
.
Specifies whether the distributed matrix
sub(
A
)
is upper or lower triangular:
if
uplo
=
'U'
or
'u'
, then the matrix is upper triangular;
if
uplo
=
'L'
or
'l'
, then the matrix is low triangular.
transa
(global)
CHARACTER*1
.
Specifies the form of
op(sub(
A
))
used in the matrix multiplication:
if
transa
= '
N
'
or
'
n
'
, then
op(sub(
A
)) = sub(
A
)
;
if
transa
= '
T
'
or
'
t
'
, then
op(sub(
A
)) = sub(
A
)'
;
if
transa
= '
C
'
or
'
c
'
, then
op(sub(
A
)) = conjg(sub(
A
)')
.
diag
(global)
CHARACTER*1
.
Specifies whether the matrix
sub(
A
)
is unit triangular:
if
diag
=
'U'
or
'u'
then the matrix is unit triangular;
if