Diff for /rpl/lapack/blas/drotmg.f between versions 1.1 and 1.6

version 1.1, 2010/01/26 15:22:45 version 1.6, 2010/08/13 21:03:40
Line 33 Line 33
 *  *
 *  DD1    (input/output) DOUBLE PRECISION  *  DD1    (input/output) DOUBLE PRECISION
 *  *
 *  DD2    (input/output) DOUBLE PRECISION   *  DD2    (input/output) DOUBLE PRECISION
 *  *
 *  DX1    (input/output) DOUBLE PRECISION   *  DX1    (input/output) DOUBLE PRECISION
 *  *
 *  DY1    (input) DOUBLE PRECISION  *  DY1    (input) DOUBLE PRECISION
 *  *
Line 71 Line 71
       IF (.NOT.DP2.EQ.ZERO) GO TO 20        IF (.NOT.DP2.EQ.ZERO) GO TO 20
       DFLAG = -TWO        DFLAG = -TWO
       GO TO 260        GO TO 260
 *     REGULAR-CASE..  
    20 CONTINUE     20 CONTINUE
   *     REGULAR-CASE..
       DP1 = DD1*DX1        DP1 = DD1*DX1
       DQ2 = DP2*DY1        DQ2 = DP2*DY1
       DQ1 = DP1*DX1        DQ1 = DP1*DX1
Line 108 Line 108
       DX1 = DY1*DU        DX1 = DY1*DU
 *         GO SCALE-CHECK  *         GO SCALE-CHECK
       GO TO 100        GO TO 100
 *     PROCEDURE..ZERO-H-D-AND-DX1..  
    60 CONTINUE     60 CONTINUE
   *     PROCEDURE..ZERO-H-D-AND-DX1..
       DFLAG = -ONE        DFLAG = -ONE
       DH11 = ZERO        DH11 = ZERO
       DH12 = ZERO        DH12 = ZERO
Line 121 Line 121
       DX1 = ZERO        DX1 = ZERO
 *         RETURN..  *         RETURN..
       GO TO 220        GO TO 220
 *     PROCEDURE..FIX-H..  
    70 CONTINUE     70 CONTINUE
   *     PROCEDURE..FIX-H..
       IF (.NOT.DFLAG.GE.ZERO) GO TO 90        IF (.NOT.DFLAG.GE.ZERO) GO TO 90
 *  *
       IF (.NOT.DFLAG.EQ.ZERO) GO TO 80        IF (.NOT.DFLAG.EQ.ZERO) GO TO 80
Line 135 Line 135
       DH12 = ONE        DH12 = ONE
       DFLAG = -ONE        DFLAG = -ONE
    90 CONTINUE     90 CONTINUE
       GO TO IGO(120,150,180,210)        GO TO (150,180,210) IGO
 *     PROCEDURE..SCALE-CHECK        GO TO 120
   100 CONTINUE    100 CONTINUE
   *     PROCEDURE..SCALE-CHECK
   110 CONTINUE    110 CONTINUE
       IF (.NOT.DD1.LE.RGAMSQ) GO TO 130        IF (.NOT.DD1.LE.RGAMSQ) GO TO 130
       IF (DD1.EQ.ZERO) GO TO 160        IF (DD1.EQ.ZERO) GO TO 160
       ASSIGN 120 TO IGO        IGO = 0
 *              FIX-H..  *              FIX-H..
       GO TO 70        GO TO 70
   120 CONTINUE    120 CONTINUE
Line 153 Line 154
   130 CONTINUE    130 CONTINUE
   140 CONTINUE    140 CONTINUE
       IF (.NOT.DD1.GE.GAMSQ) GO TO 160        IF (.NOT.DD1.GE.GAMSQ) GO TO 160
       ASSIGN 150 TO IGO        IGO = 1
 *              FIX-H..  *              FIX-H..
       GO TO 70        GO TO 70
   150 CONTINUE    150 CONTINUE
Line 166 Line 167
   170 CONTINUE    170 CONTINUE
       IF (.NOT.DABS(DD2).LE.RGAMSQ) GO TO 190        IF (.NOT.DABS(DD2).LE.RGAMSQ) GO TO 190
       IF (DD2.EQ.ZERO) GO TO 220        IF (DD2.EQ.ZERO) GO TO 220
       ASSIGN 180 TO IGO        IGO = 2
 *              FIX-H..  *              FIX-H..
       GO TO 70        GO TO 70
   180 CONTINUE    180 CONTINUE
Line 177 Line 178
   190 CONTINUE    190 CONTINUE
   200 CONTINUE    200 CONTINUE
       IF (.NOT.DABS(DD2).GE.GAMSQ) GO TO 220        IF (.NOT.DABS(DD2).GE.GAMSQ) GO TO 220
       ASSIGN 210 TO IGO        IGO = 3
 *              FIX-H..  *              FIX-H..
       GO TO 70        GO TO 70
   210 CONTINUE    210 CONTINUE
Line 186 Line 187
       DH22 = DH22*GAM        DH22 = DH22*GAM
       GO TO 200        GO TO 200
   220 CONTINUE    220 CONTINUE
       IF (DFLAG) 250,230,240        IF (DFLAG.LT.ZERO) THEN
            GO TO 250
         ELSE IF (DFLAG.EQ.ZERO) THEN
            GO TO 230 
         ELSE
            GO TO 240
         END IF
   230 CONTINUE    230 CONTINUE
       DPARAM(3) = DH21        DPARAM(3) = DH21
       DPARAM(4) = DH12        DPARAM(4) = DH12

Removed from v.1.1  
changed lines
  Added in v.1.6


CVSweb interface <joel.bertrand@systella.fr>