diff options
Diffstat (limited to 'Carpet/LoopControl/src/loopcontrol_fortran.h')
-rw-r--r-- | Carpet/LoopControl/src/loopcontrol_fortran.h | 64 |
1 files changed, 64 insertions, 0 deletions
diff --git a/Carpet/LoopControl/src/loopcontrol_fortran.h b/Carpet/LoopControl/src/loopcontrol_fortran.h new file mode 100644 index 000000000..e80bf358f --- /dev/null +++ b/Carpet/LoopControl/src/loopcontrol_fortran.h @@ -0,0 +1,64 @@ +/* -*-f90-mode-*- */ + +#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/**/_imax, name/**/_jmax, name/**/_kmax &&\ +integer :: i, j, k + + + +#define LC_PRIVATE3(name) \ +name/**/_lc, \ +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/**/_kmax = min (name/**/_kk - 1 + name/**/_lc%kkstep, name/**/_lc%kkmax) &&\ + do name/**/_jj = name/**/_lc%jjmin + 1, name/**/_lc%jjmax, name/**/_lc%jjstep &&\ + name/**/_jmax = min (name/**/_jj - 1 + name/**/_lc%jjstep, name/**/_lc%jjmax) &&\ + do name/**/_ii = name/**/_lc%iimin + 1, name/**/_lc%iimax, name/**/_lc%iistep &&\ + name/**/_imax = min (name/**/_ii - 1 + name/**/_lc%iistep, name/**/_lc%iimax) &&\ + &&\ + /* Fine loop */ &&\ + do k = name/**/_kk, name/**/_kmax &&\ + do j = name/**/_jj, name/**/_jmax &&\ + do i = name/**/_ii, 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 */ |