URL https://opencores.org/ocsvn/openrisc/openrisc/trunk

Subversion Repositories openrisc

[/] [openrisc/] [trunk/] [gnu-dev/] [or1k-gcc/] [gcc/] [testsuite/] [gfortran.fortran-torture/] [compile/] [pr39937.f] - Blame information for rev 801

Go to most recent revision | Details | Compare with Previous | View Log


      SUBROUTINE DTREVC( SIDE, HOWMNY, SELECT, N, T, LDT, VL, LDVL, VR,
     $                   LDVR, MM, M, WORK, INFO )
      DOUBLE PRECISION   T( LDT, * ), VL( LDVL, * ), VR( LDVR, * ),
     $                   WORK( * )
      DOUBLE PRECISION   X( 2, 2 )
      CALL DLALN2( .FALSE., 1, 1, SMIN, ONE, T( J, J ),
     $                            ZERO, X, 2, SCALE, XNORM, IERR )
      CALL DSCAL( KI, SCALE, WORK( 1+N ), 1 )
      DO 90 J = KI - 2, 1, -1
      IF( J.GT.JNXT )
     $               GO TO 90
      JNXT = J - 1
      IF( J.GT.1 ) THEN
          IF( T( J, J-1 ).NE.ZERO ) THEN
              IF( WORK( J ).GT.BIGNUM / XNORM ) THEN
                  X( 1, 1 ) = X( 1, 1 ) / XNORM
              END IF
          END IF
          CALL DLALN2( .FALSE., 2, 2, SMIN, ONE,
     $                            T( J-1, J-1 ), LDT, ONE, ONE,
     $                            XNORM, IERR )
          CALL DAXPY( J-2, -X( 1, 1 ), T( 1, J-1 ), 1,
     $                           WORK( 1+N ), 1 )
          CALL DAXPY( J-2, -X( 2, 2 ), T( 1, J ), 1,
     $                           WORK( 1+N2 ), 1 )
      END IF
   90          CONTINUE
      END

Line No.	Rev	Author	Line
1	695	jeremybenn	`SUBROUTINE DTREVC( SIDE, HOWMNY, SELECT, N, T, LDT, VL, LDVL, VR,`
2			`$ LDVR, MM, M, WORK, INFO )`
3			`DOUBLE PRECISION T( LDT, * ), VL( LDVL, * ), VR( LDVR, * ),`
4			`$ WORK( * )`
5			`DOUBLE PRECISION X( 2, 2 )`
6			`CALL DLALN2( .FALSE., 1, 1, SMIN, ONE, T( J, J ),`
7			`$ ZERO, X, 2, SCALE, XNORM, IERR )`
8			`CALL DSCAL( KI, SCALE, WORK( 1+N ), 1 )`
9			`DO 90 J = KI - 2, 1, -1`
10			`IF( J.GT.JNXT )`
11			`$ GO TO 90`
12			`JNXT = J - 1`
13			`IF( J.GT.1 ) THEN`
14			`IF( T( J, J-1 ).NE.ZERO ) THEN`
15			`IF( WORK( J ).GT.BIGNUM / XNORM ) THEN`
16			`X( 1, 1 ) = X( 1, 1 ) / XNORM`
17			`END IF`
18			`END IF`
19			`CALL DLALN2( .FALSE., 2, 2, SMIN, ONE,`
20			`$ T( J-1, J-1 ), LDT, ONE, ONE,`
21			`$ XNORM, IERR )`
22			`CALL DAXPY( J-2, -X( 1, 1 ), T( 1, J-1 ), 1,`
23			`$ WORK( 1+N ), 1 )`
24			`CALL DAXPY( J-2, -X( 2, 2 ), T( 1, J ), 1,`
25			`$ WORK( 1+N2 ), 1 )`
26			`END IF`
27			`90 CONTINUE`
28			`END`

Browse

Tools

Subversion Repositories openrisc

[/] [openrisc/] [trunk/] [gnu-dev/] [or1k-gcc/] [gcc/] [testsuite/] [gfortran.fortran-torture/] [compile/] [pr39937.f] - Blame information for rev 801