SUBROUTINE ZTRSMF ( SIDE, UPLO, TRANSA, DIAG, M, N, ALPHA, A, LDA,

     $                   B, LDB )

*     .. Scalar Arguments ..

      IMPLICIT NONE

      CHARACTER*1        SIDE, UPLO, TRANSA, DIAG

      INTEGER            M, N, LDA, LDB

      COMPLEX*16         ALPHA

*     .. Array Arguments ..

      COMPLEX*16         A( LDA, * ), B( LDB, * )

*     ..

*

*  Purpose

*  =======

*

*  ZTRSM  solves one of the matrix equations

*

*     op( A )*X = alpha*B,   or   X*op( A ) = alpha*B,

*

*  where alpha is a scalar, X and B are m by n matrices, 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' ).

*

*  The matrix X is overwritten on B.

*

*  Parameters

*  ==========

*

*  SIDE   - CHARACTER*1.

*           On entry, SIDE specifies whether op( A ) appears on the left

*           or right of X as follows:

*

*              SIDE = 'L' or 'l'   op( A )*X = alpha*B.

*

*              SIDE = 'R' or 'r'   X*op( A ) = alpha*B.

*

*           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*16      .

*           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*16       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*16       array of DIMENSION ( LDB, n ).

*           Before entry,  the leading  m by n part of the array  B must

*           contain  the  right-hand  side  matrix  B,  and  on exit  is

*           overwritten by the solution matrix  X.

*

*  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          DCONJG, MAX

*     .. Local Scalars ..

      LOGICAL            LSIDE, NOCONJ, NOUNIT, UPPER

      INTEGER            I, INFO, J, K, NROWA

      COMPLEX*16         TEMP

*     .. Parameters ..

      COMPLEX*16         ONE

      PARAMETER        ( ONE  = ( 1.0D+0, 0.0D+0 ) )

      COMPLEX*16         ZERO

      PARAMETER        ( ZERO = ( 0.0D+0, 0.0D+0 ) )

*     ..

*     .. Executable Statements ..

*

*     Test the input parameters.

*

      LSIDE  = LSAME( SIDE  , 'L' )

      IF( LSIDE )THEN

         NROWA = M

      ELSE

         NROWA = N

      END IF

      NOCONJ = (LSAME( TRANSA, 'N' ) .OR. LSAME( TRANSA, 'T' ))

      NOUNIT = LSAME( DIAG  , 'N' )

      UPPER  = LSAME( UPLO  , 'U' )

*

      INFO   = 0

      IF(      ( .NOT.LSIDE                ).AND.

     $         ( .NOT.LSAME( SIDE  , 'R' ) )      )THEN

         INFO = 1

      ELSE IF( ( .NOT.UPPER                ).AND.

     $         ( .NOT.LSAME( UPLO  , 'L' ) )      )THEN

         INFO = 2

      ELSE IF( ( .NOT.LSAME( TRANSA, 'N' ) ).AND.

     $         ( .NOT.LSAME( TRANSA, 'T' ) ).AND.

     $         ( .NOT.LSAME( TRANSA, 'R' ) ).AND.

     $         ( .NOT.LSAME( TRANSA, 'C' ) )      )THEN

         INFO = 3

      ELSE IF( ( .NOT.LSAME( DIAG  , 'U' ) ).AND.

     $         ( .NOT.LSAME( DIAG  , 'N' ) )      )THEN

         INFO = 4

      ELSE IF( M  .LT.0               )THEN

         INFO = 5

      ELSE IF( N  .LT.0               )THEN

         INFO = 6

      ELSE IF( LDA.LT.MAX( 1, NROWA ) )THEN

         INFO = 9

      ELSE IF( LDB.LT.MAX( 1, M     ) )THEN

         INFO = 11

      END IF

      IF( INFO.NE.0 )THEN

         CALL XERBLA( 'ZTRSM ', INFO )

         RETURN

      END IF

*

*     Quick return if possible.

*

      IF( N.EQ.0 )

     $   RETURN

*

*     And when  alpha.eq.zero.

*

      IF( ALPHA.EQ.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' ) .OR. LSAME( TRANSA, 'R' ) )THEN

*

*           Form  B := alpha*inv( A )*B.

*

            IF( UPPER )THEN

               DO 60, J = 1, N

                  IF( ALPHA.NE.ONE )THEN

                     DO 30, I = 1, M

                        B( I, J ) = ALPHA*B( I, J )

   30                CONTINUE

                  END IF

                  DO 50, K = M, 1, -1

                     IF( B( K, J ).NE.ZERO )THEN

                        IF( NOUNIT ) THEN

                           IF (NOCONJ) THEN

                              B( K, J ) = B( K, J )/A( K, K )

                           ELSE

                              B( K, J ) = B( K, J )/DCONJG(A( K, K ))

                           ENDIF

                        ENDIF

                        IF (NOCONJ) THEN

                           DO 40, I = 1, K - 1

                           B( I, J ) = B( I, J ) - B( K, J )*A( I, K )

   40                      CONTINUE

                        ELSE

                           DO 45, I = 1, K - 1

                    B( I, J ) = B( I, J ) - B( K, J )*DCONJG(A( I, K ))

   45                      CONTINUE

                        ENDIF

                     END IF

   50             CONTINUE

   60          CONTINUE

            ELSE

               DO 100, J = 1, N

                  IF( ALPHA.NE.ONE )THEN

                     DO 70, I = 1, M

                        B( I, J ) = ALPHA*B( I, J )

   70                CONTINUE

                  END IF

                  DO 90 K = 1, M

                     IF (NOCONJ) THEN

                     IF( B( K, J ).NE.ZERO )THEN

                        IF( NOUNIT )

     $                     B( K, J ) = B( K, J )/A( K, K )

                        DO 80, I = K + 1, M

                           B( I, J ) = B( I, J ) - B( K, J )*A( I, K )

   80                   CONTINUE

                     END IF

                     ELSE

                     IF( B( K, J ).NE.ZERO )THEN

                        IF( NOUNIT )

     $                     B( K, J ) = B( K, J )/DCONJG(A( K, K ))

                        DO 85, I = K + 1, M

                     B( I, J ) = B( I, J ) - B( K, J )*DCONJG(A( I, K ))

 85                   CONTINUE

                     END IF

                     ENDIF

   90             CONTINUE

  100          CONTINUE

            END IF

         ELSE

*

*           Form  B := alpha*inv( A' )*B

*           or    B := alpha*inv( conjg( A' ) )*B.

*

            IF( UPPER )THEN

               DO 140, J = 1, N

                  DO 130, I = 1, M

                     TEMP = ALPHA*B( I, J )

                     IF( NOCONJ )THEN

                        DO 110, K = 1, I - 1

                           TEMP = TEMP - A( K, I )*B( K, J )

  110                   CONTINUE

                        IF( NOUNIT )

     $                     TEMP = TEMP/A( I, I )

                     ELSE

                        DO 120, K = 1, I - 1

                           TEMP = TEMP - DCONJG( A( K, I ) )*B( K, J )

  120                   CONTINUE

                        IF( NOUNIT )

     $                     TEMP = TEMP/DCONJG( A( I, I ) )

                     END IF

                     B( I, J ) = TEMP

  130             CONTINUE

  140          CONTINUE

            ELSE

               DO 180, J = 1, N

                  DO 170, I = M, 1, -1

                     TEMP = ALPHA*B( I, J )

                     IF( NOCONJ )THEN

                        DO 150, K = I + 1, M

                           TEMP = TEMP - A( K, I )*B( K, J )

  150                   CONTINUE

                        IF( NOUNIT )

     $                     TEMP = TEMP/A( I, I )

                     ELSE

                        DO 160, K = I + 1, M

                           TEMP = TEMP - DCONJG( A( K, I ) )*B( K, J )

  160                   CONTINUE

                        IF( NOUNIT )

     $                     TEMP = TEMP/DCONJG( A( I, I ) )

                     END IF

                     B( I, J ) = TEMP

  170             CONTINUE

  180          CONTINUE

            END IF

         END IF

      ELSE

         IF( LSAME( TRANSA, 'N' ) .OR. LSAME( TRANSA, 'R' ) )THEN

*

*           Form  B := alpha*B*inv( A ).

*

            IF( UPPER )THEN

               DO 230, J = 1, N

                  IF( ALPHA.NE.ONE )THEN

                     DO 190, I = 1, M

                        B( I, J ) = ALPHA*B( I, J )

  190                CONTINUE

                  END IF

                  DO 210, K = 1, J - 1

                     IF( A( K, J ).NE.ZERO )THEN

                        IF (NOCONJ) THEN

                        DO 200, I = 1, M

                           B( I, J ) = B( I, J ) - A( K, J )*B( I, K )

  200                   CONTINUE

                        ELSE

                        DO 205, I = 1, M

                     B( I, J ) = B( I, J ) - DCONJG(A( K, J ))*B( I, K )

  205                   CONTINUE

                        ENDIF

                     END IF

  210             CONTINUE

                  IF( NOUNIT )THEN

                     IF (NOCONJ) THEN

                     TEMP = ONE/A( J, J )

                     ELSE

                     TEMP = ONE/DCONJG(A( J, J ))

                     ENDIF

                     DO 220, I = 1, M

                        B( I, J ) = TEMP*B( I, J )

  220                CONTINUE

                  END IF

  230          CONTINUE

            ELSE

               DO 280, J = N, 1, -1

                  IF( ALPHA.NE.ONE )THEN

                     DO 240, I = 1, M

                        B( I, J ) = ALPHA*B( I, J )

  240                CONTINUE

                  END IF

                  DO 260, K = J + 1, N

                     IF( A( K, J ).NE.ZERO )THEN

                        IF (NOCONJ) THEN

                        DO 250, I = 1, M

                           B( I, J ) = B( I, J ) - A( K, J )*B( I, K )

  250                   CONTINUE

                        ELSE

                        DO 255, I = 1, M

                    B( I, J ) = B( I, J ) - DCONJG(A( K, J ))*B( I, K )

  255                   CONTINUE

                        ENDIF

                     END IF

  260             CONTINUE

                  IF( NOUNIT )THEN

                     IF (NOCONJ) THEN

                     TEMP = ONE/A( J, J )

                     ELSE

                     TEMP = ONE/DCONJG(A( J, J ))

                     ENDIF

                     DO 270, I = 1, M

                       B( I, J ) = TEMP*B( I, J )

  270                CONTINUE

                  END IF

  280          CONTINUE

            END IF

         ELSE

*

*           Form  B := alpha*B*inv( A' )

*           or    B := alpha*B*inv( conjg( A' ) ).

*

            IF( UPPER )THEN

               DO 330, K = N, 1, -1

                  IF( NOUNIT )THEN

                     IF( NOCONJ )THEN

                        TEMP = ONE/A( K, K )

                     ELSE

                        TEMP = ONE/DCONJG( A( K, K ) )

                     END IF

                     DO 290, I = 1, M

                        B( I, K ) = TEMP*B( I, K )

  290                CONTINUE

                  END IF

                  DO 310, J = 1, K - 1

                     IF( A( J, K ).NE.ZERO )THEN

                        IF( NOCONJ )THEN

                           TEMP = A( J, K )

                        ELSE

                           TEMP = DCONJG( A( J, K ) )

                        END IF

                        DO 300, I = 1, M

                           B( I, J ) = B( I, J ) - TEMP*B( I, K )

  300                   CONTINUE

                     END IF

  310             CONTINUE

                  IF( ALPHA.NE.ONE )THEN

                     DO 320, I = 1, M

                        B( I, K ) = ALPHA*B( I, K )

  320                CONTINUE

                  END IF

  330          CONTINUE

            ELSE

               DO 380, K = 1, N

                  IF( NOUNIT )THEN

                     IF( NOCONJ )THEN

                        TEMP = ONE/A( K, K )

                     ELSE

                        TEMP = ONE/DCONJG( A( K, K ) )

                     END IF

                     DO 340, I = 1, M

                        B( I, K ) = TEMP*B( I, K )

  340                CONTINUE

                  END IF

                  DO 360, J = K + 1, N

                     IF( A( J, K ).NE.ZERO )THEN

                        IF( NOCONJ )THEN

                           TEMP = A( J, K )

                        ELSE

                           TEMP = DCONJG( A( J, K ) )

                        END IF

                        DO 350, I = 1, M

                           B( I, J ) = B( I, J ) - TEMP*B( I, K )

  350                   CONTINUE

                     END IF

  360             CONTINUE

                  IF( ALPHA.NE.ONE )THEN

                     DO 370, I = 1, M

                        B( I, K ) = ALPHA*B( I, K )

  370                CONTINUE

                  END IF

  380          CONTINUE

            END IF

         END IF

      END IF

*

      RETURN

*

*     End of ZTRSM .

*

      END