1 |
|
module simulation |
2 |
|
use definitions, ONLY :dp |
3 |
+ |
use force_wrappers, ONLY : alloc_force_wrappers |
4 |
+ |
|
5 |
|
implicit none |
6 |
|
PRIVATE |
7 |
|
|
8 |
|
|
9 |
+ |
|
10 |
+ |
real ( kind = dp ), public, dimension(3) :: box |
11 |
|
|
8 |
– |
real ( kind = dp ), dimension(3) :: box |
12 |
|
|
13 |
+ |
real ( kind = dp ), public :: rlist |
14 |
+ |
real ( kind = dp ), public :: rcut |
15 |
+ |
real ( kind = dp ), public :: rlistsq |
16 |
+ |
real ( kind = dp ), public :: rcutsq |
17 |
|
|
18 |
+ |
integer,public, allocatable, dimension(:) :: point |
19 |
+ |
integer,public, allocatable, dimension(:) :: list |
20 |
|
|
21 |
|
|
22 |
+ |
!MPI dependent routines |
23 |
|
|
24 |
+ |
#ifdef IS_MPI |
25 |
+ |
! Universal routines: All types of force calculations will need these arrays |
26 |
+ |
! Arrays specific to a type of force calculation should be declared in that module. |
27 |
+ |
real( kind = dp ), allocatable, dimension(:,:) :: qRow |
28 |
+ |
real( kind = dp ), allocatable, dimension(:,:) :: qColumn |
29 |
|
|
30 |
+ |
real( kind = dp ), allocatable, dimension(:,:) :: fRow |
31 |
+ |
real( kind = dp ), allocatable, dimension(:,:) :: fColumn |
32 |
|
|
33 |
+ |
real( kind = dp ), allocatable, dimension(:,:) :: tRow |
34 |
+ |
real( kind = dp ), allocatable, dimension(:,:) :: tColumn |
35 |
+ |
|
36 |
+ |
#endif |
37 |
+ |
|
38 |
+ |
|
39 |
+ |
|
40 |
+ |
|
41 |
+ |
|
42 |
+ |
|
43 |
|
contains |
44 |
|
|
45 |
+ |
#ifdef MPI |
46 |
+ |
! Allocated work arrays for MPI |
47 |
+ |
subroutine allocate_mpi_arrays(nDimensions,numComponents) |
48 |
+ |
integer, intent(in) :: nDimensions |
49 |
+ |
integer, intent(in) :: numComponents |
50 |
|
|
51 |
+ |
|
52 |
+ |
|
53 |
+ |
|
54 |
+ |
|
55 |
+ |
end subroutine allocate_mpi_arrays |
56 |
+ |
#endif |
57 |
+ |
|
58 |
|
subroutine set_simulation(box,rlist,rcut) |
59 |
|
|
60 |
|
|
70 |
|
|
71 |
|
end subroutine change_box_size |
72 |
|
|
73 |
< |
|
73 |
> |
|
74 |
|
end module simulation |