43 |
|
!! Calculates Long Range forces Lennard-Jones interactions. |
44 |
|
!! @author Charles F. Vardeman II |
45 |
|
!! @author Matthew Meineke |
46 |
< |
!! @version $Id: LJ.F90,v 1.19 2005-11-21 22:59:01 gezelter Exp $, $Date: 2005-11-21 22:59:01 $, $Name: not supported by cvs2svn $, $Revision: 1.19 $ |
46 |
> |
!! @version $Id: LJ.F90,v 1.20 2005-12-07 19:58:18 chuckv Exp $, $Date: 2005-12-07 19:58:18 $, $Name: not supported by cvs2svn $, $Revision: 1.20 $ |
47 |
|
|
48 |
|
|
49 |
|
module lj |
51 |
|
use vector_class |
52 |
|
use simulation |
53 |
|
use status |
54 |
+ |
use fForceOptions |
55 |
|
#ifdef IS_MPI |
56 |
|
use mpiSimulation |
57 |
|
#endif |
105 |
|
public :: setLJDefaultCutoff |
106 |
|
public :: getSigma |
107 |
|
public :: getEpsilon |
107 |
– |
public :: useGeometricMixing |
108 |
|
public :: do_lj_pair |
109 |
|
public :: destroyLJtypes |
110 |
|
|
193 |
|
|
194 |
|
end function getEpsilon |
195 |
|
|
196 |
– |
subroutine useGeometricMixing() |
197 |
– |
useGeometricDistanceMixing = .true. |
198 |
– |
haveMixingMap = .false. |
199 |
– |
return |
200 |
– |
end subroutine useGeometricMixing |
201 |
– |
|
196 |
|
subroutine createMixingMap() |
197 |
|
integer :: nLJtypes, i, j |
198 |
|
real ( kind = dp ) :: s1, s2, e1, e2 |
210 |
|
allocate(MixingMap(nLJtypes, nLJtypes)) |
211 |
|
endif |
212 |
|
|
213 |
+ |
useGeometricDistanceMixing = usesGeometricDistanceMixing() |
214 |
|
do i = 1, nLJtypes |
215 |
|
|
216 |
|
s1 = LJMap%LJtypes(i)%sigma |