Go to main content
Oracle Developer Studio 12.5 Man Pages

Exit Print View

Updated: June 2017
 
 

sellmm (3p)

Name

sellmm - matrix multiply

Synopsis

SUBROUTINE SELLMM( TRANSA, M, N, K, ALPHA, DESCRA,
*           VAL, INDX, LDA, MAXNZ,
*           B, LDB, BETA, C, LDC, WORK, LWORK )
INTEGER    TRANSA, M, N, K, DESCRA(5), LDA, MAXNZ,
*           LDB, LDC, LWORK
INTEGER    INDX(LDA,MAXNZ)
REAL       ALPHA, BETA
REAL       VAL(LDA,MAXNZ), B(LDB,*), C(LDC,*), WORK(LWORK)

SUBROUTINE SELLMM_64( TRANSA, M, N, K, ALPHA, DESCRA,
*           VAL, INDX, LDA, MAXNZ,
*           B, LDB, BETA, C, LDC, WORK, LWORK )
INTEGER*8  TRANSA, M, N, K, DESCRA(5), LDA, MAXNZ,
*           LDB, LDC, LWORK
INTEGER*8  INDX(LDA,MAXNZ)
REAL       ALPHA, BETA
REAL       VAL(LDA,MAXNZ), B(LDB,*), C(LDC,*), WORK(LWORK)


F95 INTERFACE
SUBROUTINE ELLMM( TRANSA, M, N, K, ALPHA, DESCRA, VAL, INDX,
*        LDA, MAXNZ, B, LDB, BETA, C, LDC, WORK, LWORK)
INTEGER    TRANSA, M, K, MAXNZ
INTEGER, DIMENSION(:) ::  DESCRA
INTEGER, DIMENSION(:, :) ::    INDX
REAL    ALPHA, BETA
REAL, DIMENSION(:, :) ::  VAL, B, C

SUBROUTINE ELLMM_64( TRANSA, M, N, K, ALPHA, DESCRA, VAL, INDX,
*        LDA, MAXNZ, B, LDB, BETA, C, LDC, WORK, LWORK)
INTEGER*8    TRANSA, M, K, MAXNZ
INTEGER*8, DIMENSION(:) ::  DESCRA
INTEGER*8, DIMENSION(:, :) ::    INDX
REAL    ALPHA, BETA
REAL, DIMENSION(:, :) ::  VAL, B, C





C INTERFACE
#include <sunperf.h>

void sellmm (const int transa, const int m, const int n, const int k,
const float alpha, const int* descra, const float* val, const
int* indx, const int lda, const int maxnz, const float* b,
const int ldb, const float beta, float* c, const int ldc);

void sellmm_64 (const long transa, const long m, const long n, const
long k, const float alpha, const long* descra, const float*
val, const long* indx, const long lda, const long maxnz,
const float* b, const long ldb, const float beta, float* c,
const long ldc);

Description

Oracle Solaris Studio Performance Library                           sellmm(3P)



NAME
       sellmm - Ellpack format matrix-matrix multiply

SYNOPSIS
        SUBROUTINE SELLMM( TRANSA, M, N, K, ALPHA, DESCRA,
       *           VAL, INDX, LDA, MAXNZ,
       *           B, LDB, BETA, C, LDC, WORK, LWORK )
        INTEGER    TRANSA, M, N, K, DESCRA(5), LDA, MAXNZ,
       *           LDB, LDC, LWORK
        INTEGER    INDX(LDA,MAXNZ)
        REAL       ALPHA, BETA
        REAL       VAL(LDA,MAXNZ), B(LDB,*), C(LDC,*), WORK(LWORK)

        SUBROUTINE SELLMM_64( TRANSA, M, N, K, ALPHA, DESCRA,
       *           VAL, INDX, LDA, MAXNZ,
       *           B, LDB, BETA, C, LDC, WORK, LWORK )
        INTEGER*8  TRANSA, M, N, K, DESCRA(5), LDA, MAXNZ,
       *           LDB, LDC, LWORK
        INTEGER*8  INDX(LDA,MAXNZ)
        REAL       ALPHA, BETA
        REAL       VAL(LDA,MAXNZ), B(LDB,*), C(LDC,*), WORK(LWORK)


   F95 INTERFACE
        SUBROUTINE ELLMM( TRANSA, M, N, K, ALPHA, DESCRA, VAL, INDX,
       *        LDA, MAXNZ, B, LDB, BETA, C, LDC, WORK, LWORK)
        INTEGER    TRANSA, M, K, MAXNZ
        INTEGER, DIMENSION(:) ::  DESCRA
        INTEGER, DIMENSION(:, :) ::    INDX
        REAL    ALPHA, BETA
        REAL, DIMENSION(:, :) ::  VAL, B, C

        SUBROUTINE ELLMM_64( TRANSA, M, N, K, ALPHA, DESCRA, VAL, INDX,
       *        LDA, MAXNZ, B, LDB, BETA, C, LDC, WORK, LWORK)
        INTEGER*8    TRANSA, M, K, MAXNZ
        INTEGER*8, DIMENSION(:) ::  DESCRA
        INTEGER*8, DIMENSION(:, :) ::    INDX
        REAL    ALPHA, BETA
        REAL, DIMENSION(:, :) ::  VAL, B, C





   C INTERFACE
       #include <sunperf.h>

       void sellmm (const int transa, const int m, const int n, const int k,
                 const float alpha, const int* descra, const float* val, const
                 int* indx, const int lda, const int maxnz, const float* b,
                 const int ldb, const float beta, float* c, const int ldc);

       void sellmm_64 (const long transa, const long m, const long n, const
                 long k, const float alpha, const long* descra, const float*
                 val, const long* indx, const long lda, const long maxnz,
                 const float* b, const long ldb, const float beta, float* c,
                 const long ldc);




DESCRIPTION
       sellmm performs one of the matrix-matrix operations

                C <- alpha op(A) B + beta C

       where op( A )  is one  of

       op( A ) = A   or   op( A ) = A'   or   op( A ) = conjg( A' )
                                          ( ' indicates matrix transpose),
       A is an M-by-K sparse matrix represented in the ellpack format,
       alpha and beta  are scalars, C and B are dense matrices.


ARGUMENTS
       TRANSA(input)   TRANSA specifies the form of op( A ) to be used in
                       the matrix multiplication as follows:
                         0 : operate with matrix
                         1 : operate with transpose matrix
                         2 : operate with the conjugate transpose of matrix.
                           2 is equivalent to 1 if matrix is real.
                       Unchanged on exit.

       M(input)        On entry,  M  specifies the number of rows in
                       the matrix A. Unchanged on exit.

       N(input)        On entry,  N specifies the number of columns in
                       the matrix C. Unchanged on exit.

       K(input)        On entry,  K specifies the number of columns
                       in  the matrix A. Unchanged on exit.

       ALPHA(input)    On entry, ALPHA specifies the scalar alpha. Unchanged on exit.

       DESCRA (input)  Descriptor argument.  Five element integer array.
                       DESCRA(1) matrix structure
                         0 : general
                         1 : symmetric (A=A')
                         2 : Hermitian (A= CONJG(A'))
                         3 : Triangular
                         4 : Skew(Anti)-Symmetric (A=-A')
                         5 : Diagonal
                         6 : Skew-Hermitian (A= -CONJG(A'))
                       DESCRA(2) upper/lower triangular indicator
                         1 : lower
                         2 : upper
                       DESCRA(3) main diagonal type
                         0 : non-unit
                         1 : unit
                       DESCRA(4) Array base (NOT IMPLEMENTED)
                         0 : C/C++ compatible
                         1 : Fortran compatible
                       DESCRA(5) repeated indices? (NOT IMPLEMENTED)
                         0 : unknown
                         1 : no repeated indices

       VAL(input)      On entry, VAL is a two-dimensional LDA-by-MAXNZ array
                       such that VAL(I,:) consists of non-zero elements
                       in row I of A, padded by zero values if the row
                       contains less than MAXNZ. Unchanged on exit.

       INDX(input)     On entry, INDX  is an integer two-dimensional
                       LDA-by-MAXNZ array such that INDX(I,:)
                       consists of the column indices of the
                       nonzero elements in row I, padded by the integer
                       value I if the number of nonzeros is less than
                       MAXNZ. Unchanged on exit.

       LDA(input)      On entry,  LDA specifies the leading dimension of VAL
                       and INDX.  Unchanged on exit.

       MAXNZ(input)    On entry, MAXNZ specifies the  max number of
                       nonzeros elements per row. Unchanged on exit.

       B (input)       Array of DIMENSION ( LDB, N ).
                       Before entry with  TRANSA = 0,  the leading  k by n
                       part of the array  B  must contain the matrix  B,  otherwise
                       the leading  m by n  part of the array  B  must contain  the
                       matrix B. Unchanged on exit.

       LDB (input)     On entry, LDB specifies the first dimension of B as declared
                       in the calling (sub) program. Unchanged on exit.

       BETA (input)    On entry, BETA specifies the scalar beta. Unchanged on exit.

       C(input/output) Array of DIMENSION ( LDC, N ).
                       Before entry with  TRANSA = 0,  the leading  m by n
                       part of the array  C  must contain the matrix C,  otherwise
                       the leading  k by n  part of the array  C must contain  the
                       matrix C. On exit, the array  C  is overwritten by the  matrix
                       ( alpha*op( A )* B  + beta*C ).

       LDC (input)     On entry, LDC specifies the first dimension of C as declared
                       in the calling (sub) program. Unchanged on exit.

       WORK (is not referenced in the current version)

       LWORK (is not referenced in the current version)


SEE ALSO
       Libsunperf  SPARSE BLAS is fully parallel and compatible with NIST FOR-
       TRAN Sparse Blas but the sources are different.  Libsunperf SPARSE BLAS
       is free of bugs found in NIST FORTRAN Sparse Blas.  Besides several new
       features and routines are implemented.

       NIST FORTRAN Sparse Blas User's Guide available at:

       http://math.nist.gov/mcsd/Staff/KRemington/fspblas/

       Based on the standard proposed in

       "Document for the Basic Linear Algebra Subprograms (BLAS) Standard",
       University of Tennessee, Knoxville, Tennessee, 1996:

       http://www.netlib.org/utk/papers/sparse.ps

       The routine is designed so that it provides a possibility to use just
       one sparse matrix representation of a general matrix A for computing
       matrix-matrix multiply for another sparse matrix composed  by  trian-
       gles and/or the main diagonal of A. The full description of the feature
       for point entry formats is given in section NOTES/BUGS for the scoomm
       manpage.




3rd Berkeley Distribution         7 Nov 2015                        sellmm(3P)