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...