4 |
|
|
5 |
|
!! @author Charles F. Vardeman II |
6 |
|
!! @author Matthew Meineke |
7 |
< |
!! @version $Id: do_Forces.F90,v 1.16 2003-07-01 22:39:53 gezelter Exp $, $Date: 2003-07-01 22:39:53 $, $Name: not supported by cvs2svn $, $Revision: 1.16 $ |
7 |
> |
!! @version $Id: do_Forces.F90,v 1.19 2003-07-14 23:06:09 gezelter Exp $, $Date: 2003-07-14 23:06:09 $, $Name: not supported by cvs2svn $, $Revision: 1.19 $ |
8 |
|
|
9 |
|
module do_Forces |
10 |
|
use force_globals |
524 |
|
integer :: me_i, me_j |
525 |
|
|
526 |
|
r = sqrt(rijsq) |
527 |
– |
|
528 |
– |
|
527 |
|
|
528 |
|
#ifdef IS_MPI |
529 |
|
if (tagRow(i) .eq. tagColumn(j)) then |
587 |
|
endif |
588 |
|
endif |
589 |
|
|
590 |
+ |
|
591 |
+ |
|
592 |
|
end subroutine do_pair |
593 |
|
|
594 |
|
|
608 |
|
if( .not.boxIsOrthorhombic ) then |
609 |
|
! calc the scaled coordinates. |
610 |
|
|
611 |
< |
scaled = matmul(d, HmatInv) |
611 |
> |
scaled = matmul(HmatInv, d) |
612 |
|
|
613 |
|
! wrap the scaled coordinates |
614 |
|
|
615 |
< |
do i = 1, 3 |
616 |
< |
scaled(i) = scaled(i) - anint(scaled(i)) |
617 |
< |
enddo |
615 |
> |
scaled = scaled - anint(scaled) |
616 |
> |
|
617 |
|
|
618 |
|
! calc the wrapped real coordinates from the wrapped scaled |
619 |
|
! coordinates |
620 |
|
|
621 |
< |
d = matmul(scaled,Hmat) |
621 |
> |
d = matmul(Hmat,scaled) |
622 |
|
|
623 |
|
else |
624 |
|
! calc the scaled coordinates. |