version 1.15, 2018/05/29 07:19:43
|
version 1.16, 2023/08/07 08:38:45
|
Line 149
|
Line 149
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date December 2016 |
|
* |
|
*> \ingroup complex16_blas_level3 |
*> \ingroup complex16_blas_level3 |
* |
* |
*> \par Further Details: |
*> \par Further Details: |
Line 173
|
Line 171
|
* ===================================================================== |
* ===================================================================== |
SUBROUTINE ZHERK(UPLO,TRANS,N,K,ALPHA,A,LDA,BETA,C,LDC) |
SUBROUTINE ZHERK(UPLO,TRANS,N,K,ALPHA,A,LDA,BETA,C,LDC) |
* |
* |
* -- Reference BLAS level3 routine (version 3.7.0) -- |
* -- Reference BLAS level3 routine -- |
* -- Reference BLAS is a software package provided by Univ. of Tennessee, -- |
* -- Reference BLAS is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* December 2016 |
|
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
DOUBLE PRECISION ALPHA,BETA |
DOUBLE PRECISION ALPHA,BETA |
Line 355
|
Line 352
|
200 CONTINUE |
200 CONTINUE |
RTEMP = ZERO |
RTEMP = ZERO |
DO 210 L = 1,K |
DO 210 L = 1,K |
RTEMP = RTEMP + DCONJG(A(L,J))*A(L,J) |
RTEMP = RTEMP + DBLE(DCONJG(A(L,J))*A(L,J)) |
210 CONTINUE |
210 CONTINUE |
IF (BETA.EQ.ZERO) THEN |
IF (BETA.EQ.ZERO) THEN |
C(J,J) = ALPHA*RTEMP |
C(J,J) = ALPHA*RTEMP |
Line 367
|
Line 364
|
DO 260 J = 1,N |
DO 260 J = 1,N |
RTEMP = ZERO |
RTEMP = ZERO |
DO 230 L = 1,K |
DO 230 L = 1,K |
RTEMP = RTEMP + DCONJG(A(L,J))*A(L,J) |
RTEMP = RTEMP + DBLE(DCONJG(A(L,J))*A(L,J)) |
230 CONTINUE |
230 CONTINUE |
IF (BETA.EQ.ZERO) THEN |
IF (BETA.EQ.ZERO) THEN |
C(J,J) = ALPHA*RTEMP |
C(J,J) = ALPHA*RTEMP |
Line 391
|
Line 388
|
* |
* |
RETURN |
RETURN |
* |
* |
* End of ZHERK . |
* End of ZHERK |
* |
* |
END |
END |