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
SNAP bik flag #3142
SNAP bik flag #3142
Conversation
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.
This looks nice and clean. The only problem is that the treatment of the last column (reference potential) is not clearly defined in the code or in the documentation. I suggest making a clear choice e.g. 1) snap all[irow][lastcol] = 0;
for all b_ik rows 2) same, but with reference energy in first b_ik row 3) Put per-atom reference energy in every b_ik row. Whatever choice is made should be explicitly documented.
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.
Still need to describe final column in documentation.
doc/src/compute_sna_atom.rst
Outdated
descriptors :math:`B_{i,k}` summed over all atoms *i* to produce | ||
:math:`B_k`, becomes bispectrum rows equal to the number of atoms. Thus, | ||
the resulting bispectrum rows are :math:`B_{i,k}` instead of just | ||
:math:`B_k`. |
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.
Add "In this case, the entries in the final column for these rows are set to zero."
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.
Done
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.
@athomps @charlessievers as a reminder. If you both consider this to be reading for merging, you need to assign the PR to me.
Summary
Added in flag for snap compute to provide per atom bispectrum descriptors instead of summed bispectrum descriptors.
Author(s)
Charlie Sievers - Sandia National Labs - charliesievers@cox.net
Licensing
By submitting this pull request, I agree, that my contribution will be included in LAMMPS and redistributed under either the GNU General Public License version 2 (GPL v2) or the GNU Lesser General Public License version 2.1 (LGPL v2.1).
Implementation Notes
Tested file examples/snap/in.snap.compute and produced the same output.
Post Submission Checklist