PyAutoFEP icon indicating copy to clipboard operation
PyAutoFEP copied to clipboard

Analyzing incomplete runs

Open xy21hb opened this issue 1 year ago • 4 comments

If some pairs are broken, could analyze_results.py still work, suppose e.g. if only FXR_12_FXR_88 pair of the tutorial is broken, e.g. system is unstable, only intermediate pdbs are written and gromacs fails to collect energy files. I wonder how to obtain the results from other pairs, considering the generated pertubation_map.

xy21hb avatar Sep 11 '23 03:09 xy21hb

Error msg is below, some steps of structures are not stable, and GMX fails to collect edr energy files

image

image

xy21hb avatar Sep 11 '23 03:09 xy21hb

I'm theory, running with --no_checks should allow for such.

However, your error message suggests that there are some edr files in the rerun dir, which could cause errors anyway. Would you mind retrying with --no_checks, and, in case it doesn't work, post either the .tgz or the directory structure (should your data be sensitive)?

Thank you very much.

luancarvalhomartins avatar Sep 12 '23 08:09 luancarvalhomartins

Though I have not checked the structure, I believe it is due to some unstable ligs since the system is a large membrane protein. Running with --no_checks would give the same error msg. So, to analyze incomplete runs, simply excluding those runs in pack.sh would work.

xy21hb avatar Sep 13 '23 01:09 xy21hb

Ok, so --no_checks wasn't enough. Let me take a look at the code to fix this.

Thank you again for reporting.

luancarvalhomartins avatar Sep 13 '23 08:09 luancarvalhomartins