-
Notifications
You must be signed in to change notification settings - Fork 0
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.
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
Concave particle pressure validation #61
Conversation
pre-commit.ci autofix |
for more information, see https://pre-commit.ci
pre-commit.ci autofix |
for more information, see https://pre-commit.ci
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Thanks! Please post the analysis results including the nvt_gpu
jobs (run on Great Lakes or Delta).
hoomd_validation/init.py
Outdated
|
||
subprojects = [alj_2d, lj_fluid, lj_union, hard_disk, hard_sphere] | ||
subprojects = [ | ||
alj_2d, lj_fluid, lj_union, hard_disk, hard_sphere, simple_polygon |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
alj_2d, lj_fluid, lj_union, hard_disk, hard_sphere, simple_polygon | |
alj_2d, lj_fluid, lj_union, hard_disk, hard_sphere, simple_polygon, |
And run pre-commit
to format one item per line.
hoomd_validation/simple_polygon.py
Outdated
# Copyright (c) 2022-2023 The Regents of the University of Michigan. | ||
# Part of HOOMD-blue, released under the BSD 3-Clause License. | ||
|
||
"""Hard disk equation of state validation test.""" |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
"""Hard disk equation of state validation test.""" | |
"""Simple polygon equation of state validation test.""" |
hoomd_validation/simple_polygon.py
Outdated
# reference statepoint from: http://dx.doi.org/10.1016/j.jcp.2013.07.023 | ||
# Assume the same pressure for a 128**2 system as this is in the fluid | ||
# from first test: mean +- std_err_means = | ||
# 6.205344838518146 +- 0.0011239315686705335 | ||
# based on nvt sims at density = 0.75 over 32 replicas |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Not all of these comment lines can be true at the same time. Remove the invalid ones.
hoomd_validation/simple_polygon.py
Outdated
RANDOMIZE_STEPS = 20_000 | ||
EQUILIBRATE_STEPS = 100_000 | ||
RUN_STEPS = 2_500_000 | ||
RESTART_STEPS = RUN_STEPS // 10 |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
RESTART_STEPS = RUN_STEPS // 10 | |
RESTART_STEPS = RUN_STEPS // 50 |
Shorten the restart period from ~30 minutes down to ~6 (tested on cheme-hodges).
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good! Thanks!
Description
Performs cross validation using NVT and NPT simulations of a hard particle with a concave shape. Computes the pressure in NVT simulations, and then runs NPT simulations and calculates the density the mean pressure.
Motivation and context
Validates the implementation of the SDF-based pressure calculation where the expansive contribution to the pressure is nonzero (i.e., for non-convex hard particles).
Resolves #43
How has this been tested?
Ran on PSC-Bridges2. Results:
Checklist:
CONTRIBUTORS.md
) in the pull request source branch.