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

Issue with Handling Topology Files with Three Columns in ligand_alchemy.py #50

Open
h28363231 opened this issue Nov 20, 2024 · 1 comment

Comments

@h28363231
Copy link

Dear Developer Team,

I am currently working with the develop branch of pmx and encountered an issue in the file pmx/src/pmx/ligand_alchemy.py, specifically within the method _get_ff_entry. The issue arises when handling topology files that contain only three columns, such as the example below:

[ bonds ]
; ai aj funct b0 Kb
1 2 1
1 26 1
1 27 1
1 28 1

In this method, there is a line of code: out = cp.deepcopy(b), and later, the code attempts to access out[3:][0]. This results in a problem because there are only three columns in the topology file, hence, accessing out[3:][0] does not exist and causes an error.

Could you please advise on how to modify the script to handle this scenario appropriately or consider updating the code to handle such cases in the develop branch?

Thank you for your assistance and looking forward to your guidance.

Best regards,

Jack Liu

@h28363231 h28363231 reopened this Nov 20, 2024
@vgapsys
Copy link
Member

vgapsys commented Nov 20, 2024

This function is used when constructing hybrid topologies. For that it is necessary to have parameters explicitly defined in the input topology files. In this case, bond length and force constant are missing

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants