Skip to content
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

7.2.0 fails gradients of bare ecps which <=7.0.2 works #801

Closed
chemsh opened this issue Jun 19, 2023 · 2 comments
Closed

7.2.0 fails gradients of bare ecps which <=7.0.2 works #801

chemsh opened this issue Jun 19, 2023 · 2 comments
Labels

Comments

@chemsh
Copy link

chemsh commented Jun 19, 2023

Dear NWChem Team,

we find a fundamental bug when testing NWChem 7.2.0 for ChemShell. Please see the attached example input and script used for compiling NWChem (both renamed as .txt). This is a test case of gradients on the background point charges which has been stable for many years. With NWChem 7.2.0 it fails with the message, either in serial or parallel:

0:cannot locate region: density matrix [0:-1 ,1:15 ]:Received an Error in Communication
Abort(-999) on node 0 (rank 0 in comm 496): application called MPI_Abort(comm=0x84000002, -999) - process 0

I've traced back to function hnd_elfcon_0 around line 262 in property/hnd_elfcon.F:
call ga_get(g_dens, bf_lo_i, bf_hi_i, bf_lo_j, where bf_lo_i and bf_hi_i became 0 and -1, respectively. I think it's a bug in the code because these lines are new in 7.2.0. @dmejiar

As metal oxide is one of the most common types of simulation in ChemShell, we appreciate if you could take a look at this issue. Thank you very much.

Kind regards,

You

compile_nwchem_on_ubuntu.txt
mgo_nwchem_bug.txt

@chemsh chemsh changed the title 7.2.0 fails calculating electric field contribution for which it worked in 7.0.2 7.2.0 fails calculating electric field contribution for which 7.0.2 works Jun 19, 2023
@chemsh chemsh changed the title 7.2.0 fails calculating electric field contribution for which 7.0.2 works 7.2.0 fails calculating electric field contribution for which <=7.0.2 works Jun 19, 2023
@edoapra edoapra changed the title 7.2.0 fails calculating electric field contribution for which <=7.0.2 works 7.2.0 fails gradients of bare ecps which <=7.0.2 works Jun 20, 2023
edoapra added a commit to edoapra/nwchem that referenced this issue Jun 20, 2023
@edoapra edoapra added the bug label Jun 20, 2023
edoapra added a commit that referenced this issue Jun 20, 2023
@edoapra
Copy link
Collaborator

edoapra commented Jun 21, 2023

A fix now now available both in the master and hotfix/release-7-2-0 branch
5a45dbf

@chemsh
Copy link
Author

chemsh commented Jun 23, 2023

Dear Edo,
@edoapra

thank you so much for the prompt response and hotfix! I can confirm hotfix/release-7-2-0 works for our tests again. I am closing the current Issue now.

Kind regards,

You

@chemsh chemsh closed this as completed Jun 23, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

2 participants