Subroutine */
| 15644 | } /* zdscal_ */ |
| 15645 | |
| 15646 | /* Subroutine */ int zgemm_(char *transa, char *transb, integer *m, integer * |
| 15647 | n, integer *k, doublecomplex *alpha, doublecomplex *a, integer *lda, |
| 15648 | doublecomplex *b, integer *ldb, doublecomplex *beta, doublecomplex * |
| 15649 | c__, integer *ldc) |
| 15650 | { |
| 15651 | /* System generated locals */ |
| 15652 | integer a_dim1, a_offset, b_dim1, b_offset, c_dim1, c_offset, i__1, i__2, |
| 15653 | i__3, i__4, i__5, i__6; |
| 15654 | doublecomplex z__1, z__2, z__3, z__4; |
| 15655 | |
| 15656 | /* Local variables */ |
| 15657 | static integer i__, j, l, info; |
| 15658 | static logical nota, notb; |
| 15659 | static doublecomplex temp; |
| 15660 | static logical conja, conjb; |
| 15661 | extern logical lsame_(char *, char *); |
| 15662 | static integer nrowa, nrowb; |
| 15663 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 15664 | |
| 15665 | |
| 15666 | /* |
| 15667 | Purpose |
| 15668 | ======= |
| 15669 | |
| 15670 | ZGEMM performs one of the matrix-matrix operations |
| 15671 | |
| 15672 | C := alpha*op( A )*op( B ) + beta*C, |
| 15673 | |
| 15674 | where op( X ) is one of |
| 15675 | |
| 15676 | op( X ) = X or op( X ) = X' or op( X ) = conjg( X' ), |
| 15677 | |
| 15678 | alpha and beta are scalars, and A, B and C are matrices, with op( A ) |
| 15679 | an m by k matrix, op( B ) a k by n matrix and C an m by n matrix. |
| 15680 | |
| 15681 | Arguments |
| 15682 | ========== |
| 15683 | |
| 15684 | TRANSA - CHARACTER*1. |
| 15685 | On entry, TRANSA specifies the form of op( A ) to be used in |
| 15686 | the matrix multiplication as follows: |
| 15687 | |
| 15688 | TRANSA = 'N' or 'n', op( A ) = A. |
| 15689 | |
| 15690 | TRANSA = 'T' or 't', op( A ) = A'. |
| 15691 | |
| 15692 | TRANSA = 'C' or 'c', op( A ) = conjg( A' ). |
| 15693 | |
| 15694 | Unchanged on exit. |
| 15695 | |
| 15696 | TRANSB - CHARACTER*1. |
| 15697 | On entry, TRANSB specifies the form of op( B ) to be used in |
| 15698 | the matrix multiplication as follows: |
| 15699 | |
| 15700 | TRANSB = 'N' or 'n', op( B ) = B. |
| 15701 | |
| 15702 | TRANSB = 'T' or 't', op( B ) = B'. |
| 15703 |