NAME

     ccoomm - coordinate matrix-matrix multiply

SYNOPSIS

       SUBROUTINE CCOOMM( TRANSA, M, N, K, ALPHA, DESCRA,
      *           VAL, INDX, JNDX, NNZ,
      *           B, LDB, BETA, C, LDC, WORK, LWORK )
       INTEGER    TRANSA, M, N, K, DESCRA(5), NNZ
      *           LDB, LDC, LWORK
       INTEGER    INDX(NNZ), JNDX(NNZ)
       COMPLEX    ALPHA, BETA
       COMPLEX    VAL(NNZ), B(LDB,*), C(LDC,*), WORK(LWORK)

       SUBROUTINE CCOOMM_64( TRANSA, M, N, K, ALPHA, DESCRA,
      *           VAL, INDX, JNDX, NNZ,
      *           B, LDB, BETA, C, LDC, WORK, LWORK )
       INTEGER*8  TRANSA, M, N, K, DESCRA(5), NNZ
      *           LDB, LDC, LWORK
       INTEGER*8  INDX(NNZ), JNDX(NNZ)
       COMPLEX    ALPHA, BETA
       COMPLEX    VAL(NNZ), B(LDB,*), C(LDC,*), WORK(LWORK)

     F95 INTERFACE

       SUBROUTINE COOMM( TRANSA, M, [N], K, ALPHA, DESCRA,
      *           VAL, INDX, JNDX, NNZ, B, [LDB], BETA, C, [LDC],
      *           [WORK], [LWORK] )
       INTEGER TRANSA, M, K,  NNZ
       INTEGER, DIMENSION(:) ::   DESCRA, INDX, JNDX
       COMPLEX    ALPHA, BETA
       COMPLEX, DIMENSION(:) :: VAL
       COMPLEX, DIMENSION(:, :) ::  B, C

       SUBROUTINE COOMM_64( TRANSA, M, [N], K, ALPHA, DESCRA,
      *           VAL, INDX, JNDX, NNZ, B, [LDB], BETA, C, [LDC],
      *           [WORK], [LWORK] )
       INTEGER*8 TRANSA, M, K,  NNZ
       INTEGER*8, DIMENSION(:) ::   DESCRA, INDX, JNDX
       COMPLEX    ALPHA, BETA
       COMPLEX, DIMENSION(:) :: VAL
       COMPLEX, DIMENSION(:, :) ::  B, C

DESCRIPTION

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

      where ALPHA and BETA are scalar, C and B are dense matrices,
      A is a matrix represented in coordinate format and
      op( A )  is one  of
      op( A ) = A   or   op( A ) = A'   or   op( A ) = conjg( A' ).
                                         ( ' indicates matrix transpose)

ARGUMENTS

      TRANSA        Indicates how to operate with the sparse matrix
                      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.

      M             Number of rows in matrix A

      N             Number of columns in matrix C

      K             Number of columns in matrix A

      ALPHA         Scalar parameter

      DESCRA()      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()         scalar array of length NNZ consisting of the
                    non-zero entries of A, in any order.

      INDX()        integer array of length NNZ consisting of the
                    corresponding row indices of the entries of A.

      JNDX()        integer array of length NNZ consisting of the
                    corresponding column indices of the entries of A.

      NNZ           number of non-zero elements in A.
      B()           rectangular array with first dimension LDB.

      LDB           leading dimension of B

      BETA          Scalar parameter

      C()           rectangular array with first dimension LDC.

      LDC           leading dimension of C

      WORK()        scratch array of length LWORK. WORK is not
                    referenced in the current version.

      LWORK         length of WORK array. LWORK is not referenced
                    in the current version.

Contents

NAME

SYNOPSIS

DESCRIPTION

ARGUMENTS

SEE ALSO