Goodman-lab / DP5

Python workflow for DP5 and DP4 analysis of organic molecules
Other
173 stars 99 forks source link

NWChem.py line 242 typo - Populating NMROutputFiles instead of OptOutputFiles. #40

Closed danny305 closed 4 years ago

danny305 commented 4 years ago

In the RunOptCalc() function in NWChem.py I think there is a typo at line 242:

The for loop should populate isoOptOutputFiles instead of isoNMROutputFiles?

So change:

for iso in Isomers:
        iso.NMROutputFiles.extend([x[:-3] + '.nwo' for x in iso.OptInputFiles if (x[:-3] + '.nwo') in Completed])

To

for iso in Isomers:
        iso.OptOutputFiles.extend([x[:-3] + '.nwo' for x in iso.OptInputFiles if (x[:-3] + '.nwo') in Completed])

After running the RunOptCalc() and trying to ReadGeometries(), I am getting the following error: NWChem.py, ReadGeometries: No geometry optimisation output files found, geometries could not be read. Quitting.

I believe this is because of the typo ^^.

danny305 commented 4 years ago

Also, I believe that line 236 meant to print the OptInputFiles and not the NMRInputFiles. This is just an observation from comparing RunOptCalc() to RunECalcs() and RunNMRCalcs().

zjp commented 4 years ago

Are you trying to use precalculated data? If so, can you try my unmerged patch in #38 and see if it works for you? Possibly the second part is more relevant.

danny305 commented 4 years ago

No I was trying to generate data. By making the changes I suggested above, NWChem optimization works.

KristapsE commented 4 years ago

Yup, that definitely looks like a typo on our part, I'll try to fix it tonight and merge it. Good catch!

KristapsE commented 4 years ago

So I have made the changes you suggested - they totally make sense to me, but could you please test the current version and let me know if it works as expected for you?

danny305 commented 4 years ago

Works on my end. Thanks.