File: pr39937.f

package info (click to toggle)
gcc-arm-none-eabi 15%3A12.2.rel1-1
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 959,712 kB
  • sloc: cpp: 3,275,382; ansic: 2,061,766; ada: 840,956; f90: 208,513; makefile: 76,132; asm: 73,433; xml: 50,448; exp: 34,146; sh: 32,436; objc: 15,637; fortran: 14,012; python: 11,991; pascal: 6,787; awk: 4,779; perl: 3,054; yacc: 338; ml: 285; lex: 201; haskell: 122
file content (30 lines) | stat: -rw-r--r-- 1,304 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
 
C { dg-do compile }
C { dg-options "-std=legacy" }
      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 ) ! { dg-warning "Type mismatch" }
      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,  ! { dg-warning "Type mismatch" }
     $                            XNORM, IERR )  ! { dg-warning "Type mismatch" }
          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