Diff for /rpl/lapack/blas/drotmg.f between versions 1.13 and 1.14

version 1.13, 2016/08/27 15:37:52 version 1.14, 2017/06/17 10:53:43
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 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 December 2016
 *  *
 *> \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.7.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  *     December 2016
 *  *
 *     .. 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.13  
changed lines
  Added in v.1.14


CVSweb interface <joel.bertrand@systella.fr>