35 |
|
end type lj_atype |
36 |
|
|
37 |
|
!! Pointer type for atype ident array |
38 |
< |
type :: lj_atypePtr |
38 |
> |
type, public :: lj_atypePtr |
39 |
|
type (lj_atype), pointer :: this => null() |
40 |
|
end type lj_atypePtr |
41 |
|
|
292 |
|
real ( kind = dp ), dimension(ndim,) :: q |
293 |
|
real ( kind = dp ), dimension(ndim,nLRparticles) :: f |
294 |
|
real ( kind = dp ) :: potE |
295 |
< |
|
295 |
> |
logical ( kind = 2) :: do_pot |
296 |
|
#ifdef MPI |
297 |
|
real( kind = DP ), dimension(3,ncol) :: efr |
298 |
|
real( kind = DP ) :: pot_local |
302 |
|
|
303 |
|
real( kind = DP ) :: pe |
304 |
|
logical, :: update_nlist |
305 |
– |
logical :: do_pot |
305 |
|
|
306 |
+ |
|
307 |
|
integer :: i, j, jbeg, jend, jnab, idim, jdim, idim2, jdim2, dim, dim2 |
308 |
|
integer :: nlist |
309 |
|
integer :: j_start |
325 |
|
|
326 |
|
call check(update_nlist) |
327 |
|
|
328 |
– |
do_pot = .false. |
329 |
– |
if (present(pe)) do_pot = .true. |
330 |
– |
|
328 |
|
#ifndef IS_MPI |
329 |
|
nloops = nloops + 1 |
330 |
|
pot = 0.0E0_DP |
546 |
|
|
547 |
|
if (do_pot) then |
548 |
|
! scatter/gather pot_row into the members of my column |
549 |
< |
call scatter(e_row,e_tmp,plan_row) |
549 |
> |
call scatter(e_row,e_tmp,plan_row) |
550 |
|
|
551 |
|
! scatter/gather pot_local into all other procs |
552 |
|
! add resultant to get total pot |