dgefmm.l

C++编写的高性能矩阵乘法的Stranssen算法

.TH DGEFMM l "12 November 1996" "PRISM version 1.0" "PRISM routine (version 1.0)"
.TH DGEFMM l "12 November 1996" "PRISM routine" "PRISM routine"
.SH NAME
DGEFMM - Winograd's variant of Strassen's algorithm to perform one of the matrix-matrix operations
.br
C := alpha*op( A )*op( B ) + beta*C.
.SH SYNOPSIS
.TP 20
SUBROUTINE DGEFMM ( TRANSA, TRANSB, M, N, K, ALPHA, A, LDA, 
B, LDB, BETA, C, LDC )
.TP 22
.ti +4
CHARACTER*1
TRANSA, TRANSB
.TP 22
.ti +4
INTEGER
M, N, K, LDA, LDB, LDC
.TP 22
.ti +4
DOUBLE PRECISION 
ALPHA, BETA
.TP 22
.ti +4
DOUBLE PRECISION 
A( LDA, * ), B( LDB, * ), C( LDC, * )
.SH PURPOSE
DGEFMM  performs one of the matrix-matrix operations

C := alpha*op( A )*op( B ) + beta*C,

where  op( X ) is one of
.br

   op( X ) = X   or   op( X ) = X' (transpose of X),
.br

alpha and beta are scalars, and A, B and C are matrices, with op( A )
an m by k matrix,  op( B )  a  k by n matrix and  C an m by n matrix.

DGEFMM is designed to be plug-and-play with the Level 3 BLAS routine
DGEMM.  It uses Winograd's variant of Strassen's algorithm to reduce
the operations count relative to conventional matrix multiplication
(DGEMM).  It has a FORTRAN interface.  All arrays are stored in column
major format.  The necessary temporary memory is dynamically malloced
inside the routine.

.SH PARAMETERS
.TP 8
TRANSA  - CHARACTER*1.
On entry, TRANSA specifies the form of op( A ) to be used in
the matrix multiplication as follows:

TRANSA = 'N' or 'n',  op( A ) = A.

TRANSA = 'T' or 't',  op( A ) = A'.

TRANSA = 'C' or 'c',  op( A ) = A'.

Unchanged on exit.
.TP 8
TRANSB  - CHARACTER*1.
On entry, TRANSB specifies the form of op( B ) to be used in
the matrix multiplication as follows:

TRANSB = 'N' or 'n',  op( B ) = B.

TRANSB = 'T' or 't',  op( B ) = B'.

TRANSB = 'C' or 'c',  op( B ) = B'.

Unchanged on exit.
.TP 8
M       - INTEGER.
On entry,  M  specifies  the number  of rows  of the  matrix
op( A )  and of the  matrix  C.  M  must  be at least  zero.
Unchanged on exit.
.TP 8
N       - INTEGER.
On entry,  N  specifies the number  of columns of the matrix
op( B ) and the number of columns of the matrix C. N must be
at least zero.
Unchanged on exit.
.TP 8
K       - INTEGER.
On entry,  K  specifies  the number of columns of the matrix
op( A ) and the number of rows of the matrix op( B ). K must
be at least  zero.
Unchanged on exit.
.TP 8
ALPHA   - DOUBLE PRECISION.
On entry, ALPHA specifies the scalar alpha.
Unchanged on exit.
.TP 8
A       - DOUBLE PRECISION array of DIMENSION ( LDA, ka ), 
where ka is
k  when  TRANSA = 'N' or 'n',  and is  m  otherwise.
Before entry with  TRANSA = 'N' or 'n',  the leading  m by k
part of the array  A  must contain the matrix  A,  otherwise
the leading  k by m  part of the array  A  must contain  the
matrix A.
Unchanged on exit.
.TP 8
LDA     - INTEGER.
On entry, LDA specifies the first dimension of A as declared
in the calling (sub) program. When  TRANSA = 'N' or 'n' then
LDA must be at least  max( 1, m ), otherwise  LDA must be at
least  max( 1, k ).
Unchanged on exit.
.TP 8
B       - DOUBLE PRECISION array of DIMENSION ( LDB, kb ), 
where kb is
n  when  TRANSB = 'N' or 'n',  and is  k  otherwise.
Before entry with  TRANSB = 'N' or 'n',  the leading  k by n
part of the array  B  must contain the matrix  B,  otherwise
the leading  n by k  part of the array  B  must contain  the
matrix B.
Unchanged on exit.
.TP 8
LDB     - INTEGER.
On entry, LDB specifies the first dimension of B as declared
in the calling (sub) program. When  TRANSB = 'N' or 'n' then
LDB must be at least  max( 1, k ), otherwise  LDB must be at
least  max( 1, n ).
Unchanged on exit.
.TP 8
BETA    - DOUBLE PRECISION.
On entry,  BETA  specifies the scalar  beta.  When  BETA  is
supplied as zero then C need not be set on input.
Unchanged on exit.
.TP 8
C       - DOUBLE PRECISION array of DIMENSION ( LDC, n ).
Before entry, the leading  m by n  part of the array  C must
contain the matrix  C,  except when  beta  is zero, in which
case C need not be set on entry.
On exit, the array  C  is overwritten by the  m by n  matrix
( alpha*op( A )*op( B ) + beta*C ).
.TP 8
LDC     - INTEGER.
On entry, LDC specifies the first dimension of C as declared
in  the  calling  (sub)  program.   LDC  must  be  at  least
max( 1, m ).
Unchanged on exit.

.SH SEE ALSO
.BR DGEFMM_MEM (l),
.BR fmm (l),
.BR TMP_DGEFMM_MEM (l),
.BR tmp_fmm (l).

.SH FURTHER INFORMATION
.LP
See the DGEFMM User's Guide provided in the doc/ directory of the distribution.

.SH BUGS
.LP
There are no known bugs.

.SH COPYRIGHT
See COPYRIGHT file provided with the software distribution.

.SH AUTHORS
.IP
Steven Huss-Lederman, University of Wisconsin
.br
Elaine M. Jacobson, Center for Computing Sciences
.br
Jeremy R. Johnson, Drexel University
.br
Anna Tsao, Center for Computing Sciences
.br
Thomas Turnbull, Center for Computing Sciences