File: RunTest.sh

package info (click to toggle)
cpptraj 5.1.0%2Bdfsg-5
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 139,188 kB
  • sloc: cpp: 145,622; ansic: 34,635; sh: 11,365; f90: 971; makefile: 770; awk: 242
file content (94 lines) | stat: -rwxr-xr-x 2,485 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
#!/bin/bash

. ../MasterTest.sh

# Clean
CleanFiles dummy.rst7 dummy.pdb.1 strip.in *.tz2.truncoct.parm7 Complex.crd \
           Receptor.crd Ligand.crd res1.tz2.crd strip.pdb.chignolin.parm7
INPUT="-i strip.in"

# NOTE: strip is also tested in Test_Center
UNITNAME='One frame strip command test'
CheckFor maxthreads 1
if [ $? -eq 0 ] ; then
  cat > strip.in <<EOF
noprogress
parm ../tz2.truncoct.parm7
trajin ../tz2.truncoct.crd 1 1 parm ../tz2.truncoct.parm7
strip :14-16,18-99999 outprefix strip
trajout dummy.pdb pdb multi parm ../tz2.truncoct.parm7 chainid X nobox
trajout dummy.rst7 restart parm ../tz2.truncoct.parm7
EOF
  RunCpptraj "$UNITNAME"
  DoTest dummy.pdb.save dummy.pdb.1
  DoTest dummy.rst7.save dummy.rst7
  # Tell diff to ignore the VERSION line
  DoTest strip.tz2.truncoct.parm7.save strip.tz2.truncoct.parm7 -I %VERSION
fi

UNITNAME='Unstrip (Lig/Rec/Complex) command test'
CheckFor maxthreads 1
if [ $? -eq 0 ] ; then
  # Unstrip Test
  cat > strip.in <<EOF
noprogress
parm ../tz2.truncoct.parm7
trajin ../tz2.truncoct.crd 1 1 parm ../tz2.truncoct.parm7

# Count res 17 as ligand, 1-13 as receptor, 1-13,17 as complex

# Complex
strip !(:1-13,17) outprefix complex
outtraj Complex.crd
unstrip

# Receptor
strip :WAT outprefix receptor
outtraj Receptor.crd
unstrip

# Ligand
strip !(:17) outprefix ligand parmopts nopdbinfo
outtraj Ligand.crd
EOF
  RunCpptraj "$UNITNAME"
  # Tell diff to ignore the VERSION line
  DoTest strip.tz2.truncoct.parm7.save complex.tz2.truncoct.parm7 -I %VERSION
  DoTest receptor.tz2.truncoct.parm7.save receptor.tz2.truncoct.parm7 -I %VERSION
  DoTest ligand.tz2.truncoct.parm7.save ligand.tz2.truncoct.parm7 -I %VERSION
  DoTest Ligand.crd.save Ligand.crd
  DoTest Receptor.crd.save Receptor.crd
  DoTest Complex.crd.save Complex.crd
fi

# Strip test that will work in parallel
UNITNAME='Multiple frame strip command test'
CheckFor netcdf
if [ $? -eq 0 ] ; then
  cat > strip.in <<EOF
noprogress
parm ../tz2.parm7
trajin ../tz2.nc
strip !(:1) nobox
trajout res1.tz2.crd
EOF
  RunCpptraj "Multi frame strip command test."
  DoTest res1.tz2.crd.save res1.tz2.crd
fi

# Check that stripping topology with PDB info works
UNITNAME='Strip topology with PDB info'
CheckFor maxthreads 1
if [ $? -eq 0 ] ; then
  cat > strip.in <<EOF
parm pdb.chignolin.parm7
parmstrip @/H
parmwrite out strip.pdb.chignolin.parm7
EOF
  RunCpptraj "$UNITNAME"
  DoTest strip.pdb.chignolin.parm7.save strip.pdb.chignolin.parm7 -I %VERSION
fi

EndTest

exit 0