BLAS / ctrmm.f

Fortran project BLAS, source module ctrmm.f.

Source module last modified on Thu, 2 Jul 1998, 23:17;
HTML image of Fortran source automatically generated by for2html on Sun, 23 Jun 2002, 15:10.


      SUBROUTINE CTRMM ( SIDE, UPLO, TRANSA, DIAG, M, N, ALPHA, A, LDA,
     $                   B, LDB )
#     .. Scalar Arguments ..
      CHARACTER*1        SIDE, UPLO, TRANSA, DIAG
      INTEGER            M, N, LDA, LDB
      COMPLEX            ALPHA
#     .. Array Arguments ..
      COMPLEX            A( LDA, * ), B( LDB, * )
#     ..
#
#  Purpose
#  =======
#
#  CTRMM  performs one of the matrix-matrix operations
#
#     B := alpha*op( A )*B,   or   B := alpha*B*op( A )
#
#  where  alpha  is a scalar,  B  is an m by n matrix,  A  is a unit, or
#  non-unit,  upper or lower triangular matrix  and  op( A )  is one  of
#
#     op( A ) = A   or   op( A ) = A'   or   op( A ) = conjg( A' ).
#
#  Parameters
#  ==========
#
#  SIDE   - CHARACTER*1.
#           On entry,  SIDE specifies whether  op( A ) multiplies B from
#           the left or right as follows:
#
#              SIDE = 'L' or 'l'   B := alpha*op( A )*B.
#
#              SIDE = 'R' or 'r'   B := alpha*B*op( A ).
#
#           Unchanged on exit.
#
#  UPLO   - CHARACTER*1.
#           On entry, UPLO specifies whether the matrix A is an upper or
#           lower triangular matrix as follows:
#
#              UPLO = 'U' or 'u'   A is an upper triangular matrix.
#
#              UPLO = 'L' or 'l'   A is a lower triangular matrix.
#
#           Unchanged on exit.
#
#  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 ) = conjg( A' ).
#
#           Unchanged on exit.
#
#  DIAG   - CHARACTER*1.
#           On entry, DIAG specifies whether or not A is unit triangular
#           as follows:
#
#              DIAG = 'U' or 'u'   A is assumed to be unit triangular.
#
#              DIAG = 'N' or 'n'   A is not assumed to be unit
#                                  triangular.
#
#           Unchanged on exit.
#
#  M      - INTEGER.
#           On entry, M specifies the number of rows of B. M must be at
#           least zero.
#           Unchanged on exit.
#
#  N      - INTEGER.
#           On entry, N specifies the number of columns of B.  N must be
#           at least zero.
#           Unchanged on exit.
#
#  ALPHA  - COMPLEX         .
#           On entry,  ALPHA specifies the scalar  alpha. When  alpha is
#           zero then  A is not referenced and  B need not be set before
#           entry.
#           Unchanged on exit.
#
#  A      - COMPLEX          array of DIMENSION ( LDA, k ), where k is m
#           when  SIDE = 'L' or 'l'  and is  n  when  SIDE = 'R' or 'r'.
#           Before entry  with  UPLO = 'U' or 'u',  the  leading  k by k
#           upper triangular part of the array  A must contain the upper
#           triangular matrix  and the strictly lower triangular part of
#           A is not referenced.
#           Before entry  with  UPLO = 'L' or 'l',  the  leading  k by k
#           lower triangular part of the array  A must contain the lower
#           triangular matrix  and the strictly upper triangular part of
#           A is not referenced.
#           Note that when  DIAG = 'U' or 'u',  the diagonal elements of
#           A  are not referenced either,  but are assumed to be  unity.
#           Unchanged on exit.
#
#  LDA    - INTEGER.
#           On entry, LDA specifies the first dimension of A as declared
#           in the calling (sub) program.  When  SIDE = 'L' or 'l'  then
#           LDA  must be at least  max( 1, m ),  when  SIDE = 'R' or 'r'
#           then LDA must be at least max( 1, n ).
#           Unchanged on exit.
#
#  B      - COMPLEX          array of DIMENSION ( LDB, n ).
#           Before entry,  the leading  m by n part of the array  B must
#           contain the matrix  B,  and  on exit  is overwritten  by the
#           transformed matrix.
#
#  LDB    - INTEGER.
#           On entry, LDB specifies the first dimension of B as declared
#           in  the  calling  (sub)  program.   LDB  must  be  at  least
#           max( 1, m ).
#           Unchanged on exit.
#
#
#  Level 3 Blas routine.
#
#  -- Written on 8-February-1989.
#     Jack Dongarra, Argonne National Laboratory.
#     Iain Duff, AERE Harwell.
#     Jeremy Du Croz, Numerical Algorithms Group Ltd.
#     Sven Hammarling, Numerical Algorithms Group Ltd.
#
#
#     .. External Functions ..
      LOGICAL            LSAME
      EXTERNAL           LSAME
#     .. External Subroutines ..
      EXTERNAL           XERBLA
#     .. Intrinsic Functions ..
      INTRINSIC          CONJG, MAX
#     .. Local Scalars ..
      LOGICAL            LSIDE, NOCONJ, NOUNIT, UPPER
      INTEGER            I, INFO, J, K, NROWA
      COMPLEX            TEMP
#     .. Parameters ..
      COMPLEX            ONE
      PARAMETER        ( ONE  = ( 1.0E+0, 0.0E+0 ) )
      COMPLEX            ZERO
      PARAMETER        ( ZERO = ( 0.0E+0, 0.0E+0 ) )
#     ..
#     .. Executable Statements ..
#
#     Test the input parameters.
#
      LSIDE  = LSAME( SIDE  , 'L' )
      IF( LSIDE )THEN
         NROWA = M
      ELSE
         NROWA = N
      END IF
      NOCONJ = LSAME( TRANSA, 'T' )
      NOUNIT = LSAME( DIAG  , 'N' )
      UPPER  = LSAME( UPLO  , 'U' )
#
      INFO   = 0
      IF(      ( ! LSIDE                )&&
     $         ( ! LSAME( SIDE  , 'R' ) )      )THEN
         INFO = 1
      ELSE IF( ( ! UPPER                )&&
     $         ( ! LSAME( UPLO  , 'L' ) )      )THEN
         INFO = 2
      ELSE IF( ( ! LSAME( TRANSA, 'N' ) )&&
     $         ( ! LSAME( TRANSA, 'T' ) )&&
     $         ( ! LSAME( TRANSA, 'C' ) )      )THEN
         INFO = 3
      ELSE IF( ( ! LSAME( DIAG  , 'U' ) )&&
     $         ( ! LSAME( DIAG  , 'N' ) )      )THEN
         INFO = 4
      ELSE IF( M  <0               )THEN
         INFO = 5
      ELSE IF( N  <0               )THEN
         INFO = 6
      ELSE IF( LDA<MAX( 1, NROWA ) )THEN
         INFO = 9
      ELSE IF( LDB<MAX( 1, M     ) )THEN
         INFO = 11
      END IF
      IF( INFO!=0 )THEN
         CALL XERBLA( 'CTRMM ', INFO )
         RETURN
      END IF
#
#     Quick return if possible.
#
      IF( N==0 )
     $   RETURN
#
#     And when  alpha.eq.zero.
#
      IF( ALPHA==ZERO )THEN
         DO 20, J = 1, N
            DO 10, I = 1, M
               B( I, J ) = ZERO
   10       CONTINUE
   20    CONTINUE
         RETURN
      END IF
#
#     Start the operations.
#
      IF( LSIDE )THEN
         IF( LSAME( TRANSA, 'N' ) )THEN
#
#           Form  B := alpha*A*B.
#
            IF( UPPER )THEN
               DO 50, J = 1, N
                  DO 40, K = 1, M
                     IF( B( K, J )!=ZERO )THEN
                        TEMP = ALPHA*B( K, J )
                        DO 30, I = 1, K - 1
                           B( I, J ) = B( I, J ) + TEMP*A( I, K )
   30                   CONTINUE
                        IF( NOUNIT )
     $                     TEMP = TEMP*A( K, K )
                        B( K, J ) = TEMP
                     END IF
   40             CONTINUE
   50          CONTINUE
            ELSE
               DO 80, J = 1, N
                  DO 70 K = M, 1, -1
                     IF( B( K, J )!=ZERO )THEN
                        TEMP      = ALPHA*B( K, J )
                        B( K, J ) = TEMP
                        IF( NOUNIT )
     $                     B( K, J ) = B( K, J )*A( K, K )
                        DO 60, I = K + 1, M
                           B( I, J ) = B( I, J ) + TEMP*A( I, K )
   60                   CONTINUE
                     END IF
   70             CONTINUE
   80          CONTINUE
            END IF
         ELSE
#
#           Form  B := alpha*A'*B   or   B := alpha*conjg( A' )*B.
#
            IF( UPPER )THEN
               DO 120, J = 1, N
                  DO 110, I = M, 1, -1
                     TEMP = B( I, J )
                     IF( NOCONJ )THEN
                        IF( NOUNIT )
     $                     TEMP = TEMP*A( I, I )
                        DO 90, K = 1, I - 1
                           TEMP = TEMP + A( K, I )*B( K, J )
   90                   CONTINUE
                     ELSE
                        IF( NOUNIT )
     $                     TEMP = TEMP*CONJG( A( I, I ) )
                        DO 100, K = 1, I - 1
                           TEMP = TEMP + CONJG( A( K, I ) )*B( K, J )
  100                   CONTINUE
                     END IF
                     B( I, J ) = ALPHA*TEMP
  110             CONTINUE
  120          CONTINUE
            ELSE
               DO 160, J = 1, N
                  DO 150, I = 1, M
                     TEMP = B( I, J )
                     IF( NOCONJ )THEN
                        IF( NOUNIT )
     $                     TEMP = TEMP*A( I, I )
                        DO 130, K = I + 1, M
                           TEMP = TEMP + A( K, I )*B( K, J )
  130                   CONTINUE
                     ELSE
                        IF( NOUNIT )
     $                     TEMP = TEMP*CONJG( A( I, I ) )
                        DO 140, K = I + 1, M
                           TEMP = TEMP + CONJG( A( K, I ) )*B( K, J )
  140                   CONTINUE
                     END IF
                     B( I, J ) = ALPHA*TEMP
  150             CONTINUE
  160          CONTINUE
            END IF
         END IF
      ELSE
         IF( LSAME( TRANSA, 'N' ) )THEN
#
#           Form  B := alpha*B*A.
#
            IF( UPPER )THEN
               DO 200, J = N, 1, -1
                  TEMP = ALPHA
                  IF( NOUNIT )
     $               TEMP = TEMP*A( J, J )
                  DO 170, I = 1, M
                     B( I, J ) = TEMP*B( I, J )
  170             CONTINUE
                  DO 190, K = 1, J - 1
                     IF( A( K, J )!=ZERO )THEN
                        TEMP = ALPHA*A( K, J )
                        DO 180, I = 1, M
                           B( I, J ) = B( I, J ) + TEMP*B( I, K )
  180                   CONTINUE
                     END IF
  190             CONTINUE
  200          CONTINUE
            ELSE
               DO 240, J = 1, N
                  TEMP = ALPHA
                  IF( NOUNIT )
     $               TEMP = TEMP*A( J, J )
                  DO 210, I = 1, M
                     B( I, J ) = TEMP*B( I, J )
  210             CONTINUE
                  DO 230, K = J + 1, N
                     IF( A( K, J )!=ZERO )THEN
                        TEMP = ALPHA*A( K, J )
                        DO 220, I = 1, M
                           B( I, J ) = B( I, J ) + TEMP*B( I, K )
  220                   CONTINUE
                     END IF
  230             CONTINUE
  240          CONTINUE
            END IF
         ELSE
#
#           Form  B := alpha*B*A'   or   B := alpha*B*conjg( A' ).
#
            IF( UPPER )THEN
               DO 280, K = 1, N
                  DO 260, J = 1, K - 1
                     IF( A( J, K )!=ZERO )THEN
                        IF( NOCONJ )THEN
                           TEMP = ALPHA*A( J, K )
                        ELSE
                           TEMP = ALPHA*CONJG( A( J, K ) )
                        END IF
                        DO 250, I = 1, M
                           B( I, J ) = B( I, J ) + TEMP*B( I, K )
  250                   CONTINUE
                     END IF
  260             CONTINUE
                  TEMP = ALPHA
                  IF( NOUNIT )THEN
                     IF( NOCONJ )THEN
                        TEMP = TEMP*A( K, K )
                     ELSE
                        TEMP = TEMP*CONJG( A( K, K ) )
                     END IF
                  END IF
                  IF( TEMP!=ONE )THEN
                     DO 270, I = 1, M
                        B( I, K ) = TEMP*B( I, K )
  270                CONTINUE
                  END IF
  280          CONTINUE
            ELSE
               DO 320, K = N, 1, -1
                  DO 300, J = K + 1, N
                     IF( A( J, K )!=ZERO )THEN
                        IF( NOCONJ )THEN
                           TEMP = ALPHA*A( J, K )
                        ELSE
                           TEMP = ALPHA*CONJG( A( J, K ) )
                        END IF
                        DO 290, I = 1, M
                           B( I, J ) = B( I, J ) + TEMP*B( I, K )
  290                   CONTINUE
                     END IF
  300             CONTINUE
                  TEMP = ALPHA
                  IF( NOUNIT )THEN
                     IF( NOCONJ )THEN
                        TEMP = TEMP*A( K, K )
                     ELSE
                        TEMP = TEMP*CONJG( A( K, K ) )
                     END IF
                  END IF
                  IF( TEMP!=ONE )THEN
                     DO 310, I = 1, M
                        B( I, K ) = TEMP*B( I, K )
  310                CONTINUE
                  END IF
  320          CONTINUE
            END IF
         END IF
      END IF
#
      RETURN
#
#     End of CTRMM .
#
      END