add example for bondlist option validation

This commit is contained in:
jrgissing
2022-02-06 16:14:30 -05:00
parent caafe2ff26
commit ac7db5041f
4 changed files with 296 additions and 0 deletions

View File

@ -0,0 +1,32 @@
# three orthogonal periodic CNTs
# demo for replicating triply looped system
# infinite loops in x, y, z
# includes bonded interactions across box corners
# includes bonds, angles, dihedrals, impropers (class2)
units real
boundary p p p
atom_style full
pair_style lj/class2 10
angle_style class2
bond_style class2
dihedral_style class2
improper_style class2
read_data three_periodic_CNTs.data.gz
replicate 2 2 2 bondlist
fix 1 all nve
run 100
# write_restart replicate.restart
# write_data replicate.data

View File

@ -0,0 +1,132 @@
LAMMPS (7 Jan 2022)
# three orthogonal periodic CNTs
# demo for replicating triply looped system
# infinite loops in x, y, z
# includes bonded interactions across box corners
# includes bonds, angles, dihedrals, impropers (class2)
units real
boundary p p p
atom_style full
pair_style lj/class2 10
angle_style class2
bond_style class2
dihedral_style class2
improper_style class2
read_data three_periodic_CNTs.data.gz
Reading data file ...
orthogonal box = (0 0 0) to (80.96 80.96 80.96)
1 by 1 by 1 MPI processor grid
reading atoms ...
3168 atoms
reading velocities ...
3168 velocities
scanning bonds ...
3 = max bonds/atom
scanning angles ...
3 = max angles/atom
scanning dihedrals ...
12 = max dihedrals/atom
scanning impropers ...
1 = max impropers/atom
reading bonds ...
4752 bonds
reading angles ...
9504 angles
reading dihedrals ...
19008 dihedrals
reading impropers ...
3168 impropers
Finding 1-2 1-3 1-4 neighbors ...
special bond factors lj: 0 0 0
special bond factors coul: 0 0 0
3 = max # of 1-2 neighbors
6 = max # of 1-3 neighbors
18 = max # of 1-4 neighbors
18 = max # of special neighbors
special bonds CPU = 0.003 seconds
read_data CPU = 0.115 seconds
replicate 2 2 2 bondlist
Replicating atoms ...
orthogonal box = (0 0 0) to (161.92 161.92 161.92)
1 by 1 by 1 MPI processor grid
bounding box image = (0 0 0) to (0 0 0)
bounding box extra memory = 1.62 MB
average # of replicas added to proc = 8.00 out of 8 (100.00%)
25344 atoms
38016 bonds
76032 angles
152064 dihedrals
25344 impropers
Finding 1-2 1-3 1-4 neighbors ...
special bond factors lj: 0 0 0
special bond factors coul: 0 0 0
3 = max # of 1-2 neighbors
6 = max # of 1-3 neighbors
18 = max # of 1-4 neighbors
18 = max # of special neighbors
special bonds CPU = 0.021 seconds
replicate CPU = 0.052 seconds
fix 1 all nve
run 100
generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
update every 1 steps, delay 10 steps, check yes
max neighbors/atom: 2000, page size: 100000
master list distance cutoff = 12
ghost atom cutoff = 12
binsize = 6, bins = 27 27 27
1 neighbor lists, perpetual/occasional/extra = 1 0 0
(1) pair lj/class2, perpetual
attributes: half, newton on
pair build: half/bin/newton
stencil: half/bin/3d
bin: standard
WARNING: Inconsistent image flags (../domain.cpp:814)
Per MPI rank memory allocation (min/avg/max) = 51.87 | 51.87 | 51.87 Mbytes
Step Temp E_pair E_mol TotEng Press
0 0 -14266.189 1466925.5 1452659.3 -29908.753
100 2155.9128 -17224.188 1306769.8 1452409 1985.2082
Loop time of 15.0972 on 1 procs for 100 steps with 25344 atoms
Performance: 0.572 ns/day, 41.937 hours/ns, 6.624 timesteps/s
100.0% CPU use with 1 MPI tasks x no OpenMP threads
MPI task timing breakdown:
Section | min time | avg time | max time |%varavg| %total
---------------------------------------------------------------
Pair | 3.7175 | 3.7175 | 3.7175 | 0.0 | 24.62
Bond | 11.222 | 11.222 | 11.222 | 0.0 | 74.33
Neigh | 0.1032 | 0.1032 | 0.1032 | 0.0 | 0.68
Comm | 0.014487 | 0.014487 | 0.014487 | 0.0 | 0.10
Output | 0.00010894 | 0.00010894 | 0.00010894 | 0.0 | 0.00
Modify | 0.027112 | 0.027112 | 0.027112 | 0.0 | 0.18
Other | | 0.01274 | | | 0.08
Nlocal: 25344 ave 25344 max 25344 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost: 22737 ave 22737 max 22737 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs: 2.89358e+06 ave 2.89358e+06 max 2.89358e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Total # of neighbors = 2893576
Ave neighs/atom = 114.17203
Ave special neighs/atom = 18
Neighbor list builds = 1
Dangerous builds = 0
# write_restart replicate.restart
# write_data replicate.data
Total wall time: 0:00:15

View File

@ -0,0 +1,132 @@
LAMMPS (7 Jan 2022)
# three orthogonal periodic CNTs
# demo for replicating triply looped system
# infinite loops in x, y, z
# includes bonded interactions across box corners
# includes bonds, angles, dihedrals, impropers (class2)
units real
boundary p p p
atom_style full
pair_style lj/class2 10
angle_style class2
bond_style class2
dihedral_style class2
improper_style class2
read_data three_periodic_CNTs.data.gz
Reading data file ...
orthogonal box = (0 0 0) to (80.96 80.96 80.96)
1 by 2 by 2 MPI processor grid
reading atoms ...
3168 atoms
reading velocities ...
3168 velocities
scanning bonds ...
3 = max bonds/atom
scanning angles ...
3 = max angles/atom
scanning dihedrals ...
12 = max dihedrals/atom
scanning impropers ...
1 = max impropers/atom
reading bonds ...
4752 bonds
reading angles ...
9504 angles
reading dihedrals ...
19008 dihedrals
reading impropers ...
3168 impropers
Finding 1-2 1-3 1-4 neighbors ...
special bond factors lj: 0 0 0
special bond factors coul: 0 0 0
3 = max # of 1-2 neighbors
6 = max # of 1-3 neighbors
18 = max # of 1-4 neighbors
18 = max # of special neighbors
special bonds CPU = 0.001 seconds
read_data CPU = 0.315 seconds
replicate 2 2 2 bondlist
Replicating atoms ...
orthogonal box = (0 0 0) to (161.92 161.92 161.92)
1 by 2 by 2 MPI processor grid
bounding box image = (0 0 0) to (0 0 0)
bounding box extra memory = 1.62 MB
average # of replicas added to proc = 4.50 out of 8 (56.25%)
25344 atoms
38016 bonds
76032 angles
152064 dihedrals
25344 impropers
Finding 1-2 1-3 1-4 neighbors ...
special bond factors lj: 0 0 0
special bond factors coul: 0 0 0
3 = max # of 1-2 neighbors
6 = max # of 1-3 neighbors
18 = max # of 1-4 neighbors
18 = max # of special neighbors
special bonds CPU = 0.007 seconds
replicate CPU = 0.020 seconds
fix 1 all nve
run 100
generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
update every 1 steps, delay 10 steps, check yes
max neighbors/atom: 2000, page size: 100000
master list distance cutoff = 12
ghost atom cutoff = 12
binsize = 6, bins = 27 27 27
1 neighbor lists, perpetual/occasional/extra = 1 0 0
(1) pair lj/class2, perpetual
attributes: half, newton on
pair build: half/bin/newton
stencil: half/bin/3d
bin: standard
WARNING: Inconsistent image flags (../domain.cpp:814)
Per MPI rank memory allocation (min/avg/max) = 28.69 | 28.69 | 28.69 Mbytes
Step Temp E_pair E_mol TotEng Press
0 0 -14266.189 1466925.5 1452659.3 -29908.753
100 2155.9128 -17224.188 1306769.8 1452409 1985.2082
Loop time of 4.13195 on 4 procs for 100 steps with 25344 atoms
Performance: 2.091 ns/day, 11.478 hours/ns, 24.202 timesteps/s
100.0% CPU use with 4 MPI tasks x no OpenMP threads
MPI task timing breakdown:
Section | min time | avg time | max time |%varavg| %total
---------------------------------------------------------------
Pair | 0.99372 | 1.014 | 1.0672 | 3.1 | 24.54
Bond | 2.9449 | 2.9795 | 3.0136 | 2.0 | 72.11
Neigh | 0.026695 | 0.026706 | 0.026716 | 0.0 | 0.65
Comm | 0.012404 | 0.099546 | 0.15425 | 18.7 | 2.41
Output | 4.3822e-05 | 5.0693e-05 | 7.0516e-05 | 0.0 | 0.00
Modify | 0.0074219 | 0.0074887 | 0.0075411 | 0.1 | 0.18
Other | | 0.004616 | | | 0.11
Nlocal: 6336 ave 6336 max 6336 min
Histogram: 4 0 0 0 0 0 0 0 0 0
Nghost: 10558 ave 10558 max 10558 min
Histogram: 4 0 0 0 0 0 0 0 0 0
Neighs: 723394 ave 723394 max 723394 min
Histogram: 4 0 0 0 0 0 0 0 0 0
Total # of neighbors = 2893576
Ave neighs/atom = 114.17203
Ave special neighs/atom = 18
Neighbor list builds = 1
Dangerous builds = 0
# write_restart replicate.restart
# write_data replicate.data
Total wall time: 0:00:04

Binary file not shown.