Diff for /rpl/lapack/blas/drotmg.f between versions 1.11 and 1.15

version 1.11, 2012/12/14 14:22:02 version 1.15, 2018/05/29 06:55:14
Line 2 Line 2
 *  *
 *  =========== DOCUMENTATION ===========  *  =========== DOCUMENTATION ===========
 *  *
 * Online html documentation available at   * Online html documentation available at
 *            http://www.netlib.org/lapack/explore-html/   *            http://www.netlib.org/lapack/explore-html/
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
 *  *
 *       SUBROUTINE DROTMG(DD1,DD2,DX1,DY1,DPARAM)  *       SUBROUTINE DROTMG(DD1,DD2,DX1,DY1,DPARAM)
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       DOUBLE PRECISION DD1,DD2,DX1,DY1  *       DOUBLE PRECISION DD1,DD2,DX1,DY1
 *       ..  *       ..
 *       .. Array Arguments ..  *       .. Array Arguments ..
 *       DOUBLE PRECISION DPARAM(5)  *       DOUBLE PRECISION DPARAM(5)
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
Line 65 Line 65
 *>          DY1 is DOUBLE PRECISION  *>          DY1 is DOUBLE PRECISION
 *> \endverbatim  *> \endverbatim
 *>  *>
 *> \param[in,out] DPARAM  *> \param[out] DPARAM
 *> \verbatim  *> \verbatim
 *>          DPARAM is DOUBLE PRECISION array, dimension 5  *>          DPARAM is DOUBLE PRECISION array, dimension (5)
 *>     DPARAM(1)=DFLAG  *>     DPARAM(1)=DFLAG
 *>     DPARAM(2)=DH11  *>     DPARAM(2)=DH11
 *>     DPARAM(3)=DH21  *>     DPARAM(3)=DH21
Line 78 Line 78
 *  Authors:  *  Authors:
 *  ========  *  ========
 *  *
 *> \author Univ. of Tennessee   *> \author Univ. of Tennessee
 *> \author Univ. of California Berkeley   *> \author Univ. of California Berkeley
 *> \author Univ. of Colorado Denver   *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.   *> \author NAG Ltd.
 *  *
 *> \date November 2011  *> \date November 2017
 *  *
 *> \ingroup double_blas_level1  *> \ingroup double_blas_level1
 *  *
 *  =====================================================================  *  =====================================================================
       SUBROUTINE DROTMG(DD1,DD2,DX1,DY1,DPARAM)        SUBROUTINE DROTMG(DD1,DD2,DX1,DY1,DPARAM)
 *  *
 *  -- Reference BLAS level1 routine (version 3.4.0) --  *  -- Reference BLAS level1 routine (version 3.8.0) --
 *  -- 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..--
 *     November 2011  *     November 2017
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       DOUBLE PRECISION DD1,DD2,DX1,DY1        DOUBLE PRECISION DD1,DD2,DX1,DY1
Line 135 Line 135
             DFLAG = -TWO              DFLAG = -TWO
             DPARAM(1) = DFLAG              DPARAM(1) = DFLAG
             RETURN              RETURN
          END IF            END IF
 *        REGULAR-CASE..  *        REGULAR-CASE..
          DP1 = DD1*DX1           DP1 = DD1*DX1
          DQ2 = DP2*DY1           DQ2 = DP2*DY1
Line 203 Line 203
                END IF                 END IF
             ENDDO              ENDDO
          END IF           END IF
     
          IF (DD2.NE.ZERO) THEN           IF (DD2.NE.ZERO) THEN
             DO WHILE ( (DABS(DD2).LE.RGAMSQ) .OR. (DABS(DD2).GE.GAMSQ) )              DO WHILE ( (DABS(DD2).LE.RGAMSQ) .OR. (DABS(DD2).GE.GAMSQ) )
                IF (DFLAG.EQ.ZERO) THEN                 IF (DFLAG.EQ.ZERO) THEN
Line 223 Line 223
                   DD2 = DD2/GAM**2                    DD2 = DD2/GAM**2
                   DH21 = DH21*GAM                    DH21 = DH21*GAM
                   DH22 = DH22*GAM                    DH22 = DH22*GAM
                END IF                       END IF
             END DO              END DO
          END IF           END IF
        
       END IF        END IF
   
       IF (DFLAG.LT.ZERO) THEN        IF (DFLAG.LT.ZERO) THEN
Line 236 Line 236
          DPARAM(5) = DH22           DPARAM(5) = DH22
       ELSE IF (DFLAG.EQ.ZERO) THEN        ELSE IF (DFLAG.EQ.ZERO) THEN
          DPARAM(3) = DH21           DPARAM(3) = DH21
          DPARAM(4) = DH12            DPARAM(4) = DH12
       ELSE        ELSE
          DPARAM(2) = DH11           DPARAM(2) = DH11
          DPARAM(5) = DH22           DPARAM(5) = DH22
Line 245 Line 245
       DPARAM(1) = DFLAG        DPARAM(1) = DFLAG
       RETURN        RETURN
       END        END
         
        
        
        

Removed from v.1.11  
changed lines
  Added in v.1.15


CVSweb interface <joel.bertrand@systella.fr>