| 6 |
|
!! @author Charles F. Vardeman II |
| 7 |
|
!! @author Matthew Meineke |
| 8 |
|
!! @author J. Daniel Gezelter |
| 9 |
< |
!! @version $Id: neighborLists.F90,v 1.1.1.1 2003-03-21 17:42:12 mmeineke Exp $, |
| 9 |
> |
!! @version $Id: neighborLists.F90,v 1.3 2003-04-08 17:16:22 chuckv Exp $, |
| 10 |
|
|
| 11 |
|
module neighborLists |
| 12 |
|
|
| 32 |
|
integer, dimension(:),public, pointer :: list => null() |
| 33 |
|
!! Position array of previous positions for check. Allocated first time |
| 34 |
|
!! into saveNeighborList. |
| 35 |
< |
real( kind = dp ), dimension(:,:), allocatable , save :: q0 |
| 35 |
> |
real( kind = dp ), dimension(:,:), allocatable, save :: q0 |
| 36 |
|
!! Current list size |
| 37 |
|
integer, save :: listSize |
| 38 |
|
!--------------MODULE ACCESS--------------------------> |
| 64 |
|
.not. associated(list) ) then |
| 65 |
|
allocate(point(natoms),stat=alloc_error) |
| 66 |
|
if (alloc_error /= 0) then |
| 67 |
+ |
write(default_error,*) & |
| 68 |
+ |
"ExpandNeighborLists: Error in allocating point" |
| 69 |
|
error = -1 |
| 70 |
|
return |
| 71 |
|
end if |
| 72 |
|
allocate(list(listMultiplier * natoms),stat=alloc_error) |
| 73 |
|
if (alloc_error /= 0) then |
| 74 |
+ |
write(default_error,*) & |
| 75 |
+ |
"ExpandNeighborLists: Error in allocating list" |
| 76 |
|
error = -1 |
| 77 |
|
return |
| 78 |
|
end if |
| 85 |
|
.not. associated(list) ) then |
| 86 |
|
allocate(point(getNRow(plan_row)),stat=alloc_error) |
| 87 |
|
if (alloc_error /= 0) then |
| 88 |
< |
error = -1 |
| 88 |
> |
write(default_error,*) & |
| 89 |
> |
"ExpandNeighborLists: Error in allocating MPI point" |
| 90 |
> |
error = -1 |
| 91 |
|
return |
| 92 |
|
end if |
| 93 |
|
allocate(list(listMultiplier * getNCol(plan_col)),stat=alloc_error) |
| 94 |
|
if (alloc_error /= 0) then |
| 95 |
+ |
write(default_error,*) & |
| 96 |
+ |
"ExpandNeighborLists: Error in allocating MPI list" |
| 97 |
|
error = -1 |
| 98 |
|
return |
| 99 |
|
end if |
| 107 |
|
|
| 108 |
|
! Check to see if we have exceeded the maximum number of allocations. |
| 109 |
|
if (nAllocations > maxAllocations) then |
| 110 |
+ |
write(default_error,*) & |
| 111 |
+ |
"ExpandNeighborList: exceeded maximum number of re-allocations" |
| 112 |
|
error = -1 |
| 113 |
|
return |
| 114 |
|
else !! Expand the list. |
| 115 |
|
oldSize = size(list) |
| 116 |
|
|
| 117 |
+ |
|
| 118 |
|
#ifdef IS_MPI !! MPI |
| 119 |
|
newSize = listMultiplier * getNCol(plan_col) + oldSize |
| 120 |
|
#else |
| 162 |
|
real( kind = DP ) :: dispmx_tmp |
| 163 |
|
real( kind = dp ) :: skin_thickness |
| 164 |
|
integer :: nlocal |
| 165 |
< |
|
| 165 |
> |
|
| 166 |
> |
|
| 167 |
|
nlocal = natoms |
| 168 |
|
skin_thickness = rlist - rcut |
| 169 |
|
dispmx = 0.0E0_DP |
| 174 |
|
update_nlist = .true. |
| 175 |
|
return |
| 176 |
|
end if |
| 177 |
< |
|
| 177 |
> |
|
| 178 |
> |
|
| 179 |
|
#ifdef MPI |
| 180 |
|
|
| 181 |
|
dispmx_tmp = 0.0E0_DP |
| 209 |
|
!! Saves neighbor list for comparison in check. |
| 210 |
|
!! Save_neighborList will work even if the number of |
| 211 |
|
!! local atoms has changed. |
| 212 |
< |
subroutine saveNeighborList(q) |
| 213 |
< |
real(kind = dp ), dimension(:,:), intent(in) :: q |
| 212 |
> |
subroutine saveNeighborList(natoms, q) |
| 213 |
> |
|
| 214 |
> |
integer :: natoms |
| 215 |
> |
real(kind = dp ), dimension(3,natoms), intent(in) :: q |
| 216 |
|
integer :: list_size |
| 217 |
+ |
|
| 218 |
|
|
| 219 |
|
!! get size of list |
| 220 |
< |
list_size = size(q) |
| 220 |
> |
list_size = natoms |
| 221 |
|
|
| 222 |
|
if (.not. allocated(q0)) then |
| 223 |
|
allocate(q0(3,list_size)) |
