biod.pnpi.spb.ru / alexxy / gromacs.git / commitdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | patch | inline | side by side (parent: 56c92cb)
Unit tests for insolidangle selections.
authorTeemu Murtola <teemu.murtola@gmail.com>
Fri, 14 Sep 2012 17:24:48 +0000 (20:24 +0300)
committerGerrit Code Review <gerrit@gerrit.gromacs.org>
Thu, 4 Oct 2012 12:20:37 +0000 (14:20 +0200)
The current input data is not particularly elaborate, and more careful
selection could possibly test more corner cases.  But coverage even with
the current 500+500 random test points should be relatively good.

In addition to the tests themselves and input data for them,
added a Python script that was used to generate the test data.

Part of #651.

Change-Id: I80acd431b2963bf65fbeed53fe6be4f885d07ec5

src/gromacs/selection/tests/gensphere.py [new file with mode: 0755] patch | blob
src/gromacs/selection/tests/refdata/SelectionCollectionDataTest_HandlesInSolidAngleKeyword.xml [new file with mode: 0644] patch | blob
src/gromacs/selection/tests/selectioncollection.cpp patch | blob | history
src/gromacs/selection/tests/sphere.gro [new file with mode: 0644] patch | blob

diff --git a/src/gromacs/selection/tests/gensphere.py b/src/gromacs/selection/tests/gensphere.py
new file mode 100755 (executable)
index 0000000..bea3df5
--- /dev/null
+++ b/src/gromacs/selection/tests/gensphere.py
@@ -0,0 +1,110 @@
+#!/usr/bin/python
+"""Script for generating spherical test configurations."""
+
+import math
+import random
+import sys
+
+def dot(a, b):
+    return sum(x*y for (x, y) in zip(a, b))
+
+def norm(a):
+    return math.sqrt(dot(a, a))
+
+def angle(a, b):
+    return math.degrees(math.acos(dot(a, b) / (norm(a) * norm(b))))
+
+def minangle(a, list):
+    minangle = 180
+    minindex = -1
+    for index, x in enumerate(list):
+        xangle = angle(a, x)
+        if xangle < minangle:
+            minangle = xangle
+            minindex = index
+    return (minangle, minindex)
+
+def get_single_vec():
+    while True:
+        x = random.randint(-1000, 1000) / 1000.0
+        y = random.randint(-1000, 1000) / 1000.0
+        z = random.randint(-1000, 1000) / 1000.0
+        pos = (x, y, z)
+        dist = norm(pos)
+        if dist <= 1.0 and dist > 0.25:
+            return pos
+
+def write_gro_title(fp, title, atomcount):
+    fp.write(title + '\n')
+    fp.write('%5d\n' % (atomcount))
+
+def write_gro_atom(fp, resnr, resname, atomname, index, x):
+    fp.write('%5d%-5s%5s%5d%8.3f%8.3f%8.3f\n' %
+            (resnr, resname, atomname, index, x[0], x[1], x[2]))
+
+def write_gro_box(fp, box):
+    fp.write('%10.5f%10.5f%10.5f\n' % box)
+
+random.seed(1097)
+center = (0, 0, 0)
+cutoff = 20
+possamples = 500
+negsamples = 500
+totsamples = 10000
+
+sys.stderr.write("Generating reference points\n")
+refpoints = []
+refpoints.append((0, 0, -1))
+refpoints.append((-0.5, 0.6, 0.1))
+refpoints.append((-0.5, -0.5, 0.25))
+while len(refpoints) < 30:
+    pos = get_single_vec()
+    if pos[0] > 0 and pos[1] > pos[0] and pos[2] > 0:
+        refpoints.append(pos)
+
+sys.stderr.write("Generating test points\n")
+postestpoints = []
+negtestpoints = []
+hits = 0
+samplecount = 0
+while samplecount < totsamples or len(postestpoints) < possamples or len(negtestpoints) < negsamples:
+    pos = get_single_vec()
+    (pangle, index) = minangle(pos, refpoints)
+    if pangle < cutoff:
+        hits += 1
+        if len(postestpoints) < possamples:
+            postestpoints.append(pos)
+    if pangle > cutoff:
+        if len(negtestpoints) < negsamples:
+            negtestpoints.append(pos)
+    samplecount += 1
+
+cfrac = float(hits) / samplecount
+errest = math.sqrt((cfrac - cfrac*cfrac) / samplecount)
+sys.stderr.write('Cutoff: %f angles\n' % (cutoff))
+sys.stderr.write('Estimated covered fraction: %f +- %f\n' % (cfrac, errest))
+
+debugfp = open('debug.txt', 'w')
+fp = sys.stdout
+count = 1 + len(refpoints) + len(postestpoints) + len(negtestpoints)
+write_gro_title(fp, 'Spherical test case, cutoff %f, cfrac %f +- %f' %
+        (cutoff, cfrac, errest) , count)
+n = 1
+write_gro_atom(fp, 1, 'C', 'C', n, center)
+n += 1
+for i in range(len(refpoints)):
+    write_gro_atom(fp, 2, 'R', 'R', n, refpoints[i])
+    n += 1
+for i in range(len(postestpoints)):
+    write_gro_atom(fp, 3, 'TP', 'TP', n, postestpoints[i])
+    x = postestpoints[i]
+    xangle, index = minangle(x, refpoints)
+    refx = refpoints[index]
+    debugfp.write('%3d%8.3f%8.3f%8.3f  %4.1f  %2d%8.3f%8.3f%8.3f\n' %
+            (n-1, x[0], x[1], x[2], xangle, index, refx[0], refx[1], refx[2]))
+    n += 1
+for i in range(len(negtestpoints)):
+    write_gro_atom(fp, 4, 'TN', 'TN', n, negtestpoints[i])
+    n += 1
+write_gro_box(fp, (10, 10, 10))
+debugfp.close()
diff --git a/src/gromacs/selection/tests/refdata/SelectionCollectionDataTest_HandlesInSolidAngleKeyword.xml b/src/gromacs/selection/tests/refdata/SelectionCollectionDataTest_HandlesInSolidAngleKeyword.xml
new file mode 100644 (file)
index 0000000..14469a4
--- /dev/null
+++ b/src/gromacs/selection/tests/refdata/SelectionCollectionDataTest_HandlesInSolidAngleKeyword.xml
@@ -0,0 +1,34 @@
+<?xml version="1.0"?>
+<?xml-stylesheet type="text/xsl" href="referencedata.xsl"?>
+<ReferenceData>
+  <ParsedSelections Name="Parsed">
+    <ParsedSelection Name="Selection1">
+      <String Name="Input">resname TP and not insolidangle center cog of resname C span resname R cutoff 20</String>
+      <String Name="Name">resname TP and not insolidangle center cog of resname C span resname R cutoff 20</String>
+      <String Name="Text">resname TP and not insolidangle center cog of resname C span resname R cutoff 20</String>
+      <Bool Name="Dynamic">true</Bool>
+    </ParsedSelection>
+    <ParsedSelection Name="Selection2">
+      <String Name="Input">resname TN and insolidangle center cog of resname C span resname R cutoff 20</String>
+      <String Name="Name">resname TN and insolidangle center cog of resname C span resname R cutoff 20</String>
+      <String Name="Text">resname TN and insolidangle center cog of resname C span resname R cutoff 20</String>
+      <Bool Name="Dynamic">true</Bool>
+    </ParsedSelection>
+  </ParsedSelections>
+  <CompiledSelections Name="Compiled">
+    <Selection Name="Selection1"/>
+    <Selection Name="Selection2"/>
+  </CompiledSelections>
+  <EvaluatedSelections Name="Frame1">
+    <Selection Name="Selection1">
+      <Sequence Name="Atoms">
+        <Int Name="Length">0</Int>
+      </Sequence>
+    </Selection>
+    <Selection Name="Selection2">
+      <Sequence Name="Atoms">
+        <Int Name="Length">0</Int>
+      </Sequence>
+    </Selection>
+  </EvaluatedSelections>
+</ReferenceData>
diff --git a/src/gromacs/selection/tests/selectioncollection.cpp b/src/gromacs/selection/tests/selectioncollection.cpp
index e51f4ddd3fdc8dcdf6fe8b20378898c00a1e62b8..b56f1a7e32b65e7e239def758f72a041a1d10ec3 100644 (file)
--- a/src/gromacs/selection/tests/selectioncollection.cpp
+++ b/src/gromacs/selection/tests/selectioncollection.cpp
@@ -159,7 +159,8 @@ class SelectionCollectionDataTest : public SelectionCollectionTest
             efTestEvaluation            = 1<<0,
             efTestPositionAtoms         = 1<<1,
             efTestPositionCoordinates   = 1<<2,
-            efTestPositionMapping       = 1<<3
+            efTestPositionMapping       = 1<<3,
+            efDontTestCompiledAtoms     = 1<<8
         };
         typedef gmx::FlagsTemplate<TestFlag> TestFlags;
 
@@ -292,7 +293,10 @@ SelectionCollectionDataTest::checkCompiled()
         std::string id = gmx::formatString("Selection%d", static_cast<int>(i + 1));
         TestReferenceChecker selcompound(
                 compound.checkCompound("Selection", id.c_str()));
-        checkSelection(&selcompound, sel_[i], flags_ & mask);
+        if (!flags_.test(efDontTestCompiledAtoms))
+        {
+            checkSelection(&selcompound, sel_[i], flags_ & mask);
+        }
     }
 }
 
@@ -690,7 +694,17 @@ TEST_F(SelectionCollectionDataTest, HandlesWithinKeyword)
 }
 
 
-// TODO: Add test for "insolidangle"
+TEST_F(SelectionCollectionDataTest, HandlesInSolidAngleKeyword)
+{
+    // Both of these should evaluate to empty on a correct implementation.
+    static const char * const selections[] = {
+        "resname TP and not insolidangle center cog of resname C span resname R cutoff 20",
+        "resname TN and insolidangle center cog of resname C span resname R cutoff 20",
+        NULL
+    };
+    setFlags(TestFlags() | efDontTestCompiledAtoms | efTestEvaluation);
+    runTest("sphere.gro", selections);
+}
 
 
 TEST_F(SelectionCollectionDataTest, HandlesPermuteModifier)
diff --git a/src/gromacs/selection/tests/sphere.gro b/src/gromacs/selection/tests/sphere.gro
new file mode 100644 (file)
index 0000000..5550191
--- /dev/null
+++ b/src/gromacs/selection/tests/sphere.gro
@@ -0,0 +1,1034 @@
+Spherical test case, cutoff 20.000000, cfrac 0.242800 +- 0.004288
+ 1031
+    1C        C    1   0.000   0.000   0.000
+    2R        R    2   0.000   0.000  -1.000
+    2R        R    3  -0.500   0.600   0.100
+    2R        R    4  -0.500  -0.500   0.250
+    2R        R    5   0.473   0.609   0.105
+    2R        R    6   0.385   0.849   0.162
+    2R        R    7   0.123   0.660   0.652
+    2R        R    8   0.194   0.467   0.571
+    2R        R    9   0.469   0.698   0.036
+    2R        R   10   0.544   0.653   0.164
+    2R        R   11   0.474   0.706   0.521
+    2R        R   12   0.285   0.303   0.580
+    2R        R   13   0.332   0.432   0.751
+    2R        R   14   0.101   0.372   0.589
+    2R        R   15   0.018   0.820   0.437
+    2R        R   16   0.134   0.679   0.403
+    2R        R   17   0.190   0.310   0.760
+    2R        R   18   0.155   0.577   0.011
+    2R        R   19   0.449   0.487   0.409
+    2R        R   20   0.133   0.963   0.112
+    2R        R   21   0.078   0.685   0.221
+    2R        R   22   0.154   0.221   0.178
+    2R        R   23   0.171   0.330   0.560
+    2R        R   24   0.503   0.753   0.351
+    2R        R   25   0.398   0.490   0.692
+    2R        R   26   0.090   0.095   0.453
+    2R        R   27   0.586   0.744   0.042
+    2R        R   28   0.449   0.644   0.606
+    2R        R   29   0.433   0.479   0.238
+    2R        R   30   0.179   0.860   0.067
+    2R        R   31   0.171   0.199   0.898
+    3TP      TP   32   0.255   0.382   0.306
+    3TP      TP   33   0.060  -0.012   0.969
+    3TP      TP   34   0.303   0.425  -0.078
+    3TP      TP   35  -0.319   0.761   0.117
+    3TP      TP   36   0.481   0.174   0.742
+    3TP      TP   37   0.106   0.488   0.608
+    3TP      TP   38  -0.283  -0.236   0.032
+    3TP      TP   39   0.229   0.713   0.346
+    3TP      TP   40   0.259   0.174  -0.026
+    3TP      TP   41   0.393   0.265   0.474
+    3TP      TP   42   0.793   0.424   0.412
+    3TP      TP   43  -0.097   0.505   0.082
+    3TP      TP   44  -0.122   0.182  -0.684
+    3TP      TP   45   0.195   0.033  -0.736
+    3TP      TP   46  -0.034  -0.012   0.938
+    3TP      TP   47   0.154   0.496   0.163
+    3TP      TP   48   0.543   0.698   0.300
+    3TP      TP   49   0.347   0.442   0.174
+    3TP      TP   50   0.322   0.462   0.709
+    3TP      TP   51   0.114  -0.072  -0.742
+    3TP      TP   52  -0.063   0.000   0.978
+    3TP      TP   53   0.069   0.459   0.844
+    3TP      TP   54  -0.022   0.832  -0.150
+    3TP      TP   55   0.337   0.407   0.108
+    3TP      TP   56   0.196  -0.105  -0.934
+    3TP      TP   57   0.221   0.616   0.627
+    3TP      TP   58   0.310   0.469  -0.109
+    3TP      TP   59  -0.021   0.492   0.270
+    3TP      TP   60  -0.052   0.032   0.782
+    3TP      TP   61   0.159   0.063  -0.842
+    3TP      TP   62   0.308   0.111   0.474
+    3TP      TP   63  -0.628   0.712   0.129
+    3TP      TP   64   0.095   0.104   0.353
+    3TP      TP   65  -0.247   0.779   0.277
+    3TP      TP   66   0.059   0.620   0.000
+    3TP      TP   67   0.423   0.385   0.792
+    3TP      TP   68   0.111   0.216   0.322
+    3TP      TP   69   0.480   0.427  -0.093
+    3TP      TP   70  -0.484  -0.404   0.234
+    3TP      TP   71   0.139   0.173   0.970
+    3TP      TP   72  -0.268   0.144  -0.932
+    3TP      TP   73  -0.308  -0.226   0.132
+    3TP      TP   74   0.044   0.165  -0.685
+    3TP      TP   75  -0.790   0.527  -0.058
+    3TP      TP   76   0.403   0.396   0.699
+    3TP      TP   77   0.535   0.704   0.157
+    3TP      TP   78   0.635   0.243   0.715
+    3TP      TP   79   0.445   0.193   0.676
+    3TP      TP   80   0.399   0.450   0.703
+    3TP      TP   81  -0.441   0.810   0.269
+    3TP      TP   82   0.045  -0.064   0.789
+    3TP      TP   83   0.574   0.461  -0.165
+    3TP      TP   84   0.037  -0.114  -0.601
+    3TP      TP   85  -0.112  -0.166  -0.667
+    3TP      TP   86   0.137   0.459   0.258
+    3TP      TP   87   0.273   0.409   0.692
+    3TP      TP   88   0.375   0.423   0.573
+    3TP      TP   89   0.470   0.512   0.204
+    3TP      TP   90  -0.466   0.304   0.184
+    3TP      TP   91  -0.504   0.698   0.113
+    3TP      TP   92   0.387   0.413   0.614
+    3TP      TP   93   0.027   0.066   0.286
+    3TP      TP   94  -0.359  -0.777   0.258
+    3TP      TP   95  -0.019   0.251  -0.818
+    3TP      TP   96   0.114   0.366   0.481
+    3TP      TP   97   0.395   0.150   0.566
+    3TP      TP   98   0.242   0.337   0.467
+    3TP      TP   99   0.547   0.487  -0.097
+    3TP      TP  100  -0.040   0.433   0.804
+    3TP      TP  101   0.305   0.155   0.724
+    3TP      TP  102   0.465   0.092   0.774
+    3TP      TP  103  -0.127   0.190  -0.864
+    3TP      TP  104   0.009  -0.028   0.409
+    3TP      TP  105   0.103  -0.055   0.678
+    3TP      TP  106   0.404   0.316   0.485
+    3TP      TP  107   0.181   0.531  -0.043
+    3TP      TP  108  -0.513  -0.664   0.053
+    3TP      TP  109  -0.442  -0.413   0.265
+    3TP      TP  110  -0.348  -0.619   0.238
+    3TP      TP  111  -0.736   0.637   0.111
+    3TP      TP  112  -0.371  -0.575   0.500
+    3TP      TP  113  -0.365  -0.756   0.320
+    3TP      TP  114   0.570   0.357   0.420
+    3TP      TP  115   0.136  -0.071   0.754
+    3TP      TP  116   0.169  -0.127  -0.670
+    3TP      TP  117   0.272   0.432  -0.025
+    3TP      TP  118  -0.646   0.524   0.361
+    3TP      TP  119   0.166   0.004  -0.838
+    3TP      TP  120   0.303   0.017  -0.928
+    3TP      TP  121  -0.719   0.673   0.144
+    3TP      TP  122  -0.010   0.840   0.017
+    3TP      TP  123   0.260   0.105  -0.842
+    3TP      TP  124  -0.104   0.593   0.349
+    3TP      TP  125   0.311   0.397   0.170
+    3TP      TP  126   0.022   0.640   0.173
+    3TP      TP  127  -0.177   0.601   0.187
+    3TP      TP  128  -0.043   0.903  -0.044
+    3TP      TP  129  -0.051   0.276   0.587
+    3TP      TP  130  -0.113   0.907   0.375
+    3TP      TP  131  -0.357   0.568   0.237
+    3TP      TP  132  -0.210   0.824   0.283
+    3TP      TP  133  -0.386   0.479   0.144
+    3TP      TP  134  -0.252   0.411  -0.088
+    3TP      TP  135   0.183   0.506   0.629
+    3TP      TP  136   0.046  -0.196  -0.856
+    3TP      TP  137   0.416   0.506  -0.023
+    3TP      TP  138   0.247   0.674  -0.080
+    3TP      TP  139  -0.811  -0.521   0.112
+    3TP      TP  140   0.216   0.184   0.426
+    3TP      TP  141   0.050   0.589   0.292
+    3TP      TP  142   0.114   0.330   0.309
+    3TP      TP  143   0.261   0.332   0.023
+    3TP      TP  144   0.393   0.610   0.424
+    3TP      TP  145  -0.057   0.435   0.095
+    3TP      TP  146  -0.003   0.839  -0.095
+    3TP      TP  147  -0.236   0.798   0.317
+    3TP      TP  148   0.648   0.330   0.493
+    3TP      TP  149   0.560   0.549   0.084
+    3TP      TP  150   0.163   0.298  -0.100
+    3TP      TP  151  -0.773  -0.420   0.268
+    3TP      TP  152   0.473   0.217   0.315
+    3TP      TP  153  -0.532  -0.317   0.424
+    3TP      TP  154   0.130   0.213   0.742
+    3TP      TP  155   0.179   0.176   0.142
+    3TP      TP  156   0.440   0.889  -0.111
+    3TP      TP  157  -0.539   0.818   0.146
+    3TP      TP  158   0.119  -0.017  -0.664
+    3TP      TP  159   0.072   0.350   0.344
+    3TP      TP  160  -0.690   0.499   0.161
+    3TP      TP  161  -0.006   0.725   0.436
+    3TP      TP  162  -0.067  -0.052  -0.861
+    3TP      TP  163  -0.174   0.521   0.344
+    3TP      TP  164   0.199   0.613  -0.178
+    3TP      TP  165   0.342   0.725  -0.194
+    3TP      TP  166  -0.488  -0.556   0.365
+    3TP      TP  167  -0.726  -0.590   0.264
+    3TP      TP  168  -0.266   0.637   0.023
+    3TP      TP  169   0.009   0.412   0.431
+    3TP      TP  170  -0.002   0.024  -0.714
+    3TP      TP  171   0.210   0.274   0.344
+    3TP      TP  172   0.769   0.426   0.427
+    3TP      TP  173   0.137   0.731   0.579
+    3TP      TP  174  -0.529  -0.282   0.382
+    3TP      TP  175  -0.132   0.699   0.099
+    3TP      TP  176  -0.727  -0.608   0.079
+    3TP      TP  177   0.036  -0.222  -0.958
+    3TP      TP  178  -0.431   0.515   0.184
+    3TP      TP  179  -0.443   0.301  -0.020
+    3TP      TP  180   0.241   0.027   0.773
+    3TP      TP  181   0.264   0.488   0.496
+    3TP      TP  182  -0.719   0.675   0.025
+    3TP      TP  183   0.572   0.602  -0.235
+    3TP      TP  184  -0.351   0.368   0.102
+    3TP      TP  185   0.731   0.321   0.545
+    3TP      TP  186   0.180   0.142   0.360
+    3TP      TP  187   0.048   0.320   0.361
+    3TP      TP  188   0.460   0.365   0.491
+    3TP      TP  189  -0.382   0.913   0.019
+    3TP      TP  190  -0.039   0.658  -0.083
+    3TP      TP  191   0.371   0.701   0.533
+    3TP      TP  192  -0.092   0.462   0.281
+    3TP      TP  193   0.547   0.395   0.143
+    3TP      TP  194  -0.164   0.664   0.564
+    3TP      TP  195   0.308   0.885   0.198
+    3TP      TP  196   0.238   0.275   0.553
+    3TP      TP  197  -0.455   0.813  -0.028
+    3TP      TP  198  -0.593  -0.448   0.217
+    3TP      TP  199  -0.242   0.773   0.356
+    3TP      TP  200   0.055   0.076   0.288
+    3TP      TP  201  -0.233  -0.243   0.270
+    3TP      TP  202   0.113   0.221   0.149
+    3TP      TP  203   0.140   0.250   0.799
+    3TP      TP  204   0.214  -0.055  -0.945
+    3TP      TP  205   0.158   0.067  -0.681
+    3TP      TP  206   0.069   0.285   0.901
+    3TP      TP  207  -0.426  -0.489   0.016
+    3TP      TP  208  -0.575   0.419   0.043
+    3TP      TP  209   0.683   0.677  -0.248
+    3TP      TP  210  -0.414  -0.598   0.528
+    3TP      TP  211   0.309   0.377   0.483
+    3TP      TP  212   0.507   0.598   0.154
+    3TP      TP  213   0.084  -0.035   0.544
+    3TP      TP  214   0.783   0.564   0.057
+    3TP      TP  215  -0.057  -0.049  -0.804
+    3TP      TP  216  -0.702   0.620   0.037
+    3TP      TP  217  -0.125   0.327   0.885
+    3TP      TP  218   0.032   0.157  -0.817
+    3TP      TP  219   0.112   0.890   0.413
+    3TP      TP  220  -0.421  -0.430   0.471
+    3TP      TP  221  -0.073   0.938  -0.145
+    3TP      TP  222   0.252   0.926  -0.009
+    3TP      TP  223  -0.423   0.509   0.010
+    3TP      TP  224   0.230  -0.076  -0.815
+    3TP      TP  225   0.183   0.888  -0.105
+    3TP      TP  226   0.332   0.755   0.356
+    3TP      TP  227   0.518   0.632   0.245
+    3TP      TP  228  -0.074   0.189   0.509
+    3TP      TP  229  -0.518  -0.749   0.239
+    3TP      TP  230  -0.026   0.109  -0.603
+    3TP      TP  231  -0.603  -0.306   0.263
+    3TP      TP  232  -0.417   0.486   0.161
+    3TP      TP  233   0.026   0.831   0.056
+    3TP      TP  234   0.345   0.839   0.398
+    3TP      TP  235  -0.690  -0.535   0.234
+    3TP      TP  236  -0.422  -0.472   0.394
+    3TP      TP  237  -0.559   0.434   0.275
+    3TP      TP  238  -0.539   0.492   0.105
+    3TP      TP  239  -0.493  -0.466   0.280
+    3TP      TP  240  -0.236   0.166  -0.934
+    3TP      TP  241   0.324   0.700   0.336
+    3TP      TP  242  -0.781  -0.447   0.147
+    3TP      TP  243   0.071   0.577  -0.068
+    3TP      TP  244   0.525   0.818  -0.100
+    3TP      TP  245   0.066  -0.028  -0.515
+    3TP      TP  246  -0.032   0.210   0.893
+    3TP      TP  247   0.680   0.581  -0.045
+    3TP      TP  248  -0.396   0.326   0.078
+    3TP      TP  249   0.160   0.150   0.967
+    3TP      TP  250  -0.104   0.968  -0.131
+    3TP      TP  251  -0.129  -0.001  -0.649
+    3TP      TP  252  -0.269   0.805   0.504
+    3TP      TP  253  -0.646   0.749  -0.010
+    3TP      TP  254   0.627   0.538   0.016
+    3TP      TP  255   0.535   0.355   0.179
+    3TP      TP  256   0.125   0.165   0.635
+    3TP      TP  257   0.240   0.265   0.454
+    3TP      TP  258   0.138   0.233   0.932
+    3TP      TP  259   0.068   0.124   0.561
+    3TP      TP  260  -0.243  -0.393   0.202
+    3TP      TP  261  -0.043   0.117   0.567
+    3TP      TP  262   0.522   0.692  -0.236
+    3TP      TP  263  -0.619  -0.685   0.150
+    3TP      TP  264  -0.194   0.562   0.446
+    3TP      TP  265   0.394   0.345   0.502
+    3TP      TP  266   0.123   0.405   0.371
+    3TP      TP  267  -0.199   0.715   0.234
+    3TP      TP  268   0.410   0.518   0.175
+    3TP      TP  269  -0.073   0.204   0.925
+    3TP      TP  270  -0.075   0.340   0.667
+    3TP      TP  271  -0.055   0.595   0.714
+    3TP      TP  272   0.354   0.245   0.397
+    3TP      TP  273  -0.002   0.136  -0.759
+    3TP      TP  274   0.463   0.476   0.547
+    3TP      TP  275   0.593   0.376   0.257
+    3TP      TP  276   0.255   0.229   0.856
+    3TP      TP  277  -0.658   0.616  -0.020
+    3TP      TP  278  -0.292  -0.500   0.123
+    3TP      TP  279   0.051   0.662   0.568
+    3TP      TP  280   0.489   0.077   0.861
+    3TP      TP  281   0.168   0.255   0.644
+    3TP      TP  282  -0.320  -0.589   0.393
+    3TP      TP  283   0.189   0.267   0.762
+    3TP      TP  284  -0.535   0.364   0.202
+    3TP      TP  285   0.003   0.743  -0.024
+    3TP      TP  286  -0.130   0.703   0.342
+    3TP      TP  287  -0.104   0.899   0.064
+    3TP      TP  288  -0.288  -0.452   0.055
+    3TP      TP  289   0.327   0.791   0.209
+    3TP      TP  290   0.174  -0.065  -0.616
+    3TP      TP  291   0.545   0.334   0.414
+    3TP      TP  292   0.254   0.646   0.207
+    3TP      TP  293   0.424   0.676   0.399
+    3TP      TP  294   0.197   0.967  -0.066
+    3TP      TP  295  -0.483  -0.286   0.329
+    3TP      TP  296  -0.726  -0.387   0.493
+    3TP      TP  297  -0.542  -0.753   0.300
+    3TP      TP  298  -0.632  -0.450   0.356
+    3TP      TP  299  -0.119   0.441   0.302
+    3TP      TP  300   0.163   0.331   0.418
+    3TP      TP  301  -0.013   0.036  -0.520
+    3TP      TP  302  -0.650  -0.550   0.158
+    3TP      TP  303   0.265   0.240   0.799
+    3TP      TP  304  -0.013   0.839   0.338
+    3TP      TP  305   0.067   0.155  -0.753
+    3TP      TP  306  -0.132   0.020  -0.872
+    3TP      TP  307  -0.533   0.759  -0.097
+    3TP      TP  308   0.029   0.198   0.814
+    3TP      TP  309   0.568   0.303   0.740
+    3TP      TP  310  -0.066   0.431   0.649
+    3TP      TP  311   0.423   0.820   0.106
+    3TP      TP  312  -0.484   0.693  -0.123
+    3TP      TP  313  -0.072   0.394   0.251
+    3TP      TP  314  -0.232   0.172  -0.834
+    3TP      TP  315  -0.562   0.522   0.274
+    3TP      TP  316  -0.511   0.483   0.272
+    3TP      TP  317   0.039   0.685   0.569
+    3TP      TP  318  -0.155  -0.177   0.179
+    3TP      TP  319  -0.399   0.429   0.151
+    3TP      TP  320   0.354   0.299   0.713
+    3TP      TP  321  -0.269   0.729   0.539
+    3TP      TP  322   0.483   0.224   0.679
+    3TP      TP  323  -0.005   0.189   0.425
+    3TP      TP  324   0.179   0.807   0.513
+    3TP      TP  325   0.348   0.408   0.459
+    3TP      TP  326  -0.070   0.990   0.089
+    3TP      TP  327   0.133   0.210   0.881
+    3TP      TP  328  -0.386   0.777   0.177
+    3TP      TP  329   0.263   0.058   0.862
+    3TP      TP  330   0.058   0.295   0.047
+    3TP      TP  331  -0.163  -0.139  -0.857
+    3TP      TP  332   0.071   0.309  -0.945
+    3TP      TP  333  -0.006   0.399   0.554
+    3TP      TP  334  -0.185  -0.087  -0.777
+    3TP      TP  335  -0.036   0.221  -0.753
+    3TP      TP  336  -0.041   0.638   0.560
+    3TP      TP  337  -0.514  -0.423   0.168
+    3TP      TP  338   0.131  -0.044   0.851
+    3TP      TP  339   0.049   0.127   0.959
+    3TP      TP  340  -0.392   0.543   0.071
+    3TP      TP  341  -0.199   0.051  -0.752
+    3TP      TP  342  -0.010   0.240   0.913
+    3TP      TP  343   0.334   0.561   0.511
+    3TP      TP  344  -0.133   0.657   0.379
+    3TP      TP  345   0.401   0.313   0.342
+    3TP      TP  346  -0.178  -0.285   0.183
+    3TP      TP  347   0.434   0.887   0.058
+    3TP      TP  348  -0.413   0.713   0.177
+    3TP      TP  349   0.055   0.150  -0.460
+    3TP      TP  350  -0.638  -0.610   0.416
+    3TP      TP  351  -0.038   0.836   0.082
+    3TP      TP  352  -0.646  -0.524   0.349
+    3TP      TP  353  -0.252   0.169   0.107
+    3TP      TP  354  -0.455  -0.445   0.125
+    3TP      TP  355   0.049   0.461   0.617
+    3TP      TP  356   0.442   0.393   0.534
+    3TP      TP  357   0.313   0.498  -0.074
+    3TP      TP  358  -0.050  -0.110  -0.805
+    3TP      TP  359  -0.698  -0.397   0.337
+    3TP      TP  360  -0.204   0.267   0.044
+    3TP      TP  361   0.502   0.123   0.687
+    3TP      TP  362  -0.012  -0.127  -0.688
+    3TP      TP  363  -0.138   0.081  -0.672
+    3TP      TP  364  -0.352  -0.561   0.184
+    3TP      TP  365   0.722   0.479   0.015
+    3TP      TP  366  -0.684   0.629   0.167
+    3TP      TP  367   0.265   0.955   0.042
+    3TP      TP  368   0.109   0.165   0.973
+    3TP      TP  369   0.660   0.446   0.195
+    3TP      TP  370  -0.026   0.676   0.551
+    3TP      TP  371  -0.610  -0.463   0.016
+    3TP      TP  372   0.380   0.136   0.758
+    3TP      TP  373   0.231   0.823  -0.029
+    3TP      TP  374   0.355   0.792  -0.067
+    3TP      TP  375  -0.137   0.213   0.900
+    3TP      TP  376   0.021   0.850   0.143
+    3TP      TP  377   0.240   0.040  -0.686
+    3TP      TP  378   0.329   0.223   0.480
+    3TP      TP  379   0.170   0.623   0.552
+    3TP      TP  380   0.061   0.939   0.114
+    3TP      TP  381  -0.028   0.855   0.471
+    3TP      TP  382   0.229   0.739   0.294
+    3TP      TP  383   0.107   0.860   0.357
+    3TP      TP  384  -0.677  -0.398   0.440
+    3TP      TP  385   0.180   0.324   0.792
+    3TP      TP  386  -0.040   0.355   0.317
+    3TP      TP  387   0.138   0.801   0.158
+    3TP      TP  388  -0.105   0.706   0.148
+    3TP      TP  389   0.118  -0.222  -0.706
+    3TP      TP  390  -0.632  -0.457   0.364
+    3TP      TP  391   0.307   0.426   0.466
+    3TP      TP  392   0.000   0.708   0.073
+    3TP      TP  393  -0.477  -0.335   0.182
+    3TP      TP  394   0.078   0.964  -0.190
+    3TP      TP  395   0.031   0.025  -0.818
+    3TP      TP  396   0.334  -0.052   0.939
+    3TP      TP  397   0.232   0.580   0.558
+    3TP      TP  398   0.103   0.165   0.600
+    3TP      TP  399   0.401   0.557   0.037
+    3TP      TP  400   0.101   0.889   0.325
+    3TP      TP  401  -0.091   0.969  -0.014
+    3TP      TP  402   0.051   0.904   0.376
+    3TP      TP  403  -0.317   0.599   0.002
+    3TP      TP  404   0.315   0.065   0.744
+    3TP      TP  405   0.122   0.572   0.629
+    3TP      TP  406  -0.254  -0.063  -0.869
+    3TP      TP  407  -0.046  -0.078  -0.559
+    3TP      TP  408   0.187   0.144   0.206
+    3TP      TP  409   0.236   0.030   0.437
+    3TP      TP  410   0.191   0.648   0.035
+    3TP      TP  411   0.035   0.546  -0.062
+    3TP      TP  412   0.669   0.544   0.493
+    3TP      TP  413  -0.456   0.702   0.119
+    3TP      TP  414   0.167  -0.021   0.928
+    3TP      TP  415  -0.117   0.926   0.094
+    3TP      TP  416  -0.358  -0.732   0.413
+    3TP      TP  417  -0.579   0.787   0.028
+    3TP      TP  418   0.344   0.518   0.336
+    3TP      TP  419  -0.029   0.711   0.097
+    3TP      TP  420  -0.092   0.855   0.021
+    3TP      TP  421  -0.627   0.386   0.102
+    3TP      TP  422  -0.414   0.788   0.358
+    3TP      TP  423  -0.470   0.695   0.302
+    3TP      TP  424   0.469   0.415   0.536
+    3TP      TP  425   0.431   0.516   0.691
+    3TP      TP  426  -0.314   0.531   0.280
+    3TP      TP  427  -0.732   0.611   0.206
+    3TP      TP  428   0.197  -0.064  -0.667
+    3TP      TP  429  -0.368  -0.252   0.170
+    3TP      TP  430   0.110  -0.029   0.664
+    3TP      TP  431   0.551   0.128   0.761
+    3TP      TP  432  -0.129   0.373   0.230
+    3TP      TP  433   0.064   0.862   0.230
+    3TP      TP  434   0.436   0.738  -0.095
+    3TP      TP  435  -0.699  -0.546   0.058
+    3TP      TP  436  -0.659  -0.313   0.349
+    3TP      TP  437   0.064   0.073   0.728
+    3TP      TP  438   0.582   0.728   0.341
+    3TP      TP  439  -0.077   0.174  -0.913
+    3TP      TP  440   0.006   0.648   0.000
+    3TP      TP  441   0.138   0.485   0.779
+    3TP      TP  442   0.596   0.363   0.109
+    3TP      TP  443   0.460   0.500   0.686
+    3TP      TP  444  -0.043   0.132  -0.865
+    3TP      TP  445   0.083   0.219   0.313
+    3TP      TP  446   0.041   0.761   0.441
+    3TP      TP  447   0.095   0.801  -0.195
+    3TP      TP  448  -0.277   0.627  -0.048
+    3TP      TP  449   0.042   0.368  -0.057
+    3TP      TP  450  -0.469   0.803  -0.107
+    3TP      TP  451  -0.393   0.841   0.182
+    3TP      TP  452  -0.568   0.677  -0.185
+    3TP      TP  453   0.007   0.126   0.721
+    3TP      TP  454   0.333   0.282   0.199
+    3TP      TP  455  -0.278   0.676   0.075
+    3TP      TP  456  -0.367  -0.578   0.368
+    3TP      TP  457  -0.532  -0.667   0.119
+    3TP      TP  458   0.628   0.721   0.225
+    3TP      TP  459   0.555   0.761  -0.010
+    3TP      TP  460  -0.066   0.137   0.882
+    3TP      TP  461   0.105   0.674  -0.014
+    3TP      TP  462  -0.222   0.695   0.572
+    3TP      TP  463   0.391   0.089   0.643
+    3TP      TP  464   0.589   0.300   0.272
+    3TP      TP  465   0.134  -0.072  -0.658
+    3TP      TP  466   0.724   0.676   0.047
+    3TP      TP  467   0.594   0.345   0.701
+    3TP      TP  468  -0.470   0.339   0.008
+    3TP      TP  469   0.455   0.823   0.117
+    3TP      TP  470   0.415   0.381   0.550
+    3TP      TP  471  -0.158   0.174  -0.818
+    3TP      TP  472  -0.250   0.123  -0.852
+    3TP      TP  473  -0.406  -0.794   0.241
+    3TP      TP  474   0.748   0.660   0.035
+    3TP      TP  475   0.276   0.701   0.056
+    3TP      TP  476   0.448   0.702   0.170
+    3TP      TP  477  -0.577  -0.682   0.417
+    3TP      TP  478   0.400   0.769   0.205
+    3TP      TP  479  -0.443  -0.691   0.076
+    3TP      TP  480   0.512   0.809   0.272
+    3TP      TP  481  -0.520   0.438  -0.077
+    3TP      TP  482   0.371   0.378   0.109
+    3TP      TP  483   0.183   0.691   0.519
+    3TP      TP  484   0.244   0.592  -0.132
+    3TP      TP  485  -0.129   0.330   0.670
+    3TP      TP  486  -0.049  -0.018   0.757
+    3TP      TP  487  -0.212   0.430   0.021
+    3TP      TP  488   0.355   0.241   0.734
+    3TP      TP  489  -0.026   0.252   0.436
+    3TP      TP  490  -0.479  -0.275   0.253
+    3TP      TP  491   0.257   0.895   0.344
+    3TP      TP  492  -0.009   0.479   0.410
+    3TP      TP  493  -0.180   0.013  -0.845
+    3TP      TP  494   0.081   0.881   0.036
+    3TP      TP  495   0.022   0.277  -0.864
+    3TP      TP  496   0.210   0.357   0.000
+    3TP      TP  497  -0.037   0.184  -0.562
+    3TP      TP  498   0.080   0.779   0.205
+    3TP      TP  499  -0.149   0.745   0.440
+    3TP      TP  500  -0.343  -0.418   0.003
+    3TP      TP  501   0.034   0.613   0.259
+    3TP      TP  502  -0.001   0.770   0.387
+    3TP      TP  503   0.762   0.431   0.303
+    3TP      TP  504  -0.385  -0.635   0.419
+    3TP      TP  505  -0.225  -0.274   0.024
+    3TP      TP  506  -0.065   0.742   0.069
+    3TP      TP  507   0.132  -0.082   0.831
+    3TP      TP  508   0.338   0.736   0.488
+    3TP      TP  509  -0.390   0.723   0.348
+    3TP      TP  510  -0.174   0.707   0.195
+    3TP      TP  511   0.348   0.685   0.132
+    3TP      TP  512   0.639   0.742   0.120
+    3TP      TP  513   0.236   0.894  -0.192
+    3TP      TP  514  -0.396   0.624  -0.106
+    3TP      TP  515  -0.059   0.041   0.416
+    3TP      TP  516   0.237   0.422  -0.108
+    3TP      TP  517  -0.621  -0.438   0.500
+    3TP      TP  518  -0.356  -0.438   0.408
+    3TP      TP  519   0.710   0.481   0.155
+    3TP      TP  520   0.314   0.337  -0.015
+    3TP      TP  521   0.459   0.498   0.456
+    3TP      TP  522   0.031  -0.103  -0.354
+    3TP      TP  523   0.013   0.736   0.264
+    3TP      TP  524   0.132   0.629   0.712
+    3TP      TP  525   0.030   0.049   0.829
+    3TP      TP  526   0.086  -0.065  -0.771
+    3TP      TP  527  -0.638  -0.565   0.136
+    3TP      TP  528   0.628   0.592   0.250
+    3TP      TP  529   0.277   0.481   0.750
+    3TP      TP  530  -0.242  -0.247   0.103
+    3TP      TP  531  -0.606  -0.531   0.486
+    4TN      TN  532  -0.288   0.785   0.014
+    4TN      TN  533   0.089  -0.279   0.940
+    4TN      TN  534   0.668  -0.566   0.371
+    4TN      TN  535  -0.104   0.120   0.307
+    4TN      TN  536   0.053  -0.982   0.131
+    4TN      TN  537   0.024  -0.890  -0.231
+    4TN      TN  538   0.070  -0.333   0.787
+    4TN      TN  539  -0.280   0.094  -0.348
+    4TN      TN  540   0.477  -0.649   0.265
+    4TN      TN  541   0.792   0.249  -0.334
+    4TN      TN  542  -0.870   0.097   0.419
+    4TN      TN  543  -0.538  -0.002  -0.163
+    4TN      TN  544   0.146   0.330  -0.738
+    4TN      TN  545   0.401  -0.051  -0.007
+    4TN      TN  546   0.040   0.097  -0.249
+    4TN      TN  547  -0.782  -0.484  -0.255
+    4TN      TN  548   0.870  -0.253   0.297
+    4TN      TN  549   0.447  -0.108  -0.879
+    4TN      TN  550  -0.767  -0.027  -0.514
+    4TN      TN  551   0.761  -0.270  -0.425
+    4TN      TN  552   0.523  -0.501   0.352
+    4TN      TN  553  -0.507   0.363  -0.335
+    4TN      TN  554  -0.350  -0.201  -0.584
+    4TN      TN  555   0.217  -0.540   0.711
+    4TN      TN  556   0.441   0.226  -0.604
+    4TN      TN  557   0.191  -0.694   0.512
+    4TN      TN  558  -0.061  -0.831   0.039
+    4TN      TN  559  -0.364   0.347  -0.793
+    4TN      TN  560  -0.778  -0.323  -0.491
+    4TN      TN  561   0.197   0.614  -0.712
+    4TN      TN  562  -0.172   0.380   0.822
+    4TN      TN  563  -0.429   0.286  -0.704
+    4TN      TN  564  -0.159  -0.529  -0.073
+    4TN      TN  565  -0.191   0.694   0.136
+    4TN      TN  566   0.040   0.455  -0.521
+    4TN      TN  567   0.188   0.254  -0.161
+    4TN      TN  568  -0.030  -0.392  -0.713
+    4TN      TN  569   0.294  -0.768  -0.525
+    4TN      TN  570  -0.113   0.216  -0.502
+    4TN      TN  571   0.235   0.130  -0.108
+    4TN      TN  572  -0.742   0.515   0.378
+    4TN      TN  573  -0.123  -0.625   0.331
+    4TN      TN  574   0.884   0.293  -0.290
+    4TN      TN  575   0.069   0.549  -0.804
+    4TN      TN  576  -0.617  -0.110  -0.544
+    4TN      TN  577   0.559  -0.240  -0.173
+    4TN      TN  578  -0.880   0.108  -0.300
+    4TN      TN  579   0.885   0.209   0.130
+    4TN      TN  580  -0.501   0.598  -0.289
+    4TN      TN  581  -0.533  -0.037   0.126
+    4TN      TN  582   0.392  -0.391  -0.313
+    4TN      TN  583   0.055   0.262  -0.724
+    4TN      TN  584   0.961   0.165  -0.145
+    4TN      TN  585   0.335   0.218  -0.438
+    4TN      TN  586   0.380   0.579  -0.644
+    4TN      TN  587  -0.232  -0.183  -0.687
+    4TN      TN  588   0.268  -0.765   0.314
+    4TN      TN  589  -0.358  -0.448  -0.556
+    4TN      TN  590   0.398  -0.599   0.151
+    4TN      TN  591   0.063  -0.717  -0.128
+    4TN      TN  592  -0.236   0.074   0.168
+    4TN      TN  593  -0.882   0.109  -0.127
+    4TN      TN  594  -0.267   0.546  -0.388
+    4TN      TN  595   0.340   0.046  -0.314
+    4TN      TN  596   0.587  -0.023  -0.599
+    4TN      TN  597   0.616  -0.106  -0.235
+    4TN      TN  598   0.183   0.733  -0.327
+    4TN      TN  599  -0.905   0.073  -0.196
+    4TN      TN  600  -0.469   0.057   0.275
+    4TN      TN  601  -0.743   0.443  -0.103
+    4TN      TN  602  -0.803  -0.075  -0.124
+    4TN      TN  603  -0.308  -0.135   0.204
+    4TN      TN  604   0.574   0.568  -0.549
+    4TN      TN  605  -0.216   0.175  -0.444
+    4TN      TN  606  -0.095   0.255  -0.041
+    4TN      TN  607   0.802  -0.257   0.361
+    4TN      TN  608   0.680   0.072   0.202
+    4TN      TN  609   0.334   0.254  -0.900
+    4TN      TN  610  -0.425  -0.312   0.431
+    4TN      TN  611  -0.197   0.306   0.788
+    4TN      TN  612  -0.239   0.658   0.160
+    4TN      TN  613  -0.065  -0.325  -0.417
+    4TN      TN  614  -0.768  -0.293  -0.339
+    4TN      TN  615   0.407  -0.010   0.720
+    4TN      TN  616   0.439  -0.058   0.641
+    4TN      TN  617  -0.007   0.244  -0.579
+    4TN      TN  618  -0.495   0.626  -0.223
+    4TN      TN  619   0.423   0.103  -0.324
+    4TN      TN  620  -0.449  -0.141   0.419
+    4TN      TN  621  -0.797  -0.018  -0.173
+    4TN      TN  622   0.012  -0.294  -0.175
+    4TN      TN  623  -0.507  -0.700  -0.142
+    4TN      TN  624  -0.388  -0.092   0.386
+    4TN      TN  625  -0.554   0.064  -0.750
+    4TN      TN  626   0.617  -0.453  -0.396
+    4TN      TN  627   0.835   0.146  -0.159
+    4TN      TN  628  -0.754  -0.318  -0.468
+    4TN      TN  629  -0.264  -0.685  -0.233
+    4TN      TN  630  -0.382   0.213  -0.731
+    4TN      TN  631   0.582  -0.272  -0.240
+    4TN      TN  632   0.324   0.626  -0.475
+    4TN      TN  633  -0.723   0.280   0.104
+    4TN      TN  634  -0.236  -0.408  -0.222
+    4TN      TN  635   0.491  -0.767   0.359
+    4TN      TN  636  -0.820  -0.425   0.117
+    4TN      TN  637   0.739  -0.353  -0.169
+    4TN      TN  638  -0.395  -0.295   0.812
+    4TN      TN  639  -0.522  -0.016  -0.339
+    4TN      TN  640   0.664  -0.377   0.116
+    4TN      TN  641   0.377   0.130   0.152
+    4TN      TN  642   0.048  -0.580   0.768
+    4TN      TN  643  -0.241   0.914  -0.065
+    4TN      TN  644   0.195  -0.572  -0.148
+    4TN      TN  645   0.510   0.337  -0.165
+    4TN      TN  646  -0.316  -0.049   0.420
+    4TN      TN  647   0.095  -0.509  -0.119
+    4TN      TN  648  -0.618  -0.026   0.682
+    4TN      TN  649  -0.633  -0.519  -0.330
+    4TN      TN  650   0.504  -0.331  -0.479
+    4TN      TN  651  -0.031  -0.105   0.624
+    4TN      TN  652  -0.290  -0.637   0.683
+    4TN      TN  653   0.373  -0.070  -0.619
+    4TN      TN  654   0.592  -0.298   0.207
+    4TN      TN  655   0.501  -0.285   0.222
+    4TN      TN  656   0.651  -0.651   0.255
+    4TN      TN  657   0.522   0.059  -0.043
+    4TN      TN  658   0.280  -0.278  -0.023
+    4TN      TN  659  -0.759  -0.005  -0.260
+    4TN      TN  660  -0.691  -0.233  -0.243
+    4TN      TN  661   0.045  -0.690  -0.275
+    4TN      TN  662   0.743  -0.391  -0.186
+    4TN      TN  663   0.589  -0.404   0.619
+    4TN      TN  664  -0.666   0.523  -0.511
+    4TN      TN  665  -0.292   0.606  -0.655
+    4TN      TN  666  -0.078  -0.247  -0.386
+    4TN      TN  667  -0.423  -0.425   0.508
+    4TN      TN  668   0.375  -0.086  -0.847
+    4TN      TN  669   0.088  -0.260   0.004
+    4TN      TN  670   0.137   0.373  -0.450
+    4TN      TN  671   0.439  -0.252   0.221
+    4TN      TN  672  -0.233   0.552   0.514
+    4TN      TN  673   0.268  -0.252  -0.802
+    4TN      TN  674   0.774   0.084  -0.438
+    4TN      TN  675   0.677  -0.548  -0.031
+    4TN      TN  676  -0.253   0.163  -0.391
+    4TN      TN  677  -0.198  -0.627  -0.687
+    4TN      TN  678  -0.294   0.581   0.659
+    4TN      TN  679  -0.267  -0.379   0.850
+    4TN      TN  680  -0.197  -0.430  -0.313
+    4TN      TN  681  -0.406  -0.492  -0.310
+    4TN      TN  682  -0.008  -0.668   0.168
+    4TN      TN  683   0.141  -0.472  -0.486
+    4TN      TN  684  -0.539   0.701  -0.429
+    4TN      TN  685  -0.378  -0.597  -0.007
+    4TN      TN  686   0.180  -0.135   0.171
+    4TN      TN  687   0.331  -0.376  -0.710
+    4TN      TN  688   0.358  -0.158  -0.040
+    4TN      TN  689  -0.310  -0.855  -0.039
+    4TN      TN  690  -0.385  -0.450  -0.242
+    4TN      TN  691  -0.724  -0.480  -0.054
+    4TN      TN  692  -0.763  -0.487  -0.307
+    4TN      TN  693  -0.263   0.829  -0.066
+    4TN      TN  694  -0.225   0.134   0.747
+    4TN      TN  695  -0.308   0.228  -0.517
+    4TN      TN  696   0.797   0.341   0.290
+    4TN      TN  697  -0.134  -0.767  -0.163
+    4TN      TN  698  -0.524   0.175   0.193
+    4TN      TN  699  -0.083  -0.753   0.583
+    4TN      TN  700   0.491   0.688  -0.458
+    4TN      TN  701   0.388  -0.469   0.225
+    4TN      TN  702   0.725   0.005  -0.231
+    4TN      TN  703  -0.278  -0.238  -0.561
+    4TN      TN  704  -0.697   0.255  -0.485
+    4TN      TN  705  -0.044   0.421  -0.768
+    4TN      TN  706  -0.355   0.016  -0.762
+    4TN      TN  707   0.060  -0.283  -0.460
+    4TN      TN  708  -0.379  -0.441  -0.097
+    4TN      TN  709  -0.649   0.359  -0.629
+    4TN      TN  710   0.892   0.073  -0.366
+    4TN      TN  711  -0.581   0.061   0.776
+    4TN      TN  712  -0.355   0.343  -0.271
+    4TN      TN  713  -0.553  -0.016  -0.064
+    4TN      TN  714   0.467   0.237  -0.836
+    4TN      TN  715   0.759   0.053   0.504
+    4TN      TN  716  -0.223   0.116   0.877
+    4TN      TN  717   0.760  -0.558  -0.313
+    4TN      TN  718   0.254  -0.618   0.056
+    4TN      TN  719   0.380  -0.673  -0.195
+    4TN      TN  720   0.570  -0.006  -0.064
+    4TN      TN  721   0.454   0.074  -0.234
+    4TN      TN  722   0.561  -0.185  -0.419
+    4TN      TN  723  -0.679   0.070   0.621
+    4TN      TN  724   0.655   0.114  -0.379
+    4TN      TN  725   0.315  -0.868   0.075
+    4TN      TN  726   0.878  -0.357  -0.302
+    4TN      TN  727  -0.692   0.136   0.224
+    4TN      TN  728   0.793  -0.584   0.031
+    4TN      TN  729   0.306  -0.319   0.711
+    4TN      TN  730   0.615   0.325  -0.341
+    4TN      TN  731  -0.841   0.492  -0.086
+    4TN      TN  732   0.731   0.121  -0.100
+    4TN      TN  733  -0.679   0.000   0.228
+    4TN      TN  734  -0.201   0.558   0.187
+    4TN      TN  735   0.301  -0.207  -0.486
+    4TN      TN  736  -0.577   0.130  -0.121
+    4TN      TN  737   0.411  -0.808  -0.287
+    4TN      TN  738   0.578   0.551  -0.298
+    4TN      TN  739  -0.077  -0.776  -0.547
+    4TN      TN  740   0.457  -0.547  -0.136
+    4TN      TN  741  -0.476   0.009   0.165
+    4TN      TN  742  -0.543  -0.128  -0.275
+    4TN      TN  743  -0.362  -0.519  -0.355
+    4TN      TN  744  -0.084   0.328  -0.886
+    4TN      TN  745   0.391  -0.308   0.602
+    4TN      TN  746  -0.184  -0.806   0.409
+    4TN      TN  747   0.039  -0.458  -0.201
+    4TN      TN  748  -0.009  -0.396  -0.166
+    4TN      TN  749   0.065   0.418  -0.765
+    4TN      TN  750   0.160  -0.611  -0.065
+    4TN      TN  751   0.297  -0.724   0.038
+    4TN      TN  752  -0.797  -0.224   0.463
+    4TN      TN  753  -0.433  -0.815  -0.284
+    4TN      TN  754   0.595   0.316   0.066
+    4TN      TN  755   0.534   0.119  -0.321
+    4TN      TN  756  -0.404  -0.143  -0.682
+    4TN      TN  757   0.699  -0.320  -0.431
+    4TN      TN  758  -0.256  -0.163  -0.117
+    4TN      TN  759  -0.407  -0.209   0.653
+    4TN      TN  760   0.506   0.083   0.304
+    4TN      TN  761   0.365  -0.271  -0.741
+    4TN      TN  762  -0.486  -0.130  -0.027
+    4TN      TN  763  -0.352   0.805  -0.120
+    4TN      TN  764  -0.802  -0.140  -0.060
+    4TN      TN  765  -0.309  -0.168   0.832
+    4TN      TN  766  -0.390  -0.460  -0.243
+    4TN      TN  767   0.568   0.123   0.364
+    4TN      TN  768   0.625  -0.081  -0.471
+    4TN      TN  769   0.031  -0.804   0.165
+    4TN      TN  770  -0.160   0.517   0.151
+    4TN      TN  771  -0.915   0.075  -0.097
+    4TN      TN  772   0.464  -0.668  -0.390
+    4TN      TN  773  -0.370   0.478  -0.621
+    4TN      TN  774   0.606   0.094  -0.339
+    4TN      TN  775   0.546  -0.680   0.300
+    4TN      TN  776   0.399   0.151   0.171
+    4TN      TN  777   0.716  -0.430  -0.115
+    4TN      TN  778  -0.012   0.480  -0.398
+    4TN      TN  779  -0.112  -0.052  -0.272
+    4TN      TN  780  -0.771   0.179  -0.276
+    4TN      TN  781   0.595   0.110   0.595
+    4TN      TN  782  -0.557  -0.357   0.664
+    4TN      TN  783  -0.889   0.081  -0.147
+    4TN      TN  784   0.343  -0.112   0.195
+    4TN      TN  785   0.719  -0.111   0.187
+    4TN      TN  786   0.499  -0.255   0.521
+    4TN      TN  787  -0.375  -0.707   0.527
+    4TN      TN  788   0.890  -0.030  -0.287
+    4TN      TN  789   0.048  -0.566  -0.465
+    4TN      TN  790   0.313  -0.453  -0.748
+    4TN      TN  791   0.098  -0.247  -0.391
+    4TN      TN  792   0.487   0.151   0.072
+    4TN      TN  793   0.169  -0.585  -0.415
+    4TN      TN  794  -0.040   0.374  -0.325
+    4TN      TN  795   0.244  -0.326   0.865
+    4TN      TN  796   0.091  -0.235  -0.452
+    4TN      TN  797  -0.004  -0.871   0.403
+    4TN      TN  798   0.689  -0.162  -0.242
+    4TN      TN  799   0.346  -0.023  -0.081
+    4TN      TN  800   0.090   0.489  -0.757
+    4TN      TN  801   0.378  -0.349  -0.758
+    4TN      TN  802   0.354   0.301  -0.204
+    4TN      TN  803   0.370   0.219  -0.875
+    4TN      TN  804  -0.370   0.278  -0.095
+    4TN      TN  805   0.852  -0.328  -0.147
+    4TN      TN  806  -0.314  -0.251  -0.615
+    4TN      TN  807  -0.556   0.366   0.641
+    4TN      TN  808   0.419  -0.053   0.201
+    4TN      TN  809   0.530   0.431  -0.392
+    4TN      TN  810   0.411   0.197  -0.759
+    4TN      TN  811  -0.898   0.370   0.128
+    4TN      TN  812  -0.310   0.009  -0.534
+    4TN      TN  813  -0.383  -0.117   0.017
+    4TN      TN  814   0.102  -0.904  -0.235
+    4TN      TN  815   0.152  -0.542  -0.140
+    4TN      TN  816   0.771   0.280  -0.527
+    4TN      TN  817  -0.302   0.146   0.176
+    4TN      TN  818   0.602   0.389  -0.451
+    4TN      TN  819   0.414   0.236  -0.017
+    4TN      TN  820  -0.448   0.367   0.556
+    4TN      TN  821  -0.294  -0.032  -0.794
+    4TN      TN  822  -0.403   0.004   0.480
+    4TN      TN  823   0.183   0.281  -0.133
+    4TN      TN  824  -0.566  -0.532  -0.604
+    4TN      TN  825   0.407   0.128   0.182
+    4TN      TN  826   0.592   0.130   0.526
+    4TN      TN  827   0.731  -0.388  -0.026
+    4TN      TN  828  -0.398  -0.801  -0.433
+    4TN      TN  829   0.329  -0.368  -0.457
+    4TN      TN  830  -0.562   0.119   0.283
+    4TN      TN  831  -0.499  -0.092  -0.322
+    4TN      TN  832  -0.051  -0.425  -0.281
+    4TN      TN  833  -0.660  -0.319  -0.676
+    4TN      TN  834  -0.792   0.101   0.468
+    4TN      TN  835   0.204  -0.179  -0.417
+    4TN      TN  836  -0.798   0.191  -0.515
+    4TN      TN  837   0.377  -0.864   0.001
+    4TN      TN  838   0.155  -0.047   0.305
+    4TN      TN  839   0.132   0.193  -0.234
+    4TN      TN  840   0.574  -0.654   0.344
+    4TN      TN  841   0.236  -0.767   0.476
+    4TN      TN  842   0.370   0.854  -0.311
+    4TN      TN  843   0.685  -0.046   0.581
+    4TN      TN  844  -0.002  -0.397  -0.655
+    4TN      TN  845   0.197  -0.905   0.121
+    4TN      TN  846   0.351  -0.114   0.722
+    4TN      TN  847  -0.263   0.150   0.147
+    4TN      TN  848   0.519  -0.837  -0.155
+    4TN      TN  849  -0.365  -0.102  -0.641
+    4TN      TN  850   0.074  -0.334  -0.857
+    4TN      TN  851   0.518  -0.089   0.724
+    4TN      TN  852  -0.121   0.494  -0.808
+    4TN      TN  853  -0.451   0.114   0.249
+    4TN      TN  854   0.160  -0.308   0.891
+    4TN      TN  855  -0.850  -0.020   0.125
+    4TN      TN  856   0.683   0.258   0.244
+    4TN      TN  857  -0.210   0.445  -0.500
+    4TN      TN  858   0.254  -0.285  -0.409
+    4TN      TN  859  -0.872   0.109   0.327
+    4TN      TN  860   0.043  -0.403   0.527
+    4TN      TN  861  -0.138   0.694  -0.125
+    4TN      TN  862   0.752   0.399  -0.270
+    4TN      TN  863  -0.144   0.648  -0.499
+    4TN      TN  864  -0.542   0.195   0.150
+    4TN      TN  865  -0.428   0.081  -0.454
+    4TN      TN  866  -0.211   0.297  -0.344
+    4TN      TN  867  -0.739   0.413   0.276
+    4TN      TN  868   0.865   0.077  -0.158
+    4TN      TN  869  -0.328  -0.394  -0.753
+    4TN      TN  870   0.703  -0.256  -0.498
+    4TN      TN  871   0.263   0.695  -0.600
+    4TN      TN  872  -0.196   0.631  -0.603
+    4TN      TN  873   0.198  -0.313   0.347
+    4TN      TN  874   0.866  -0.280   0.282
+    4TN      TN  875   0.047  -0.312  -0.484
+    4TN      TN  876  -0.214   0.491  -0.388
+    4TN      TN  877  -0.207   0.238   0.769
+    4TN      TN  878   0.635  -0.518  -0.375
+    4TN      TN  879  -0.725  -0.051   0.559
+    4TN      TN  880  -0.240   0.717   0.090
+    4TN      TN  881  -0.796   0.250   0.507
+    4TN      TN  882  -0.521   0.014  -0.160
+    4TN      TN  883   0.603   0.209   0.347
+    4TN      TN  884  -0.319   0.171   0.544
+    4TN      TN  885  -0.630  -0.461  -0.028
+    4TN      TN  886  -0.536   0.795  -0.273
+    4TN      TN  887   0.572   0.180  -0.222
+    4TN      TN  888  -0.544   0.280   0.386
+    4TN      TN  889   0.866   0.238   0.156
+    4TN      TN  890   0.350   0.148  -0.265
+    4TN      TN  891  -0.677  -0.051   0.022
+    4TN      TN  892   0.058  -0.269  -0.367
+    4TN      TN  893  -0.133   0.658  -0.340
+    4TN      TN  894   0.306  -0.495  -0.638
+    4TN      TN  895  -0.360   0.454  -0.711
+    4TN      TN  896  -0.384  -0.203  -0.122
+    4TN      TN  897   0.169  -0.475   0.021
+    4TN      TN  898  -0.605   0.411  -0.518
+    4TN      TN  899   0.567  -0.761   0.089
+    4TN      TN  900   0.417   0.105  -0.086
+    4TN      TN  901   0.253  -0.854   0.131
+    4TN      TN  902  -0.472   0.736   0.476
+    4TN      TN  903   0.051  -0.250  -0.512
+    4TN      TN  904   0.223   0.239  -0.159
+    4TN      TN  905  -0.055   0.450  -0.354
+    4TN      TN  906  -0.007   0.934  -0.198
+    4TN      TN  907  -0.515  -0.332  -0.423
+    4TN      TN  908   0.322  -0.322   0.456
+    4TN      TN  909   0.771  -0.213  -0.197
+    4TN      TN  910   0.628   0.150   0.321
+    4TN      TN  911   0.663  -0.224  -0.473
+    4TN      TN  912   0.167   0.141  -0.201
+    4TN      TN  913   0.834   0.072   0.217
+    4TN      TN  914  -0.774  -0.166  -0.421
+    4TN      TN  915  -0.253   0.885   0.211
+    4TN      TN  916  -0.324   0.065  -0.292
+    4TN      TN  917  -0.206   0.371  -0.703
+    4TN      TN  918  -0.568   0.285  -0.233
+    4TN      TN  919  -0.646  -0.079  -0.087
+    4TN      TN  920   0.274  -0.003  -0.078
+    4TN      TN  921   0.505  -0.070  -0.315
+    4TN      TN  922   0.608   0.297   0.098
+    4TN      TN  923  -0.263  -0.173   0.935
+    4TN      TN  924  -0.269  -0.734  -0.351
+    4TN      TN  925   0.412   0.291  -0.425
+    4TN      TN  926  -0.375  -0.090  -0.437
+    4TN      TN  927  -0.584   0.339   0.402
+    4TN      TN  928  -0.071  -0.761   0.083
+    4TN      TN  929   0.052   0.379  -0.678
+    4TN      TN  930   0.794  -0.234  -0.309
+    4TN      TN  931  -0.388  -0.090  -0.793
+    4TN      TN  932  -0.336  -0.835  -0.120
+    4TN      TN  933  -0.437  -0.258  -0.509
+    4TN      TN  934   0.686  -0.503  -0.228
+    4TN      TN  935   0.270   0.295  -0.803
+    4TN      TN  936  -0.440  -0.228   0.361
+    4TN      TN  937   0.668   0.228   0.595
+    4TN      TN  938   0.150  -0.894   0.399
+    4TN      TN  939   0.348   0.058   0.378
+    4TN      TN  940   0.542   0.015  -0.162
+    4TN      TN  941  -0.053   0.501  -0.094
+    4TN      TN  942  -0.205   0.183  -0.450
+    4TN      TN  943  -0.434   0.211  -0.615
+    4TN      TN  944   0.901  -0.015  -0.212
+    4TN      TN  945   0.402  -0.153   0.025
+    4TN      TN  946   0.835   0.128   0.025
+    4TN      TN  947   0.886   0.066  -0.289
+    4TN      TN  948  -0.393   0.561   0.447
+    4TN      TN  949   0.027   0.813  -0.426
+    4TN      TN  950  -0.103   0.759  -0.546
+    4TN      TN  951  -0.523   0.141   0.257
+    4TN      TN  952   0.083  -0.456  -0.303
+    4TN      TN  953  -0.362   0.031  -0.482
+    4TN      TN  954  -0.383  -0.361   0.460
+    4TN      TN  955  -0.620  -0.165   0.072
+    4TN      TN  956   0.715  -0.382   0.105
+    4TN      TN  957  -0.710  -0.076  -0.067
+    4TN      TN  958   0.165  -0.300  -0.559
+    4TN      TN  959   0.686  -0.387  -0.391
+    4TN      TN  960   0.632  -0.042   0.606
+    4TN      TN  961   0.035  -0.198   0.481
+    4TN      TN  962  -0.025  -0.844   0.437
+    4TN      TN  963   0.191   0.068  -0.462
+    4TN      TN  964  -0.233  -0.177  -0.773
+    4TN      TN  965  -0.551  -0.690  -0.299
+    4TN      TN  966  -0.249   0.931   0.213
+    4TN      TN  967  -0.728   0.488  -0.424
+    4TN      TN  968   0.501   0.057   0.747
+    4TN      TN  969   0.421  -0.583  -0.557
+    4TN      TN  970   0.004   0.509  -0.618
+    4TN      TN  971  -0.177  -0.044   0.874
+    4TN      TN  972  -0.531   0.545  -0.443
+    4TN      TN  973  -0.435   0.006  -0.350
+    4TN      TN  974   0.962  -0.077   0.134
+    4TN      TN  975   0.727  -0.279   0.465
+    4TN      TN  976   0.154  -0.848  -0.227
+    4TN      TN  977  -0.331  -0.433  -0.737
+    4TN      TN  978  -0.537  -0.204  -0.262
+    4TN      TN  979   0.009  -0.301   0.698
+    4TN      TN  980   0.082  -0.598   0.235
+    4TN      TN  981  -0.214  -0.468  -0.246
+    4TN      TN  982   0.079  -0.431   0.284
+    4TN      TN  983   0.188   0.174  -0.354
+    4TN      TN  984  -0.612   0.059   0.183
+    4TN      TN  985  -0.602  -0.090  -0.613
+    4TN      TN  986   0.858   0.046  -0.425
+    4TN      TN  987  -0.743  -0.119   0.106
+    4TN      TN  988   0.282  -0.239   0.692
+    4TN      TN  989   0.019   0.208  -0.145
+    4TN      TN  990  -0.651  -0.058  -0.672
+    4TN      TN  991  -0.356  -0.840   0.373
+    4TN      TN  992  -0.083  -0.487   0.152
+    4TN      TN  993   0.151  -0.164   0.897
+    4TN      TN  994  -0.247  -0.680  -0.373
+    4TN      TN  995  -0.432   0.083   0.542
+    4TN      TN  996   0.734  -0.217   0.015
+    4TN      TN  997   0.139  -0.945  -0.191
+    4TN      TN  998  -0.022  -0.692  -0.319
+    4TN      TN  999  -0.221   0.354  -0.835
+    4TN      TN 1000  -0.334  -0.051   0.533
+    4TN      TN 1001   0.394   0.496  -0.512
+    4TN      TN 1002   0.144   0.728  -0.534
+    4TN      TN 1003  -0.795   0.033  -0.072
+    4TN      TN 1004  -0.455   0.389  -0.575
+    4TN      TN 1005  -0.496  -0.601  -0.587
+    4TN      TN 1006  -0.280  -0.248   0.762
+    4TN      TN 1007  -0.480  -0.274  -0.592
+    4TN      TN 1008   0.437   0.602  -0.667
+    4TN      TN 1009  -0.849  -0.309  -0.377
+    4TN      TN 1010  -0.394   0.158   0.886
+    4TN      TN 1011  -0.286   0.322  -0.359
+    4TN      TN 1012  -0.802  -0.041   0.030
+    4TN      TN 1013   0.010  -0.616   0.363
+    4TN      TN 1014  -0.288  -0.636   0.480
+    4TN      TN 1015  -0.236  -0.386  -0.425
+    4TN      TN 1016   0.314   0.393  -0.296
+    4TN      TN 1017   0.666  -0.458   0.493
+    4TN      TN 1018  -0.607   0.453   0.484
+    4TN      TN 1019  -0.423   0.011   0.258
+    4TN      TN 1020   0.681  -0.023   0.377
+    4TN      TN 1021   0.324  -0.328   0.805
+    4TN      TN 1022  -0.264   0.762  -0.099
+    4TN      TN 1023  -0.261   0.052   0.727
+    4TN      TN 1024   0.643  -0.114   0.450
+    4TN      TN 1025   0.514   0.733  -0.330
+    4TN      TN 1026  -0.474   0.085   0.415
+    4TN      TN 1027  -0.344   0.349   0.323
+    4TN      TN 1028  -0.391   0.170   0.513
+    4TN      TN 1029   0.230  -0.372  -0.206
+    4TN      TN 1030   0.821  -0.074   0.257
+    4TN      TN 1031  -0.369   0.030  -0.602
+  10.00000  10.00000  10.00000
Local GROMACS mirror with custom stuff