Treecode Test Runs: procedure and sample outputs.
=================================================

1.  Generate test data:

	*   Compile and run program testdata.f,

	*   select test 3 (off-axis parabolic collision),

	*   specify 4096 particles, and

	*   use random number seed 12345.

	Because testdata.f uses a portable random number generator,
	it should yield nearly identical results on any computer with
	integers at least 32 bits long.

	Instead of copying and running testdata.f, you can copy the
        files it produces when run with the above parameters; these
	are in testdata/treebodi.Z and testdata/treepars.

2.  Compile simulation code.

	Results are available for

	*   code_c, compiled using gcc,

	*   code_f77, compiled using the Sun f77 compiler, and

	*   code_f77vs, compiled using the Sun f77 compiler.

	The results given below were obtained with the compilation
	flags in the supplied makefiles.  When using the Sun f77
	compiler, a significant speed-up can be obtained with the
	flags -fast -O3.  These flags were not used in these test
	runs because they apparently result in floating-point math
	which does not exactly conform to the IEEE standard.  In
	production calculations the benefit of this speed-up far
	outweighs the slight loss of standardized behavior.

3.  Run simulation code.

	For the C code, use the treebodi file generated by testdata.f
	as the input file.  This simulation should be run for 24 time
	units; default values should be used for all other parameters.

	For either of the F77 codes, the treebodi and treepars files
	generated by testdata.f provide the necessary input data.

4.  Compare output files.

	Listed below are excerpts from the log files for all three
	versions of the code tested.  Results are given every two
	or three time units to keep this file to a reasonable size.

	Copies of the final particle outputs for all three codes are
	available in testdata/treebodf_c.Z, testdata/treebodf_f77.Z,
	and testdata/treebodf_f77vs.Z.

Output file for c code:
=======================


Hierarchical code: input file treebodi

     nbody     dtime       eps     theta   usequad     dtout     tstop
      4096   0.03125    0.0250      1.00     false    0.2500   24.0000


      tnow       T+U       T/U     nttot     nbavg     ncavg   cputime
     0.000   -0.4952   -0.6007   1280049       101       211      0.42

              cm pos   -0.0000   -0.0000    0.0000
              cm vel   -0.0000    0.0000   -0.0000
              am vec    0.0137   -0.0112   -0.5906

      tnow       T+U       T/U     nttot     nbavg     ncavg   cputime
     2.000   -0.4948   -0.6720   1356785       113       217     19.75

              cm pos    0.0002   -0.0001    0.0001
              cm vel    0.0001   -0.0001    0.0001
              am vec    0.0139   -0.0106   -0.5914

      tnow       T+U       T/U     nttot     nbavg     ncavg   cputime
     4.000   -0.4937   -0.7302   1325922       107       216     39.68

              cm pos    0.0006    0.0001    0.0005
              cm vel    0.0002    0.0003    0.0002
              am vec    0.0139   -0.0114   -0.5921

      tnow       T+U       T/U     nttot     nbavg     ncavg   cputime
     6.000   -0.4908   -0.5527   1157641        93       189     57.82

              cm pos    0.0004    0.0012    0.0010
              cm vel   -0.0001    0.0005    0.0002
              am vec    0.0139   -0.0118   -0.5917

      tnow       T+U       T/U     nttot     nbavg     ncavg   cputime
     8.000   -0.4901   -0.4782   1135117        80       196     74.74

              cm pos   -0.0001    0.0019    0.0013
              cm vel   -0.0003    0.0004    0.0001
              am vec    0.0138   -0.0120   -0.5913

      tnow       T+U       T/U     nttot     nbavg     ncavg   cputime
    10.000   -0.4904   -0.4767   1206863        78       216     92.25

              cm pos   -0.0005    0.0027    0.0016
              cm vel   -0.0002    0.0005    0.0002
              am vec    0.0139   -0.0119   -0.5907

      tnow       T+U       T/U     nttot     nbavg     ncavg   cputime
    12.000   -0.4912   -0.5847   1261317        84       222    110.17

              cm pos   -0.0014    0.0038    0.0020
              cm vel   -0.0009    0.0006    0.0002
              am vec    0.0143   -0.0116   -0.5903

      tnow       T+U       T/U     nttot     nbavg     ncavg   cputime
    14.000   -0.4899   -0.5475   1274093        79       231    128.89

              cm pos   -0.0035    0.0045    0.0028
              cm vel   -0.0011    0.0003    0.0005
              am vec    0.0142   -0.0116   -0.5898

      tnow       T+U       T/U     nttot     nbavg     ncavg   cputime
    16.000   -0.4894   -0.5101   1287971        73       240    147.64

              cm pos   -0.0056    0.0048    0.0038
              cm vel   -0.0010   -0.0001    0.0007
              am vec    0.0139   -0.0116   -0.5893

      tnow       T+U       T/U     nttot     nbavg     ncavg   cputime
    18.000   -0.4891   -0.5096   1340893        78       248    167.29

              cm pos   -0.0076    0.0043    0.0054
              cm vel   -0.0013   -0.0002    0.0008
              am vec    0.0140   -0.0116   -0.5886

      tnow       T+U       T/U     nttot     nbavg     ncavg   cputime
    20.000   -0.4895   -0.5054   1326785        65       258    186.89

              cm pos   -0.0104    0.0038    0.0072
              cm vel   -0.0017   -0.0002    0.0009
              am vec    0.0135   -0.0117   -0.5883

      tnow       T+U       T/U     nttot     nbavg     ncavg   cputime
    22.000   -0.4893   -0.5143   1344037        65       262    206.47

              cm pos   -0.0141    0.0034    0.0087
              cm vel   -0.0020   -0.0002    0.0007
              am vec    0.0132   -0.0119   -0.5882

      tnow       T+U       T/U     nttot     nbavg     ncavg   cputime
    24.000   -0.4887   -0.4972   1363251        77       255    226.30

              cm pos   -0.0187    0.0033    0.0102
              cm vel   -0.0026   -0.0001    0.0009
              am vec    0.0131   -0.0119   -0.5880

Output file for f77 code:
=========================


 ************************************************************************
 *                                                                      *
 *                   OFF-AXIS PARABOLIC COLLISION                       *
 *                                                                      *
 ************************************************************************
 *                                                                      *
 *    NBODY =   4096    NSTEPS = 768    NOUT =   8    DTIME =  0.031    *
 *                                                                      *
 *    USEQUAD = F       EPS = 0.025     THETA = 1.00                    *
 *                                                                      *
 *    OPTIONS:  IMPROVED ACCEPTANCE CRITERION                           *
 *                                                                      *
 ************************************************************************

      CPU TIME USED:          MKTREE           ACCEL           TOTAL
        (SECONDS)               0.53           25.85           26.38



      TIME    NCELL     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
     0.000     2005    415633    864416     431     312     0.49


             ETOT            EKIN            EPOT            MASS
      -0.49519456      0.74482475     -1.24001931      2.00000000


                                X               Y               Z
             JTOT    0.136990E-01   -0.112055E-01   -0.590601E+00
            CMPOS   -0.232959E-08   -0.645760E-08    0.605388E-08
            CMVEL   -0.275682E-08    0.201258E-08   -0.534061E-08



      TIME    NCELL     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
     3.000     1980    469673    908766     452     336    45.25


             ETOT            EKIN            EPOT            MASS
      -0.49588402      1.49206731     -1.98795133      2.00000000


                                X               Y               Z
             JTOT    0.140775E-01   -0.110357E-01   -0.591705E+00
            CMPOS    0.395534E-03   -0.877114E-04    0.274383E-03
            CMVEL    0.267150E-03    0.566982E-04    0.148868E-03



      TIME    NCELL     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
     6.000     1992    383134    774507     415     282    88.36


             ETOT            EKIN            EPOT            MASS
      -0.49076807      0.60635389     -1.09712196      2.00000000


                                X               Y               Z
             JTOT    0.138876E-01   -0.118438E-01   -0.591716E+00
            CMPOS    0.416759E-03    0.116280E-02    0.987755E-03
            CMVEL   -0.115034E-03    0.484350E-03    0.151270E-03



      TIME    NCELL     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
     9.000     2045    343172    837296     464     288   127.09


             ETOT            EKIN            EPOT            MASS
      -0.49027172      0.44109990     -0.93137162      2.00000000


                                X               Y               Z
             JTOT    0.139822E-01   -0.116935E-01   -0.590993E+00
            CMPOS   -0.273633E-03    0.231222E-02    0.145751E-02
            CMVEL   -0.162535E-03    0.341463E-03    0.181434E-03



      TIME    NCELL     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
    12.000     1997    348158    913159     466     307   167.91


             ETOT            EKIN            EPOT            MASS
      -0.49123695      0.69161398     -1.18285094      2.00000000


                                X               Y               Z
             JTOT    0.142662E-01   -0.115921E-01   -0.590334E+00
            CMPOS   -0.143990E-02    0.375875E-02    0.199977E-02
            CMVEL   -0.904272E-03    0.605516E-03    0.208400E-03



      TIME    NCELL     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
    15.000     2031    312716    956343     463     309   210.95


             ETOT            EKIN            EPOT            MASS
      -0.48977183      0.53115676     -1.02092858      2.00000000


                                X               Y               Z
             JTOT    0.140191E-01   -0.116358E-01   -0.589584E+00
            CMPOS   -0.459956E-02    0.483009E-02    0.327451E-02
            CMVEL   -0.975260E-03    0.112764E-03    0.420869E-03



      TIME    NCELL     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
    18.000     2030    321554   1019339     502     327   254.55


             ETOT            EKIN            EPOT            MASS
      -0.48906365      0.50813019     -0.99719384      2.00000000


                                X               Y               Z
             JTOT    0.139505E-01   -0.116020E-01   -0.588638E+00
            CMPOS   -0.762772E-02    0.432624E-02    0.535228E-02
            CMVEL   -0.126711E-02   -0.246431E-03    0.849218E-03



      TIME    NCELL     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
    21.000     2012    271448   1070537     504     327   297.90


             ETOT            EKIN            EPOT            MASS
      -0.48922213      0.50960476     -0.99882689      2.00000000


                                X               Y               Z
             JTOT    0.133806E-01   -0.117888E-01   -0.588226E+00
            CMPOS   -0.121959E-01    0.357830E-02    0.806369E-02
            CMVEL   -0.180108E-02   -0.219636E-03    0.700039E-03



      TIME    NCELL     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
    24.000     2002    317569   1045682     497     332   341.86


             ETOT            EKIN            EPOT            MASS
      -0.48874727      0.48337979     -0.97212706      2.00000000


                                X               Y               Z
             JTOT    0.131010E-01   -0.119031E-01   -0.587990E+00
            CMPOS   -0.186730E-01    0.327638E-02    0.101806E-01
            CMVEL   -0.260708E-02   -0.539074E-04    0.858967E-03


 ELAPSED CPU TIME (SECONDS) :         20519.78


 CALCULATION TERMINATED

Output file for f77vs code:
===========================


 ************************************************************************
 *                                                                      *
 *                   OFF-AXIS PARABOLIC COLLISION                       *
 *                                                                      *
 ************************************************************************
 *                                                                      *
 *    NBODY =   4096    NSTEPS = 768    NOUT =   8    DTIME =  0.031    *
 *                                                                      *
 *    USEQUAD = F       EPS = 0.025     THETA = 1.00  NCRIT =  8        *
 *                                                                      *
 *    OPTIONS:  IMPROVED ACCEPTANCE CRITERION                           *
 *              VECTORIZED FORCE CALCULATION                            *
 *              SHARED INTERACTION LISTS                                *
 *                                                                      *
 ************************************************************************

      CPU TIME USED:          MKTREE           ACCEL           TOTAL
        (SECONDS)               0.80           32.78           33.58



      TIME    NCELL    NWALK     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
     0.000     2005     1543    695695   1001744     600     414     0.60


             ETOT            EKIN            EPOT            MASS
      -0.49525674      0.74482475     -1.24008149      2.00000000


                                X               Y               Z
             JTOT    0.136990E-01   -0.112055E-01   -0.590601E+00
            CMPOS   -0.232959E-08   -0.645760E-08    0.605388E-08
            CMVEL   -0.275682E-08    0.201258E-08   -0.534061E-08



      TIME    NCELL    NWALK     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
     3.000     1972     1537    760353   1092690     650     452    59.60


             ETOT            EKIN            EPOT            MASS
      -0.49587231      1.49389705     -1.98976936      2.00000000


                                X               Y               Z
             JTOT    0.139691E-01   -0.109249E-01   -0.590967E+00
            CMPOS    0.320845E-03    0.706815E-04    0.261589E-03
            CMVEL    0.184032E-03    0.121970E-03    0.941926E-04



      TIME    NCELL    NWALK     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
     6.000     2040     1418    663677    941193     576     391   116.84


             ETOT            EKIN            EPOT            MASS
      -0.49143437      0.60582622     -1.09726059      2.00000000


                                X               Y               Z
             JTOT    0.136362E-01   -0.116420E-01   -0.590923E+00
            CMPOS    0.298247E-03    0.122615E-02    0.854659E-03
            CMVEL   -0.103009E-03    0.397338E-03    0.104080E-03



      TIME    NCELL    NWALK     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
     9.000     2056     1528    615092    985772     591     390   169.99


             ETOT            EKIN            EPOT            MASS
      -0.49131201      0.44395392     -0.93526592      2.00000000


                                X               Y               Z
             JTOT    0.133558E-01   -0.117876E-01   -0.590557E+00
            CMPOS   -0.111103E-03    0.239759E-02    0.116498E-02
            CMVEL    0.787582E-04    0.306918E-03    0.743887E-04



      TIME    NCELL    NWALK     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
    12.000     2021     1452    619593   1073098     646     413   225.18


             ETOT            EKIN            EPOT            MASS
      -0.49222251      0.72787289     -1.22009540      2.00000000


                                X               Y               Z
             JTOT    0.131380E-01   -0.120884E-01   -0.590118E+00
            CMPOS   -0.383021E-03    0.331688E-02    0.122757E-02
            CMVEL   -0.391726E-03    0.244461E-03   -0.489363E-04



      TIME    NCELL    NWALK     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
    15.000     2016     1513    587015   1102080     609     412   280.72


             ETOT            EKIN            EPOT            MASS
      -0.49115105      0.53816310     -1.02931415      2.00000000


                                X               Y               Z
             JTOT    0.132909E-01   -0.122790E-01   -0.589711E+00
            CMPOS   -0.200323E-02    0.365158E-02    0.156386E-02
            CMVEL   -0.687920E-03    0.205579E-03    0.281508E-03



      TIME    NCELL    NWALK     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
    18.000     2026     1562    575112   1193974     658     431   337.57


             ETOT            EKIN            EPOT            MASS
      -0.49071363      0.49310492     -0.98381855      2.00000000


                                X               Y               Z
             JTOT    0.133793E-01   -0.122721E-01   -0.589400E+00
            CMPOS   -0.413122E-02    0.460767E-02    0.266036E-02
            CMVEL   -0.768484E-03    0.465876E-03    0.312542E-03



      TIME    NCELL    NWALK     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
    21.000     2030     1503    551725   1202457     654     428   393.88


             ETOT            EKIN            EPOT            MASS
      -0.49080282      0.49854121     -0.98934403      2.00000000


                                X               Y               Z
             JTOT    0.133969E-01   -0.122925E-01   -0.589172E+00
            CMPOS   -0.620984E-02    0.566706E-02    0.399186E-02
            CMVEL   -0.668664E-03    0.189953E-03    0.486575E-03



      TIME    NCELL    NWALK     NBTOT     NCTOT   NTMAX   NTAVG  CPUTIME
    24.000     1998     1494    592505   1226799     675     444   450.20


             ETOT            EKIN            EPOT            MASS
      -0.49047028      0.50807083     -0.99854111      2.00000000


                                X               Y               Z
             JTOT    0.135297E-01   -0.123606E-01   -0.589019E+00
            CMPOS   -0.879267E-02    0.577577E-02    0.604823E-02
            CMVEL   -0.102523E-02    0.360163E-04    0.976902E-03


 ELAPSED CPU TIME (SECONDS) :         27019.93


 CALCULATION TERMINATED
