Added analysis functionalities for landmark selection and dissimilarity matrix calculation

680d30d6 · Gareth Tribello · f8047308 · 680d30d6 · 680d30d6 · 680d30d6
Commit 680d30d6 authored 9 years ago by Gareth Tribello
--- a/regtest/analysis/rt-calc-dissims/Makefile
+++ b/regtest/analysis/rt-calc-dissims/Makefile
+include ../../scripts/test.make
--- a/regtest/analysis/rt-calc-dissims/analysis.0.mymatrix3.dat.reference
+++ b/regtest/analysis/rt-calc-dissims/analysis.0.mymatrix3.dat.reference
+   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000
+   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000
+   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000
+   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000
+  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000
+  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000
--- a/regtest/analysis/rt-calc-dissims/colv_in
+++ b/regtest/analysis/rt-calc-dissims/colv_in
+#! FIELDS time data
+0 0
+1 0
+2 1
+3 2
+4 3
+5 4
+6 5 
+7 6
+8 7 
+9 8 
+10 9
+11 10
+12 11
--- a/regtest/analysis/rt-calc-dissims/config
+++ b/regtest/analysis/rt-calc-dissims/config
+type=driver
+arg="--noatoms --plumed plumed.dat"
+
--- a/regtest/analysis/rt-calc-dissims/mymatrix.dat.reference
+++ b/regtest/analysis/rt-calc-dissims/mymatrix.dat.reference
+   0.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000  81.0000 100.0000 121.0000
+   0.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000  81.0000 100.0000 121.0000
+   1.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000  81.0000 100.0000
+   4.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000  81.0000
+   9.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000
+  16.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000
+  25.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000
+  36.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000
+  49.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000
+  64.0000  64.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000
+  81.0000  81.0000  64.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000
+ 100.0000 100.0000  81.0000  64.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000
+ 121.0000 121.0000 100.0000  81.0000  64.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000
--- a/regtest/analysis/rt-calc-dissims/mymatrix2.dat.reference
+++ b/regtest/analysis/rt-calc-dissims/mymatrix2.dat.reference
+   0.0000   1.0000   9.0000  25.0000  49.0000
+   1.0000   0.0000   4.0000  16.0000  36.0000
+   9.0000   4.0000   0.0000   4.0000  16.0000
+  25.0000  16.0000   4.0000   0.0000   4.0000
+  49.0000  36.0000  16.0000   4.0000   0.0000
--- a/regtest/analysis/rt-calc-dissims/mymatrix3.dat.reference
+++ b/regtest/analysis/rt-calc-dissims/mymatrix3.dat.reference
+   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000
+   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000
+   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000
+   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000
+  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000
+  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000
--- a/regtest/analysis/rt-calc-dissims/output-fps.pdb.reference
+++ b/regtest/analysis/rt-calc-dissims/output-fps.pdb.reference
+DESCRIPTION: analysis data from calculation done at time 12.000000 
+REMARK WEIGHT=6.000000
+REMARK ARG=d1
+REMARK d1=7.000000 
+END
+DESCRIPTION: analysis data from calculation done at time 12.000000 
+REMARK WEIGHT=5.000000
+REMARK ARG=d1
+REMARK d1=0.000000 
+END
+DESCRIPTION: analysis data from calculation done at time 12.000000 
+REMARK WEIGHT=2.000000
+REMARK ARG=d1
+REMARK d1=11.000000 
+END
--- a/regtest/analysis/rt-calc-dissims/output-stride.pdb.reference
+++ b/regtest/analysis/rt-calc-dissims/output-stride.pdb.reference
+DESCRIPTION: analysis data from calculation done at time 12.000000 
+REMARK WEIGHT=2.000000
+REMARK ARG=d1
+REMARK d1=0.000000 
+END
+DESCRIPTION: analysis data from calculation done at time 12.000000 
+REMARK WEIGHT=2.000000
+REMARK ARG=d1
+REMARK d1=1.000000 
+END
+DESCRIPTION: analysis data from calculation done at time 12.000000 
+REMARK WEIGHT=2.000000
+REMARK ARG=d1
+REMARK d1=3.000000 
+END
+DESCRIPTION: analysis data from calculation done at time 12.000000 
+REMARK WEIGHT=2.000000
+REMARK ARG=d1
+REMARK d1=5.000000 
+END
+DESCRIPTION: analysis data from calculation done at time 12.000000 
+REMARK WEIGHT=5.000000
+REMARK ARG=d1
+REMARK d1=7.000000 
+END
--- a/regtest/analysis/rt-calc-dissims/plumed.dat
+++ b/regtest/analysis/rt-calc-dissims/plumed.dat
+d1: READ FILE=colv_in VALUES=data
+
+ss1: EUCLIDEAN_DISSIMILARITIES STRIDE=1 USE_ALL_DATA ARG=d1 
+PRINT_DISSIMILARITY_MATRIX USE_OUTPUT_DATA_FROM=ss1 FILE=mymatrix.dat FMT=%8.4f
+
+ll1: LANDMARK_SELECT_STRIDE ARG=d1 USE_OUTPUT_DATA_FROM=ss1 NLANDMARKS=5 
+ss2: EUCLIDEAN_DISSIMILARITIES USE_OUTPUT_DATA_FROM=ll1 
+OUTPUT_ANALYSIS_DATA_TO_PDB USE_OUTPUT_DATA_FROM=ll1 FILE=output-stride.pdb
+PRINT_DISSIMILARITY_MATRIX USE_OUTPUT_DATA_FROM=ll1 FILE=mymatrix2.dat FMT=%8.4f
+
+ll2: LANDMARK_SELECT_FPS USE_OUTPUT_DATA_FROM=ss1 NLANDMARKS=3
+OUTPUT_ANALYSIS_DATA_TO_PDB USE_OUTPUT_DATA_FROM=ll2 FILE=output-fps.pdb
+
+ss3: EUCLIDEAN_DISSIMILARITIES STRIDE=1 RUN=6 ARG=d1
+oo: PRINT_DISSIMILARITY_MATRIX USE_OUTPUT_DATA_FROM=ss3 FILE=mymatrix3.dat FMT=%8.4f
--- a/regtest/analysis/rt-read-dissims/Makefile
+++ b/regtest/analysis/rt-read-dissims/Makefile
+include ../../scripts/test.make
--- a/regtest/analysis/rt-read-dissims/config
+++ b/regtest/analysis/rt-read-dissims/config
+type=driver
+arg="--noatoms --plumed plumed.dat"
+
--- a/regtest/analysis/rt-read-dissims/mymatrix.dat
+++ b/regtest/analysis/rt-read-dissims/mymatrix.dat
+   0.0000   1.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000  81.0000 100.0000 121.0000
+   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000  81.0000 100.0000 121.0000
+   1.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000  81.0000 100.0000
+   4.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000  81.0000
+   9.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000
+  16.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000
+  25.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000
+  36.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000
+  49.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000
+  64.0000  64.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000
+  81.0000  81.0000  64.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000
+ 100.0000 100.0000  81.0000  64.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000
+ 121.0000 121.0000 100.0000  81.0000  64.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000
--- a/regtest/analysis/rt-read-dissims/mymatrix_out.dat.reference
+++ b/regtest/analysis/rt-read-dissims/mymatrix_out.dat.reference
+   0.0000   1.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000  81.0000 100.0000 121.0000
+   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000  81.0000 100.0000 121.0000
+   1.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000  81.0000 100.0000
+   4.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000  81.0000
+   9.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000  64.0000
+  16.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000  49.0000
+  25.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000  36.0000
+  36.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000  25.0000
+  49.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000  16.0000
+  64.0000  64.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000   9.0000
+  81.0000  81.0000  64.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000   4.0000
+ 100.0000 100.0000  81.0000  64.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000   1.0000
+ 121.0000 121.0000 100.0000  81.0000  64.0000  49.0000  36.0000  25.0000  16.0000   9.0000   4.0000   1.0000   0.0000
--- a/regtest/analysis/rt-read-dissims/mymatrix_stride.dat.reference
+++ b/regtest/analysis/rt-read-dissims/mymatrix_stride.dat.reference
+   0.0000   1.0000   9.0000  25.0000  49.0000
+   1.0000   0.0000   4.0000  16.0000  36.0000
+   9.0000   4.0000   0.0000   4.0000  16.0000
+  25.0000  16.0000   4.0000   0.0000   4.0000
+  49.0000  36.0000  16.0000   4.0000   0.0000
--- a/regtest/analysis/rt-read-dissims/plumed.dat
+++ b/regtest/analysis/rt-read-dissims/plumed.dat
+r1: READ_DISSIMILARITY_MATRIX FILE=mymatrix.dat 
+PRINT_DISSIMILARITY_MATRIX USE_OUTPUT_DATA_FROM=r1 FILE=mymatrix_out.dat FMT=%8.4f
+
+l1: LANDMARK_SELECT_STRIDE USE_OUTPUT_DATA_FROM=r1 NLANDMARKS=5
+PRINT_DISSIMILARITY_MATRIX USE_OUTPUT_DATA_FROM=l1 FILE=mymatrix_stride.dat FMT=%8.4f
--- a/regtest/analysis/rt-store-print/Makefile
+++ b/regtest/analysis/rt-store-print/Makefile
+include ../../scripts/test.make
--- a/regtest/analysis/rt-store-print/analysis.0.output-half.pdb.reference
+++ b/regtest/analysis/rt-store-print/analysis.0.output-half.pdb.reference
+#! FIELDS d1
+0.000000 
+1.000000 
+2.000000 
+3.000000 
+4.000000 
+5.000000 
--- a/regtest/analysis/rt-store-print/analysis.0.output-stride2.pdb.reference
+++ b/regtest/analysis/rt-store-print/analysis.0.output-stride2.pdb.reference
+#! FIELDS d1
+0.000000 
+2.000000 
+4.000000 
--- a/regtest/analysis/rt-store-print/colv_in
+++ b/regtest/analysis/rt-store-print/colv_in
+#! FIELDS time data
+0 0
+1 0
+2 1
+3 2
+4 3
+5 4
+6 5 
+7 6
+8 7 
+9 8 
+10 9
+11 10
+12 11