File: forall1.f90

package info (click to toggle)
lfortran 0.45.0-1
  • links: PTS, VCS
  • area: main
  • in suites: sid, trixie
  • size: 46,332 kB
  • sloc: cpp: 137,068; f90: 51,260; python: 6,444; ansic: 4,277; yacc: 2,285; fortran: 806; sh: 524; makefile: 30; javascript: 15
file content (40 lines) | stat: -rw-r--r-- 1,054 bytes parent folder | download
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
31
32
33
34
35
36
37
38
39
40
program forall1
implicit none
! Only for Syntax check

!>>>>>>>>>>>>> Single line statments <<<<<<<<<<<<<<
    FORALL (i=1:n) A(i,i) = B(i)
    forall (i=1:n:3, j=1:n:5) A(i,j) = SIN(A(j,i))
    forall (i=1:n, j=1:n, i/=j) A(i,j) = REAL(i+j)
    forall (i=1:n, j=1:m, A(i,j).NE.0) &
           A(i,j) = 1/A(i,j)
    forall (i=1:1000, j=1:1000, i /= j) A(i,j) = A(j,i)

!>>>>>>>>>>>>> Multipe line statements <<<<<<<<<<<<<<
    forall (j = 1:n) shared(i) local(x) default(none)
        forall (i=1:j) A(i,j) = B(i)
    end forall

    forall (i = 1:N) reduce(*: s)
        s = s + a(i)
    end forall

    FORALL(i = 3:N + 1, j = 3:N + 1)
        C(i, j) = C(i, j + 2) + C(i, j - 2)
        D(i, j) = C(i, j) + C(i + 2, j) + C(i - 2, j)
    end FORALL

    forall (x=1:100, J(x)>0)
        where (I(x,:)<0)
            I(x,:)=0
        elsewhere
            I(x,:)=1
        end where
    end forall

    outer: forall (i=1:100)
        inner: forall (j=1:100, i.NE.j)
            A(i,j) = A(j,i)
        end forall inner
    end forall outer
end program