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

true/false properties in .mae file not working #41

Open
arosale4 opened this issue Mar 14, 2017 · 4 comments
Open

true/false properties in .mae file not working #41

arosale4 opened this issue Mar 14, 2017 · 4 comments
Labels

Comments

@arosale4
Copy link
Contributor

Whether I use 0 or 1 for b_cs_first_match_only, or b_cs_both_enantiomers I can not seem to reduce the number of structures produced. I have even gone into merge.py and set first_match_only to True, and it would still match every instance of a pattern. I'm trying every version of schrodinger that we have available to see if that is a problem. I have also noticed that the pattern PD-C2=C2 will also match PD.C2=C2.

@ericchansen
Copy link
Owner

Any updates on this?

@ericchansen
Copy link
Owner

@arosale4 Do your recent changes solve this issue?

@arosale4
Copy link
Contributor Author

Thanks for drawing my attention to this. The recent changes may help with it, but indirectly. I do think we fixed this in 68da398?

@arosale4
Copy link
Contributor Author

Or at least we addressed the issue of reducing matches from first_match_only() to one.

ericchansen added a commit that referenced this issue Oct 31, 2018
See Q2MM#40 (I hope that the magic GitHub Markdown turns that into a link).
Also, I got rid of some blank spaces.

Anyway, it turns out that Schrodinger's [getBond()](http://content.schrodinger.com/Docs/r2015-4/python_api/api/schrodinger.structure.Structure-class.html#getBond)
function can return None. Personally, I would have made that raise an exception
of some sort, but I mean, I didn't write it. Anyway, here's an example.

Ex.) Don't take this syntax too literally. It's just an illustration.

N-P (single bond)
bond = merge.getBond(N, P)
print(bond)
>>> bond(1, 2)

N P (no bond)
bond = merge.getBond(N, P)
print(bond)
>>> None

The code now checks for this behavior and adds the bond if it's missing.
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