Skip to content

Initial addition of angle constraints - #8

Open
rsdefever wants to merge 5 commits into
mainfrom
add/angle-constraints
Open

Initial addition of angle constraints#8
rsdefever wants to merge 5 commits into
mainfrom
add/angle-constraints

Conversation

@rsdefever

Copy link
Copy Markdown
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
Copy Markdown
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
Copy Markdown
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
Copy Markdown
Owner Author

/azp run

@azure-pipelines

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

@rsdefever

rsdefever commented Nov 18, 2020

Copy link
Copy Markdown
Owner Author

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