[BACK]Return to drotm.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / blas

Diff for /rpl/lapack/blas/drotm.f between versions 1.6 and 1.17

version 1.6, 2010/08/13 21:03:40 version 1.17, 2023/08/07 08:38:43
Line 1 Line 1
       SUBROUTINE DROTM(N,DX,INCX,DY,INCY,DPARAM)  *> \brief \b DROTM
 *     .. Scalar Arguments ..  
       INTEGER INCX,INCY,N  
 *     ..  
 *     .. Array Arguments ..  
       DOUBLE PRECISION DPARAM(5),DX(*),DY(*)  
 *     ..  
 *  
 *  Purpose  
 *  =======  
 *  
 *     APPLY THE MODIFIED GIVENS TRANSFORMATION, H, TO THE 2 BY N MATRIX  
 *  
 *     (DX**T) , WHERE **T INDICATES TRANSPOSE. THE ELEMENTS OF DX ARE IN  
 *     (DY**T)  
 *  
 *     DX(LX+I*INCX), I = 0 TO N-1, WHERE LX = 1 IF INCX .GE. 0, ELSE  
 *     LX = (-INCX)*N, AND SIMILARLY FOR SY USING LY AND INCY.  
 *     WITH DPARAM(1)=DFLAG, H HAS ONE OF THE FOLLOWING FORMS..  
 *  
 *     DFLAG=-1.D0     DFLAG=0.D0        DFLAG=1.D0     DFLAG=-2.D0  
 *  
 *       (DH11  DH12)    (1.D0  DH12)    (DH11  1.D0)    (1.D0  0.D0)  
 *     H=(          )    (          )    (          )    (          )  
 *       (DH21  DH22),   (DH21  1.D0),   (-1.D0 DH22),   (0.D0  1.D0).  
 *     SEE DROTMG FOR A DESCRIPTION OF DATA STORAGE IN DPARAM.  
 *  *
 *  Arguments  *  =========== DOCUMENTATION ===========
 *  =========  
 *  *
 *  N      (input) INTEGER  * Online html documentation available at
 *         number of elements in input vector(s)  *            http://www.netlib.org/lapack/explore-html/
 *  *
 *  DX     (input/output) DOUBLE PRECISION array, dimension N  *  Definition:
 *         double precision vector with N elements  *  ===========
   *
   *       SUBROUTINE DROTM(N,DX,INCX,DY,INCY,DPARAM)
   *
   *       .. Scalar Arguments ..
   *       INTEGER INCX,INCY,N
   *       ..
   *       .. Array Arguments ..
   *       DOUBLE PRECISION DPARAM(5),DX(*),DY(*)
   *       ..
   *
   *
   *> \par Purpose:
   *  =============
   *>
   *> \verbatim
   *>
   *>    APPLY THE MODIFIED GIVENS TRANSFORMATION, H, TO THE 2 BY N MATRIX
   *>
   *>    (DX**T) , WHERE **T INDICATES TRANSPOSE. THE ELEMENTS OF DX ARE IN
   *>    (DY**T)
   *>
   *>    DX(LX+I*INCX), I = 0 TO N-1, WHERE LX = 1 IF INCX .GE. 0, ELSE
   *>    LX = (-INCX)*N, AND SIMILARLY FOR SY USING LY AND INCY.
   *>    WITH DPARAM(1)=DFLAG, H HAS ONE OF THE FOLLOWING FORMS..
   *>
   *>    DFLAG=-1.D0     DFLAG=0.D0        DFLAG=1.D0     DFLAG=-2.D0
   *>
   *>      (DH11  DH12)    (1.D0  DH12)    (DH11  1.D0)    (1.D0  0.D0)
   *>    H=(          )    (          )    (          )    (          )
   *>      (DH21  DH22),   (DH21  1.D0),   (-1.D0 DH22),   (0.D0  1.D0).
   *>    SEE DROTMG FOR A DESCRIPTION OF DATA STORAGE IN DPARAM.
   *> \endverbatim
   *
   *  Arguments:
   *  ==========
   *
   *> \param[in] N
   *> \verbatim
   *>          N is INTEGER
   *>         number of elements in input vector(s)
   *> \endverbatim
   *>
   *> \param[in,out] DX
   *> \verbatim
   *>          DX is DOUBLE PRECISION array, dimension ( 1 + ( N - 1 )*abs( INCX ) )
   *> \endverbatim
   *>
   *> \param[in] INCX
   *> \verbatim
   *>          INCX is INTEGER
   *>         storage spacing between elements of DX
   *> \endverbatim
   *>
   *> \param[in,out] DY
   *> \verbatim
   *>          DY is DOUBLE PRECISION array, dimension ( 1 + ( N - 1 )*abs( INCY ) )
   *> \endverbatim
   *>
   *> \param[in] INCY
   *> \verbatim
   *>          INCY is INTEGER
   *>         storage spacing between elements of DY
   *> \endverbatim
   *>
   *> \param[in] DPARAM
   *> \verbatim
   *>          DPARAM is DOUBLE PRECISION array, dimension (5)
   *>     DPARAM(1)=DFLAG
   *>     DPARAM(2)=DH11
   *>     DPARAM(3)=DH21
   *>     DPARAM(4)=DH12
   *>     DPARAM(5)=DH22
   *> \endverbatim
   *
   *  Authors:
   *  ========
   *
   *> \author Univ. of Tennessee
   *> \author Univ. of California Berkeley
   *> \author Univ. of Colorado Denver
   *> \author NAG Ltd.
 *  *
 *  INCX   (input) INTEGER  *> \ingroup double_blas_level1
 *         storage spacing between elements of DX  
 *  *
 *  DY     (input/output) DOUBLE PRECISION array, dimension N  *  =====================================================================
 *         double precision vector with N elements        SUBROUTINE DROTM(N,DX,INCX,DY,INCY,DPARAM)
 *  *
 *  INCY   (input) INTEGER  *  -- Reference BLAS level1 routine --
 *         storage spacing between elements of DY  *  -- Reference BLAS is a software package provided by Univ. of Tennessee,    --
   *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
 *  *
 *  DPARAM (input/output)  DOUBLE PRECISION array, dimension 5   *     .. Scalar Arguments ..
 *     DPARAM(1)=DFLAG        INTEGER INCX,INCY,N
 *     DPARAM(2)=DH11  *     ..
 *     DPARAM(3)=DH21  *     .. Array Arguments ..
 *     DPARAM(4)=DH12        DOUBLE PRECISION DPARAM(5),DX(*),DY(*)
 *     DPARAM(5)=DH22  *     ..
 *  *
 *  =====================================================================  *  =====================================================================
 *  *
Line 61 Line 116
 *     ..  *     ..
 *  *
       DFLAG = DPARAM(1)        DFLAG = DPARAM(1)
       IF (N.LE.0 .OR. (DFLAG+TWO.EQ.ZERO)) GO TO 140        IF (N.LE.0 .OR. (DFLAG+TWO.EQ.ZERO)) RETURN
       IF (.NOT. (INCX.EQ.INCY.AND.INCX.GT.0)) GO TO 70        IF (INCX.EQ.INCY.AND.INCX.GT.0) THEN
 *  *
       NSTEPS = N*INCX           NSTEPS = N*INCX
       IF (DFLAG.LT.ZERO) THEN           IF (DFLAG.LT.ZERO) THEN
          GO TO 50              DH11 = DPARAM(2)
       ELSE IF (DFLAG.EQ.ZERO) THEN              DH12 = DPARAM(4)
          GO TO 10               DH21 = DPARAM(3)
       ELSE              DH22 = DPARAM(5)
          GO TO 30              DO I = 1,NSTEPS,INCX
       END IF                 W = DX(I)
    10 CONTINUE                 Z = DY(I)
       DH12 = DPARAM(4)                 DX(I) = W*DH11 + Z*DH12
       DH21 = DPARAM(3)                 DY(I) = W*DH21 + Z*DH22
       DO 20 I = 1,NSTEPS,INCX              END DO
           W = DX(I)           ELSE IF (DFLAG.EQ.ZERO) THEN
           Z = DY(I)              DH12 = DPARAM(4)
           DX(I) = W + Z*DH12              DH21 = DPARAM(3)
           DY(I) = W*DH21 + Z              DO I = 1,NSTEPS,INCX
    20 CONTINUE                 W = DX(I)
       GO TO 140                 Z = DY(I)
    30 CONTINUE                 DX(I) = W + Z*DH12
       DH11 = DPARAM(2)                 DY(I) = W*DH21 + Z
       DH22 = DPARAM(5)              END DO
       DO 40 I = 1,NSTEPS,INCX           ELSE
           W = DX(I)              DH11 = DPARAM(2)
           Z = DY(I)              DH22 = DPARAM(5)
           DX(I) = W*DH11 + Z              DO I = 1,NSTEPS,INCX
           DY(I) = -W + DH22*Z                 W = DX(I)
    40 CONTINUE                 Z = DY(I)
       GO TO 140                 DX(I) = W*DH11 + Z
    50 CONTINUE                 DY(I) = -W + DH22*Z
       DH11 = DPARAM(2)              END DO
       DH12 = DPARAM(4)           END IF
       DH21 = DPARAM(3)  
       DH22 = DPARAM(5)  
       DO 60 I = 1,NSTEPS,INCX  
           W = DX(I)  
           Z = DY(I)  
           DX(I) = W*DH11 + Z*DH12  
           DY(I) = W*DH21 + Z*DH22  
    60 CONTINUE  
       GO TO 140  
    70 CONTINUE  
       KX = 1  
       KY = 1  
       IF (INCX.LT.0) KX = 1 + (1-N)*INCX  
       IF (INCY.LT.0) KY = 1 + (1-N)*INCY  
 *  
       IF (DFLAG.LT.ZERO) THEN  
          GO TO 120  
       ELSE IF (DFLAG.EQ.ZERO) THEN  
          GO TO 80   
       ELSE        ELSE
          GO TO 100           KX = 1
            KY = 1
            IF (INCX.LT.0) KX = 1 + (1-N)*INCX
            IF (INCY.LT.0) KY = 1 + (1-N)*INCY
   *
            IF (DFLAG.LT.ZERO) THEN
               DH11 = DPARAM(2)
               DH12 = DPARAM(4)
               DH21 = DPARAM(3)
               DH22 = DPARAM(5)
               DO I = 1,N
                  W = DX(KX)
                  Z = DY(KY)
                  DX(KX) = W*DH11 + Z*DH12
                  DY(KY) = W*DH21 + Z*DH22
                  KX = KX + INCX
                  KY = KY + INCY
               END DO
            ELSE IF (DFLAG.EQ.ZERO) THEN
               DH12 = DPARAM(4)
               DH21 = DPARAM(3)
               DO I = 1,N
                  W = DX(KX)
                  Z = DY(KY)
                  DX(KX) = W + Z*DH12
                  DY(KY) = W*DH21 + Z
                  KX = KX + INCX
                  KY = KY + INCY
               END DO
            ELSE
                DH11 = DPARAM(2)
                DH22 = DPARAM(5)
                DO I = 1,N
                   W = DX(KX)
                   Z = DY(KY)
                   DX(KX) = W*DH11 + Z
                   DY(KY) = -W + DH22*Z
                   KX = KX + INCX
                   KY = KY + INCY
               END DO
            END IF
       END IF        END IF
    80 CONTINUE  
       DH12 = DPARAM(4)  
       DH21 = DPARAM(3)  
       DO 90 I = 1,N  
           W = DX(KX)  
           Z = DY(KY)  
           DX(KX) = W + Z*DH12  
           DY(KY) = W*DH21 + Z  
           KX = KX + INCX  
           KY = KY + INCY  
    90 CONTINUE  
       GO TO 140  
   100 CONTINUE  
       DH11 = DPARAM(2)  
       DH22 = DPARAM(5)  
       DO 110 I = 1,N  
           W = DX(KX)  
           Z = DY(KY)  
           DX(KX) = W*DH11 + Z  
           DY(KY) = -W + DH22*Z  
           KX = KX + INCX  
           KY = KY + INCY  
   110 CONTINUE  
       GO TO 140  
   120 CONTINUE  
       DH11 = DPARAM(2)  
       DH12 = DPARAM(4)  
       DH21 = DPARAM(3)  
       DH22 = DPARAM(5)  
       DO 130 I = 1,N  
           W = DX(KX)  
           Z = DY(KY)  
           DX(KX) = W*DH11 + Z*DH12  
           DY(KY) = W*DH21 + Z*DH22  
           KX = KX + INCX  
           KY = KY + INCY  
   130 CONTINUE  
   140 CONTINUE  
       RETURN        RETURN
   *
   *     End of DROTM
   *
       END        END
Removed from v.1.6  
changed lines
  Added in v.1.17

CVSweb interface <joel.bertrand@systella.fr>