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
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
|
/* -*-f90-*- */
#ifndef LOOPCONTROL_FORTRAN_H
#define LOOPCONTROL_FORTRAN_H
#define LC_DECLARE3(name, i,j,k) &&\
type (lc_statmap_t), save :: name/**/_lm &&\
logical, save :: name/**/_initialised = .false. &&\
type (lc_control_t) :: name/**/_lc &&\
integer :: name/**/_ii, name/**/_jj, name/**/_kk &&\
integer :: name/**/_imin, name/**/_jmin, name/**/_kmin &&\
integer :: name/**/_imax, name/**/_jmax, name/**/_kmax &&\
integer :: i, j, k
#define LC_PRIVATE3(name) \
name/**/_lc, \
name/**/_imin, name/**/_jmin, name/**/_kmin, \
name/**/_imax, name/**/_jmax, name/**/_kmax
#define LC_LOOP3(name, i,j,k, imin,jmin,kmin, imax,jmax,kmax, ilsh,jlsh,klsh) &&\
if (.not. name/**/_initialised) then &&\
!$omp single &&\
call lc_statmap_init (name/**/_lm, "name") &&\
!$omp end single &&\
!$omp single &&\
/* Set this flag only after initialising */ &&\
name/**/_initialised = .true. &&\
!$omp end single &&\
end if &&\
call lc_control_init (name/**/_lc, name/**/_lm, imin,jmin,kmin, imax,jmax,kmax, ilsh,jlsh,klsh) &&\
&&\
/* Coarse loop */ &&\
do name/**/_kk = name/**/_lc%kkmin + 1, name/**/_lc%kkmax, name/**/_lc%kkstep &&\
name/**/_kmin = name/**/_kk + name/**/_lc%kkkkmax &&\
name/**/_kmax = min (name/**/_kmin - 1 + name/**/_lc%kkkkstep, name/**/_lc%kkmax) &&\
do name/**/_jj = name/**/_lc%jjmin + 1, name/**/_lc%jjmax, name/**/_lc%jjstep &&\
name/**/_jmin = name/**/_jj + name/**/_lc%jjjjmax &&\
name/**/_jmax = min (name/**/_jj - 1 + name/**/_lc%jjjjmax, name/**/_lc%jjmax) &&\
do name/**/_ii = name/**/_lc%iimin + 1, name/**/_lc%iimax, name/**/_lc%iistep &&\
name/**/_imin = name/**/_ii + name/**/_lc%iiiimax &&\
name/**/_imax = min (name/**/_ii - 1 + name/**/_lc%iiiimax, name/**/_lc%iimax) &&\
&&\
/* Fine loop */ &&\
do k = name/**/_kmin, name/**/_kmax &&\
do j = name/**/_jmin, name/**/_jmax &&\
do i = name/**/_imin, name/**/_imax
#define LC_ENDLOOP3(name) &&\
end do &&\
end do &&\
end do &&\
&&\
end do &&\
end do &&\
end do &&\
&&\
call lc_control_finish (name/**/_lc)
#endif /* #ifndef LOOPCONTROL_FORTRAN_H */
|