Skip to content

Commit

Permalink
New experimental autoimage mode (#1119)
Browse files Browse the repository at this point in the history
* Start preparing for new autoimage method

* Add mode variable

* Try to implement algorithm which looks at the original orientation of a
fixed molecule w.r.t. the anchor molecule and tries to move it in as
close a configuration as possible to the original orientation

* Actually image the fixed unit if needed

* Convert back to Cartesian

* Wrap mobile molecules

* Hide some debug info

* Add converttofrac command for converting coordinates to fractional space

* For testing a new imaged distance calc

* Action to test speed of distance calculations

* Test out a new imaging routine that tries to only look for images that
can potentially be close

* Create a version of the distance calc that takes wrapped fractional
coords

* Do the 'familiar' truncated octahedron wrapping.

* Add some code doc

* Hide debug

* Hide converttofrac

* Note that converttofrac is only for debug

* Note that 'moveanchor' can still be active for mode byvec for the first
frame

* Add a description for 'mode byvec' in help.

* Clean up help text

* Fix up help text

* Add missing manual entry for moveanchor keyword of autoiamge. Add manual
entry for new mode keyword of autoimage.

* 6.29.8. Revision bump for addition of 'mode byvec' keyword to autoimage.

* Add autoimage mode byvec test

* Add warning that mode byvec still needs more testing

* Revision bump for experimental autoimage mode byvec

* Revert to 6.29.8. Forgot I already bumped the revision

* 6.29.9. 6.29.8 was PR #1116
  • Loading branch information
drroe authored Dec 6, 2024
1 parent 20512f3 commit 77f5392
Show file tree
Hide file tree
Showing 17 changed files with 87,590 additions and 18 deletions.
63 changes: 58 additions & 5 deletions doc/cpptraj.lyx
Original file line number Diff line number Diff line change
Expand Up @@ -24348,15 +24348,15 @@ autoimage
\end_layout

\begin_layout LyX-Code
autoimage [<mask> | anchor <mask> [fixed <mask>] [mobile <mask>]]
autoimage [{<mask> | anchor <mask> [fixed <mask>] [mobile <mask>]}]
\end_layout

\begin_layout LyX-Code
[origin] [firstatom] [familiar | triclinic]
\begin_inset Separator latexpar
\end_inset

[origin] [firstatom] [{familiar|triclinic}] [moveanchor]
\end_layout

\begin_layout LyX-Code
[mode {bydist|byvec}]
\end_layout

\begin_deeper
Expand Down Expand Up @@ -24416,6 +24416,45 @@ anchor
[triclinic] Force general triclinic imaging.
\end_layout

\begin_layout Description
[moveanchor] When treating
\begin_inset Quotes eld
\end_inset

fixed
\begin_inset Quotes erd
\end_inset

molecules,
the anchor point will be set to the previous "fixed" molecule;
this is only expected to work well when "fixed" molecules that are sequential are also geometrically close.
Most useful in more condensed systems like those containing membranes.
\end_layout

\begin_layout Description
[mode
\begin_inset space ~
\end_inset

{bydist|byvec}] How to treat 'fixed' molecules.
\end_layout

\begin_deeper
\begin_layout Description
bydist The default.
'Fixed' molecules will use the image closest to the "anchor" molecule.
\end_layout

\begin_layout Description
byvec Fixed molecules will use the image closest to their orientation with respect to the anchor in the first frame.
May work better than
\series bold
'bydist'
\series default
for systems containing membranes.
\end_layout

\end_deeper
\end_deeper
\begin_layout Standard
Automatically center and image (by molecule) a trajectory with periodic boundaries.
Expand Down Expand Up @@ -24459,6 +24498,20 @@ The autoimage command works for the majority of systems;
it is recommended to choose as the anchor a small region which should lie near the center of your system.
For example,
in a protein dimer system one could choose a single residue that is near the center of the interface between the two monomers.
The
\series bold
'moveanchor'
\series default
and
\series bold
'mode
\series default

\series bold
byvec'
\series default
options may also help in cases where the system is more condensed,
such as those containing membranes.
\end_layout

\begin_layout Subsection
Expand Down
Loading

0 comments on commit 77f5392

Please sign in to comment.