-
Notifications
You must be signed in to change notification settings - Fork 23
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
Molecules in high pressure environments #467
base: master
Are you sure you want to change the base?
Conversation
hugoastrom
commented
Oct 11, 2023
•
edited
Loading
edited
- Confinement potential class
- Incorporate cavity function
- Test for atoms
- there seems to be problems with convergence
|
||
double ConfinementFunction::evalf(const mrcpp::Coord<3> &r) const { | ||
double f = 0.0; | ||
const mrcpp::Coord<3> center = {0.0, 0.0, 0.0}; |
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.
Is the center hard-wired to the origin? Maybe it would be wise to be able to pass this from input.
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.
The plan is to obtain the centers as it is implemented in the Cavity function but it is under construction
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.
I left a few comments for you to consider.
@@ -268,6 +269,7 @@ json driver::scf::run(const json &json_scf, Molecule &mol) { | |||
auto relativity = json_scf["scf_solver"]["relativity"]; | |||
auto environment = json_scf["scf_solver"]["environment"]; | |||
auto external_field = json_scf["scf_solver"]["external_field"]; | |||
auto confinement_potential = json_scf["scf_solver"]["confinement_potential"]; |
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.
Seems you have some indentation issues which should be fixed
src/driver.cpp
Outdated
/////////////// Confinement Potential ///////////////// | ||
/////////////////////////////////////////////////////////// | ||
if (json_fock.contains("confinement_operator")) { | ||
auto r_0 = json_fock["confinement_operator"]["r_0"]; |
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.
indentation should be 4 spaces :-)
About the indentation, I think you can set git so that it can automatically format your files before committing changes |