| 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 :: 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 |
