So far, three programs which make use of GRAPE-6 exist.
grape6 is a simple shared-timestep direct-summation program. Its source files live in /usr2/makino/src/grape. The makefile for this program is Makefile.grape6. The source file is not too easy to read, since it's the result of the evolution since the days of GRAPE-1...
Anyway, you can run the program by:
/usr2/makino/src/grape/grape6 /usr2/makino/src/grape/testparm2x
The output would look like:
(read_config_file) MAXCHIP, NCHIP = 16 16, CHIPS:
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
testing LED ...
LED test end.
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_with_fifo_2p.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_2p.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_2p.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/pb_jp/cbfin.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_with_fifo_rh.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_with_fifo_rh.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_with_fifo_rh.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_with_fifo_rh.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_with_fifo_rh.ttf
enter n, dt, tstop, dtout, dtsnapout, eps, pesample:
LONG F format read
Enter out_snap file name( ~ for no output)
enter scales for x and v:enter power index for rho = r\^\{-p\}:enter omega:CPU sec.= 7.43
T= 0.000 E = -0.266529909
DE= 0.00000000e+00 V.R. = 0.244443
initialize_grape6, retcode = 0
argc = 2, argv[1]=/usr2/makino/src/grape/testparm2x
Readinf DOUBLE parameters....
n=8192 dt= 0.0312 tstop= 5.000 dtout= 0.250 dtsnapout= 0.031
eps= 0.25000 pesample=8
output snap = ~
x, v scale factor = 1.000000 0.200000 power = 0.000000 omega=0.800000
Etot = -0.404249 1.616996 0.786404
Enter diag:CPU sec.= 7.43
T= 0.000 E = -0.266529909
DE= 0.00000000e+00 V.R. = 0.244443
CM : 1.87567e-17 -2.35949e-17 2.78504e-17
CMV:-1.39456e-17 -6.74916e-18 -8.18234e-19
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 7.44
CPU sec.= 8.52
T= 0.250 E = -0.266529875
DE=-1.28383736e-07 V.R. = 0.246534
Enter diag:CPU sec.= 8.52
T= 0.250 E = -0.266529875
DE=-1.28383736e-07 V.R. = 0.246534
CM : 2.64545e-17 -3.08998e-17 2.51061e-17
CMV:-2.38524e-18 -1.23057e-17 -1.66357e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 8.52
CPU sec.= 9.60
T= 0.500 E = -0.266529766
DE=-5.34901433e-07 V.R. = 0.253028
Enter diag:CPU sec.= 9.60
T= 0.500 E = -0.266529766
DE=-5.34901433e-07 V.R. = 0.253028
CM : 8.55164e-18 -4.96294e-17 3.81842e-17
CMV:-1.08556e-17 8.26704e-19 -5.11947e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 9.61
CPU sec.= 10.69
T= 0.750 E = -0.266529566
DE=-1.28686658e-06 V.R. = 0.264026
Enter diag:CPU sec.= 10.69
T= 0.750 E = -0.266529566
DE=-1.28686658e-06 V.R. = 0.264026
CM : 1.21295e-17 -3.98173e-17 1.91633e-17
CMV:-1.27123e-17 -1.13164e-17 4.09964e-19
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 10.69
CPU sec.= 11.78
T= 1.000 E = -0.266529238
DE=-2.51785768e-06 V.R. = 0.279704
Enter diag:CPU sec.= 11.78
T= 1.000 E = -0.266529238
DE=-2.51785768e-06 V.R. = 0.279704
CM : 1.30646e-17 -3.74321e-17 2.04914e-17
CMV: 5.57009e-18 -7.43356e-18 -1.31121e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 11.78
CPU sec.= 12.86
T= 1.250 E = -0.266528718
DE=-4.46870731e-06 V.R. = 0.300321
Enter diag:CPU sec.= 12.86
T= 1.250 E = -0.266528718
DE=-4.46870731e-06 V.R. = 0.300321
CM : 8.13152e-20 -5.15267e-17 2.80131e-17
CMV:-7.42678e-18 -7.54198e-18 2.98494e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 12.87
CPU sec.= 13.95
T= 1.500 E = -0.266527893
DE=-7.56231468e-06 V.R. = 0.326237
Enter diag:CPU sec.= 13.95
T= 1.500 E = -0.266527893
DE=-7.56231468e-06 V.R. = 0.326237
CM : 6.92873e-18 -4.38560e-17 2.48791e-17
CMV:-9.89334e-18 -5.29226e-18 -2.58514e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 13.95
CPU sec.= 15.04
T= 1.750 E = -0.266526574
DE=-1.25127809e-05 V.R. = 0.357927
Enter diag:CPU sec.= 15.04
T= 1.750 E = -0.266526574
DE=-1.25127809e-05 V.R. = 0.357927
CM : 9.88657e-18 -1.52330e-17 2.93226e-17
CMV:-4.43168e-18 -1.79097e-17 -6.33581e-19
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 15.04
CPU sec.= 16.12
T= 2.000 E = -0.266524499
DE=-2.02997490e-05 V.R. = 0.395985
Enter diag:CPU sec.= 16.12
T= 2.000 E = -0.266524499
DE=-2.02997490e-05 V.R. = 0.395985
CM : 2.45563e-17 -6.04172e-17 2.43209e-17
CMV: 1.04354e-18 2.92057e-18 -1.01712e-17
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 16.13
CPU sec.= 17.21
T= 2.250 E = -0.266521847
DE=-3.02470587e-05 V.R. = 0.440989
Enter diag:CPU sec.= 17.21
T= 2.250 E = -0.266521847
DE=-3.02470587e-05 V.R. = 0.440989
CM : 1.48875e-17 -5.50775e-17 2.72044e-17
CMV:-1.58836e-17 -3.45589e-19 -4.31648e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 17.21
CPU sec.= 18.29
T= 2.500 E = -0.266524723
DE=-1.94587790e-05 V.R. = 0.492597
Enter diag:CPU sec.= 18.29
T= 2.500 E = -0.266524723
DE=-1.94587790e-05 V.R. = 0.492597
CM :-6.35614e-18 -5.10117e-17 1.12926e-17
CMV: 4.55365e-18 -1.06049e-17 1.79571e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 18.30
CPU sec.= 19.38
T= 2.750 E = -0.266569675
DE= 1.49178145e-04 V.R. = 0.544339
Enter diag:CPU sec.= 19.38
T= 2.750 E = -0.266569675
DE= 1.49178145e-04 V.R. = 0.544339
CM :-3.93023e-18 -5.45760e-17 1.72761e-17
CMV:-1.30985e-17 5.14996e-18 3.81842e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 19.38
CPU sec.= 20.46
T= 3.000 E = -0.266596553
DE= 2.49980823e-04 V.R. = 0.576093
Enter diag:CPU sec.= 20.46
T= 3.000 E = -0.266596553
DE= 2.49980823e-04 V.R. = 0.576093
CM :-8.60585e-18 -7.16929e-17 2.26056e-17
CMV:-1.44199e-17 -6.89824e-18 -1.09132e-17
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 20.47
CPU sec.= 21.55
T= 3.250 E = -0.266566256
DE= 1.36352907e-04 V.R. = 0.586571
Enter diag:CPU sec.= 21.55
T= 3.250 E = -0.266566256
DE= 1.36352907e-04 V.R. = 0.586571
CM : 1.58971e-17 -5.20010e-17 4.50994e-17
CMV: 4.58075e-18 4.91957e-18 7.93839e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 21.55
CPU sec.= 22.63
T= 3.500 E = -0.266572612
DE= 1.60192107e-04 V.R. = 0.591606
Enter diag:CPU sec.= 22.63
T= 3.500 E = -0.266572612
DE= 1.60192107e-04 V.R. = 0.591606
CM : 2.15079e-17 -7.04122e-17 2.40998e-17
CMV:-1.62224e-17 -4.33681e-18 5.67512e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 22.64
CPU sec.= 23.72
T= 3.750 E = -0.266583967
DE= 2.02781044e-04 V.R. = 0.591086
Enter diag:CPU sec.= 23.72
T= 3.750 E = -0.266583967
DE= 2.02781044e-04 V.R. = 0.591086
CM : 5.27193e-18 -7.40493e-17 1.45080e-17
CMV:-4.45878e-18 -8.29415e-18 -5.85808e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 23.72
CPU sec.= 24.81
T= 4.000 E = -0.266573960
DE= 1.65247645e-04 V.R. = 0.580252
Enter diag:CPU sec.= 24.81
T= 4.000 E = -0.266573960
DE= 1.65247645e-04 V.R. = 0.580252
CM :-9.96111e-18 -7.44711e-17 1.76318e-17
CMV:-7.45389e-19 -6.91179e-18 -2.19890e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 24.81
CPU sec.= 25.90
T= 4.250 E = -0.266555207
DE= 9.49058810e-05 V.R. = 0.560385
Enter diag:CPU sec.= 25.90
T= 4.250 E = -0.266555207
DE= 9.49058810e-05 V.R. = 0.560385
CM : 1.63579e-17 -6.38460e-17 2.36339e-17
CMV:-2.85145e-17 -9.93400e-18 5.25499e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 25.90
CPU sec.= 26.98
T= 4.500 E = -0.266542816
DE= 4.84248779e-05 V.R. = 0.536870
Enter diag:CPU sec.= 26.98
T= 4.500 E = -0.266542816
DE= 4.84248779e-05 V.R. = 0.536870
CM :-7.83336e-18 -7.74120e-17 2.45445e-17
CMV:-9.85269e-18 -1.29020e-17 -4.65190e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 26.98
CPU sec.= 28.07
T= 4.750 E = -0.266538246
DE= 3.12788341e-05 V.R. = 0.514305
Enter diag:CPU sec.= 28.07
T= 4.750 E = -0.266538246
DE= 3.12788341e-05 V.R. = 0.514305
CM :-1.46367e-18 -6.38053e-17 1.26218e-17
CMV:-5.25838e-18 -5.43456e-18 -1.29867e-17
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 28.07
CPU sec.= 29.15
T= 5.000 E = -0.266539550
DE= 3.61696770e-05 V.R. = 0.494341
Enter diag:CPU sec.= 29.15
T= 5.000 E = -0.266539550
DE= 3.61696770e-05 V.R. = 0.494341
CM :-9.41731e-18 -6.57027e-17 1.44555e-17
CMV:-4.45878e-18 7.45389e-19 2.54449e-18
AM : 2.44037e-03 -8.75096e-04 -4.04957e-01
Exit diag:CPU sec.= 29.16
Errors: jp, ip, ecc, ecc(u), cm = 0 0 0 0 0
nbody1 is a basic individual- (block-) timestep integrator. The particular version of nbody1 with GRAPE-6 support lives in source files live in /usr2/makino/src/bhnbody1. The makefile for this program is Makefile.grape6. The source file is again not too easy to read, since it's the result of the evolution since the days of GRAPE-2...
Anyway, you can run the program by:
/usr2/makino/src/bhnbody1/nbody1_g6 < /usr2/makino/src/bhnbody1/samplein
The output would look like:
G6NCLUSTER = 1
N NBH NRAND ETA DELTAT TCRIT QE CUTOFF
200 0 42 0.02 0.1 0.5 0.00001 0.00
OPTIONS 0 2 0 2 1 0 0 0 0 1
data n = 2048
Body end
X end
V end
exit data
call g6open
(read_config_file) MAXCHIP, NCHIP = 16 16, CHIPS:
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
testing LED ...
LED test end.
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_with_fifo_2p.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_2p.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_2p.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/pb_jp/cbfin.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_with_fifo_rh.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_with_fifo_rh.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_with_fifo_rh.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_with_fifo_rh.ttf
(send_fpga_data) file to open: /usr2/makino/src/grape6board/processor_board/fo_unit_with_fifo_rh.ttf
return g6open
fpoly0, eps2 = 1.000000000000000E-006 npipe = 48
ic, nj = 1 2048
Errors: jp, ip, ecc, ecc(u), cm = 0 0 0 0 0
T = 0.000 Q = 0.50 STEPS = 0 0 0
DE= 0.000E+00 0.000E+00 E = -0.25066
CM : -0.130498E-07 -0.978288E-08 -0.265915E-07
CMV : 0.510450E-08 -0.886712E-09 -0.433698E-09
AM : -0.838052E-02 -0.543369E-03 0.823090E-02
new eta = 2.200000000000000E-002
CPU min = 0.122586
Number of pipes = 48
intgrt, 0.125244140625000 0.125000000000000
CPU min = 0.141260
Errors: jp, ip, ecc, ecc(u), cm = 0 0 0 0 0
T = 0.125 Q = 0.50 STEPS = 55332 774 0
DE= 2.535E-08 -5.050E-08 E = -0.25066
CM : 0.196896E-08 -0.416279E-08 -0.218563E-08
CMV : 0.439653E-08 -0.483714E-09 -0.118700E-08
AM : -0.838052E-02 -0.543364E-03 0.823092E-02
new eta = 2.420000000000000E-002
CPU min = 0.141276
Number of pipes = 48
intgrt, 0.250488281250000 0.250000000000000
CPU min = 0.157315
Errors: jp, ip, ecc, ecc(u), cm = 0 0 0 0 0
T = 0.250 Q = 0.50 STEPS = 103030 1600 0
DE= 1.270E-08 -7.582E-08 E = -0.25066
CM : -0.272468E-08 -0.172687E-07 0.233900E-08
CMV : 0.372973E-08 -0.493510E-08 -0.741453E-08
AM : -0.838052E-02 -0.543359E-03 0.823092E-02
new eta = 2.662000000000000E-002
CPU min = 0.157331
Number of pipes = 48
intgrt, 0.375122070312500 0.375000000000000
CPU min = 0.171825
Errors: jp, ip, ecc, ecc(u), cm = 0 0 0 0 0
T = 0.375 Q = 0.50 STEPS = 146129 2218 0
DE= 1.458E-08 -1.049E-07 E = -0.25066
CM : 0.238592E-08 -0.286509E-09 -0.960970E-08
CMV : 0.273282E-08 0.490193E-08 -0.218992E-07
AM : -0.838054E-02 -0.543361E-03 0.823092E-02
new eta = 2.928200000000001E-002
CPU min = 0.171841
Number of pipes = 48
intgrt, 0.500244140625000 0.500000000000000
CPU min = 0.185635
Errors: jp, ip, ecc, ecc(u), cm = 0 0 0 0 0
T = 0.500 Q = 0.50 STEPS = 186832 2713 0
DE= 2.113E-08 -1.471E-07 E = -0.25066
CM : 0.366066E-08 -0.332199E-08 -0.149043E-07
CMV : 0.225550E-08 0.159237E-08 -0.219382E-07
AM : -0.838053E-02 -0.543371E-03 0.823092E-02
new eta = 3.221020000000001E-002
CPU min = 0.185668
Number of pipes = 48
intgrt, 0.625122070312500 0.625000000000000
CPU min = 0.198616
Errors: jp, ip, ecc, ecc(u), cm = 0 0 0 0 0
T = 0.625 Q = 0.50 STEPS = 224995 3309 0
DE= 7.072E-09 -1.612E-07 E = -0.25066
CM : -0.215099E-08 -0.329190E-08 -0.834841E-08
CMV : -0.187221E-08 -0.532395E-08 -0.378905E-07
AM : -0.838054E-02 -0.543368E-03 0.823092E-02
new eta = 3.543122000000001E-002
TIME = 0.63 TCOMP = 0.20 KZ(1) = 0
Kira is a rather fancy
-body integration program specialized
to star clusters, with the capability to handle stellar evolution,
binary evolution, stellar collisions, galactic tidal fields and all
the ``realistic'' additional physics. This program will however need
some separate documentation...