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

Initial addition of angle constraints #8

Open
wants to merge 5 commits into
base: main
Choose a base branch
from

Conversation

rsdefever
Copy link
Owner

@adowling2 can you have a look at my attempted implementation of angle constraints. I'm getting the following error, even for a very simple test case. Any ideas?

EXIT: Converged to a point of local infeasibility. Problem may be infeasible.                                                                                                                                                                                                                             

@adowling2
Copy link
Collaborator

@rsdefever This is on my todo list. I might not get to it for a week.

@rsdefever rsdefever force-pushed the add/angle-constraints branch from e3109c0 to 5cb41a8 Compare November 18, 2020 20:22
@rsdefever rsdefever marked this pull request as ready for review November 18, 2020 20:23
@rsdefever rsdefever changed the title Initial addition of angle constraints. Tests not passing. Initial addition of angle constraints Nov 18, 2020
@rsdefever
Copy link
Owner Author

The angle constraints now seem to be working in a simple case (an important simple case!) of a water molecule (so two bond constraints with one angle constraint). A more difficult test case of a dimethyl ether molecule still reports the same error on the solver. When I check the geometry, all the bonds meet the bond-length constraints. However, there are several angles that do not meet the angle constraints. I am unsure if it is feasible for all the angles to be satisfied, although the total number of constraints is 21 and the number of variables is 27.

Here is the molecule:

vmdscene

@rsdefever
Copy link
Owner Author

/azp run

@azure-pipelines
Copy link

Azure Pipelines could not run because the pipeline triggers exclude this branch/path.

@rsdefever
Copy link
Owner Author

rsdefever commented Nov 18, 2020

Also of note is that the molecule changes configuration a fair bit more than in the case of the bond constraints alone. Blue = start; color based on atom type = end; cyan = C, red = O, white = H.

dme-start-end

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

Successfully merging this pull request may close these issues.

2 participants