-
Notifications
You must be signed in to change notification settings - Fork 5
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
Add Schwenke's gas phase H2O analytical potential #166
Conversation
Codecov ReportAttention:
Additional details and impacted files@@ Coverage Diff @@
## master #166 +/- ##
==========================================
- Coverage 92.37% 92.28% -0.09%
==========================================
Files 43 45 +2
Lines 6084 6182 +98
Branches 741 746 +5
==========================================
+ Hits 5620 5705 +85
- Misses 452 465 +13
Partials 12 12
Flags with carried forward coverage won't be shown. Click here to find out more.
|
c808f6d
to
176651a
Compare
176651a
to
cee8338
Compare
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.
@mbarnfield63 It would be great if you could go through the changes in this PR, and make a comment where things are not clear to you. Focus on file force_h2o.F90
. The actual potential is in h2o_schwenke.f
, but that has been taken essentially verbatim from the paper (and we already went through it together once).
Changes to the other files are not that important, they are mostly just hooking up the new potential trough ABIN. (and there are a couple of unrelated changes as well).
I am going to merge this to the master branch, which you can then pull-in and start your own branch on your fork for the numerical forces.
Note that I have renamed the module and function names a bit to make them consistent, so your scan script would need to be updated accrodingly.
|
||
! TODO: Implement numerical forces generally for all potentials | ||
! For now, they can be implemented here and hardcoded for a specific H2O potential | ||
subroutine numerical_forces(x, y, z, Epot, fx, fy, fz, natom, nbeads) |
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.
@mbarnfield63 Here's where you can start with your implementation.
Ultimately we'll want to have the implementation of numerical forces generalized and in a separate module, but that's not important at this stage.
No description provided.