-
Notifications
You must be signed in to change notification settings - Fork 24
Add guide for using GauXC from C and Fortran #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
base: main
Are you sure you want to change the base?
Conversation
601b733 to
0d9fda0
Compare
29716cd to
6a7fa80
Compare
1e71806 to
b8a9725
Compare
|
Some questions on the GauXC in C instructions
|
7909698 to
e38b615
Compare
Good point, now this all handled in GauXC since wavefunction91/GauXC#169 was merged. |
Since GauXC has MPI support I would like to highlight it as well. Also, it does impact the API for the runtime environment so it is something that can't be easily ignored. We are testing all examples with and without MPI to make sure they actually work correctly. |
I follow the naming used in GauXC: |
| ! Integrate exchange-correlation energy | ||
| vxc_s = gauxc_matrix_empty(status) | ||
| vxc_z = gauxc_matrix_empty(status) | ||
| call gauxc_eval_exc_vxc(status, integrator, p_s, p_z, model, exc, vxc_s, vxc_z) |
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.
my compiler complains it cannot find a matching argument pattern. Removing model fixed the problem in my case.
| call gauxc_eval_exc_vxc(status, integrator, p_s, p_z, model, exc, vxc_s, vxc_z) | |
| call gauxc_eval_exc_vxc(status, integrator, p_s, p_z, exc, vxc_s, vxc_z) |
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.
Model should be passed here to dispatch to Skala interface, if you can't find a missing interface you might not have Skala enabled in GauXC.
| integer :: error | ||
|
|
||
| #ifdef GAUXC_HAS_MPI | ||
| call MPI_Init(error) |
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 mpi communicator error is now the same variable as the error in cli. I am not sure wether this in causing conflicts, but giving separate variables didn't solve the mpi issue in my case.
See wavefunction91/GauXC#172
See wavefunction91/GauXC#174
Preview at https://awvwgk.github.io/skala/gauxc/c-library.html and https://awvwgk.github.io/skala/gauxc/fortran-library.html