Amber-MD / cpptraj

Biomolecular simulation trajectory/data analysis.
Other
139 stars 64 forks source link

Add TICA, max/min distance calc, and data projection commands #1075

Closed drroe closed 7 months ago

drroe commented 7 months ago

Version 6.24.1.

Adds new commands to cpptraj: one for performing time-independent correlation analysis (TICA), one for calculating the minimum distance between atoms or pairs of residues/molecules, one for calculating the maximum distance, and one for projecting 1D data sets along eigenmodes.

New analysis command: tica

        { crdset <COORDS set name> [mask <mask>] |
          data <input set arg1> ... }
        [lag <time lag>] [map {kinetic|commute|none}]
        [name <output set name>] [out <file>] [cumvarout <file>]

New actions: mindist / maxdist

        mask1 <mask1> [mask2 <mask2>] [{byatom|byres|bymol}]
        [noimage] [name <setname>] [out <file>] [resoffset <#>]

New analysis command: projectdata

        evecs <evecs dataset> [name <name>] [out <outfile>] [beg <beg>] [end <end>]
        {[dihedrals <dataset arg>] | [data <dataset arg> ...]}

This is a useful (often faster) alternative to project for 1D data sets when you don't have to project coordinates. The projection command has also been updated to allow projection of regular 1D data sets (on e.g. TICA modes). The matrix command will be updated to allow covariance calculation on 1D data sets in another PR.

Adds tests and updates the manual.