1- LAMMPS (4 Jan 2019 )
1+ LAMMPS (24 Aug 2020 )
22units real
33boundary p p p
44
@@ -11,7 +11,8 @@ pair_style lj/cut/coul/long 12.0 12.0
1111kspace_style pppm 1.0e-3
1212
1313read_data data.swm4-ndp
14- orthogonal box = (-12 -12 -12) to (12 12 12)
14+ Reading data file ...
15+ orthogonal box = (-12.000000 -12.000000 -12.000000) to (12.000000 12.000000 12.000000)
1516 2 by 2 by 2 MPI processor grid
1617 reading atoms ...
1718 2500 atoms
@@ -23,10 +24,15 @@ read_data data.swm4-ndp
2324 2000 bonds
2425 reading angles ...
2526 500 angles
26- 4 = max # of 1-2 neighbors
27- 3 = max # of 1-3 neighbors
28- 3 = max # of 1-4 neighbors
29- 4 = max # of special neighbors
27+ Finding 1-2 1-3 1-4 neighbors ...
28+ special bond factors lj: 0.0 0.0 0.5
29+ special bond factors coul: 0.0 0.0 0.5
30+ 4 = max # of 1-2 neighbors
31+ 3 = max # of 1-3 neighbors
32+ 3 = max # of 1-4 neighbors
33+ 4 = max # of special neighbors
34+ special bonds CPU = 0.001 seconds
35+ read_data CPU = 0.018 seconds
3036
3137pair_coeff 1 1 0.210939 3.183950 # ODw ODw
3238pair_coeff * 2* 0.000000 0.0
@@ -54,14 +60,15 @@ velocity DRUDES create ${TDK} 12345
5460velocity DRUDES create 1 12345
5561
5662delete_bonds ATOMS multi
63+ System init for delete_bonds ...
5764PPPM initialization ...
58- using 12-bit tables for long-range coulomb (.. /kspace.cpp:321 )
59- G vector (1/distance) = 0.217457
65+ using 12-bit tables for long-range coulomb (src /kspace.cpp:328 )
66+ G vector (1/distance) = 0.21745695
6067 grid = 10 10 10
6168 stencil order = 5
62- estimated absolute RMS force accuracy = 0.197288
63- estimated relative force accuracy = 0.000594128
64- using double precision FFTs
69+ estimated absolute RMS force accuracy = 0.19728846
70+ estimated relative force accuracy = 0.00059412835
71+ using double precision FFTW3
6572 3d grid and FFT values/proc = 1000 200
6673Rebuild special list taking Drude particles into account
6774Old max number of 1-2 to 1-4 neighbors: 4
@@ -78,6 +85,7 @@ Neighbor list info ...
7885 pair build: half/bin/newton
7986 stencil: half/bin/3d/newton
8087 bin: standard
88+ Deleting bonds ...
8189 2000 total bonds, 500 turned on, 1500 turned off
8290 500 total angles, 0 turned on, 500 turned off
8391 0 total dihedrals, 0 turned on, 0 turned off
@@ -93,10 +101,11 @@ fix RIGID ATOMS rigid/npt/small molecule temp 300 ${TK} 100 iso ${PBAR} ${PBAR}
93101fix RIGID ATOMS rigid/npt/small molecule temp 300 300 100 iso ${PBAR} ${PBAR} 500
94102fix RIGID ATOMS rigid/npt/small molecule temp 300 300 100 iso 1 ${PBAR} 500
95103fix RIGID ATOMS rigid/npt/small molecule temp 300 300 100 iso 1 1 500
96- 500 rigid bodies with 2000 atoms
97- 0.832099 = max distance from body owner to body atom
104+ create bodies CPU = 0.000 seconds
105+ 500 rigid bodies with 2000 atoms
106+ 0.83209864 = max distance from body owner to body atom
98107fix_modify RIGID temp TATOM press thermo_press
99- WARNING: Temperature for fix modify is not for group all (../ fix_rigid_nh_small.cpp:1338 )
108+ WARNING: Temperature for fix modify is not for group all (src/RIGID/ fix_rigid_nh_small.cpp:1336 )
100109fix NVT DRUDES nvt temp ${TDK} ${TDK} 100.0
101110fix NVT DRUDES nvt temp 1 ${TDK} 100.0
102111fix NVT DRUDES nvt temp 1 1 100.0
@@ -107,45 +116,45 @@ thermo 20
107116
108117run 20
109118PPPM initialization ...
110- using 12-bit tables for long-range coulomb (.. /kspace.cpp:321 )
111- G vector (1/distance) = 0.217457
119+ using 12-bit tables for long-range coulomb (src /kspace.cpp:328 )
120+ G vector (1/distance) = 0.21745695
112121 grid = 10 10 10
113122 stencil order = 5
114- estimated absolute RMS force accuracy = 0.197288
115- estimated relative force accuracy = 0.000594128
116- using double precision FFTs
123+ estimated absolute RMS force accuracy = 0.19728846
124+ estimated relative force accuracy = 0.00059412835
125+ using double precision FFTW3
117126 3d grid and FFT values/proc = 1000 200
118- Per MPI rank memory allocation (min/avg/max) = 15.58 | 15.89 | 15.96 Mbytes
127+ Per MPI rank memory allocation (min/avg/max) = 15.82 | 16.13 | 16.21 Mbytes
119128Step CPU TotEng KinEng PotEng E_bond E_angle E_vdwl E_coul E_long Press Volume Temp c_TEMP[1] c_TEMP[2]
120129 0 0 4568.0415 903.41148 3664.63 623.46454 0 2890.2449 157955.66 -157804.74 244591.08 13824 202.18534 295.64644 14.858785
121- 20 0.25570471 3742.8202 5338.828 -1596.0078 567.26788 0 662.20797 155098.52 -157924 31680.274 13979.575 1194.8407 1690.2355 201.6614
122- Loop time of 0.255785 on 8 procs for 20 steps with 2500 atoms
130+ 20 0.28928059 3742.8114 5338.8112 -1595.9999 567.27181 0 662.20936 155098.52 -157924 31679.868 13979.575 1194.8369 1690.2317 201.65774
131+ Loop time of 0.289343 on 8 procs for 20 steps with 2500 atoms
123132
124- Performance: 6.756 ns/day, 3.553 hours/ns, 78.191 timesteps/s
125- 98 .3% CPU use with 8 MPI tasks x no OpenMP threads
133+ Performance: 5.972 ns/day, 4.019 hours/ns, 69.122 timesteps/s
134+ 93 .3% CPU use with 8 MPI tasks x no OpenMP threads
126135
127136MPI task timing breakdown:
128137Section | min time | avg time | max time |%varavg| %total
129138---------------------------------------------------------------
130- Pair | 0.12109 | 0.13272 | 0.13908 | 1.8 | 51.89
131- Bond | 5.4672e -05 | 5.9783e -05 | 6.5615e -05 | 0.0 | 0.02
132- Kspace | 0.01205 | 0.018419 | 0.030091 | 4.8 | 7.20
133- Neigh | 0.0098117 | 0.009816 | 0.0098189 | 0.0 | 3.84
134- Comm | 0.013448 | 0.013578 | 0.013707 | 0.1 | 5.31
135- Output | 5.9536e -05 | 6.7641e-05 | 0.00010464 | 0.0 | 0.03
136- Modify | 0.080548 | 0.080682 | 0.080814 | 0.0 | 31.54
137- Other | | 0.0004443 | | | 0.17
138-
139- Nlocal: 312.5 ave 326 max 292 min
139+ Pair | 0.1288 | 0.14409 | 0.17346 | 3.3 | 49.80
140+ Bond | 7.6162e -05 | 8.3223e -05 | 9.4231e -05 | 0.0 | 0.03
141+ Kspace | 0.01032 | 0.039771 | 0.054956 | 6.3 | 13.75
142+ Neigh | 0.010001 | 0.010005 | 0.010013 | 0.0 | 3.46
143+ Comm | 0.012662 | 0.013401 | 0.01375 | 0.4 | 4.63
144+ Output | 7.2818e -05 | 0.00016643 | 0.0002415 | 0.0 | 0.06
145+ Modify | 0.080679 | 0.081196 | 0.082105 | 0.2 | 28.06
146+ Other | | 0.000631 | | | 0.22
147+
148+ Nlocal: 312.500 ave 326 max 292 min
140149Histogram: 1 0 1 0 1 0 1 1 1 2
141- Nghost: 11354.5 ave 11462 max 11209 min
150+ Nghost: 11354.5 ave 11462 max 11209 min
142151Histogram: 1 0 0 2 1 0 0 1 1 2
143- Neighs: 322375 ave 344480 max 285556 min
152+ Neighs: 322375.0 ave 344481 max 285557 min
144153Histogram: 1 0 0 0 0 2 1 2 0 2
145154
146- Total # of neighbors = 2579001
147- Ave neighs/atom = 1031.6
148- Ave special neighs/atom = 4
155+ Total # of neighbors = 2579003
156+ Ave neighs/atom = 1031.6012
157+ Ave special neighs/atom = 4.0000000
149158Neighbor list builds = 1
150159Dangerous builds = 0
151160
0 commit comments