next up previous contents
Next: 13 Known bugs and Up: grape6user Previous: 11 RUN-TIME SUPPORT

Subsections


12 Sample programs

So far, three programs which make use of GRAPE-6 exist.

12.1 grape6

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

12.2 nbody1

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

12.3 Kira

Kira is a rather fancy $N$-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...


next up previous contents
Next: 13 Known bugs and Up: grape6user Previous: 11 RUN-TIME SUPPORT
Jun Makino
平成17年1月31日