Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Collected small fixes and updates #3976

Merged
merged 11 commits into from
Nov 18, 2023
Merged

Collected small fixes and updates #3976

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
11 commits
Select commit Hold shift + click to select a range
1f509d9
make output more consistent with other LAMMPS code
akohlmey Nov 16, 2023
7f2a6b1
update h2 fix pimd/langevin example
akohlmey Nov 16, 2023
ba1c848
Merge branch 'develop' into numerial_propagation
akohlmey Nov 16, 2023
ac5f287
explicitly include `if constexpr` only when C++17 or later is selected
akohlmey Nov 16, 2023
ddd5cc1
avoid segfault on command errors in force style unit tests and print …
akohlmey Nov 16, 2023
d739faf
disable `if constexpr` for all Intel compilers
akohlmey Nov 17, 2023
ea8234a
treat immediate variables like equal style also for evaluating refere…
akohlmey Nov 17, 2023
e59a665
add tests for immediate variables
akohlmey Nov 17, 2023
a4da559
check if creating unix domain socket failed
akohlmey Nov 17, 2023
ab183eb
Merge remote-tracking branch 'github/develop' into collected-small-fixes
akohlmey Nov 17, 2023
5223f2d
Merge remote-tracking branch 'akohlmey/numerial_propagation' into col…
akohlmey Nov 17, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
28 changes: 28 additions & 0 deletions examples/PACKAGES/pimd/langevin_metal_units/in.langevin.metal
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,28 @@
variable ibead uloop 99 pad

units metal
atom_style atomic
atom_modify map yes
boundary p p p
pair_style lj/cut 9.5251
read_data data.metalnpt${ibead}

pair_coeff * * 0.00965188 3.4
pair_modify shift yes

mass 1 39.948

timestep 0.001

velocity all create 0.0 ${ibead}

fix 1 all pimd/langevin method nmpimd ensemble npt integrator obabo thermostat PILE_L 1234 tau 1.0 temp 113.15 iso 1.0 barostat BZP taup 1.0 fixcom no

thermo_style custom step temp f_1[*] vol press
thermo 100
thermo_modify norm no

# dump dcd all custom 100 ${ibead}.dcd id type xu yu zu vx vy vz ix iy iz fx fy fz
# dump_modify dcd sort id format line "%d %d %.16f %.16f %.16f %.16f %.16f %.16f %d %d %d %.16f %.16f %.16f"

run 1000
2 changes: 1 addition & 1 deletion ...PACKAGES/pimd/langevin_metal_units/in.lmp → ...ngevin_metal_units/in.pimd-langevin.metal
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,7 +16,7 @@ timestep 0.001

velocity all create 0.0 ${ibead}

fix 1 all pimd/langevin ensemble npt integrator obabo thermostat PILE_L 1234 tau 1.0 temp 113.15 iso 1.0 barostat BZP taup 1.0 fixcom no
fix 1 all pimd/langevin method pimd ensemble nvt integrator obabo thermostat PILE_L 1234 tau 1.0 temp 113.15 taup 1.0 fixcom no

thermo_style custom step temp f_1[*] vol press
thermo 100
Expand Down
2 changes: 1 addition & 1 deletion ...etal_units/log.14Jun23.langevin.metal.g++ → ...etal_units/log.16Nov23.langevin.metal.g++
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,2 +1,2 @@
LAMMPS (28 Mar 2023)
LAMMPS (3 Aug 2023)
Running on 4 partitions of processors
62 changes: 31 additions & 31 deletions ...al_units/log.14Jun23.langevin.metal.g++.0 → ...al_units/log.16Nov23.langevin.metal.g++.0
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,5 +1,6 @@
LAMMPS (28 Mar 2023)
LAMMPS (3 Aug 2023)
Processor partition = 0
OMP_NUM_THREADS environment is not set. Defaulting to 1 thread. (src/comm.cpp:98)
using 1 OpenMP thread(s) per MPI task
variable ibead uloop 99 pad

Expand Down Expand Up @@ -29,10 +30,10 @@ timestep 0.001
velocity all create 0.0 ${ibead}
velocity all create 0.0 01

fix 1 all pimd/langevin ensemble npt integrator obabo thermostat PILE_L 1234 tau 1.0 temp 113.15 iso 1.0 barostat BZP taup 1.0 fixcom no
fix 1 all pimd/langevin method nmpimd ensemble npt integrator obabo thermostat PILE_L 1234 tau 1.0 temp 113.15 iso 1.0 barostat BZP taup 1.0 fixcom no

Initializing PIMD BZP barostat...
The barostat mass is W = 2.3401256650800001e+01
The barostat mass is W = 2.3401256650800001e+01

thermo_style custom step temp f_1[*] vol press
thermo 100
Expand All @@ -43,14 +44,13 @@ thermo_modify norm no

run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule

Initializing PI Langevin equation thermostat...
Bead ID | omega | tau | c1 | c2
0 0.00000000e+00 1.00000000e+00 9.99500125e-01 3.16148726e-02
1 8.37986825e+01 5.96668092e-03 9.19616372e-01 3.92817678e-01
2 1.18509233e+02 4.21908054e-03 8.88243614e-01 4.59372705e-01
3 8.37986825e+01 5.96668092e-03 9.19616372e-01 3.92817678e-01
PILE_L thermostat successfully initialized!
Bead ID | omega | tau | c1 | c2
0 0.00000000e+00 1.00000000e+00 9.99500125e-01 3.16148726e-02
1 8.37986825e+01 5.96668092e-03 9.19616372e-01 3.92817678e-01
2 1.18509233e+02 4.21908054e-03 8.88243614e-01 4.59372705e-01
3 8.37986825e+01 5.96668092e-03 9.19616372e-01 3.92817678e-01
PILE_L thermostat successfully initialized!

Neighbor list info ...
update: every = 1 steps, delay = 0 steps, check = yes
Expand All @@ -66,31 +66,31 @@ Neighbor list info ...
bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.121 | 3.121 | 3.121 Mbytes
Step Temp f_1[1] f_1[2] f_1[3] f_1[4] f_1[5] f_1[6] f_1[7] f_1[8] f_1[9] f_1[10] f_1[11] f_1[12] f_1[13] f_1[14] f_1[15] Volume Press
0 0 0 0 -7.3046601 4.3005229 -21.877018 -8.7249482 2.9571502 -1743.5332 -698.49808 -172.07477 0 0 0.033460054 -0.37064378 4.216227 13402.228 -164.47373
100 149.95804 3.8573359 0 -7.7921375 42.886648 -23.396327 -1.980193 2.954003 -1888.0547 -1648.7118 -332.0298 -0.099139345 0.11500091 0.033044702 -0.3701566 42.83112 13235.861 -101.30374
200 245.00113 6.3021074 0 -8.2639651 41.690123 -22.521598 -4.273021 2.9600599 -1906.2904 -1609.02 -265.94404 -0.20527926 0.49305948 0.031504957 -0.36829556 41.729191 12619.125 112.22426
300 300.57486 7.7316177 0 -8.2986331 43.180131 -21.755813 -7.7032433 2.9714114 -1968.7685 290.49656 251.72564 -0.21935745 0.56300721 0.029467915 -0.36568855 43.236828 11803.2 814.45889
400 368.08438 9.4681493 0 -8.4800193 49.109699 -24.824142 2.9744597 2.9794185 -2335.993 1368.7398 570.03286 -0.028366234 0.0094148316 0.028338146 -0.36416383 49.028096 11350.678 1202.0398
500 419.32066 10.786088 0 -8.640773 45.427771 -22.825143 16.22356 2.9684828 -2113.91 -272.84753 185.53392 0.091614289 0.098205455 0.028793585 -0.36478567 45.368325 11533.101 952.59748
600 385.4127 9.9138817 0 -8.4356035 47.783726 -22.456104 6.837575 2.967236 -2023.8117 -918.27943 -2.4106994 0.093360761 0.10198539 0.029589188 -0.36584873 47.725157 11851.775 676.62913
700 360.14242 9.2638601 0 -8.2900275 42.626187 -20.571698 -5.7252564 2.9560528 -1806.9448 -1418.2247 -148.41657 0.075011202 0.065835696 0.030359455 -0.36685105 42.558523 12160.301 456.91446
800 346.92167 8.923786 0 -8.0694169 45.160336 -21.885719 -6.7745694 2.9575472 -1894.3641 -1329.3179 -136.42193 0.011114896 0.0014455064 0.030808183 -0.3674233 45.076543 12340.037 454.60123
900 364.39442 9.3732334 0 -8.0415668 45.604542 -21.816625 5.586068 2.9578604 -1890.4653 -1271.1107 -111.89061 -0.020285587 0.0048148677 0.030774258 -0.36738033 45.521594 12326.448 499.75868
1000 390.77042 10.051697 0 -8.1948009 45.264242 -22.833545 6.9260573 2.960122 -2007.6188 -1179.7125 -70.907567 -0.062733519 0.046047757 0.030329191 -0.36681215 45.191633 12148.179 572.98799
Loop time of 0.248186 on 1 procs for 1000 steps with 200 atoms

Performance: 348.126 ns/day, 0.069 hours/ns, 4029.238 timesteps/s, 805.848 katom-step/s
99.6% CPU use with 1 MPI tasks x 1 OpenMP threads
0 0 0 0 -7.3046601 4.3005229 3.3062167 -8.7249482 2.9571502 -1743.5332 -698.49808 -172.07477 0 0 0.033460054 -0.37064378 4.216227 13402.228 -164.47373
100 149.95804 3.8573359 0 -7.7921375 42.886648 2.9263894 -1.980193 2.954003 -1888.0547 -1648.7118 -332.0298 -0.099139345 0.11500091 0.033044702 -0.3701566 42.83112 13235.861 -101.30374
200 245.00113 6.3021074 0 -8.2639651 41.690123 3.1450717 -4.273021 2.9600599 -1906.2904 -1609.02 -265.94404 -0.20527926 0.49305948 0.031504957 -0.36829556 41.729191 12619.125 112.22426
300 300.57486 7.7316177 0 -8.2986331 43.180131 3.336518 -7.7032433 2.9714114 -1968.7685 290.49656 251.72564 -0.21935745 0.56300721 0.029467915 -0.36568855 43.236828 11803.2 814.45889
400 368.08438 9.4681493 0 -8.4800193 49.109699 2.5694358 2.9744597 2.9794185 -2335.993 1368.7398 570.03286 -0.028366234 0.0094148316 0.028338146 -0.36416383 49.028096 11350.678 1202.0398
500 419.32066 10.786088 0 -8.640773 45.427771 3.0691855 16.22356 2.9684828 -2113.91 -272.84753 185.53392 0.091614289 0.098205455 0.028793585 -0.36478567 45.368325 11533.101 952.59748
600 385.4127 9.9138817 0 -8.4356035 47.783726 3.1614452 6.837575 2.967236 -2023.8117 -918.27943 -2.4106994 0.093360761 0.10198539 0.029589188 -0.36584873 47.725157 11851.775 676.62913
700 360.14242 9.2638601 0 -8.2900275 42.626187 3.6325468 -5.7252564 2.9560528 -1806.9448 -1418.2247 -148.41657 0.075011202 0.065835696 0.030359455 -0.36685105 42.558523 12160.301 456.91446
800 346.92167 8.923786 0 -8.0694169 45.160336 3.3040415 -6.7745694 2.9575472 -1894.3641 -1329.3179 -136.42193 0.011114896 0.0014455064 0.030808183 -0.3674233 45.076543 12340.037 454.60123
900 364.39442 9.3732334 0 -8.0415668 45.604542 3.321315 5.586068 2.9578604 -1890.4653 -1271.1107 -111.89061 -0.020285587 0.0048148677 0.030774258 -0.36738033 45.521594 12326.448 499.75868
1000 390.77042 10.051697 0 -8.1948009 45.264242 3.0670849 6.9260573 2.960122 -2007.6188 -1179.7125 -70.907567 -0.062733519 0.046047757 0.030329191 -0.36681215 45.191633 12148.179 572.98799
Loop time of 0.218442 on 1 procs for 1000 steps with 200 atoms

Performance: 395.527 ns/day, 0.061 hours/ns, 4577.865 timesteps/s, 915.573 katom-step/s
98.9% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section | min time | avg time | max time |%varavg| %total
---------------------------------------------------------------
Pair | 0.14541 | 0.14541 | 0.14541 | 0.0 | 58.59
Neigh | 0.00099082 | 0.00099082 | 0.00099082 | 0.0 | 0.40
Comm | 0.0039966 | 0.0039966 | 0.0039966 | 0.0 | 1.61
Output | 0.00016346 | 0.00016346 | 0.00016346 | 0.0 | 0.07
Modify | 0.096205 | 0.096205 | 0.096205 | 0.0 | 38.76
Other | | 0.001425 | | | 0.57
Pair | 0.11918 | 0.11918 | 0.11918 | 0.0 | 54.56
Neigh | 0.0010314 | 0.0010314 | 0.0010314 | 0.0 | 0.47
Comm | 0.0046197 | 0.0046197 | 0.0046197 | 0.0 | 2.11
Output | 0.0001329 | 0.0001329 | 0.0001329 | 0.0 | 0.06
Modify | 0.092616 | 0.092616 | 0.092616 | 0.0 | 42.40
Other | | 0.0008616 | | | 0.39

Nlocal: 200 ave 200 max 200 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Expand Down