File: pr68565.f90

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-- 832 bytes parent folder | download | duplicates (3)
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
! { dg-do run }
! { dg-options "-floop-nest-optimize" }

MODULE test
  IMPLICIT NONE
  TYPE subset_type
     INTEGER                                    :: ncon_tot
     REAL(KIND=8),DIMENSION(:,:),ALLOCATABLE   :: coeff
  END TYPE
CONTAINS
  SUBROUTINE foo(subset)
      TYPE(subset_type)                        :: subset
      INTEGER :: icon1
      DO icon1=1,subset%ncon_tot
       subset%coeff(:,icon1)=subset%coeff(:,icon1)/&
        SQRT(DOT_PRODUCT(subset%coeff(:,icon1),subset%coeff(:,icon1)))
      END DO
  END SUBROUTINE
END MODULE

USE test
    TYPE(subset_type)                        :: subset
    INTEGER, VOLATILE :: n1=7,n2=4
    ALLOCATE(subset%coeff(n1,n2))
    CALL RANDOM_NUMBER(subset%coeff)
    subset%coeff=subset%coeff-0.5
    subset%ncon_tot=n2
    CALL foo(subset)
    WRITE(6,*) MAXVAL(subset%coeff)
END